一種巖體中油氣信息的分析方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于地質(zhì)勘探領(lǐng)域，設(shè)及一種巖體中油氣信息的分析方法。
【背景技術(shù)】
[0002]AV0屬性分析技術(shù)（Ampl;UudeVersusOffset,振幅隨偏移距的變化）是一項(xiàng)應(yīng) 用在地質(zhì)勘探，特別是地下油氣資源勘探領(lǐng)域的巖體特性分析技術(shù)，AV0技術(shù)通過在巖體上 放炮得到的反射波的縱波和橫波隨入射角變化的特征，來探討反射波和入射波的比，即反 射系數(shù)響應(yīng)隨入射角的變化，來確定反射界面上覆、下伏介質(zhì)的巖性特征及物性參數(shù)。AV0 技術(shù)是一項(xiàng)求取巖石的彈性參數(shù)、分析和識(shí)別巖性及含油氣性的地震勘探技術(shù)。
[0003] 反射系數(shù)與入射角的關(guān)系符合Zoeppritz(佐普里茲）方程，完整的Zoeppritz方 程形式非常復(fù)雜，物理意義不明確，很難直接分析介質(zhì)參數(shù)對(duì)振幅或反射系數(shù)的影響，為了 進(jìn)一步反映不同彈性參數(shù)與反射系數(shù)之間的內(nèi)在關(guān)系，前人做了大量的工作，分別從不同 的角度對(duì)Zoeppritz方程進(jìn)行了詳細(xì)的討論和簡(jiǎn)化，提出了不同假設(shè)條件下的縱、橫波反 射振幅的近似表達(dá)式。1985年，化uey(蘇耶）對(duì)前人各種近似進(jìn)行重組，并進(jìn)一步研究了 泊松比對(duì)反射系數(shù)的影響。首次提出了反射系數(shù)的AV0截距和梯度的概念，證明了相對(duì)反 射系數(shù)隨入射角（或炮檢距）的變化梯度主要由泊松比0的變化來決定，給出了用不同角 度項(xiàng)表示的反射系數(shù)近似公式：
[0006] R〇= (AVp/Vp+A p/p)/2
[0007] 該里Rpp為反射系數(shù)，i為入射角，Vp為縱波振幅，P為炮點(diǎn)巖體密度，AVp/ Vp，Ap/p分別表示縱波變化率和巖體密度變化率。
[0008] Shuey公式①式把反射系數(shù)視為小角度項(xiàng)（第一項(xiàng)）、中等角度項(xiàng)（第二項(xiàng)）和大 角度項(xiàng)（第S項(xiàng)）之和。當(dāng)入射角稍大時(shí)（0。<i<30° )，第S項(xiàng)的tg2i-sin2i<0.083,而 Aa/a又較小，第S項(xiàng)可W忽略，于是化U巧公式簡(jiǎn)化為化U巧兩項(xiàng)式；
[0009]
[0010] 利用②式對(duì)小角度和中角度入射角的各個(gè)反射系數(shù)進(jìn)行采集，通過數(shù)學(xué)分析得出 簡(jiǎn)化的化uey兩項(xiàng)式②的各項(xiàng)系數(shù)，將②式簡(jiǎn)寫為：
[0011]
[001引求出截距P和梯度GW后，可W計(jì)算得出泊松比0。
[0013] 從化uey近似公式②可W看出，只有第二項(xiàng)包含顯示含油氣特征的泊松比0，而 該式中的第二項(xiàng)除了受到泊松比O的影響外，還要受到地層速度變化率Aa/a、密度變 化率Ap/p的影響，仍然不能較為敏感地反映地層泊松比o的變化。
[0014] 常規(guī)AV0分析一般是基于化uey兩項(xiàng)式，而舍棄了大角度項(xiàng)，眾所周知，大角度道 集含有更多的橫波信息，因此利用化uey近似公式②對(duì)各系數(shù)進(jìn)行擬合，所得出的結(jié)果不 能完全反映橫波的真實(shí)信息，此外，舍棄大角度項(xiàng)意味著引入了不可能消除的理論誤差，因 此得出的化uey兩項(xiàng)式的截距P和梯度G存在誤差，得出的泊松比參數(shù)0也不夠精確。從 而影響對(duì)地底巖層油氣資源判斷的準(zhǔn)確性。

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

[0015] 為克服現(xiàn)有技術(shù)中舍棄了大角度道集，造成泊松比0求解不夠精確，影響判斷地 底巖層油氣資源的技術(shù)問題。本發(fā)明提供一種巖體中油氣信息的分析方法。
[0016] 本發(fā)明的技術(shù)方案如下：巖體中油氣信息的分析方法，包括如下步驟：
[0017] (1)對(duì)目標(biāo)巖體炮擊，在不同的入射角i利用檢波器測(cè)得對(duì)應(yīng)的反射系數(shù)化P和入 射角i，取得若干組對(duì)應(yīng)的入射角i和反射系數(shù)化P的檢測(cè)值；
[0018] (2)對(duì)入射角i大于〇的所有檢測(cè)值，構(gòu)造樣本集H(x，y);其中，樣本集每個(gè)樣本 的X=sin2i，y=Rpp*cos2i;〇為預(yù)先設(shè)定的大角度道集中入射角i的下限值；
[0019] (3)使樣本集H對(duì)方程F:y=R+Wx+Vx2進(jìn)行擬合，得出R、W、V的擬合值；
[0020] 其中，樣本集H是采用最小平方法擬合得出R、W、V的值；最小平方法擬合包括如 下步驟：
[0021] (31)設(shè)定信噪比參數(shù)M1，剔除樣本集H中信噪比大于Ml的樣本點(diǎn)，剔除后的樣本 集為H1，記錄H1與H的樣本點(diǎn)的比例P1%，如果P1%大于預(yù)先設(shè)定的樣本剔除最大系數(shù) 60%時(shí)，進(jìn)入步驟（34)，否則進(jìn)入步驟（32);
[0022] (32)利用樣本集化對(duì)方程F進(jìn)行初始擬合，得出R、W、V的值，將擬合出的R、W、V 帶入F方程，計(jì)算擬合中使用的各個(gè)樣本點(diǎn)相對(duì)擬合后的方程F的誤差，如果全部樣本點(diǎn)的 誤差都小于預(yù)先設(shè)定的誤差參數(shù)EPS,則進(jìn)入步驟（35)，否則進(jìn)入步驟（33);
[0023] (33)剔除誤差大于邸S的樣本點(diǎn)構(gòu)成樣本集肥，肥與原始樣本集H的樣本點(diǎn)比 例如果大于80 %，則進(jìn)入步驟（34)，否則進(jìn)入步驟（32);
[0024] (34)利用樣本集H1對(duì)方程F進(jìn)行初始擬合，得出R、W、V的值；
[0025] (35)保存輸出最后得出的R、W、V;
[002引 （4)根據(jù)R=P，W= - (P-G)，V= (C-G)，求得P、G、C的值，其中為截距P、G為梯 度G、C為擬橫波C，利用截距P、梯度G、擬橫波C對(duì)巖體油氣信息進(jìn)行屬性分析。
[0027] 在擬合過程中采用不斷逼近AV0標(biāo)準(zhǔn)曲線的方法，其基本思路為通過逐步擬 合-剔除-擬合的過程來求取穩(wěn)定的各AV0屬性，使得真實(shí)結(jié)果更加可靠。
[0028] 具體的，如權(quán)利要求1或2所述的巖體中油氣信息的分析方法，其特征在于：所述 步驟似中的〇為45度。
[0029] 進(jìn)一步，所述步驟（2)還包括：對(duì)入射角不大于〇的所有檢測(cè)值構(gòu)造樣本集 化(sin2i，Rpp)，利用Rpp=P+Gsin2i和樣本集化對(duì)截距P、梯度G進(jìn)行擬合計(jì)算。
[0030]采用本發(fā)明所述的巖體中油氣信息的分析方法，對(duì)入射大角度道集進(jìn)行了AV0分 析，消除了對(duì)化uey方程擬合前存在的理論誤差，得到的P、G、C系數(shù)更精確，更真實(shí)的反映 了目標(biāo)巖體的油氣信息。
【具體實(shí)施方式】
[0031] 下面結(jié)合實(shí)施例，對(duì)本發(fā)明作進(jìn)一步地詳細(xì)說明，但本發(fā)明的實(shí)施方式不限于此。
[0032] 本發(fā)明所述巖體中油氣信息的分析方法，包括如下步驟：
[0033] (1)對(duì)目標(biāo)巖體炮擊，在不同的入射角i利用檢波器測(cè)得對(duì)應(yīng)的反射系數(shù)化P和入 射角i，取得若干組對(duì)應(yīng)的入射角i和反射系數(shù)化P的檢測(cè)值；
[0034] (2)對(duì)入射角i大于〇的所有檢測(cè)值，構(gòu)造樣本集H(x，y);其中，樣本集每個(gè)樣本 的X=sin2i，y=R