基于流體體積模量avo反演的流體識(shí)別的方法及系統(tǒng)的制作方法
【專利摘要】本發(fā)明提供一種基于流體體積模量AVO反演的流體識(shí)別的方法及系統(tǒng)，所述方法包括：采集目的儲(chǔ)層的疊前地震資料、測(cè)井資料以及地質(zhì)成果資料；構(gòu)建表征流體體積模量的AVO近似公式；根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料以及所述的表征流體體積模量的AVO近似公式方程提取子波；根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料以及子波確定流體體積模量數(shù)據(jù)體；根據(jù)所述的流體體積模量數(shù)據(jù)體識(shí)別當(dāng)前儲(chǔ)層孔隙中流體的類型。利用表征流體體積模量的AVO近似公式進(jìn)行AVO反演可以估算出流體體積模量參數(shù)，避免流體識(shí)別假象，提高儲(chǔ)層孔隙流體類型預(yù)測(cè)的可靠性。
【專利說(shuō)明】基于流體體積模量AVO反演的流體識(shí)別的方法及系統(tǒng)

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明關(guān)于石油地球物理勘探領(lǐng)域，特別是關(guān)于疊前地震的反演技術(shù)，具體的講 是一種基于流體體積模量AVO反演的流體識(shí)別的方法及系統(tǒng)。

【背景技術(shù)】
[0002] 從20世紀(jì)70年代開(kāi)始，由于地震技術(shù)數(shù)字化和計(jì)算機(jī)技術(shù)的飛速發(fā)展，利用地震 波動(dòng)力學(xué)問(wèn)題進(jìn)行儲(chǔ)層預(yù)測(cè)和流體識(shí)別成為石油工業(yè)界的研究熱點(diǎn)。其中，AVCKAmplitude Versus Offset，振幅隨偏移距的變化）反演技術(shù)作為一項(xiàng)根據(jù)振幅隨偏移距的變化規(guī)律預(yù) 測(cè)地下巖性及含流體性質(zhì)的主流技術(shù)，其核心思想是通過(guò)研究反射系數(shù)隨入射角度的變化 規(guī)律提取需求的彈性參數(shù)，理論基礎(chǔ)是描述平面波在水平分界面反射和透射的Zoeppritz 方程。由于基于Zoeppritz方程的參數(shù)提取屬非線性反問(wèn)題，計(jì)算效率低和受噪音影響大 等問(wèn)題制約了其在實(shí)際生產(chǎn)中的應(yīng)用。
[0003] 為了克服由Zo印pritz方程導(dǎo)出的反射系數(shù)形式復(fù)雜及不易進(jìn)行數(shù)值計(jì)算的困 難，許多學(xué)者對(duì)Zoeppritz方程進(jìn)行了簡(jiǎn)化。Koefoed將原來(lái)7個(gè)獨(dú)立變量簡(jiǎn)化為5個(gè)獨(dú)立 變量；Bortfeld詳細(xì)論述了垂直入射的平面縱波反射系數(shù)近似計(jì)算方法，并給出了區(qū)分流 體和固體的簡(jiǎn)化方程；Aki和Richards在假設(shè)相鄰地層介質(zhì)彈性參數(shù)變化較小的情況下， 給出了包含縱橫波速度和密度相對(duì)變化量的縱波反射系數(shù)線性近似公式。在此基礎(chǔ)上，許 多學(xué)者對(duì)Aki-Richards方程進(jìn)行了重新推導(dǎo)、歸納，分別提出了以不同彈性參數(shù)表示的縱 波反射系數(shù)。其中，Shuey給出了突出泊松比的相對(duì)反射系數(shù)近似表達(dá)形式；Fatti等人給 出了以相對(duì)波阻抗變化表示的反射系數(shù)近似公式；Mallick給出了用射線參數(shù)表示的反射 系數(shù)近似形式；Goodway推導(dǎo)了以拉梅參數(shù)相對(duì)變化量表示的反射系數(shù)近似公式。
[0004] 隨著石油勘探開(kāi)發(fā)的深入，地下儲(chǔ)層的含流體識(shí)別問(wèn)題成為研究的熱點(diǎn)，迫切需 要從地震資料中提取更多可以反映地下介質(zhì)中的流體特征的參數(shù)。Russell等人結(jié)合飽和 流體多孔彈性介質(zhì)理論對(duì)Aki-Richards近似進(jìn)行重組，提出了包含Gassmann流體項(xiàng)f的 Russell反射系數(shù)近似公式，突出體現(xiàn)了巖石孔隙流體效應(yīng)對(duì)反射系數(shù)的影響，推動(dòng)了基于 AVO反演的流體識(shí)別技術(shù)的發(fā)展。但是由于地下巖石的固液耦合特性，現(xiàn)階段基于AVO反 演方法估算的彈性參數(shù)（如波阻抗，拉梅常數(shù)，Gassmann流體項(xiàng)等）均會(huì)受巖石固體效應(yīng) (巖石基質(zhì)，孔隙度等）影響產(chǎn)生流體識(shí)別假象，降低了地質(zhì)情況復(fù)雜地區(qū)儲(chǔ)層含油氣預(yù)測(cè) 的可靠性。
[0005] 因此，為了有效的提高地質(zhì)情況復(fù)雜地區(qū)的儲(chǔ)層流體表征敏感性和流體識(shí)別可靠 性，現(xiàn)有技術(shù)中急需一種包含對(duì)孔隙流體信息更為敏感彈性參數(shù)的反射系數(shù)近似公式以及 相應(yīng)的AVO反演方法來(lái)識(shí)別儲(chǔ)層流體。


【發(fā)明內(nèi)容】

[0006] 為了克服現(xiàn)有技術(shù)存在的基于常規(guī)AVO近似估算的巖石彈性參數(shù)受固液耦合特 性影響易產(chǎn)生儲(chǔ)層含流體判識(shí)假象的問(wèn)題，本發(fā)明提供了一種基于流體體積模量AVO反演 的流體識(shí)別的方法及系統(tǒng)，在反射系數(shù)近似中實(shí)現(xiàn)流體體積模量參數(shù)分離，流體體積模量 對(duì)儲(chǔ)層流體信息更為敏感，利用表征流體體積模量的AVO近似公式進(jìn)行AVO反演可以估算 出流體體積模量參數(shù)，避免流體識(shí)別假象，提高儲(chǔ)層孔隙流體類型預(yù)測(cè)的可靠性。
[0007] 本發(fā)明的目的之一是，提供一種基于流體體積模量AVO反演的流體識(shí)別的方法， 包括：采集目的儲(chǔ)層的疊前地震資料、測(cè)井資料以及地質(zhì)成果資料；構(gòu)建表征流體體積模 量的AVO近似公式；根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料以及所述的表征流 體體積模量的AVO近似公式方程提取子波；根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果 資料以及子波確定流體體積模量數(shù)據(jù)體；根據(jù)所述的流體體積模量數(shù)據(jù)體識(shí)別當(dāng)前儲(chǔ)層孔 隙中流體的類型。
[0008] 本發(fā)明的目的之一是，提供了一種基于流體體積模量AVO反演的流體識(shí)別的系 統(tǒng)，包括：資料采集裝置，用于采集目的儲(chǔ)層的疊前地震資料、測(cè)井資料以及地質(zhì)成果資料； 近似公式構(gòu)建裝置，用于構(gòu)建表征流體體積模量的AVO近似公式；子波提取裝置，用于根據(jù) 所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料以及所述的表征流體體積模量的AVO近似 公式方程提取子波；流體體積模量數(shù)據(jù)體確定裝置，用于根據(jù)所述的疊前地震資料、測(cè)井資 料、地質(zhì)成果資料以及子波確定流體體積模量數(shù)據(jù)體；流體類型識(shí)別裝置，用于根據(jù)所述的 流體體積模量數(shù)據(jù)體識(shí)別當(dāng)前儲(chǔ)層孔隙中流體的類型。
[0009] 本發(fā)明的有益效果在于，提供了一種基于流體體積模量AVO反演的流體識(shí)別的方 法及系統(tǒng)，基于常規(guī)AVO近似估算的巖石彈性參數(shù)受固液耦合特性影響易產(chǎn)生儲(chǔ)層含流體 判識(shí)假象，本發(fā)明基于孔隙彈性介質(zhì)理論以及臨界孔隙度模型，在反射系數(shù)近似中實(shí)現(xiàn)流 體體積模量參數(shù)分離，流體體積模量對(duì)儲(chǔ)層流體信息更為敏感，且最大限度的消除了儲(chǔ)層 巖石孔隙度等固體效應(yīng)造成的流體識(shí)別假象，利用表征流體體積模量的AVO近似公式進(jìn)行 AVO反演可以估算出流體體積模量參數(shù)，避免流體識(shí)別假象，提高儲(chǔ)層孔隙流體類型預(yù)測(cè)的 可靠性。
[0010] 為讓本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂，下文特舉較佳實(shí)施例， 并配合所附圖式，作詳細(xì)說(shuō)明如下。

【專利附圖】

【附圖說(shuō)明】
[0011] 為了更清楚地說(shuō)明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本 發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以 根據(jù)這些附圖獲得其他的附圖。
[0012] 圖1為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的方法的 實(shí)施方式一的流程圖；
[0013] 圖2為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的方法的 實(shí)施方式二的流程圖；
[0014] 圖3為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的方法的 實(shí)施方式三的流程圖；
[0015] 圖4為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的方法的 實(shí)施方式四的流程圖；
[0016] 圖5為圖1中的步驟S103的具體流程圖；
[0017] 圖6為圖1中的步驟S104的具體流程圖；
[0018]圖7為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的系統(tǒng)的 實(shí)施方式一的結(jié)構(gòu)框圖；
[0019] 圖8為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的系統(tǒng)的 實(shí)施方式二的結(jié)構(gòu)框圖；
[0020] 圖9為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的系統(tǒng)的 實(shí)施方式三的結(jié)構(gòu)框圖；
[0021] 圖10為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的系統(tǒng) 的實(shí)施方式四的結(jié)構(gòu)框圖；
[0022] 圖11為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的系統(tǒng) 中的子波提取裝置300的具體結(jié)構(gòu)框圖；
[0023] 圖12為本發(fā)明實(shí)施例提供的一種基于流體體積模量AVO反演的流體識(shí)別的系統(tǒng) 中的流體體積模量數(shù)據(jù)體確定裝置400的具體結(jié)構(gòu)框圖；
[0024] 圖13為孔隙流體類型為水和氣情況下，縱波速度隨孔隙度與含水飽和度的變化 趨勢(shì)圖；
[0025] 圖14為孔隙流體類型為水和氣情況下，Gassmann流體項(xiàng)隨孔隙度與含水飽和度 的變化趨勢(shì)圖；
[0026] 圖15為孔隙流體類型為水和氣情況下，流體體積模量隨孔隙度與含水飽和度的 變化趨勢(shì)圖；
[0027] 圖16為一維地層模型曲線示意圖；
[0028] 圖17為未加噪聲時(shí)的疊前CMP道集示意圖；
[0029] 圖18為利用未加噪聲時(shí)的疊前CMP道集進(jìn)行反演得到的反演結(jié)果示意圖；
[0030] 圖19為信噪比為1:1的疊前CMP道集示意圖；
[0031] 圖20為利用信噪比為1:1的疊前CMP道集進(jìn)行反演得到的反演結(jié)果示意圖；
[0032] 圖21為研究測(cè)線的疊加地震剖面示意圖；
[0033] 圖22為估算的流體體積模量Kf的剖面示意圖。

【具體實(shí)施方式】
[0034] 下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完 整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；?本發(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他 實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0035] 由于地下巖石的固液耦合特性，現(xiàn)階段基于AVO反演方法估算的彈性參數(shù)（如波 阻抗，Gassmann流體項(xiàng)等）均會(huì)受巖石固體效應(yīng)（巖石基質(zhì)，孔隙度等）影響產(chǎn)生流體識(shí) 別假象，降低了地質(zhì)情況復(fù)雜地區(qū)儲(chǔ)層含油氣預(yù)測(cè)的可靠性。本發(fā)明基于孔隙彈性介質(zhì)理 論以及臨界孔隙度模型推導(dǎo)新的AVO近似公式，在反射系數(shù)近似中實(shí)現(xiàn)流體體積模量參數(shù) 分離，基于新的近似公式進(jìn)行AVO反演可以得到對(duì)儲(chǔ)層流體信息更為敏感的彈性參數(shù)-- 流體體積模量。
[0036] 圖1為本發(fā)明提出的一種基于流體體積模量AVO反演的流體識(shí)別的方法的實(shí)施方 式一的具體流程圖，為了提高流體體積模量估算的穩(wěn)定性，本發(fā)明基于貝葉斯理論建立反 演目標(biāo)函數(shù)，利用重加權(quán)最小二乘迭代算法循環(huán)迭代求解目標(biāo)函數(shù)，實(shí)現(xiàn)地震尺度流體體 積模量的定量估算。由圖1可知，在實(shí)施方式一中，所述的方法包括：
[0037] SlOl :采集目的儲(chǔ)層的疊前地震資料、測(cè)井資料以及地質(zhì)成果資料。
[0038] 疊前地震資料的采集在具體的實(shí)施方式中可通過(guò)如下方式：在勘探工區(qū)中，根據(jù) 目的層的深度和地質(zhì)特點(diǎn)，考慮動(dòng)校正拉伸、干擾波、多次波等因素來(lái)設(shè)計(jì)寬方位角的三維 觀測(cè)系統(tǒng)，保證有足夠的偏移距、方位角。經(jīng)過(guò)激發(fā)、接收，得到可以滿足AVO分析需求的寬 方位大角度的疊前地震資料。
[0039] 測(cè)井資料的采集在具體的實(shí)施方式中可通過(guò)如下方式：在勘探工區(qū)中進(jìn)行全波列 測(cè)井，得到測(cè)井資料，主要包含縱波速度、橫波速度、密度等全波列測(cè)井曲線，孔隙度等解釋 成果曲線，測(cè)井層位，錄井資料等。
[0040] S102 :構(gòu)建表征流體體積模量的AVO近似公式。
[0041] AVO (Amplitude Versus Offset,振幅隨偏移距的變化）近似公式是進(jìn)行AVO反演 的基礎(chǔ)，本發(fā)明基于孔隙彈性介質(zhì)理論以及臨界孔隙度模型，推導(dǎo)得到的表征流體體積模 量的AVO近似公式，其為基于AVO反演的流體識(shí)別提供了理論基礎(chǔ)。
[0042] Russell等人從Aki-Richards近似公式出發(fā)，推導(dǎo)了突出儲(chǔ)層流體特征的反射特 征近似公式：

【權(quán)利要求】
1. 一種基于流體體積模量AVO反演的流體識(shí)別的方法，其特征是，所述的方法具體包 括： 采集目的儲(chǔ)層的疊前地震資料、測(cè)井資料以及地質(zhì)成果資料； 構(gòu)建表征流體體積模量的振幅隨偏移距的變化AVO近似公式； 根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料以及所述的表征流體體積模量的AVO近似公式提取子波； 根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料以及子波確定流體體積模量數(shù)據(jù) 體； 根據(jù)所述的流體體積模量數(shù)據(jù)體識(shí)別當(dāng)前儲(chǔ)層孔隙中流體的類型。
2.根據(jù)權(quán)利要求1所述的方法，其特征是，所述的方法在構(gòu)建表征流體體積模量的AVO 近似公式之前還包括： 對(duì)所述的疊前地震資料進(jìn)行保幅處理，得到共中心點(diǎn)CMP道集資料。
3.根據(jù)權(quán)利要求2所述的方法，其特征是，所述的方法在構(gòu)建表征流體體積模量的AVO 近似公式之前還包括： 對(duì)所述的測(cè)井資料進(jìn)行環(huán)境校正； 對(duì)經(jīng)過(guò)環(huán)境校正后的測(cè)井資料進(jìn)行奇異值消除； 根據(jù)經(jīng)過(guò)奇異值消除后的測(cè)井資料確定所述目的儲(chǔ)層的干巖石縱橫波速度比、飽和巖 石縱橫波速度比； 根據(jù)所述的目的儲(chǔ)層的干巖樣縱橫波速度比以及所述的測(cè)井資料確定與所述測(cè)井資 料對(duì)應(yīng)的流體體積模量曲線以及剪切模量曲線。
4.根據(jù)權(quán)利要求3所述的方法，其特征是，所述的方法在構(gòu)建表征流體體積模量的AVO 近似公式之前還包括： 根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料對(duì)目的儲(chǔ)層進(jìn)行精細(xì)層位追蹤，得 到地震層位數(shù)據(jù)。
5.根據(jù)權(quán)利要求4所述的方法，其特征是，所述的表征流體體積模量的AVO近似公式 為：
其中，fm =φμ為地下巖石的固體項(xiàng)，Θ為入射角度，Kf為地下巖石的流體體積模量， 乙為地下巖石的固體項(xiàng)，Φ為地下巖石的孔隙度，P為地下巖石的密度，μ為地下巖石的 剪切模量，ΛKf為界面兩側(cè)的流體體積模量的差值，Λfm為界面兩側(cè)固體項(xiàng)的差值，ΛP 為界面兩側(cè)的密度的差值，ΛΦ為界面兩側(cè)的孔隙度的差值，= 為干巖石縱 橫波速度比的平方，為飽和巖石縱橫波速度比的平方。
6. 根據(jù)權(quán)利要求5所述的方法，其特征是，根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì) 成果資料以及表征流體體積模量的AVO近似公式提取子波具體包括： 從所述的測(cè)井資料中提取出地下巖石的密度以及孔隙度； 根據(jù)所述的剪切模量曲線以及所述的孔隙度確定地下巖石的固體項(xiàng)； 根據(jù)所述地下巖石的密度、固體項(xiàng)、孔隙度、流體體積模量曲線、表征流體體積模量的AVO近似公式、疊前CMP道集以及地震層位數(shù)據(jù)提取子波。
7.根據(jù)權(quán)利要求6所述的方法，其特征是，根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì) 成果資料以及子波確定流體體積模量數(shù)據(jù)體具體包括： 利用所述的疊前CMP道集、測(cè)井資料、地質(zhì)成果資料以及子波，基于貝葉斯理論，使用 高斯分布函數(shù)為似然函數(shù)、柯西分布為先驗(yàn)約束建立AVO反演目標(biāo)函數(shù)； 利用重加權(quán)最小二乘迭代算法求解所述的AVO反演目標(biāo)函數(shù)，得到流體體積模量的相 對(duì)變化量； 利用道積分思想將所述的相對(duì)變化量轉(zhuǎn)化為流體體積模量數(shù)據(jù)體。
8.根據(jù)權(quán)利要求7所述的方法，其特征是，建立的AVO反演目標(biāo)函數(shù)為： (GtG+nQ+αCTC)m=GTd+aCtξ
其中，巧為CMP道集資料中噪音的方差，σ〗為彈性參數(shù)相對(duì)變化量的方差，λ、α為 加權(quán)系數(shù)，m為待反演的彈性參數(shù)相對(duì)變化量矩陣，G為正演算子矩陣，C為積分矩陣，d為CMP道集矩陣，Q為斜對(duì)角加權(quán)矩陣，M為模型約束矩陣，Mtl為M的均值。
9. 一種基于流體體積模量AVO反演的流體識(shí)別的系統(tǒng)，其特征是，所述的系統(tǒng)具體包 括： 資料采集裝置，用于采集目的儲(chǔ)層的疊前地震資料、測(cè)井資料以及地質(zhì)成果資料； 近似公式構(gòu)建裝置，用于構(gòu)建表征流體體積模量的振幅隨偏移距的變化AVO近似公 式； 子波提取裝置，用于根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料以及所述的表 征流體體積模量的AVO近似公式提取子波； 流體體積模量數(shù)據(jù)體確定裝置，用于根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資 料以及子波確定流體體積模量數(shù)據(jù)體； 流體類型識(shí)別裝置，用于根據(jù)所述的流體體積模量數(shù)據(jù)體識(shí)別當(dāng)前儲(chǔ)層孔隙中流體的 類型。
10. 根據(jù)權(quán)利要求9所述的系統(tǒng)，其特征是，所述的系統(tǒng)還包括： 保幅處理裝置，用于對(duì)所述的疊前地震資料進(jìn)行保幅處理，得到共中心點(diǎn)CMP道集資 料。
11. 根據(jù)權(quán)利要求10所述的系統(tǒng)，其特征是，所述的系統(tǒng)還包括： 環(huán)境校正裝置，用于對(duì)所述的測(cè)井資料進(jìn)行環(huán)境校正； 奇異值消除裝置，用于對(duì)經(jīng)過(guò)環(huán)境校正后的測(cè)井資料進(jìn)行奇異值消除； 縱橫速度比確定裝置，用于根據(jù)經(jīng)過(guò)奇異值消除后的測(cè)井資料確定所述目的儲(chǔ)層的干 巖石縱橫波速度比、飽和巖石縱橫波速度比； 剪切模量曲線確定裝置，用于根據(jù)所述的目的儲(chǔ)層的干巖樣縱橫波速度比以及所述的 測(cè)井資料確定與所述測(cè)井資料對(duì)應(yīng)的流體體積模量曲線以及剪切模量曲線。
12. 根據(jù)權(quán)利要求11所述的系統(tǒng)，其特征是，所述的系統(tǒng)還包括： 層位追蹤裝置，用于根據(jù)所述的疊前地震資料、測(cè)井資料、地質(zhì)成果資料對(duì)目的儲(chǔ)層進(jìn) 行精細(xì)層位追蹤，得到地震層位數(shù)據(jù)。
13.根據(jù)權(quán)利要求12所述的系統(tǒng)，其特征是，所述的表征流體體積模量的AVO近似公式 為： \ - -/
其中，fm= φμ為地下巖石的固體項(xiàng)，Θ為入射角度，Kf為地下巖石的流體體積模量， 乙為地下巖石的固體項(xiàng)，Φ為地下巖石的孔隙度，P為地下巖石的密度，μ為地下巖石的 剪切模量，ΛKf為界面兩側(cè)的流體體積模量的差值，Λfm為界面兩側(cè)固體項(xiàng)的差值，ΛP 為界面兩側(cè)的密度的差值，ΛΦ為界面兩側(cè)的孔隙度的差值，ri=[匕/C]2為干巖石縱 」drv 橫波速度比的平方，匕/C]2為飽和巖石縱橫波速度比的平方。 L -1Sat
14.根據(jù)權(quán)利要求13所述的系統(tǒng)，其特征是，所述的子波提取裝置具體包括： 孔隙度提取模塊，用于從所述的測(cè)井資料中提取出地下巖石的密度以及孔隙度； 固體項(xiàng)確定模塊，用于根據(jù)所述的剪切模量曲線以及所述的孔隙度確定地下巖石的固 體項(xiàng)； 子波提取模塊，用于根據(jù)所述地下巖石的密度、固體項(xiàng)、孔隙度、流體體積模量曲線、表 征流體體積模量的AVO近似公式、疊前CMP道集以及地震層位數(shù)據(jù)提取子波。
15.根據(jù)權(quán)利要求14所述的系統(tǒng)，其特征是，所述的流體體積模量數(shù)據(jù)體確定裝置具 體包括： 反演目標(biāo)函數(shù)建立模塊，用于利用所述的疊前CMP道集、測(cè)井資料、地質(zhì)成果資料以及 子波，基于貝葉斯理論，使用高斯分布函數(shù)為似然函數(shù)、柯西分布為先驗(yàn)約束建立AVO反演 目標(biāo)函數(shù)； 相對(duì)變化量確定模塊，用于利用重加權(quán)最小二乘迭代算法求解所述的AVO反演目標(biāo)函 數(shù)，得到流體體積模量的相對(duì)變化量； 轉(zhuǎn)化模塊，用于利用道積分思想將所述的相對(duì)變化量轉(zhuǎn)化為流體體積模量數(shù)據(jù)體。
16.根據(jù)權(quán)利要求15所述的系統(tǒng)，其特征是，建立的AVO反演目標(biāo)函數(shù)為：

其中，< 為CMP道集資料中噪音的方差，σ·〗為彈性參數(shù)相對(duì)變化量的方差，λ、α為 加權(quán)系數(shù)，m為待反演的彈性參數(shù)相對(duì)變化量矩陣，G為正演算子矩陣，C為積分矩陣，d為CMP道集矩陣，Q為斜對(duì)角加權(quán)矩陣，M為模型約束矩陣，Mtl為M的均值。
【文檔編號(hào)】G01V1/30GK104459777SQ201410725562
【公開(kāi)日】2015年3月25日 申請(qǐng)日期:2014年12月3日 優(yōu)先權(quán)日:2014年12月3日
【發(fā)明者】曹丹平, 李超, 印興耀, 張世鑫, 吳國(guó)忱, 宗兆云 申請(qǐng)人:中國(guó)石油天然氣股份有限公司, 中國(guó)石油大學(xué)(華東)