專利名稱::利用地震與井控制數(shù)據(jù)的統(tǒng)計比較優(yōu)化幅值反演的制作方法
技術(shù)領(lǐng)域:
:本發(fā)明總體上涉及地下構(gòu)造的地震探測,更具體地,涉及使用井控制數(shù)據(jù)來改進地震處理和幅值反演的質(zhì)量。
背景技術(shù):
:在最近25年期間,石油和天然氣工業(yè)尋求從地震數(shù)據(jù)得到更多的地下特性信息一尤其是關(guān)于地下孔流、孔隙率、巖相、壓力和幾何結(jié)構(gòu)的詳細信息。這是部分地通過使用地震數(shù)據(jù)與井記錄信息的組合來提高從地震數(shù)據(jù)估計出的地下特性的精度而完成的。這些努力通過改進的勘探成功率和開發(fā)井"最佳鉆井位置(sweetspot)"的選擇而得到回報。這個成功在很大程度上是基于相對聲阻(在交界面處密度與壓縮速度改變的乘積)的復(fù)原。圖1顯示對其尋找詳細地下信息的多個地下構(gòu)造10。地下構(gòu)造10可包括一個或多個儲層或生產(chǎn)區(qū)域。地震探查是通過使用地震聲源12和接收機14在這個感興趣的地下區(qū)域進行的。接收機14記錄以波前從源12行進到地下構(gòu)造10的聲波,在地下構(gòu)造10中聲波的一部分被反射回接收機14,并被記錄為地震軌跡。記錄車16常常被用來存儲所記錄的數(shù)據(jù)。另外,優(yōu)選的是,鉆井20穿過地下構(gòu)造IO,記錄工具在鉆井20中工作,得到關(guān)于地下構(gòu)造10的彈性特性的地下特性信息。圖2是包括單個源12和以2D幾何關(guān)系布置的多個接收機14的單個爆破(shot)剖面的示意圖。該示意圖顯示射線路徑幾何關(guān)系(垂直于傳播波前),主要(直接路徑)地震反射能量行進,從源位置到達反射點22,然后回到接收機位置。顯示了作為對于最大源-接收機距離的反射角的角度e〃爆破剖面是對于形成共同中點(CMP)集合所需的基局部震探查組成部分,如圖3所示。在這個圖中,顯示了對于最接近的源-接收機對的反射角&。CMP軌跡集合是一系列爆破剖面軌跡的集合,它們具有不同的爆破到接收機的距離和共同的源或中點位置。如果地質(zhì)層沒有下傾的話,這個位置對應(yīng)于共同的地下反射點24。因為地質(zhì)層常常具有復(fù)雜的結(jié)構(gòu),CMP集合軌跡通過使用各種不同的成像和噪聲抑制技術(shù)被處理成共同反射點(CPR)軌跡集合。這些集合是被輸入到幅值反演算法的軌跡的典型集合。這些集合包含關(guān)于主反射強度或幅值如何在公共地質(zhì)交界面處根據(jù)射線路徑到達角而改變的信息。對于得到所觀察的反射幅值的地下交界特性的反演需要許多到達角軌跡。圖4描繪了幅值對偏移(AVO)主反射信號的巻積模型,圖中顯示軌跡的主反射信號幅值可能與在地質(zhì)構(gòu)造10之間的交界面26有怎樣的關(guān)系。來自井記錄的數(shù)據(jù)和從鉆井20取得的巖芯顯示,地質(zhì)構(gòu)造的物理特性典型地隨不同的構(gòu)造而不同。通過這些構(gòu)造傳播的聲波對于構(gòu)造的彈性特性,包括壓縮速度Vp、剪切速度Vs、和密度特性p最敏感,而對它們的吸收特性的敏感程度較低。在兩種構(gòu)造之間的交界面處彈性特性的差別決定了反射的和透射的波前的幅值相對于向下的或到達的波前的幅值的大小。如圖4所示,反射信號幅值的大小可以用三個彈性反射率或通過交界面的彈性特性的均值歸一化的交界面彈性特性的差值的加權(quán)和來近似。對于每個反射率的加權(quán)函數(shù)是如公式(l)描述的反射角的三角函數(shù)。A(0)sRo+G*sin2e+RP*sin2Gtan2G(1)其中A(0"以角度e反射的幅值9[平均交界面角度-(eren+erefr)/2Ro[p-阻抗反射率]-Avp/2vp+Ap/2pG[梯度反射率l-Avp/2vp-K(Ap/2p+Avs/vs)K=(2vs/vp)2;以及Rp[p-反射率-Avp/2vp7由許多構(gòu)造交界面構(gòu)成的地質(zhì)部分對于傳播的地震波前起的作用將像由對于聲波在交界面之間行進所花費的時間間隔開的一系列反射率。主反射信號響應(yīng)是通過用傳播的小波的拷貝或波前擾動替代每個反射率,并被縮放反射率的幅值和正負號而被計算。表面記錄的主響應(yīng)是所有的時間移位的縮放的小波或小波的巻積與交界面反射率級數(shù)之和。在圖4的底部的兩個箭頭表明AVO幅值反演的目標是把AVO軌跡變換成它們的分量反射率,而AVO幅值分析的目標是指生成反演的反射率的、在交界面處的地質(zhì)特性的改變。在鉆到地質(zhì)構(gòu)造10中的鉆井20中記錄的井記錄數(shù)據(jù)可被用來估計地質(zhì)構(gòu)造10的彈性和吸收特性并隨后估計主反射信號幅值,作為全部記錄的地震探查響應(yīng)的一部分。井記錄數(shù)據(jù)可被用來統(tǒng)計地表征預(yù)期的反射響應(yīng)并提供處理質(zhì)量度量,如下面說明的。對于與角度相關(guān)的幅值反演的基礎(chǔ)工作是在1950年代,在那時Bortfeld,R.,1961,Approximationstothereflectionandtransmissioncoefficientsofplanelongitudinalandtransversewaves,Geophys,Prosp.,v.9,p.485-502,描述了對于Zoeppritz反射系數(shù)公式的線性化表達式。Lindseth,R.O.,1979,Synthicsoniclogs—aprocessforstratigraphicinterpretation:Geophysics,44,p3-26,實現(xiàn)了對于阻抗的軌跡幅值的反演。此后不久,由Ostrander,W.J"1984,Plane-wavereflectionforgassandsatnonnormalanglesofincidence:Geophysies,49,1637-1648,實現(xiàn)了定量AVO分析。隨后,利用定量和定性幅值對偏移(AVO)反演來估計地下地質(zhì)特性。幅值反演包括從通過增加反射角或爆破到接收機偏移而排序的公共地下反射點軌跡集合,預(yù)測一個或多個分量反射率(與已知的小波進行巻積)造成AVO反射響應(yīng)的過程。圖5顯示在兩個媒介30和32之間的單個交界面、對于具有單位幅值到達的到達脈沖小波的反射角Gren和折射角0refr、以及和值與差值彈性特性,它們在確定反射幅值脈沖的強度時是重要的。還顯示了媒介30具有特性vpl,vsl和Pl,媒介30具有特性vp2,vs2和p2。表l提供了用于確定在交界面處反射幅值的因,。表l:用于確定在交界面處反射幅值的因素交界面反射率-Af(vs'Vp,P),g(vs'Vp,p)_vp"媒介質(zhì)壓縮波速度__vs-媒介質(zhì)剪切波速度__p;樣介質(zhì)密度_Avs-vS2-vsl;vs=(vS2+vsl)/2Avp=Vp2-vp1;vp=(Vp2+vP1)/2Ap=p2_Pi;p=(p2+pi)/2表2提供了通常在AVO反演中使用的、對于反射率的定義。應(yīng)當(dāng)指出,Af(…)和g(,.,)可以是媒介的差分和平均特性的不同的函數(shù)。表2:在AVO幅值反演中使用的反射率的定義<table>tableseeoriginaldocumentpage9</column></row><table>圖6是使用完全Zoeppritz平面波反射系數(shù)公式(l)的、對于單個交界面的典型的AVO響應(yīng)的圖。圖中顯示,對于典型的地質(zhì)交界面和在小的反射角下,反射系數(shù)近似為常數(shù)。在接近30。的角度下反射系數(shù)典型地消失幾個百分點,以及在接近臨界角的角度下,反射系數(shù)接近于1。對于精確的Zoeppritz公式有許多近似。這些近似是從Zoeppritz公式的復(fù)雜構(gòu)和從地震數(shù)據(jù)獲取與處理方法可得到的與角度有關(guān)的幅值信息的類型導(dǎo)出的。Aki-Richardszi^式,這是平面波Zo印pritz公式的線性化版本,產(chǎn)生對于法線到次臨界入射角和對于小的對比的交界面的主反射幅值的完美的近似。它也是作為許多幅值反演算法的基礎(chǔ)的公式。根據(jù)圖5中顯示的幾何關(guān)系和彈性媒介參數(shù),公式可被寫為<formula>formulaseeoriginaldocumentpage10</formula>其中A(e)-以角度e反射的幅值e[平均交界面角度卜(eren+erefr)/2Ro[p-阻抗反射率]-AVp/2Vp+Ap/2pG[梯度反射率^AVp/2Vp-K(Ap/2p+AVS/VS);K=(2VS/VP)2Rp[p-反射率]-AVp/2Vp以及Vp,Vs,p,AVp,AVs和Ap是在表1中定義的。公式(l)表明,作為地下反射角的函數(shù)的AVO幅值是反射角的三角函數(shù)與三個反射率Rn,G和Rp的乘積之和。這些反射率又是在地下交界面處的壓縮波速度、剪切波速度和密度特性的差值與均值的函數(shù)。梯度項G是牽涉到三個另外的反射率的和與乘積的特別復(fù)雜的反射率。另外,在公式(l)中的角度是在交界面處入射角和折射角的平均值的函數(shù)。通常,公式(l)可以通過使用其它三角函數(shù)和具有以下的顯示的形式的反射率來按公式計算<formula>formulaseeoriginaldocumentpage10</formula>測得的幅值可以通過對于三個或更多個源到接收機偏移進行反射幅值測量和從與偏移有關(guān)的行進時間精確地測量過度的速度場以便估計反射角而被反演,產(chǎn)生三個交界面反射率。.反射率模型<formula>formulaseeoriginaldocumentpage10</formula>,其中-.A產(chǎn)在①i處測得的幅值Ro-阻抗的改變(AVp,p)R產(chǎn)梯度系數(shù)R尸Vp的改變Op反射角-反射器解—ANGLE_-iGEOMETRYSA,SIN20,R2(3X3)SAiSENP0iTANZ0i表3表3顯示用于幅值反演的非限制的、最小平方方法(L2模),典型地被用來求解用于其幅值和角度通過數(shù)據(jù)處理方法進行估計的反射的反射率。L2非限制的方法使得以下形式的誤差函數(shù)最小化E2=S(Ai-(Ro+R^sin29i+R2*sin2eitan26j))2(3)通過求解朋2/朋,=0況2/朋2=0其中Ai和G產(chǎn)事件幅值和反射角。定量幅值反演的目標是反演的屬性的信噪比在小的空間距離上是可與堆疊部分響應(yīng)的信噪比進行比較的。因為在定量地震幅值反演技術(shù)中預(yù)測的那些實踐是熟知的,很難單獨用公式(l)定量地得到有用的幅值反演結(jié)果。主要誤差源典型地包括由各種噪聲和在大反射角時入射角的不精確估計所影響的幅值。在具有有限的反射角范圍和不正確的信號失真效果的有噪聲的數(shù)據(jù)中,未限制的幅值反演將生成反演的反射率,與堆疊部分或p-阻抗部分相比,具有非常低的S/N比率。為了提高幅值反演結(jié)果的S/N比率,公式(3)通過使用其它誤差模(例如,Ll模)和補充地震幅值信息的各種不同的約束條件來表示。該約束條件可以是硬約束條件一改變反射率公式中的形式或變量的那些,或是軟約束條件一作為要被最小化的誤差函數(shù)所包括的附加信息。在1970年代使用的AVO公式,用于較短的偏移獲取幾何結(jié)構(gòu)和非幅值保留處理序列,用以下形式的公式對具有偏移的幅值進行建模A(X)sRo+G*X2(4)其中X-爆破到接收機偏移。這個公式使用相對于公式(l)的"硬"約束條件,暗示sin2etan20項是可忽略的,并且在單個標量常數(shù)內(nèi)恒定和以小的反射角,平方的偏移是對于sii^e項的良好的代理。去除了NMO的、非轉(zhuǎn)移的CMP集合是對于處理后的幅值的Ro和G的最小平方擬合的通常的輸入。R0和G項的時間平均能量或包絡(luò)差被用作為直接碳氫化合物指示,沒有采用其他約束條件數(shù)據(jù)。以后從測得的速度場通過使用導(dǎo)致下式的直射線近似來計算反射角A(est)三Ro+G*sin2est(5)其中0s產(chǎn)對入射角的直射線近似-tan"(X/(toV邵s))。公式(5)暗示公式(1)的sin2etan2e項是可忽略的,以及直射線角度是對于地下反射角的良好近似一當(dāng)幾乎沒有速度加速時的良好假設(shè)。這個公式的使用也發(fā)起用于通過交叉繪制出與它的梯度的事件交截而檢測異常事件行為的分析方法。但使用最小硬約束條件假設(shè)需要在跨越法線到臨界角的角度范圍上的幅值信息和在生成在本專利申請中使用的反演例子時使用的公式(l)的形式是A(G)sRo-2RSH*sin29+RP*tan2e(6)其中.Rsh=A(pVs2)/pVp2為了補償由于噪聲和失真造成的AVO數(shù)據(jù)中的不適應(yīng)性,關(guān)于在巖石壓縮速度、剪切速度和密度與Vp/Vs項的形式之間的關(guān)系的假設(shè)通過改變在反射公式中的變量而導(dǎo)致公式(l)的反映這些"硬"約束條件的版本。這些硬約束條件包括經(jīng)驗關(guān)系,如A+BVp-Vs,暗示△Vs/Vs=(BVp/Vs)AVP/VP(7)其中A和B是常數(shù)。Gardner法則決定速度與密度之間的關(guān)系P=CVpK(8)其中C和K是常數(shù),被用來表示Ap/p=KAVp/VP(9)當(dāng)參數(shù)KK2Vs/Vp被設(shè)置為常數(shù)時,公式(l)的其它表述可以借助于p-阻抗和剪切阻抗反射率結(jié)果表示。以數(shù)據(jù)加權(quán)的形式的軟約束條件和衰減約束條件也被用于幅值反演,以便考慮到在地震幅值數(shù)據(jù)中的噪聲和失真。加到公式(3)的軟約束條件導(dǎo)致以下形式的誤差函數(shù)E2=SWi(Ai-(Ro+R一sin26i+R2*sin2eitan2G,))2+WcoRo2+WdR+WC2R22+Wff(Ro,Ri,R2)2+...(10)其中Wi-幅值A(chǔ)i的信噪比估計;Wco,d,c2,產(chǎn)加到反射率的衰減因子;以及f(R。,RbRj=用于反射率的和值或差值的衰減因子。與利用公式(10)相關(guān)聯(lián)的缺點是建立用于選擇權(quán)重和衰減因子的標準。理想地,如果幅值數(shù)據(jù)沒有噪聲或失真的話,衰減權(quán)重應(yīng)當(dāng)是零,而幅值權(quán)重因子等于1。因為現(xiàn)場記錄的爆破剖面常常是比信號多的噪聲,以及因為爆破剖面信號受到各種失真影響,下面描述的數(shù)據(jù)處理的類型是在幅值反演之前非常想要的。地震數(shù)據(jù)的典型的現(xiàn)場記錄爆破剖面包含被獲取系統(tǒng)顯著歪曲的主反射信號,穿過地表從爆破點到接收機的傳輸效果,和爆破產(chǎn)生的、環(huán)境和獲取噪聲的影響。圖7是在通過和不通過初始噪聲抑制步驟顯示的正常移出(NMO)校正(對于聲波傳播速度的校正)后實際的爆破剖面圖,在初始噪聲抑制步驟中主反射應(yīng)當(dāng)呈現(xiàn)為并行平坦放置的相干事件。如圖7所示,因為各種附加噪聲(在這種情形下,地面滾動和倍數(shù))和將產(chǎn)生非常嘈雜的反演的反射率的信號失真影響,很難檢測真實數(shù)據(jù)主信號反射。圖8是描繪會使得傳播的地震小波和它的地下反射行為失真的某些重要的因素。影響因素包括地震源強度、源和接收機方向性、耦合、陣列特性和近表面分層、波前擴散或發(fā)散損失、非彈性吸收、層間多散射、局部地質(zhì)傾斜、反射交界面的彎曲、橫向地球非均勻性、傳輸系數(shù)損失、記錄系統(tǒng)濾波和所要確定的地下反射系數(shù)。如下面在表4中顯示的,地震數(shù)據(jù)處理序列的目的是去除來自地震數(shù)據(jù)的倍數(shù)和其他噪聲,以及通過應(yīng)用各種校正算法到地震探查數(shù)據(jù)而補償獲取和地球濾波對主反射的影響。某些步驟本身可包含子序列的步驟,而其它步驟可以在具有不同的參數(shù)的序列中重復(fù)一次以上,因為數(shù)據(jù)的信號和噪聲結(jié)構(gòu)變?yōu)楦黠@。表4:處理/反演序列分量通用處理步驟后遷移和預(yù)反演幅值調(diào)節(jié)步驟獲取幾何結(jié)構(gòu)分配預(yù)堆成像crp集合初始速度分析CMP噪聲抑制六數(shù)據(jù)規(guī)則化小波頻語均衡*爆破/cmp噪聲抑制*剩余事件對齊*信號校正地球/記錄空間幅值平衡小波均衡/去巻積*偏移/角度幅值平衡最后速度分析角度/偏移選擇靜音*CMP噪聲抑制a幅值反演/不確定性*預(yù)堆疊遷移/速度改善*=步驟可以使用井記錄統(tǒng)計幅值處理序列(在幅值反演之前的處理)的目的是對數(shù)據(jù)作出各種校正,同時還保留在空間和反射角或偏移方面主反射的相對幅值表現(xiàn)。這個特定的目的也可以定量地闡述為恢復(fù)地球平面波的單個標量常數(shù)內(nèi)的幅值,頻帶限制的地下反射率,以使得在每個交界面處的響應(yīng)對應(yīng)于具有已知的相對幅值、相位和定時的小波與局部交界面反射系數(shù)的巻積。幅值處理序列可以具有20個或更多的步驟,包括在表4中列出的那些步驟。表4的左面的列包括在產(chǎn)生遷移的CRP集合的通用處理序列中典型的步驟,而右面的列列出在后遷移、預(yù)反演幅值調(diào)節(jié)序列中的步驟。重要的通用處理步驟包括幾何結(jié)構(gòu)分配、速度分析、數(shù)據(jù)規(guī)則化、在不同的域中噪聲抑制的通道、信號失真校正、小波均衡校正、最后速度分析以及遷移或成像步驟。表4的右面列是預(yù)反演調(diào)節(jié)序列,它包括用于剩余噪聲抑制、剩余事件時序(速度)校正、小波頻i普均衡、剩余幅值平衡和部分CRP集合的數(shù)據(jù)角度/偏移靜音的步驟。在幅值反演之前的預(yù)反演調(diào)節(jié)序列是重要的,因為通用處理序列(表4的左面一列)常常不適當(dāng)?shù)匮a償各種獲取和地球傳輸影響。而且,對于表4的兩個列的處理步驟,用戶(典型地,處理地質(zhì)學(xué)家)需要選擇用于實施特定的噪聲或信號失真校正的每個處理步驟的參數(shù)和參數(shù)值。所選擇的數(shù)值可以通過該步驟決定性地改變數(shù)據(jù)的輸出幅值特性。例如,在去巻積步驟中,算子類型的選擇(產(chǎn)生間隙或加白)、時間變化的程度、自相關(guān)設(shè)計門的尺寸、空間平均的程度和去巻積算子的長度可以大大地改變?nèi)喎e的數(shù)據(jù)的特性。來自所有的處理步驟的所選擇的參數(shù)和參數(shù)值的組合將確定最后的幅值反演反射率的質(zhì)量。最后,來自幅值處理序列的輸出數(shù)據(jù)將針對幅值反演而被優(yōu)化。在處理序列中的每個步驟和每個參數(shù)可以通過在整個數(shù)據(jù)集中對于步驟和參數(shù)值進行全局搜索而被優(yōu)化。但對于20步驟序列,每個步驟具有3個參數(shù),每個參數(shù)有10個可能的數(shù)值,數(shù)據(jù)集必須被處理106°次來搜索參數(shù)值空間,這是完全不實際的和開銷不能容忍的建議。此外,即使在有可能生成10"個輸出數(shù)據(jù)集的情況下,建立用于選擇優(yōu)化的數(shù)據(jù)集的標準的問題仍舊存在。因此,由于在序列中的步驟的15數(shù)目、在序列中的步驟的次序、在每個步驟中關(guān)鍵參數(shù)的數(shù)目、和對于執(zhí)行各個步驟所需要的計算資源量,常規(guī)的實踐是測試和選擇在數(shù)據(jù)子集中的參數(shù)值,一次一個步驟。有經(jīng)驗的解譯者或處理者然后通過檢查各種質(zhì)量控制(QC)顯示而判斷對于該步驟或步驟集合哪些參數(shù)值生成最高質(zhì)量的輸出數(shù)據(jù)。在堆疊部分(在CRP集合中軌跡的和)中空間相干的程度、在CMP/CRP集合中的相干性或經(jīng)過處理的輸出與從井控制得到的理論響應(yīng)的相似性常常對于在CRP集合中的軌跡之間的幅值的小的橫向變化非常敏感,因為來自一個步驟的輸出影響參數(shù)值在后面的步驟中的效果,并且由于典型的QC產(chǎn)品不使用來自幅值反演的輸出,針對單個處理步驟的參數(shù)選擇判斷可能相對于針對幅值反演的優(yōu)化數(shù)據(jù)的順序目標而表現(xiàn)為次優(yōu)。圖9顯示在執(zhí)行包括在表4中所包括的那些步驟的各種處理步驟后的單個CMP和CRP集合。應(yīng)當(dāng)指出,S/N比率和信號(平坦的事件)的特性的顯著改變作為數(shù)據(jù)順序進到所指示的處理步驟。圖10比較初始NMO校正后的爆破剖面與最終的預(yù)反演CRP集合,顯示S/N比率、幅值分布和事件的相位特性已顯著改變。圖10動態(tài)地指示處理序列對于將被輸入到幅值反演的主反射信號的特性有多大的影響。由于在技術(shù)上經(jīng)歷的那些是已知的,圖10的左面的地震數(shù)據(jù)將生成帶有噪聲的幅值屬性。但即使圖10的右面顯示的數(shù)據(jù)也可能生成帶有噪聲的反演的幅值屬性。這是因為數(shù)據(jù)處理可能沒有足夠地補償各種剩余獲取和地球傳輸影響。事實上,當(dāng)未校正時,在由這樣的剩余影響引起的事件的相對時序(速度)、幅值和相位的小的變化可導(dǎo)致反演的屬性(例如,梯度)的RMS水平的200%到400%的誤差。即使在這樣的剩余影響是相當(dāng)大時,它們可能很難用典型的依賴于在CDP或CRP域中事件連續(xù)性的視覺估計的處理QC方法來檢測。在數(shù)據(jù)處理后,用戶必須選擇各種參數(shù),用于幅值反演步驟本身(例如,衰減參數(shù)和數(shù)據(jù)權(quán)重)。反演約束條件參數(shù)的選擇常常強烈地影響從帶有偏移的幅值的差值得出的那些反射率。用戶典型地使用兩個方法對來自幅值反演的輸出進行質(zhì)量控制檢查反演的反射率的公共深度點(CDP)空間相干,或分析反演的反射率與從現(xiàn)有的井控制計算出的反射率相比是如何接近。這些方法的任一項都不指示所選擇的反演參數(shù)和約束條件是否是最佳的,而且這些標準也不保證在遠離井控制的區(qū)域中的參數(shù)選擇是最佳的。實際上,反演的參數(shù)常常受到另外的處理、統(tǒng)計分析和"校準",以便提高它們的信噪比。這樣的后反演分析不引入附加的不確定性和誤差常常不能補償不適當(dāng)?shù)奶幚怼?傊?,如上所述,在傳統(tǒng)的地震處理和幅值反演中有幾個缺點。第一個缺點是幅值處理和反演序列因為差的參數(shù)選擇可能是欠最佳的,這導(dǎo)致處在和遠離井控制處反演的反射率的次最佳信噪比。第二個缺點是在井控制和遠離井控制處可應(yīng)用的對于感興趣的區(qū)域的處理和反演參數(shù)的選擇時缺乏用于QC(質(zhì)量控制)的標準。在當(dāng)前的幅值反演方法中的第三個缺點是對于補充被包含在處理后的地震數(shù)據(jù)中的信息的各種約束條件的使用。優(yōu)選地,最佳處理的CRP集合的相對幅值表現(xiàn)正比于與角度有關(guān)的交界面反射系數(shù),所以需要最小硬或軟約束條件信息,以便產(chǎn)生高信噪比(S/N)反演的屬性。當(dāng)處理后的地震數(shù)據(jù)具有差的質(zhì)量時,反演算法本身必須受到強烈的約束限制(經(jīng)由衰減和加權(quán)參數(shù)),以產(chǎn)生反演的屬性的現(xiàn)實數(shù)值。隨著約束限制的權(quán)重增加,反演的屬性變?yōu)楦蕾囉诩s束條件,而不太依賴于處理后的地震數(shù)據(jù)。這帶來問題,即反演結(jié)果的精度遠離其中不能應(yīng)用約束條件的位置。本發(fā)明解決這些缺點。
發(fā)明內(nèi)容參照圖21,本發(fā)明的一方面包括用于得到增強的地震數(shù)據(jù)和優(yōu)化的反演的反射率的方法。該方法包括根據(jù)對應(yīng)于感興趣的地下區(qū)域的頻帶限制的井反射率計算多個統(tǒng)計井特性(SCSwe,i)的步驟。通過把多個不同的地震處理序列應(yīng)用到對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù)集的至少一部分而計算多個處理的地震數(shù)據(jù)集。還通過反演多個經(jīng)過處理的地震數(shù)據(jù)集而計算多個反演的反射率組。多個統(tǒng)計地震特性(SCs^,)是根據(jù)多個反演的反射率組計算的。多個統(tǒng)計差值測量結(jié)果(SDM)是在SCSweU與多個SCSsei訓(xùn)ic之間進行計算的。根據(jù)計算出的SDM選擇最佳地震處理序列。通過利用最佳的地震處理序列來處理對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù),得到增強的地震數(shù)據(jù)。然后根據(jù)計算出的SDM選擇最佳地震反演算法。然后使用最佳地震反演算法來反演地震數(shù)據(jù),以便通過使用最佳地震反演算法來反演增強的地震數(shù)據(jù)而產(chǎn)生優(yōu)化的反演的反射率組。在所述方法中,不同的地震處理序列可包括參數(shù),至少參數(shù)的數(shù)值在不同的處理序列之間是不同的。在所述方法中,不同的地震處理序列的不同之處至少在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的次序是不同的。在所述方法中,不同的地震處理序列的不同之處至少也可以在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的數(shù)目是不同的。在所述方法中,在優(yōu)化方案中可以利用統(tǒng)計差值測量結(jié)果(SDM)來預(yù)測最佳地震處理序列。在所述方法中,可以不使用從記錄感興趣區(qū)域內(nèi)的鉆井得到的井記錄數(shù)據(jù)而計算頻帶限制的井反射率。在所述方法中,計算多個統(tǒng)計井特性(SCsweu)可包括以下步驟(i)獲取井記錄數(shù)據(jù);(ii)用井記錄處理參數(shù)選擇井記錄處理算法;和(iii)通過使用井記錄處理算法把井記錄數(shù)據(jù)轉(zhuǎn)換成頻帶限制的井反射率。井記錄數(shù)據(jù)可以從在感興趣區(qū)域內(nèi)的鉆井、類似地質(zhì)設(shè)置、多個類似的井的平均值和統(tǒng)計模型之一被獲取。在所述方法中,選擇最佳地震處理序列的步驟也可包括確定SDM中的任一個是否在預(yù)定的標準內(nèi)的步驟。在所述方法中,SCsweu可以根據(jù)剪切的、壓縮的和密度反射率的RMS水平以及在剪切的、壓縮的和密度反射率之間的相關(guān)系數(shù)進行計算。在所述方法中,SDM響應(yīng)于井反射率的S/(S+N)對地震數(shù)據(jù)反射率的S/(S+N)的比率。統(tǒng)計的井特性也可以在Zoe卯ritz反射率公式中用來,以計算AVO集合中地震軌跡的預(yù)期的RMS。在所述方法中,SDM可以反比于井反射率的預(yù)期的S/(S+N)比率與地震數(shù)據(jù)反射率的S/(S+N)比率之間的差值。參照圖22,本發(fā)明的再一方面是用于得到增強的地震數(shù)據(jù)的方法,增強的地震數(shù)據(jù)可被用于優(yōu)化的反演反射率或其它目的。該方法包括根據(jù)對應(yīng)于感興趣的地下區(qū)域的頻帶限制的井反射率,計算統(tǒng)計井特性(SCSweU)組的步驟。通過把多個不同的地震處理序列應(yīng)用到對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù)集的至少一部分而計算多個經(jīng)過處理的地震數(shù)據(jù)集。通過利用地震反演算法用固定的參數(shù)組來反演經(jīng)過處理的地震數(shù)據(jù)集而計算多個反演的反射率組。統(tǒng)計地震特性(SCs^,)組是根據(jù)反演的反射率組計算的。統(tǒng)計差值測量結(jié)果(SDM)是在SCS組與多個SCs^mic組之間進行計算的。根據(jù)計算出的SDM選擇最佳地震處理序列。地震數(shù)據(jù)通過使用最佳地震處理序列和參數(shù)而被處理,以產(chǎn)生增強的地震數(shù)據(jù)集。該方法還可包括確定SDM是否處在預(yù)定的標準內(nèi),以及重復(fù)利用地震處理序列的固定的參數(shù)組的計算和選擇步驟,直至SDM差值在預(yù)定的標準內(nèi)為止的步驟。該方法還可包括通過利用與處在預(yù)定的標準內(nèi)的SDM相關(guān)聯(lián)的地震處理序列來處理對應(yīng)于感興趣的地下的地震數(shù)據(jù)而得到增強的地震數(shù)據(jù)的步驟。地震處理序列可包括參數(shù),并且至少參數(shù)的數(shù)值在不同的處理序列之間可以是不同的。地震處理序列的不同之處至少可以在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的次序是不同的。地震處理序列的不同之處至少也可以在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的數(shù)目是不同的。本發(fā)明的一個目的是提供一種方法,其中數(shù)據(jù)處理步驟包括次序、參數(shù)和參數(shù)值被優(yōu)化,以使得為了從幅值反演得到統(tǒng)計上相關(guān)的結(jié)果只需要最小附加約束條件信息。另一個目的是提供用于估計幅值處理和幅值反演質(zhì)量的統(tǒng)計測量結(jié)果的方法。參照圖23,本發(fā)明的另一方面包括用于得到優(yōu)化的反演反射率的方法。該方法包括得到對應(yīng)于感興趣的地下區(qū)域的頻帶限制的井反射率的步驟。根據(jù)井反射率計算多個統(tǒng)計井特性(SCSweH)。得到對應(yīng)19于感興趣的地下區(qū)域的經(jīng)過處理的地震數(shù)據(jù)。通過利用多個不同的地震反演算法反演經(jīng)過處理的地震數(shù)據(jù)集而計算多個反演的反射率組。多個統(tǒng)計地震特性(SCSweU)是根據(jù)多個反演的反射率組計算的。多個統(tǒng)計差值測量結(jié)果(SDM)是在SCSwen與統(tǒng)計地質(zhì)特性SCs^mie之間進行計算的。根據(jù)計算出的SDM選擇最佳地震反演算法。地震數(shù)據(jù)通過使用最佳地震反演算法而被反演,以產(chǎn)生優(yōu)化反演反射率組。在所述方法中,不同的反演算法可以具有包含不同的數(shù)值的參數(shù)。該方法還可包括計算對于反演的反射率的不確定性數(shù)值的步驟。通過以下的說明、待決權(quán)利要求和附圖,將更好地理解本發(fā)明的這些和其它目的、特征與優(yōu)點,其中圖l顯示對于其尋求地震信息的地下構(gòu)造的示意性剖視圖;圖2是顯示對于每個軌跡的中點和反射角的爆破剖面軌跡集合;圖3是顯示在中點和反射點之間的關(guān)系的共同中點(CMP)集合;圖4顯示表示幅值反演與幅值分析的目的的幅值對偏移(AVO)地震數(shù)據(jù)的巻積模型;圖5是顯示確定反射幅值的因素的示意圖6顯示從法線到臨界入射角的反射幅值的典型的變化;圖7顯示在噪聲抑制處理步驟之前和之后施加的正常移出(NMO)的爆破剖面;圖8是顯示影響地震小波的幅值的各種因素的示意圖;圖9顯示處理序列的各種步驟對于地震軌跡的集合的影響;圖IO顯示包含噪聲和信號失真以及被輸入到幅值反演的、輸出處理的后成像調(diào)節(jié)的數(shù)據(jù)的輸入爆破剖面;圖11是用于確定從井記錄和地震數(shù)據(jù)得到的統(tǒng)計差值測量結(jié)果的實施例的概貌;圖12是描述在本發(fā)明的實施例中為了從井記錄數(shù)據(jù)和地震數(shù)據(jù)產(chǎn)生優(yōu)化的地震反射率組和統(tǒng)計差值測量結(jié)果(SDM)所采取的步驟的流程圖13是顯示為了利用優(yōu)化的統(tǒng)計差值測量結(jié)果(SDM)產(chǎn)生增強的地震數(shù)據(jù)所采取的步驟的實施例的流程圖14是顯示用于利用統(tǒng)計差值測量結(jié)果(SDM)將增強的地震數(shù)據(jù)處理成地震數(shù)據(jù)集和相關(guān)的不確定性值的步驟的流程圖15顯示在對于7個不同類型的反射率的4個不同的井中測量的RMS水平的變化;圖16顯示在7個不同反射率對之間在4個不同的井中測量的相關(guān)系數(shù)的變化;圖17顯示使用SDM選擇角度范圍參數(shù)和選擇反演衰減參數(shù)的例子;圖18顯示當(dāng)被優(yōu)化的衰減常數(shù)被用于反射率反演時,在反演的反射率輸出中的改進;圖19A和19B比較從未限制的和最佳限制的(衰減的)幅值反演得到的梯度部分和RMS水平;圖20顯示來自圖19的未限制的和限制的(衰減的)梯度顯示的S/(S+N)估計。圖21是描述在本發(fā)明的實施例中為了從井記錄數(shù)據(jù)和地震數(shù)據(jù)產(chǎn)生優(yōu)化的地震反射率組和統(tǒng)計差值測量結(jié)果(SDM)所采取的步驟的流程圖22是描述在本發(fā)明的另一個實施例中為了從井記錄和地震數(shù)據(jù)產(chǎn)生增強的地震數(shù)據(jù)和統(tǒng)計差值測量結(jié)果(SDM)所采取的步驟的流程圖;以及圖23是描述在本發(fā)明的另一個實施例中為了從井記錄和地震數(shù)據(jù)產(chǎn)生優(yōu)化的地震反射率組和統(tǒng)計差值測量結(jié)果(SDM)所采取的步驟的流程圖。具體實施例方式本發(fā)明克服了在背景部分中描述的許多缺點。本發(fā)明的總體概貌顯示于圖11。在從井控制與從反演的地震數(shù)據(jù)測量的統(tǒng)計結(jié)果之間的差值可以通過選擇井記錄處理、地震處理和反演參數(shù)的方法而被最小化。本發(fā)明是通過在多個處理步驟中利用在可能的參數(shù)范圍內(nèi)的地震數(shù)據(jù)的試驗反演而實現(xiàn)的。利用相關(guān)系數(shù)和均方根(RMS)水平的函數(shù)值被用來測量對于所選擇的參數(shù)的"良好性"或統(tǒng)計距離測量(SDM)(SDM是在反演的地震數(shù)據(jù)和井?dāng)?shù)據(jù)的統(tǒng)計特性中的差值的度量值)。導(dǎo)致具有與井控制統(tǒng)計特性最佳匹配的統(tǒng)計特性的反演的反射率的那些參數(shù)被選擇,其中優(yōu)化在幅值處理和反演序列中的參數(shù)選擇。圖12是為了得到優(yōu)化的反演的地震反射率組而可以采取的步驟的實施例的概貌流程圖。在步驟100,從井控制和/或地震數(shù)據(jù)中選擇感興趣的地下區(qū)域。在步驟110,得到可被用來估計地下壓縮波、剪切波和密度特性的井記錄數(shù)據(jù)。接著,在步驟120,從井記錄數(shù)據(jù)中選擇來自感興趣的探查或開發(fā)的區(qū)域的井記錄數(shù)據(jù)子集。在步驟130,選擇包括井記錄處理序列的一個或多個井記錄處理算法,每個具有相關(guān)聯(lián)的處理參數(shù)。在步驟140,通過使用井記錄處理序列把井記錄數(shù)據(jù)子集轉(zhuǎn)換成頻帶限制的井反射率。然后在步驟150對井反射率進行運算,獲得井反射率的統(tǒng)計井特性(SCweU)。SC油mie典型地包括零滯后自相關(guān)和互相關(guān)的比率和/或相關(guān)系數(shù)與在包容在步驟ioo規(guī)定的感興趣的時間/深度地區(qū)和感興趣區(qū)域的間隔上不同的反射率之間的RMS水平。SCswdl將包括可以從AVO幅值反演算法輸出的以下的反射率的各種函數(shù)定義△Vp/2Vp+Ap/2p-A(pVp)/pVp△VP/2Vp2(Vs/Vp)2(厶p/2p+AVs/Vs)=A(pVs2)/(pVP2)△VP/2VP-(2Vs/VP)2(Ap/2p+AVs/Vs)=RP-2RsH厶p/2p反射率符號p-阻抗p-速度RP混合剪切RsH梯度G密度RPS-阻抗RSIAVs/2Vs+Ap/2p=A(pVs)/pVsS-速度RsAVs/2Vs在步驟160,收集從感興趣的地下區(qū)域生成的地震數(shù)據(jù),在圖12的步驟200,通過使用多個地震處理步驟210把地震數(shù)據(jù)處理成增強的地震數(shù)據(jù)。每個步驟210的目標是選擇將對用于最終的幅值反演的增強的地震數(shù)據(jù)進行優(yōu)化的處理步驟、參數(shù)和參數(shù)值。步驟300包含步驟310-360,用于優(yōu)化幅值反演參數(shù)并用于反演增強的地震數(shù)據(jù),以輸出優(yōu)化的地震反射率組和相關(guān)聯(lián)的不確定性數(shù)值?,F(xiàn)在參照圖13,圖中更詳細地顯示步驟200和210。在步驟205選擇對應(yīng)于步驟120的井記錄數(shù)據(jù)子集的地震數(shù)據(jù)子集。在步驟215,如圖10所示的那些處理步驟組被選擇以處理在步驟205所選擇的地震數(shù)據(jù)。每個處理步驟可以具有一個或多個處理參數(shù),在許多情況下,步驟215可被用來優(yōu)化那些處理步驟和/或在處理序列中的處理步驟的次序以及處理參數(shù)值。某些處理參數(shù)具有固定的數(shù)值,而其它處理參數(shù)可以具有要被優(yōu)化的可變的數(shù)值。然后在步驟220,通過使用來自215的第一試驗步驟組和在處理序列中的參數(shù)值來處理在步驟205所選擇的地震數(shù)據(jù),以產(chǎn)生增強的地震數(shù)據(jù)。在步驟225,通過使用所選擇的固定反演參數(shù)組來反演這個增強的地震數(shù)據(jù)集,以產(chǎn)生相應(yīng)的地震反射率組。在步驟230,對這個反演的地震反射率組進行操作,達到具有與在步驟150從井記錄數(shù)據(jù)確定的形式相同的形式的統(tǒng)計地震特性(SC^mic)。在步驟240,通過使用統(tǒng)計特性(SC)與具有以下形式的公式的差值或比率,來計算統(tǒng)計差值測量結(jié)果(SDM):SDM=:E(SC附",i-SCsei加ic.,)2/N或SDM=S(SCweU,i/SC汰ismic,i-l)2/N(8)其中N-在地震數(shù)據(jù)子集中CRP的數(shù)目;SC^,,尸在第i個CRP處的統(tǒng)計井特性;SCseismie,在第i個CRP處的統(tǒng)計地震特性。在步驟245,采用各種策略來選擇新的試驗參數(shù)值,它將有效地定位SDM表面上的最小值。如果預(yù)期的SDM表面是簡單的,則對于參數(shù)值的網(wǎng)格的SDM的系統(tǒng)估計和重復(fù)步驟220-240可以有效地定位最佳的(即,最小的)SDM。作為替代,在圍繞試驗參數(shù)值的區(qū)域中計算的數(shù)值的和/或分析的推導(dǎo)可被用來創(chuàng)建SDM表面的局部模型。這個模型然后被用來選擇接近于建模的SDM表面的最小值的新的試驗參數(shù)值組。在步驟250,選擇對應(yīng)于最佳SDM的處理參數(shù)值組。優(yōu)選地,在SDM之間進行內(nèi)插的任選的步驟248可被用來進一步改善優(yōu)化的參數(shù)值組。在步驟260,用在步驟250所選擇的處理參數(shù)值組來處理整個地震數(shù)據(jù)集,以產(chǎn)生增強的地震數(shù)據(jù)作為單個處理步驟210的輸出。對于每個所選擇的處理步驟組,重復(fù)進行步驟210,以便產(chǎn)生增強的地震數(shù)據(jù)作為步驟200的總的輸出和作為加到最后的幅值反演的輸入。如圖12所示,地震數(shù)據(jù)可以通過幾個不同的地震處理步驟被處理,其中的某些步驟可以使用步驟210和圖13的通常的SDM方法,以及其中的另外的步驟可能不需要SDM方法,以便優(yōu)化用于反演的參數(shù)值。典型地,初始處理步驟可以是非-SDM,而SDM步驟對于在這里的表4的右面列出的后成像步驟是重要的。通常,步驟210用來通過使得增強的地震數(shù)據(jù)更容易反演而克服以前的處理參數(shù)選擇的缺點。這減小了對于在圖14的最后的幅值反演步驟中大的衰減權(quán)重和其它約束條件的要求。參照描述步驟300的圖14,來自步驟200的增強的地震數(shù)據(jù)的最后的輸出,通過使用步驟310-380的優(yōu)化的幅值反演流程,對于反射率進行反演。在步驟310,得到來自步驟200的數(shù)據(jù),然后在步驟320,從步驟200的增強的經(jīng)過處理的地震數(shù)據(jù)中選擇增強的地震數(shù)據(jù)子集。這個選擇的地震數(shù)據(jù)子集優(yōu)選地對應(yīng)于在步驟150選擇的井24記錄數(shù)據(jù)子集。在第一優(yōu)選的方式下,地震數(shù)據(jù)子集的這個選擇包括對應(yīng)于與井記錄數(shù)據(jù)相同的深度間隔的所有的地震軌跡。作為替代,地震數(shù)據(jù)的選擇可以是整個地震數(shù)據(jù)探查的大量篩選后的子集。在步驟330,通過使用所選擇的地震反演算法和試驗的反演參數(shù)值組,把來自320的增強的地震數(shù)據(jù)的這個子集反演成頻帶限制的反射率。正如這里討論的,幅值反演算法可以是幾個公式表示法之一。在步驟340,對反演的反射率進行運算,獲得具有與步驟150的統(tǒng)計井特性(SCweU)相同的分析形式的統(tǒng)計地震特性(SCsei,ie)。在步驟350,步驟150的統(tǒng)計井特性(SCwen)被投影到地震CRP集合位置,并與步驟340的統(tǒng)計地震特性(SCseismic)相比較,通過使用以上的/>式(8)的形式之一,產(chǎn)生對于試驗幅值反演參數(shù)值的統(tǒng)計差值測量結(jié)果(SDM)。對于步驟360,用試驗反演參數(shù)值重復(fù)進行步驟330-350,用來定位最佳的(最小的)SDM。優(yōu)選的方法使用CRP的子集,用反演參數(shù)值的網(wǎng)格對反演參數(shù)空間進行采樣,計算在網(wǎng)格上的SDM,然后在步驟365選擇對應(yīng)于最小的SDM的反演參數(shù)。在步驟370,整個增強的地震數(shù)據(jù)通過使用與步驟365的優(yōu)化的SDM相關(guān)聯(lián)的反演參數(shù)值組被反演,以產(chǎn)生最佳的反演反射率組。在步驟380,對于每個反演的反射率計算反比于反演的反射率的S/(S+N)[信號/(信號+噪聲)比率的不確定性數(shù)值。在步驟3卯,輸出最佳的反演反射率組和相關(guān)聯(lián)的不確定性,供在其它幅值分析過程中使用。例如,事件的交截和梯度屬性可以互相交叉繪出,以檢測顯著增加的碳氬化合物飽和水平的潛在區(qū)域。再次參照圖12,將更詳細地描述上述步驟中的每一個。在步驟100,選擇感興趣的地下或地質(zhì)區(qū)域。在典型的感興趣的天然氣和石油區(qū)域中,將進行鉆井,以便在各種感興趣的地質(zhì)構(gòu)造中勘探或產(chǎn)生碳氫化合物。在附加的勘探或開發(fā)鉆探的過程中,地震數(shù)據(jù)被用來產(chǎn)生關(guān)于這些感興趣的地質(zhì)區(qū)域的附加信息,由此提高鉆出成功的井的概率。選擇具有特征的彈性特性的感興趣的區(qū)域,導(dǎo)致幅值反演時的診斷地震響應(yīng)。地下構(gòu)造可包括一個或多個儲層或產(chǎn)生區(qū)域。另外,優(yōu)選地,井位于地下構(gòu)造中,記錄工具在井中運行,以得到關(guān)于地質(zhì)構(gòu)造的彈性特性的地下特性信息。同樣地,通過使用用來記錄反射的聲波的地震聲源和接收機對感興趣的地下區(qū)域進行地震探查。在步驟110,獲取可被用來校準和增強地震數(shù)據(jù)的井記錄。更優(yōu)選地,通過使用直接測量彈性特性的記錄工具來得到井記錄數(shù)據(jù)。作為例子,而不是限制,這些彈性特性包括壓縮速度(Vp)、剪切速度(Vs)和塊密度(p)。井記錄數(shù)據(jù)的其它非限制例子可被用來當(dāng)彈性記錄數(shù)據(jù)不可用時,幫助推導(dǎo)出彈性特性,包括gamma、電阻率和NMR記錄。也可以不用直接記錄鉆井而得到彈性特性或統(tǒng)計資料。例如,在類似的地質(zhì)設(shè)置值中可以使用來自類似鉆井或井記錄數(shù)據(jù)庫的數(shù)據(jù)。作為替代,可以從創(chuàng)建偽井記錄數(shù)據(jù)的統(tǒng)計沉積模型、地質(zhì)模型或其它知曉的估計推出彈性特性或統(tǒng)計資料。當(dāng)然,本領(lǐng)域技術(shù)人員將意識到,按照本發(fā)明,也可以獲得和使用其它類型的數(shù)據(jù),諸如巖芯數(shù)據(jù)。在步驟120,選擇對應(yīng)于來自步驟110的感興趣的地下區(qū)域的井記錄數(shù)據(jù)子集。作為替代,可以從適合于地下或感興趣的地質(zhì)區(qū)域的上述其它源選擇彈性特性或統(tǒng)計資料。此外,當(dāng)多個井穿透同一個感興趣的區(qū)域并且井記錄數(shù)據(jù)在這些井之間變化時,來自多個井的統(tǒng)計資料可以在空間上被內(nèi)插和外插到感興趣的區(qū)域。作為替代,感興趣的區(qū)域可以被擴大,以包括更大的地質(zhì)部分。在這種情況下,所投射的統(tǒng)計資料更可能是更大的感興趣區(qū)域的特性,并且不太可能以不利的方式偏離結(jié)果。在步驟130,選擇井處理算法和它們的參數(shù)。這些算法中的某些算法可被用來把井記錄數(shù)據(jù)子集變換成高信噪比的井記錄數(shù)據(jù)。例如,一系列算法可被用來校正鉆井環(huán)境影響,諸如鉆井偏離、鉆井大小、水泥粘結(jié)、鉆管校正、線深度校正、溫度影響、工具沖擊響應(yīng)等等。選擇可被用來把高信噪比井記錄數(shù)據(jù)變換成頻帶限制的反射率的其它算法。特別重要的井處理算法是把井記錄數(shù)據(jù)向上定標成在要被反演的地震數(shù)據(jù)中存在的可比較的主要波長。例如,基于等效媒介理論的Backus平均常常在井記錄處理中被使用。另一所選擇的算法是能夠把作為深度的函數(shù)的高信噪比井記錄數(shù)據(jù)變換成作為時間的函數(shù)的各種原始的反射率。計算出的反射率的例子包括密度、p阻抗、剪切阻抗、p速度、s速度、梯度反射率、混合剪切等等,如在表2中規(guī)定的??梢允褂昧硪痪幚硭惴▉戆言嫉姆瓷渎首儞Q成具有與要被反演的地震數(shù)據(jù)(表面地震數(shù)據(jù)、垂直地震剖面分布等等)相同的帶寬的頻帶限制的反射率。這個變換算法執(zhí)行所選擇的小波與原始的反射率的巻積。在步驟140,井記錄數(shù)據(jù)子集被變換成高信噪比的井記錄數(shù)據(jù),然后通過使用井記錄處理算法再被變換成頻帶限制的反射率。作為在步驟110使用井記錄數(shù)據(jù)的記錄處理算法130得到頻帶限制的反射率的替代,原始的或頻帶限制的反射率可以替代地直接從上述的數(shù)據(jù)庫、地質(zhì)模型、地震模型或其它知曉的估計中得到。在步驟150,然后對頻帶限制的井反射率進行運算,得到統(tǒng)計井特性(SCs^)。這些特性是從足夠大的深度間隔取得的,所以在感興趣區(qū)域中它們的自然預(yù)期的空間變化是慢的,并且沒有受到局部異常結(jié)構(gòu)的強烈影響,這使得SCSweU在處理或反演質(zhì)量控制中是更有用的。圖15和16是對于在表2中列出的某些反射率從間隔開10公里的4個井測量的統(tǒng)計結(jié)果的例子。圖15所示的RMS和圖16所示的相關(guān)系數(shù)有相當(dāng)小的變化程度,這表明它們可被用作為對這個區(qū)域的統(tǒng)計特性的輸入。除了這些統(tǒng)計測量在這個區(qū)域中很慢變化以夕卜,在不同的盆地的其它井中的測量結(jié)果對于相關(guān)系數(shù)和對于相對RMS水平也表現(xiàn)出非常相似的數(shù)值。這支持這些統(tǒng)計結(jié)果對于在非常遠離井控制的盆地中SC的質(zhì)量控制的使用。對于其中分量反射率在井之間更快速地變化的情況,必須把想要的統(tǒng)計特性在空間上內(nèi)插或外插到感興趣的CRP集合位置。因為反射率是相對幅值的函數(shù),而不是絕對幅值的函數(shù),統(tǒng)計特性(SCswic)優(yōu)選地是無單位的比率。作為具體的例子,在兩個反射率之間的統(tǒng)計井特性(SCwell)可以通過計算p阻抗的RMS與在包括在步驟100規(guī)定的感興趣區(qū)域的深度范圍上的梯度反射率的比率。這個比率被證明是AVO處理質(zhì)量的極其診斷的和可預(yù)測的指示符。另一個有用的診斷的測量值是測得的p阻抗(零偏移)反射率與測得的梯度反射率的相關(guān)系數(shù)。這些測量值可以組合在一起,得到用反射率之一的自相關(guān)進行歸一化的互相關(guān),以產(chǎn)生本身可以在空間上預(yù)測的統(tǒng)計特性。在p阻抗與梯度之間的有用的相對統(tǒng)計特性包括SC=aGMr0SC=CCro,g和SC=CCro,gCg/cJro=Oro,g他rojro其中SC-在p阻抗與梯度之間的統(tǒng)計特性O(shè)ro,og=p阻抗與4弟度的rmsCC則,(^在p阻抗與梯度之間的相關(guān)系數(shù)①ro,g-在P阻抗與梯度之間的互相關(guān)①r0,r0-P阻抗的自相關(guān)通常,來自兩個或更多個反射率的互相關(guān)和自相關(guān)可以組合成協(xié)方差矩陣,其特性包括它的特征向量,可以是不同的地質(zhì)環(huán)境的彈性特性的特征。在優(yōu)選的模式中,選擇一組兩個或更多個反射率,在反射率之間計算包含所有的零記錄自相關(guān)和互相關(guān)的協(xié)方差矩陣。這個矩陣的主特征向量的特征值產(chǎn)生統(tǒng)計特性,它自然地加權(quán)所有的反射率的統(tǒng)計效果,并可以預(yù)期它具有比起每個單個反射率的特性更小的空間和時間變化SCwell=理論特征值[協(xié)方差(Ro,Rb…Rw)](10)其中SCwe,尸統(tǒng)計井特性11(),111,,..1^=在地震幅值反演中感興趣的反射率協(xié)方差=(①r0,r0OrO,R1Or0,R2...①r0,rnOri,R0。r1'r10r1'r2...①r1,rn。rn,r0。rn,r1。rn,r2…。rn,rn)另一個替代方案涉及計算在反射率對之間的互相關(guān)與自相關(guān)的比率(即,信噪比),然后以各種方式組合和加權(quán)它們,利用空間變化的期望差值。加權(quán)因子可被用來以下述方式偏置井特性,使其趨向于所選擇的反射率組的信噪比SCwdl=SWi(O)Ro,Ri怖ro,ro)=SWj(CCR0,RioRi/aR0)(11)其中W產(chǎn)正比于統(tǒng)計測量值的信噪比的權(quán)重R。-具有高相關(guān)和信噪比的反射率R產(chǎn)第i個感興趣的反射率<5Ro,cvRi=Rq和Ri反射率的rmsCC則,r產(chǎn)在Ro與Ri反射率之間的相關(guān)系數(shù)①r0,w=在Ro與Ri反射率之間的互相關(guān)①r0,r。=Rq反射率的自相關(guān)通過這兩個公式(10)和(11),能夠根據(jù)可以是在特定的地質(zhì)設(shè)置值中感興趣的那些反射率計算統(tǒng)計井特性(SCSwe,i)。在步驟160,例如,獲取在感興趣的地下上的地震數(shù)據(jù),如以上對于步驟IOO描述的。這個地震數(shù)據(jù)可以是在各種地質(zhì)設(shè)置值下獲取的,包括在陸地、近?;蛏钏恢茫约熬哂懈鞣N現(xiàn)場類型,包括之D、3D、4D、VSP、多分量和時間滯后。優(yōu)選地,這個地震數(shù)據(jù)通過密度源和接收機釆樣并使用爆破到接收機偏移被記錄,允許在感興趣的地下中的高角度地下反射(高達60度)。另一個非常想要得到的特性是利用高帶寬源和小的接收機陣列。在步驟200,使用一系列處理步驟210來處理地震數(shù)據(jù)。在步驟200優(yōu)選地使得地震數(shù)據(jù)信噪比和帶寬最大化,以得到增強的地震數(shù)據(jù),其適合于作為對幅值反演的輸入。具體地,生成地震幅值,它優(yōu)選地處在頻帶限制的、平面波、與角度有關(guān)的地下反射率的單個標量常數(shù)內(nèi)。這樣的地震數(shù)據(jù)在幅值反演內(nèi)將達到需要最小的約束條件的所希望的目標。幅值處理序列可以具有20個或更多的步驟,包括在表4中顯示的那些步驟。主要地震處理步驟的非限制性例子典型地包括a)調(diào)整地震數(shù)據(jù);b)噪聲抑制;c)校正各種地球傳輸/獲取影響;d)去巻積;e)速度和各向異性估計;f)遷移;g)小波頻鐠均衡/去拉伸;h)小波相位均衡;i)剩余事件對準;j)空間幅值平衡;l)偏移幅值平衡;m)空間和時間濾波;n)偏移/角度選擇和靜音;以及o)反射率反演。這些處理步驟的總的目標,除了反演步驟以外,是從地震數(shù)據(jù)中去除噪聲,和校正在每個交界面處的響應(yīng),使得它對應(yīng)于具有已知的幅值、相位和時序的小波與交界面的與角度有關(guān)的反射系數(shù)的巻積。因為步驟g)到n)和被用于這些步驟的參數(shù)值會顯著影響反演步驟的結(jié)果,它們是用于步驟300的SDM方法的有竟?fàn)幜Φ暮蜻x方案,以及僅僅改變這樣的步驟的次序或執(zhí)行某些選擇的步驟而不執(zhí)行其他步驟。對于其它預(yù)成像步驟,SDM作為參數(shù)選擇過程可以是非常低效的和無能的。再次參照圖13,將更詳細地描述關(guān)于步驟210的示例性步驟組。在步驟205,從步驟160的地震數(shù)據(jù)中選擇包括對應(yīng)于感興趣區(qū)域的數(shù)據(jù)的地震數(shù)據(jù)集。優(yōu)選地,這個地震數(shù)據(jù)子集包括步驟160的整個數(shù)據(jù)集(以及選擇處理步驟來增強對應(yīng)于感興趣的地區(qū)的地震數(shù)據(jù)集)。還有可能選擇來自步驟160的小的適當(dāng)?shù)牡卣饠?shù)據(jù)子集,在另外的包括步驟210的步驟中處理。在步驟215,選擇地震處理步驟組。地震處理步驟包括處理要在一個或多個處理步驟中被優(yōu)化的參數(shù)和參數(shù)值的步驟組。用于必須被選擇的那些參數(shù)地震處理步驟的例子包括但不限于以上對于步驟200和在表4中列出的所有的處理步驟。選擇還可以涉及步驟本身的選擇,以及處理序列中的步驟的次序也可以變化。如上所述,一次優(yōu)化所有的步驟中的所有的參數(shù)在計算上是很難的。從以上的處理算法組,選擇一個算法或算法的組合,使得它的參數(shù)在優(yōu)化期間變化。例如,可以選擇在去巻積算法中的參數(shù)用于優(yōu)化,而用于干預(yù)處理步驟,即諸如噪聲抑制和用于速度反演的參數(shù)保持不變。因此,在每個重要的步驟內(nèi),參數(shù)組被優(yōu)化,而在其余步驟中的參數(shù)^f皮設(shè)置和保持在固定數(shù)值。接著在步驟220,通過使用所選擇的地震處理步驟和至少一個可變參數(shù)值,而其它參數(shù)值是固定的,把地震數(shù)據(jù)子集處理成增強的地震數(shù)據(jù)。在步驟225,通過使用所選擇的幅值反演算法來反演增強的地震數(shù)據(jù),它利用公式(2)或(6)或者AVO分析研究中使用的其它公式,和固定的反演參數(shù)組來生成反演的反射率。例如,反演算法可以是二項算法,它通過使用恒定的衰減參數(shù)來反演交截和梯度反射率。作為替代,三項反演算法可被用來產(chǎn)生三個所選擇的反射率的組。本領(lǐng)域技術(shù)人員將意識到,在這個步驟中可以使用許多其它替代反演算法。在步驟230,對反演的反射率進行運算,得到對于一組試驗參數(shù)值的統(tǒng)計地震特性。優(yōu)選地,被用于井?dāng)?shù)據(jù)的同一個特性算法被用于地震統(tǒng)計特性。被用于這個計算的反演反射率組是從感興趣的同一區(qū)31域選擇的,作為步驟100。來自2D或3D地震數(shù)據(jù)集所有的公共反射點集合(CRP)或它的子集是從感興趣的地區(qū)選擇的。對于每個CRP,由在該CRP處的反演反射率計算統(tǒng)計特性。這些統(tǒng)計特性被平均,以產(chǎn)生對于特定的試驗參數(shù)值組的單個統(tǒng)計特性。在來自步驟150的統(tǒng)計井特性與步驟230的統(tǒng)計地質(zhì)特性(SCseis^)之間的差值是通過使用公式(8)計算的,以產(chǎn)生統(tǒng)計差值測量結(jié)果值(SDM),即,對于特定的試驗處理參數(shù)值組的"統(tǒng)計差值測量結(jié)果"。對于本說明書,術(shù)語"良好性度量"指示統(tǒng)計地質(zhì)特性(SCsei,ic)與統(tǒng)計地質(zhì)特性(SCwel,)的相關(guān)。接著在步驟245選擇要在試驗處理算法中使用的新的試驗處理參數(shù)值組,而其余的處理算法的參數(shù)值再次保持不變。重復(fù)步驟220-240的目的是找出對于試驗參數(shù)值組的最佳SDM。如果期望的SDM表面是簡單的,則對于參數(shù)值的網(wǎng)格的SDM的系統(tǒng)估計可以有效地定位最小的SDM。每個參數(shù)具有潛在的數(shù)值的預(yù)期范圍,它們可被用于規(guī)定可能的參數(shù)值的網(wǎng)格??赡艿膮?shù)值的網(wǎng)格被用來通過重復(fù)執(zhí)行步驟220-240以計算在網(wǎng)格位置處的SDM值而創(chuàng)建SDM曲線或表面。如果預(yù)期SDM表面是非常復(fù)雜的,或當(dāng)SDM表面是多維時,在圍繞初始試驗參數(shù)值的區(qū)域中計算的數(shù)值的和/或分析的推導(dǎo)結(jié)果可被用來創(chuàng)建SDM表面的局部模型。這個模型然后被用來選擇接近于建模的SDM表面的最小值的參數(shù)值的新的試驗組。本領(lǐng)域技術(shù)人員將意識到,可以使用許多其他技術(shù)來采樣SDM表面,得出參數(shù)的數(shù)組。在步驟245,采用各種策略來選擇新的試驗參數(shù)值,它們將有效地定位在SDM表面上的最小值。如上所述,如果預(yù)期的SDM表面是簡單的,對于參數(shù)值的網(wǎng)格的SDM的系統(tǒng)估計和重復(fù)執(zhí)行步驟220-240可以有效地定位最佳的(最小的)SDM。在圖17中的SDM曲線和表面都通過使用參數(shù)值的搜索網(wǎng)格而被生成。作為替代,在圍繞初始試驗參數(shù)值的區(qū)域中計算的數(shù)值的和/或分析的推導(dǎo)結(jié)果可被用來創(chuàng)建SDM表面的局部模型。這個模型然后被用來選擇接近于建模的SDM表面的最小值的參數(shù)值的新的試驗組。在步驟250,選擇處理對應(yīng)于最佳SDM的參數(shù)值的組。圖17是SDM如何可被用來選擇對于特定的感興趣的地區(qū)的數(shù)據(jù)選擇(角度靜音)參數(shù)的一個例子,該圖顯示SDM值作為不同角度靜音值的函數(shù)的曲線。它規(guī)定了最大反射角,通過使得計算出的SDM最小化而積極地對幅值反演做出貢獻。來自大于這個角度的反射角的幅值數(shù)據(jù)然后從CRP集合中被去除(被靜音)。在步驟260,整個地震數(shù)據(jù)集用在步驟250中所選擇的處理參數(shù)值組進行處理,產(chǎn)生增強的地震數(shù)據(jù)作為來自單個處理步驟210的輸出。整個數(shù)據(jù)集僅僅通過所選擇的步驟210被處理,它是對下一個所選擇的步驟210的輸入數(shù)據(jù)。再次參照圖14,該圖中描述步驟300,來自步驟200的增強的地震數(shù)據(jù)的最終輸出通過使用步驟310-380的優(yōu)化的幅值反演流程而對于反射率進行反演。在步驟310得到來自步驟200的數(shù)據(jù),以及在步驟320從步驟200的增強的經(jīng)過處理的地震數(shù)據(jù)中選擇增強的地震數(shù)據(jù)子集。這個所選擇的地震數(shù)據(jù)子集應(yīng)當(dāng)對應(yīng)于在步驟150所選擇的井記錄子集。按第一優(yōu)選的方式,地震數(shù)據(jù)子集的這種選擇包括對應(yīng)于與井記錄數(shù)據(jù)相同的深度間隔的所有的地震軌跡。作為替代,地震數(shù)據(jù)的選擇可以是整個地震數(shù)據(jù)勘查的經(jīng)過大量篩選后的子集。在步驟330,增強的地震數(shù)據(jù)子集通過使用所選擇的地震反演算法和試驗的反演參數(shù)值組而被反演成頻帶限制的反演的反射率。在優(yōu)選的例子中,反演算法可以是對于方程(6)的誤差函數(shù)的受到約束的最小平方解,它需要衰減參數(shù)和數(shù)據(jù)加權(quán)參數(shù)的數(shù)值。通常,為了減小在地震數(shù)據(jù)中剩余噪聲和信號失真的影響,通過選擇衰減參數(shù)值和數(shù)據(jù)權(quán)重(即,在方程(6)中的Wi,Wc和Wf)而使反射率互相約束。因為衰減趨于遠離幅值數(shù)據(jù)的偏置結(jié)果,大的衰減權(quán)重是不希望的,許可地震幅值數(shù)據(jù)本身對于最終得到的反射率具有的影響少于全部影響。對于非常大的衰減因子,某些屬性可以有效地被約束為其它屬性的縮放的版本。作為替代,不帶有R2項的方程(6)可被用來確定反演參數(shù)值。本領(lǐng)域技術(shù)人員將意識到,也可以使用其它反演算法,而不受這些例子限制。在步驟340,然后對頻帶限制的地震反射率進行運算,得到等同于從井控制生成的和以上在步驟150中描述的各種統(tǒng)計地震特性。優(yōu)選地,由于它的良好的橫向可預(yù)測性,所選擇的SC將是從反演地震反射率計算出的協(xié)方差矩陣的主要特征向量。在步驟350,將統(tǒng)計地震特性與統(tǒng)計井特性進行比較,通過使用公式(8)的形式之一產(chǎn)生統(tǒng)計差值測量結(jié)果(SDM)。在步驟360,重復(fù)步驟330-350,定位幅值反演參數(shù)值,它使得SDM優(yōu)化(最小化)。用于這樣做的策略類似于在步驟245中的策略。優(yōu)選地,這個步驟是通過使用公式(6)對于可能的衰減參數(shù)值的網(wǎng)格進行試驗幅值反演而實現(xiàn)的。SDM是對于每個衰減參數(shù)值的輸出,并分析數(shù)值的網(wǎng)格找出最小值。作為替代,可以分析SDM的推導(dǎo)結(jié)果作為參數(shù)值的函數(shù),并且識別局部最小值,得到參數(shù)值的優(yōu)化的組合。還可以檢驗SDM,看它們是否對應(yīng)于令人滿意的衰減或處理參數(shù)。極其重要的標準是與最小SDM相關(guān)聯(lián)的衰減參數(shù)的大小。在幾乎相等的SDM值的兩個局部最小值被定位的情形下,選擇與最小衰減參數(shù)相關(guān)聯(lián)的那個最小值。還可以使用的另一個標準是反演的反射率的空間相干性。這個標準必須結(jié)合其他標準一起使用,因為很可能反射率像梯度那樣同時是非??臻g相干的和非常錯誤的。再一個標準可以是SDM的大小。如果這些預(yù)定的標準中的任一標準不滿足,則可以用新選擇的處理算法和/或地震參數(shù)重復(fù)執(zhí)行步驟210-370,直至預(yù)定的標準滿足為止。如果SDM足夠低,則反演的地震反射率是令人滿意的??墒褂玫钠渌鼉?yōu)化參數(shù)搜索程序的例子,作為例子而不是限制,包括對誤差表面采用泰勒級數(shù)表達式或牛頓近似的方法。本領(lǐng)域技術(shù)人員將意識到,可以利用許多其它方法來有效地搜索最佳SDM,包括使用以前確定的、在感興趣的區(qū)域中的反演參數(shù)。圖17的右面顯示對應(yīng)于梯度和密度衰減參數(shù)值的網(wǎng)格的SDM數(shù)值的輪廓面。選擇最小的SDM值和產(chǎn)生該數(shù)值的相關(guān)聯(lián)的參數(shù)。最小值是圖17中對應(yīng)于0.14的密度衰減值和0.06的梯度衰減的最后閉合的輪廓,這些數(shù)值粗略地正比于如在井控制中測量的這些參數(shù)的相對RMS水平。SDM表面表現(xiàn)出對于衰減參數(shù)值的小的改變有大的變化的事實,反映了試圖不采用試驗參數(shù)值而選擇衰減參數(shù)的困難。事實上,因為根據(jù)使用空間相干性或絕對RMS水平的現(xiàn)有的質(zhì)量控制方法很難選擇衰減參數(shù)值,在許多反演算法中常常甚至不使用約束條件(在/>式(6)中Wc-O)。在步驟370,通過使用在步驟360中選擇的反演參數(shù)值組來反演增強的地震數(shù)據(jù),以產(chǎn)生優(yōu)化反演的反射率組。圖18顯示對于合成數(shù)據(jù)集的未限制的(未衰減的)和限制的(衰減的)反演的p阻抗、混合剪切和p速度反射率。所顯示的兩種情況具有被添加到CRP集合信號軌跡的50%和100%噪聲(相干和不相干)。因為p阻抗、混合剪切和p速度反射率信號在每個DRP處是相同的,并且互相高度相關(guān),正確的結(jié)果應(yīng)當(dāng)是幾乎相同的和空間上相干的反射率。如圖18的頂部所示,對于p阻抗的未限制的反演產(chǎn)生良好的結(jié)果,而混合剪切和p速度反射率具有低信噪比,正如通過空間非相干和相當(dāng)高的RMS數(shù)值所證實的那樣。在圖18的下面部分的顯示表明當(dāng)通過使用SDM方法選擇最佳衰減約束條件時,混合剪切和p速度反射率的空間相干和相對RMS大大地改善。圖19A和19B是真實地震數(shù)據(jù)例子,與在左面的未衰減的梯度部分與在右面的最佳衰減的梯度部分比較,其中最佳衰減的參數(shù)值是通過使用本專利的SDM方法被選擇的。在頂部顯示的梯度RMS水平的大的下降和空間相干性的大的增加都是SDM衰減的梯度部分比起未衰減的梯度部分具有高得多的信噪比的指示。不僅僅需要反演的反射率,而且也需要它們的不確定性((N+S)/S)作為定量幅值分析技術(shù)的輸入。在步驟380,針對成對的反射率計算的SDM可被用來估計對于每個反演反射率的相對S/(S+N)數(shù)值(不確定性"。這些可以通過使用以下的公式計算針對地震反演反射率計算的統(tǒng)計特性(SC)與從井控制投影的SC之間的差值而得到Uri=SCweU,RiR2/l((SCweii,rir2,-scseismic,RlR2)i+SCwe",RlR2))(9)其中Ur尸與反射率Rj目關(guān)聯(lián)的S/(S+N)比率(不確定性")來自井控制的、對于R1,R2的0R1,R2/€>R1,R1互相關(guān)和自相關(guān)SCseismic,Rl,R2—來自地震的、對于R1,R2的0Ri,R2/0R1,R1互相關(guān)和自相關(guān)同樣地,可以通過比較在受到限制的(衰減的,加權(quán)的)與未受限制的幅值反演之間在輸出上的差值而估計它們。反射率之一的SDM,優(yōu)選地為p阻抗或堆疊反射率,相對于其S/(S+N)估計的想要的SDM具有高的S/(S+N)。它們利用在成對的軌跡之間的相同種類的統(tǒng)計特性來量化在每個反演的反射率中的不確定性。圖20顯示對于圖19所示的梯度部分的S/(S+N)的改變作為CPR和時間的函數(shù)。當(dāng)SDM選擇的衰減參數(shù)被選擇用于反演時,S/(S+N)有大的增加。在步驟390,輸出優(yōu)化的反演的地震反射率組和相關(guān)聯(lián)的不確定性值。這些輸出在各種幅值分析技術(shù)中被用來減小與地震數(shù)據(jù)驅(qū)動勘測和現(xiàn)場輪廓決定相關(guān)聯(lián)的風(fēng)險。雖然在上述的說明書中本發(fā)明是相對于本發(fā)明的某些優(yōu)選實施例描述的,以及為了說明起見闡述了許多細節(jié),但本領(lǐng)域技術(shù)人員將會看到,本發(fā)明易于作出改變,并且這里描述的某些其它細節(jié)可以做出很多變化而不背離本發(fā)明的基本原理。權(quán)利要求1.一種用于得到增強的地震數(shù)據(jù)和優(yōu)化的反演的反射率的方法,該方法包括(a)根據(jù)對應(yīng)于感興趣的地下區(qū)域的頻帶限制的井反射率計算多個統(tǒng)計井特性(SCswell);(b)通過把多個不同的地震處理序列應(yīng)用到對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù)集合的至少一部分而計算多個經(jīng)過處理的地震數(shù)據(jù)集;(c)通過反演多個經(jīng)過處理的地震數(shù)據(jù)集而計算多個反演的反射率組;(d)根據(jù)多個反演的反射率組計算多個統(tǒng)計地震特性(SCsseismic);(e)計算在統(tǒng)計井特性(SCswell)與多個統(tǒng)計地震特性(SCsseismic)之間的多個統(tǒng)計差值測量結(jié)果(SDM);(f)根據(jù)計算出的統(tǒng)計差值測量結(jié)果(SDM)選擇最佳地震處理序列;(g)通過利用最佳的地震處理序列來處理對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù)而得到增強的地震數(shù)據(jù);(h)根據(jù)計算出的統(tǒng)計差值測量結(jié)果(SDM)選擇最佳地震反演算法;和(i)使用最佳地震反演算法來反演地震數(shù)據(jù),以產(chǎn)生優(yōu)化的反演的反射率組。2.權(quán)利要求1的方法,其中不同的地震處理序列包括參數(shù),以及在不同的處理序列之間至少參數(shù)的數(shù)值是不同的。3.權(quán)利要求1的方法,其中不同的地震處理序列的不同之處至少在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的次序不同。4.權(quán)利要求1的方法,其中不同的地震處理序列的不同之處至少在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的數(shù)目不同。5.權(quán)利要求1的方法,其中在優(yōu)化方案中利用統(tǒng)計差值測量結(jié)果(SDM),以預(yù)測最佳地震處理序列。6.權(quán)利要求1的方法,其中不使用從記錄感興趣的區(qū)域中的鉆井所得到的井記錄數(shù)據(jù)而計算頻帶限制的井反射率。7.權(quán)利要求l的方法,其中步驟(a)還包括以下步驟(i)獲取井記錄數(shù)據(jù);(ii)用井記錄處理參數(shù)選擇井記錄處理算法;和(iii)通過使用井記錄處理算法把井記錄數(shù)據(jù)變換成頻帶限制的井反射率。8.權(quán)利要求7的方法,其中井記錄數(shù)據(jù)是從在感興趣區(qū)域中的鉆井、類似地質(zhì)設(shè)置、多個類似的井的平均值和統(tǒng)計模型中的一項獲取的。9.權(quán)利要求1的方法,其中步驟(f)還包括確定統(tǒng)計差值測量結(jié)果(SDM)中的任一個是否在預(yù)定的標準內(nèi)的步驟。10.權(quán)利要求l的方法,其中根據(jù)剪切的、壓縮的和密度反射率的RMS水平與在剪切的、壓縮的和密度反射率之間的相關(guān)系數(shù)計算統(tǒng)計井特性(SC、en)。11.權(quán)利要求l的方法,其中統(tǒng)計差值測量結(jié)果(SDM)響應(yīng)于井反射率的S/(S+N)對地震數(shù)據(jù)反射率的S/(S+N)的比率。12.權(quán)利要求8的方法,其中統(tǒng)計的井特性在Zo印pritz反射率公式中被用來計算在AVO集合中地震軌跡的預(yù)期的RMS。13.權(quán)利要求l的方法,其中統(tǒng)計差值測量結(jié)果(SDM)反比于井反射率的預(yù)期的S/(S+N)比率與地震數(shù)據(jù)反射率的S/(S+N)比率之間的差值。14.一種用于得到增強的地震數(shù)據(jù)的方法,該方法包括(a)根據(jù)對應(yīng)于感興趣的地下區(qū)域的頻帶限制的井反射率計算統(tǒng)計井特性(SCSw一組;(b)通過把多個不同的地震處理序列應(yīng)用到對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù)集的至少一部分而計算多個經(jīng)過處理的地震數(shù)據(jù)集;(c)通過利用地震反演算法用固定的參數(shù)組來反演多個經(jīng)過處理的地震數(shù)據(jù)集而計算多個反演的反射率組;(d)根據(jù)多個反演的反射率組計算統(tǒng)計地震特性(SCSsei,ic)組;(e)計算在統(tǒng)計井特性(SCSw一與多個統(tǒng)計地震特性(SCSsei,ic)之間的多個統(tǒng)計差值測量結(jié)果(SDM);(f)根據(jù)計算出的統(tǒng)計差值測量結(jié)果(SDM)選擇最佳地震處理序列;和(g)使用最佳的地震處理序列和參數(shù)來處理地震數(shù)據(jù),以得到增強的地震數(shù)據(jù)集。15.權(quán)利要求14的方法,其中地震反演算法的固定的參數(shù)組包含具有不同數(shù)值的參數(shù)。16.權(quán)利要求14的方法,還包括(h)確定統(tǒng)計差值測量結(jié)果(SDM)是否在預(yù)定的標準內(nèi);和(i)利用在步驟(b)中的不同的地震處理序列重復(fù)步驟(b)-(f),直至統(tǒng)計差值測量結(jié)果(SDM)在預(yù)定的標準內(nèi)為止。17.權(quán)利要求16的方法,還包括(j)通過利用與處在預(yù)定的標準內(nèi)的統(tǒng)計差值測量結(jié)果(SDM)相關(guān)聯(lián)的地震處理序列來處理對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù),得到增強的地震數(shù)據(jù)。18.權(quán)利要求17的方法,其中地震處理序列的固定的參數(shù)組的不同之處至少在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的次序不同。19.權(quán)利要求17的方法,其中地震處理序列的固定的參數(shù)組的不同之處至少在于,在每個相應(yīng)的地震處理序列內(nèi)的處理步驟的數(shù)目不同。20.權(quán)利要求17的方法,其中地震處理序列的固定的參數(shù)組的不同之處至少在于,在不同的處理序列之間的系數(shù)值不同。21.—種用于得到優(yōu)化的反演的反射率的方法,該方法包括(a)得到對應(yīng)于感興趣的地下區(qū)域的頻帶限制的井反射率;(b)根據(jù)對應(yīng)于感興趣的地下區(qū)域的頻帶限制的井反射率計算多個統(tǒng)計井特性(SCSwe");(c)得到對應(yīng)于感興趣的地下區(qū)域的經(jīng)過處理的地震數(shù)據(jù);(d)通過利用多個不同的地震反演算法來反演經(jīng)過處理的地震數(shù)據(jù)集而計算多個反演的反射率組;(e)根據(jù)反演的反射率組計算多個統(tǒng)計地震特性(SCXei,ic);(f)計算在統(tǒng)計井特性(SCswe)與統(tǒng)計地震特性(SCSsei訓(xùn)"之間的多個統(tǒng)計差值測量結(jié)果(SDM);(g)根據(jù)計算出的統(tǒng)計差值測量結(jié)果(SDM)選擇最佳地震反演算法;和(h)使用最佳地震反演算法來反演地震數(shù)據(jù),以產(chǎn)生優(yōu)化的反演的反射率組。22.權(quán)利要求21的方法,還包括以下步驟(i)計算對于反演的反射率的不確定性數(shù)值。全文摘要一種用于得到增強的地震數(shù)據(jù)和優(yōu)化的反演的反射率的方法包括根據(jù)對于感興趣的地下區(qū)域的頻帶限制的井反射率計算統(tǒng)計井特性。通過把不同的地震處理序列應(yīng)用到對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù)而計算經(jīng)過處理的地震數(shù)據(jù)集。通過反演經(jīng)過處理的地震數(shù)據(jù)集而計算反演的反射率組。根據(jù)反演的反射率組計算統(tǒng)計地震特性。統(tǒng)計差值測量結(jié)果(SDM)是在統(tǒng)計井特性與統(tǒng)計地震特性之間進行計算的。根據(jù)計算出的SDM選擇最佳處理序列。通過處理對應(yīng)于感興趣的地下區(qū)域的地震數(shù)據(jù)而得到增強的地震數(shù)據(jù)。根據(jù)計算出的SDM得到最佳地震反演算法。使用最佳地震反演算法來反演地震數(shù)據(jù),以產(chǎn)生優(yōu)化的反演的反射率組。文檔編號G01V1/28GK101680960SQ200880020397公開日2010年3月24日申請日期2008年6月12日優(yōu)先權(quán)日2007年6月15日發(fā)明者E·F·赫爾肯霍弗,J·沃施伯爾尼申請人:雪佛龍美國公司