基于逆算子估計(jì)的AVO反演方法研究

印興耀, 鄧煒, 宗兆云

中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院, 山東青島　266580

基于逆算子估計(jì)的AVO反演方法研究

印興耀, 鄧煒, 宗兆云

中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院, 山東青島266580

摘要傳統(tǒng)反演算法以?xún)?yōu)化算法為主，而基于逆算子估計(jì)的AVO反演算法則利用了直接求逆的思路.算法的關(guān)鍵在于尋找存在逆函數(shù)的子域，進(jìn)而可以在子域內(nèi)直接求逆，這種解決反問(wèn)題的思路不同于一般的優(yōu)化類(lèi)算法所采用的直接搜索解的方式，具有更高的效率.AVO反演利用了振幅隨著偏移距的變化特征，反演的精度受到地震資料質(zhì)量的影響，通過(guò)加入L1范數(shù)約束以及合理的初始模型有助于提高反演的穩(wěn)定性以及準(zhǔn)確度.模型測(cè)算和實(shí)際應(yīng)用表明，基于逆算子估計(jì)的AVO反演方法具有較高的精確程度和可靠性.

關(guān)鍵詞AVO反演； 逆算子估計(jì)； L1范數(shù)； 初始模型

Aki&Richards approximation is utilized to establish the inversion objective function. L1 norm constraint is considered on the basis of reasonable initial model in order to improve effciency and stability during AVO inversion process. Inverse operator estimation algorithm utilizes numerical approximation of objective function in empirically constrained subspaces. The existence of inverse mapping inside these subspaces is supposed, and thus the existence of its numerical approximation is also supposed. This numerical approximation is used for the prediction of the solution. It can be divided into the following four steps: (1) initialization, (2) selection, (3)prediction,(4) solution space correction. Since this method is approximate, the algorithm is arranged into iterative cycles and the solution is gained successively.

Model test in which noise immunity and stability based on L1 norm constraint are considered lays the foundation for the actual data inversion. In the model test, the inversion results based on synthetic seismogram with different SNR noise remain good consistency with the model value. L1 constraint is necessary, a relatively stable result can be obtained in multiple inversion processes and it helps reduce the number of singular values. The seismic data is from China′s land, inversion quality can be evaluated by well data and seismic data.Inversion results of actual seismic data maintain good consistency with logging curves.

Reliability of the solution, noise immunity and convergence speed can be improved by increasing the constraint requirements and considering reasonable initial models. Poor lateral continuity when dealing with actual data can be reduced by optimizing the objective function and the initial model. In general, AVO inversion based on inverse operator estimation is effective and potential.KeywordsAVO inversion； Inverse operator estimation; L1 norm; Initial model

1引言

疊前彈性參數(shù)反演是利用疊前地震信息、測(cè)井資料數(shù)據(jù)以及地質(zhì)信息反演地下的相關(guān)參數(shù)、預(yù)測(cè)儲(chǔ)層以及流體(楊培杰和印興耀，2008;印興耀等，2010,2014)的一種方法.相比于疊后反演方法，利用AVO反演方法可以更為精確地提取反映巖石含油氣性以及巖性的屬性參數(shù)(Goodway et al.,1997;Gray and Andersen,2001; Xingyao Yin and Shixin Zhang，2014；Zong et al.,2015;Yin et al.,2015)，從而提高流體識(shí)別和儲(chǔ)層預(yù)測(cè)的精度.

基于逆算子估計(jì)的AVO反演方法是根據(jù)Zoeppritz方程近似式(Bortfeld，1961；Aki and Richards,1980; Shuey,1985;鄭曉東,1991；Zong et al.,2012，2013a )直接求取三參數(shù)的方法，利用Aki&Richards近似方程建立反演方程，目標(biāo)函數(shù)為正演數(shù)據(jù)與觀測(cè)數(shù)據(jù)殘差的L2范數(shù)，由于反演問(wèn)題的不適定性，較小的噪聲即可對(duì)反演結(jié)果帶來(lái)較大影響，可以通過(guò)正則化改善反演結(jié)果(Tikhonov,1963)，因此本文在L2范數(shù)基礎(chǔ)上加入L1范數(shù)以及初始模型(宗兆云等，2012,2013b;Yin et al.,2013)進(jìn)行約束.求解過(guò)程中逆算子估計(jì)算法的精確程度和穩(wěn)定性至關(guān)重要，通過(guò)試驗(yàn)表明算法可以穩(wěn)定地提取三參數(shù)信息.逆算子估計(jì)算法在實(shí)現(xiàn)過(guò)程中最重要的是找到存在逆函數(shù)的區(qū)域，繼而在這些區(qū)域內(nèi)直接求逆，并且在每一次求解中循環(huán)利用線(xiàn)性回歸、徑向基函數(shù)網(wǎng)絡(luò)和克里金3種不同方法，增強(qiáng)了對(duì)于反問(wèn)題的適應(yīng)性.

本文首先闡述了基于逆算子估計(jì)的AVO反演方法的理論基礎(chǔ)，包括反演目標(biāo)函數(shù)構(gòu)建以及目標(biāo)參數(shù)求解過(guò)程，然后分別利用模型測(cè)算和實(shí)際資料應(yīng)用驗(yàn)證方法的可行性以及有效性.在模型測(cè)試部分除了進(jìn)行抗噪性測(cè)試，還考慮L1范數(shù)約束以及初始模型對(duì)反演的影響.實(shí)際應(yīng)用部分則利用實(shí)際資料直接反演得到三參數(shù)信息，并驗(yàn)證其合理性.

2理論基礎(chǔ)

2.1反演方程

1980年Aki和Richard整理了Frasier和Richards的研究成果，提出了Zoeppritz方程的近似式，即用縱橫波反射系數(shù)以及密度反射系數(shù)表示的線(xiàn)性近似式：

(1)

其中，

(2)

其中，Rp,Rs,Rd分別代表縱波速度反射系數(shù)、橫波速度反射系數(shù)以及密度反射系數(shù).利用正演模型我們可以得到某一個(gè)采樣點(diǎn)不同角度的反射系數(shù)，原則上我們只需要3個(gè)角度的信息就可以達(dá)到目的.建立如下所示的目標(biāo)函數(shù)：

(3)

(4)

(5)

(6)

參數(shù)向量代表待求參數(shù)，數(shù)據(jù)向量代表觀測(cè)數(shù)據(jù)，誤差代表精度，精度表示形式為

算法實(shí)現(xiàn)包括4個(gè)步驟，初始化、選擇、預(yù)測(cè)、校正.初始化階段需要定義求解范圍、模型向量和數(shù)據(jù)向量的維數(shù)；選擇最為關(guān)鍵，目的在于在一個(gè)較大的空間內(nèi)尋找到存在逆函數(shù)諸多的小空間；預(yù)測(cè)則是在這些空間內(nèi)求解，可以通過(guò)很多方法求解，本文中采用線(xiàn)性回歸、徑向基函數(shù)網(wǎng)絡(luò)(閻平凡和張長(zhǎng)水,2005)、克里金3種方法.實(shí)驗(yàn)證明3種方法在算法循環(huán)的過(guò)程中的循環(huán)利用可以有效地提高算法的精度和適應(yīng)性；在校正階段則使求得的解更加合理，即使得目標(biāo)函數(shù)最小.算法在上述步驟的不斷循環(huán)中實(shí)現(xiàn).

(1) 初始化

在這個(gè)階段，對(duì)各個(gè)必要的參數(shù)進(jìn)行初始化，下一次的計(jì)算參數(shù)要被限制在已定義的空間內(nèi)部.

初始模型的集合是按照統(tǒng)一的概率分布在既定的空間內(nèi).所有的模型按照誤差的大小排序，誤差最小的稱(chēng)為Mb:

Mb={pb,db,errb},errb=min(erri),

i=1,2…q.

(7)

在每次的迭代中產(chǎn)生的每個(gè)模型集合都要按照特定的幾何規(guī)則分布：先選擇一個(gè)模型并且作為集合的中心，其他的q-1個(gè)模型隨機(jī)地分布在以其為圓心、R為半徑的球體或者超球體面上.

(2) 選擇

第一次循環(huán)前，隨機(jī)選擇集合中心pc,且R=1.假設(shè){pc,R}已知，將按照下列步驟產(chǎn)生新的預(yù)測(cè)模型：

1) 集合的中心模型{pc,dc,errc}作為第一個(gè)模型；

2) 張量Cm以Choleski分解Cm=LLT；

3) 候選模型pg=pc+R Lg；

g為單位向量.根據(jù)合理的條件判斷該模型是否合理，為了提高算法的效率，在選擇預(yù)測(cè)模型的時(shí)候如果產(chǎn)生了與之前的模型非?？拷那闆r，那么就不應(yīng)該被選擇，這樣就減少了需要正演評(píng)估的計(jì)算量.

(3) 預(yù)測(cè)

被選擇的預(yù)測(cè)模型在這個(gè)階段{pi,di,erri},i=1,2…q才是可以利用的：預(yù)測(cè)模型集合{pi,di,erri}→預(yù)測(cè)算法→候選解p0.

實(shí)驗(yàn)表明，在復(fù)雜多樣的環(huán)境下只使用相同的預(yù)測(cè)算法不利于提高算法的適用性，最好的方式是：即使是在一個(gè)單一的求逆過(guò)程中都要利用多種預(yù)測(cè)算法.現(xiàn)在逆算子估計(jì)是利用3種預(yù)測(cè)算法：線(xiàn)性回歸；利用RBFN(徑向基函數(shù)網(wǎng)絡(luò))預(yù)測(cè)和“克里金”.

通過(guò)選擇，可以保證在新的模型空間內(nèi)逆函數(shù)存在，3種算法的循環(huán)使用對(duì)一個(gè)未知的反問(wèn)題有更強(qiáng)的適應(yīng)性.

(4) 校正

預(yù)測(cè)的質(zhì)量原則上可以通過(guò)對(duì)下列事項(xiàng)的考慮來(lái)控制：1) 組成集合的個(gè)體的個(gè)數(shù)q；2) 集合的大小，即直徑R；3) 集合的位置，即集合的中心pc.

模型個(gè)數(shù)q會(huì)影響計(jì)算量，q值越大，需要選擇的模型越多，計(jì)算量就會(huì)增大，但是對(duì)結(jié)果的影響較小，q值是隨著迭代次數(shù)的增加不斷自動(dòng)發(fā)生改變的，如果要求更高的精度選取較大的q值是有利的，一般選取范圍為模型參數(shù)個(gè)數(shù)的3～5倍.校正的關(guān)鍵在于集合的大小以及集合的位置，具體的校正方法如表1所示，其中nwait表示在循環(huán)過(guò)程中誤差沒(méi)有改善的循環(huán)次數(shù).可以看到目的就在于使得誤差err越來(lái)越小.

表1　校正方法表

3模型反演

下面是一個(gè)實(shí)際的井資料模型，包括橫波速度、縱波速度以及密度反射系數(shù).三參數(shù)模型值如圖1所示.

圖1　井模型參數(shù)圖(a) 縱波速度反射系數(shù)； (b) 橫波速度反射系數(shù)； (c) 密度反射系數(shù).Fig.1　Well model parameters(a) P-wave velocity reflection coefficient; (b) Shear wave velocity reflection coefficient; (c) Density reflection coefficient.

利用模型數(shù)據(jù)以及主頻為40 Hz的雷克子波合成地震記錄，反演得到如圖2所示的結(jié)果，其中藍(lán)色為反演值，紅色為模型值，青綠色為初始值，可以看到反演值與模型值保持很好的一致性，初始值為模型值的100次平滑，初始值與模型值的相關(guān)性不高，但是反演結(jié)果仍然理想.

給合成記錄加信噪比為3的噪聲，合成記錄加噪聲前后對(duì)比如圖3所示((a)為未加噪聲合成記錄，(b)為加噪聲記錄)，對(duì)密度項(xiàng)的反演結(jié)果較其他兩項(xiàng)較差，初始值較不加噪聲反演的時(shí)候與模型值的相關(guān)性提高一點(diǎn)，但反演結(jié)果與模型值也保持了很好的一致性.圖4為加信噪比時(shí)反演結(jié)果與模型值對(duì)比.

給合成記錄加上信噪比為1的噪聲，合成記錄加噪聲前后對(duì)比如圖5所示((a)未加噪聲，(b)為加噪聲)，圖6為反演結(jié)果.從圖4和圖6可以觀察到，隨著噪聲增加，反演對(duì)初始模型的依賴(lài)性相應(yīng)增加.

初始模型是對(duì)模型的不同程度的多次平滑，約束解的范圍是在初始模型基礎(chǔ)上加減一個(gè)該參數(shù)的最大值或者最小值.從模型測(cè)試的結(jié)果可以看出，無(wú)噪聲時(shí)，反演結(jié)果與模型值有很高的吻合度，對(duì)初始模型的依賴(lài)較低，加入不同程度信噪比噪聲后反演結(jié)果與模型值也保持了較好的吻合度，但是對(duì)初始模型的要求也相應(yīng)有所增加，更好的初始模型有助于得到更好的反演結(jié)果.從3個(gè)測(cè)試的結(jié)果可以看到，對(duì)于密度的反演結(jié)果相對(duì)較差，這也是目前線(xiàn)性反演研究工作的一個(gè)難點(diǎn)，并非單純只是本文所使用的算法的原因.從圖2可以看出在資料信噪比較高的情況下對(duì)初始模型的依賴(lài)性不強(qiáng).如果資料的信噪比較低，除了尋找更加合適的初始模型之外，還可以通過(guò)提高精度要求和循環(huán)次數(shù)來(lái)改善反演結(jié)果.總體來(lái)說(shuō)，算法有較高的準(zhǔn)確度和穩(wěn)定性，但是對(duì)于密度項(xiàng)的反演結(jié)果較其他兩項(xiàng)差一點(diǎn).

由于本文采用的是一種新的非優(yōu)化算法，所以對(duì)于加入L1范數(shù)以及初始模型是否對(duì)于反演有幫助是值得商榷的，因此筆者比較了不加范數(shù)和模型與加范數(shù)和模型兩種情況的反演結(jié)果.

考慮到的影響主要分為兩個(gè)方面：穩(wěn)定性和反演中奇異值數(shù)量.

為了驗(yàn)證加約束項(xiàng)對(duì)于結(jié)果穩(wěn)定性的影響，在信噪比為2的情況下，對(duì)上述的反演過(guò)程獨(dú)立進(jìn)行10次，將3參數(shù)信息利用道積分方法轉(zhuǎn)換成密度、縱橫波速度信息，得到10個(gè)反演結(jié)果.對(duì)于每一個(gè)采樣點(diǎn)，計(jì)算這10次反演結(jié)果與模型值的方差以及所有采樣點(diǎn)方差的總和的平均，得到表2，然后分別作歸一化處理得到圖7，從左至右分別是縱波速度方差、橫波速度方差和密度方差，為了得到一個(gè)趨勢(shì)筆者做了一定的平滑處理.值得注意的是，由于做了歸一化，圖7上的大小只是相對(duì)的，不能說(shuō)明兩者的絕對(duì)大小，絕對(duì)大小如表2所示.

圖2　反演結(jié)果與模型值對(duì)比圖紅色曲線(xiàn)代表模型值，藍(lán)色曲線(xiàn)為反演值，青綠色曲線(xiàn)為初始值. (a) 縱波速度反射系數(shù)反演結(jié)果比對(duì)； (b) 橫波速度反射系數(shù)反演結(jié)果比對(duì)； (c) 密度反射系數(shù)反演結(jié)果比對(duì).Fig.2　The comparison between inversion results and the model valueThe model value is colored with red.The inversion value is colored with blue and the initial value is colored with green.(a) Inversion result of P-wave velocity reflection coefficient; (b) Inversion result of shear wave velocity reflection coefficient; (c) Inversion result of density reflection coefficient.

圖3　加噪聲(信噪比為3)前(a)后(b)合成記錄對(duì)比圖Fig.3　The comparison of synthetic seismogram (S/N=3)

圖4　合成記錄加信噪比為3時(shí)噪聲反演結(jié)果與模型值對(duì)比圖紅色曲線(xiàn)代表模型值，藍(lán)色曲線(xiàn)為反演值，青綠色曲線(xiàn)為初始值. (a) 縱波速度反射系數(shù)反演結(jié)果比對(duì)； (b) 橫波速度反射系數(shù)反演結(jié)果比對(duì)； (c) 密度反射系數(shù)反演結(jié)果比對(duì).Fig.4　The comparison between inversion results and the model value with S/N=3The model value is colored with red.The inversion value is colored with blue and the initial value is colored with green. (a) Inversion result of P-wave velocity reflection coefficient; (b) Inversion result of shear wave velocity reflection coefficient; (c) Inversion result of density reflection coefficient.

圖5　加噪聲(信噪比為1)前(a)后(b)合成記錄對(duì)比圖Fig.5　The comparison of synthetic seismogram (S/N=1)

圖6　合成記錄加信噪比為1噪聲反演結(jié)果與模型值對(duì)比圖紅色曲線(xiàn)代表模型值，藍(lán)色曲線(xiàn)為反演值，青綠色曲線(xiàn)為初始值. (a) 縱波速度反射系數(shù)反演結(jié)果比對(duì)； (b) 橫波速度反射系數(shù)反演結(jié)果比對(duì)； (c) 密度反射系數(shù)反演結(jié)果比對(duì).Fig.6　The comparison between inversion results and the model value with S/N=1 The model value is colored with red. The inversion value is colored with blue and the initial value is colored with green. (a) Inversion result of P-wave velocity reflection coefficient; (b) Inversion result of shear wave velocity reflection coefficient; (c) Inversion result of density reflection coefficient.

圖7　加約束項(xiàng)與初始模型前后穩(wěn)定性對(duì)比紅色曲線(xiàn)為不加約束和初始模型，藍(lán)色代表加入約束項(xiàng)以及初始模型. (a) 縱波速度誤差； (b) 橫波速度誤差； (c) 密度誤差.Fig.7　Stability with constraint term and the initial model compared to the condition without it The error with constraint and initial model is colored with blue. The error without them is colored with red. (a) Error of P-velocity; (b) Error of shear wave velocity; (c) Error of density.

圖8　合成10°、20°、30°地震記錄以及反演剖面Fig.8　Synthetic seismograms and P-wave velocity, shear velocity, density inversion results

由表2可以得到，加了約束項(xiàng)和模型后方差總和整體上明顯小于不加約束項(xiàng)，說(shuō)明整體上穩(wěn)定性有所提高.平均誤差為方差，反應(yīng)了與模型值的差值絕對(duì)值，可以看到加L1約束以及模型后誤差明顯減小.再觀察圖7，紅色曲線(xiàn)代表未加約束項(xiàng)與模型，藍(lán)色代表加約束項(xiàng)與模型，對(duì)于每一個(gè)采樣點(diǎn)來(lái)說(shuō)，未加約束項(xiàng)和模型的縱波速度的方差變化幅度為0.3～1，而加了約束項(xiàng)和模型的則基本維持在0.5附近，結(jié)合表2不難看出，加了約束項(xiàng)和模型后對(duì)于每一個(gè)采樣點(diǎn)而言誤差是相對(duì)較小的，并且從整體上來(lái)講誤差維持在一個(gè)更加穩(wěn)定的水平.同樣對(duì)于密度有類(lèi)似的現(xiàn)象.綜上所述，加入合適的約束項(xiàng)和模型后有利于提高結(jié)果的穩(wěn)定性.

表2　方差總和數(shù)值表

其次是所耗時(shí)間，對(duì)10次反演過(guò)程經(jīng)過(guò)統(tǒng)計(jì)，300個(gè)采樣點(diǎn)平均在反演過(guò)程中由于奇異值存在而使得算法無(wú)法在既定的有限循環(huán)次數(shù)內(nèi)達(dá)到求解精度的平均次數(shù)為12，而加入約束和初始模型后平均次數(shù)為3，所以時(shí)間有大幅度減少.

綜上所述，算法可以通過(guò)加入L1范數(shù)與初始模型來(lái)提高穩(wěn)定性以及計(jì)算效率.

對(duì)于二維剖面模型，在不加噪聲的情況下，合成10°、20°、30°的地震記錄(圖8a).選取的初始模型是對(duì)每一道數(shù)據(jù)進(jìn)行100次平滑，反演得到如下縱波速度、橫波速度以及密度剖面如圖8b所示(從左自右).

二維模型反演的目的在于驗(yàn)證橫向連續(xù)性.從二維模型的反演結(jié)果看，算法具有較高的分辨率，橫向連續(xù)性較好，說(shuō)明在得到高準(zhǔn)確度的反射系數(shù)或者子波后逆算子估計(jì)算法具有很高的準(zhǔn)確性，這也為處理實(shí)際資料打下理論基礎(chǔ).

4實(shí)際應(yīng)用

實(shí)際地震數(shù)據(jù)來(lái)自中國(guó)大陸陸地某勘探工區(qū).圖9為3個(gè)部分角度疊加剖面，圖9依次分別是小角度(中心角度10°)、中角度(中心角度20°)、大角度(中心角度30°)的剖面，一共70道，采樣間隔為2 ms.在該工區(qū)提取3個(gè)角度的地震子波，加入L1范數(shù)約束以及適當(dāng)?shù)某跏寄Ｐ陀兄谔岣呓Y(jié)果的準(zhǔn)確性，反演得到如圖10所示的結(jié)果，在圖10中，左半圖是初始模型，從上到下分別為縱波速度反射系數(shù)、橫波速度反射系數(shù)以及密度反射系數(shù).初始模型與地質(zhì)先驗(yàn)信息以及井?dāng)?shù)據(jù)有關(guān)，是對(duì)工區(qū)的一個(gè)粗略的估計(jì)，所選擇解的范圍是在初始模型的基礎(chǔ)上加減一個(gè)工區(qū)該參數(shù)的最大值，這樣的一個(gè)范圍被認(rèn)為是比較合適的，因?yàn)樗x取的范圍至少要能包含合理的解，但是范圍太大又會(huì)增加反演的多解性.另外一方面，如果得到了較好的初始模型，可以縮小解的約束范圍，從而提高反演的效率和質(zhì)量，但本文并非重點(diǎn)尋找更好的初始模型，而關(guān)心的是算法本身.圖10的右半圖是反演結(jié)果，從上而下分別是縱波速度反射系數(shù)、橫波速度反射系數(shù)以及密度反射系數(shù)剖面，由初始模型和反演結(jié)果對(duì)比，初始模型可以確定解的范圍以及初次迭代值，兩者的相關(guān)性并不高，但是反演結(jié)果仍然比較理想.

圖10　縱波速度反射系數(shù)(上)橫波速度反射系數(shù)(中)密度反射系數(shù)(下)剖面Fig.10　P-wave velocity reflection coefficient (top) wave velocity reflection coefficient (mid) density reflection coefficient (below) profiles

圖11　反演值(第70道)(紅色)與軟件反演值(藍(lán)色)對(duì)比Fig.11　Inversion value (70th trace) (red) and software inversion value (blue) contrast

反演結(jié)果整體效果體現(xiàn)出了地震有限頻帶內(nèi)的信息，地震剖面上反射能量的強(qiáng)弱與反演結(jié)果是一致的，圖10中紅色代表高值，縱波速度整體上變化較大，而密度的變化較小，整體的橫向連續(xù)性有待提高.為了驗(yàn)證結(jié)果的正確性，用某反演軟件得到較為理想的反演剖面，隨機(jī)抽出第70道與結(jié)果對(duì)比如圖11所示，圖中從左至右分別是縱波速度反射系數(shù)、橫波速度反射系數(shù)以及密度反射系數(shù)，可以發(fā)現(xiàn)用逆算子估計(jì)算法反演得到的結(jié)果與用軟件反演結(jié)果保持較好的一致性，兩者在局部極值位置的反演結(jié)果有較好的統(tǒng)一.

5結(jié)論

本文采用Aki&Richards線(xiàn)性方程建立AVO反演方程，目標(biāo)函數(shù)加入初始模型以及L1范數(shù)約束，采用非優(yōu)化的逆算子估計(jì)算法進(jìn)行求解，得到縱橫波速度反射系數(shù)以及密度反射系數(shù).模型試算證明算法可以增加對(duì)解的約束條件以及加入合理初始模型，加快收斂速度且提高解的可靠性，并且提高抗噪性.在處理實(shí)際資料時(shí)橫向連續(xù)性和穩(wěn)定性較差，可以通過(guò)優(yōu)化目標(biāo)函數(shù)以及優(yōu)化初始模型的方法減小這個(gè)問(wèn)題帶來(lái)的影響.試驗(yàn)證明基于逆算子估計(jì)的AVO反演方法是有效并且具有較大潛力的.

References

Aki K, Richards P G. 1980. Quantitative Seismology: Theory and Methods. USA: W H Freeman & Co.

Bortfeld R. 1961. Approximations to the reflection and transmission coefficients of plane longitudinal and transverse WAVES.GeophysicalProspecting, 9(4): 485-502. Chen J J, Yin X Y. 2007. Three-parameter AVO waveform inversion based on Bayesian theorem.ChineseJournalofGeophysics(in Chinese), 50(4): 1251-1260.Downton J E, Lines L R. 2001. Constrained three parameter AVO inversion and uncertainty analysis.∥CSEG 2001 expanded abstracts. 251-254.

Downton J E, Lines L R. 2004. Three term AVO waveform inversion.∥SEG Technical Program Expanded Abstracts. 23(1): 215-218.

Goodway W, Chen T, Downton J. 1997. Improved AVO fluid detection and lithology discrimination using Lamé Petrophysical parameters “λρμρλμfluid stack” from P and S inversions.∥ SEG Annual Meeting Expanded Abstracts.

Gray D, Andersen E. 2001. The application of AVO and inversion to the estimation of rock properties.∥ SEG Annual Meeting Expanded Abstracts. 70: 549-552.

Kuzma H A, Rector J W. 2004. Non-linear AVO inversion using support vector machines.∥ 74th SEG Annual Meeting Expanded Abstracts. Ma X Q. 2001. Technical article: Global joint inversion for the estimation of acoustic and shear impedances from AVO derived P- and S-wave reflectivity data.FirstBreak, 19(10): 557-566.

Mogensen S, Link C. 2001. Artificial neural network solutions to AVO inversion problems.∥ 2001 SEG Annual Meeting Expanded Abstracts. Press F. 1968. Earth models obtained by Monte Carlo inversion.JournalofGeophysicalResearch, 73(16): 5223-5234.

Shuey R T. 1985. A simplification of the Zoeppritz equations.Geophysics, 50(4): 609-614. Simmons J L Jr, Backus M M. 1996. Waveform-based AVO inversion and AVO prediction-error.Geophysics, 61(6): 1575-1588. Tikhonov A N .1963. Solution of incorrectly formulated problems and the regularization method.SovietMathematicalDoklady, 4:1035-1038.

Yin X Y, Zhang S X. 2014. Bayesian inversion for effective pore-fluid bulk modulus based on fluid-matrix decoupled amplitude variation with offset approximation.Geophysics,79(5):R221-R232Yan P F, Zhang C S. 2005. Artificial Neural Networks and Evolutionary Computing (in Chinese). Beijing: Tsinghua University Press. Yang P J, Yin X Y. 2008. Non-linear quadratic programming bayesian prestack inversion.ChineseJournalofGeophysics(in Chinese), 51(6): 1876-1882.

Yin X Y, Wu G C, Zhang H Z. 1994 .The application of neural networks in the reservior prediction.JournaloftheUniversityofPetroleum(in Chinese), 18(5):20-26.

Yin X Y, Zhang S X, Zhang F C, Hao Q Y. 2010 .Elastic impedance inversion for reservoir description and fluid identification based on Russell approximation.OilGeophysicalProspecting(in Chinese),45(3): 373- 380.Yin X Y, Zong Z Y, Wu G C. 2013. Improving seismic interpretation: a high-contrast approximation to the reflection coefficient of a plane longitudinal wave.PetroleumScience, 10(4): 466-476. Yin X Y, Cao D P, Wang B L, Zong Z Y. 2014 . Progress on fluid identification method based on pre-stack seismic inversion.OilGeophysicalProspecting(in Chinese), 49(1)：22-34,46.

Yin X Y, Cui W, Zong Z Y, et al. 2014. Petrophysical property inversion of reservoirs based on elastic impedance.ChineseJ.Geophys.(in Chinese), 57(12): 4132-4140.

Yin X Y, Zong Z Y, Wu G C. 2015. Research on seismic fluid identification driven by rock physics.ScienceChinaEarthSciences, 58(2): 159-171. Zheng X D. 1991. The approximate equations of Zoeppritz equation and its applications.OilGeophysicalProspecting(in Chinese), 26(2): 129-144.

Zong Z Y, Yin X Y, Zhang F, et al. 2012. Reflection coefficient equation and pre-stack seismic inversion with Young′s modulus and Poisson ratio.ChineseJ.Geophys. (in Chinese), 55(11): 3786-3794, doi: 10.6038/j.issn.001-5733.2012.11.025. Zong Z Y, Yin X Y, Wu G C. 2012. AVO inversion and poroelasticity with P- and S-wave moduli.Geophysics, 77(6): N17-N24.

Zong Z Y, Yin X Y, Wu G C. 2013a. Elastic impedance parameterization and inversion with Young′s modulus and Poisson′s ratio.Geophysics, 78(6): N35-N42. Zong Z Y, Yin X Y, Wu G C. 2013b. Multi-parameter nonlinear inversion with exact reflection coefficient equation.JournalofAppliedGeophysics, 98: 21-32.

Zong Z Y, Yin X Y, Wu G C, et al. 2015. Elastic inverse scattering for fluid variation with time-lapse seismic data.Geophysics, 80(2): WA61-WA67.

附中文參考文獻(xiàn)

陳建江, 印興耀. 2007. 基于貝葉斯理論的AVO三參數(shù)波形反演. 地球物理學(xué)報(bào), 50(4): 1251-1260.

閻平凡, 張長(zhǎng)水. 2005. 人工神經(jīng)網(wǎng)絡(luò)與模擬進(jìn)化計(jì)算. 北京: 清華大學(xué)出版.

楊培杰, 印興耀. 2008. 非線(xiàn)性二次規(guī)劃貝葉斯疊前反演. 地球物理學(xué)報(bào), 51(6): 1876-1882.

印興耀, 吳國(guó)忱, 張洪宙. 1994 .神經(jīng)網(wǎng)絡(luò)在儲(chǔ)層橫向預(yù)測(cè)中的應(yīng)用.石油大學(xué)學(xué)報(bào):自然科學(xué)版, 18( 5): 20 -26.

印興耀, 張世鑫, 張繁昌, 郝前勇. 2010 .利用基于Russell 近似的彈性波阻抗反演進(jìn)行儲(chǔ)層描述和流體識(shí)別. 石油地球物理勘探, 45(3): 373-380.

印興耀，曹丹平，王保麗，宗兆云. 2014 .基于疊前地震反演的流體識(shí)別方法研究進(jìn)展.石油地球物理勘探，49(1): 22-34,46.

印興耀, 崔維, 宗兆云, 劉曉晶. 2014.基于彈性阻抗的儲(chǔ)層物性參數(shù)預(yù)測(cè)方法研究. 地球物理學(xué)報(bào), 57(12): 4132-4140.

鄭曉東. 1991. Zoeppritz方程的近似及其應(yīng)用. 石油地球物理勘探, 26(2): 129-144.

宗兆云, 印興耀, 張峰等. 2012. 楊氏模量和泊松比反射系數(shù)近似方程及疊前地震反演. 地球物理學(xué)報(bào), 55(11): 3786-3794.

(本文編輯胡素芳)

AVO inversion based on inverse operator estimation

YIN Xing-Yao, DENG Wei, ZONG Zhao-Yun

Geo-ScienceandTechnologyFaculty,ChinaUniversityofPetroleum,QingdaoShandong266580,China

AbstractAmplitude variation with amplitude or angle (AVO/AVA) inversion has been widely utilized in exploration geophysics to estimate the formation of elastic parameters underground. The traditional inversion algorithms are mainly optimization, while the AVO inversion based on inverse operator estimation is to inverse directly. The key is to find the subspaces which exist inverse mapping instead of searching for the solution directly as optimization algorithms do.

基金項(xiàng)目國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2013CB228604)，國(guó)家油氣重大專(zhuān)項(xiàng)(2011ZX05030-004-002)，中國(guó)博士后科學(xué)基金(2014M550379),山東省博士后基金(201401018)和山東省優(yōu)秀中青年科學(xué)科研獎(jiǎng)勵(lì)基金(2014BSE28009)聯(lián)合資助.

作者簡(jiǎn)介印興耀，男，中國(guó)石油大學(xué)(華東)教授.主要從事地球物理理論與方法方面的研究. E-mail： xyyin@upc.edu.cn

doi:10.6038/cjg20160426 中圖分類(lèi)號(hào)P631

收稿日期2015-03-20，2015-10-08收修定稿

印興耀, 鄧煒, 宗兆云. 2016. 基于逆算子估計(jì)的AVO反演方法研究.地球物理學(xué)報(bào)，59(4):1457-1468,doi:10.6038/cjg20160426.

Yin X Y, Deng W, Zong Z Y. 2016. AVO inversion based on inverse operator estimation.ChineseJ.Geophys. (in Chinese),59(4):1457-1468,doi:10.6038/cjg20160426.