本發(fā)明涉及地震勘探領(lǐng)域，特別涉及一種裂縫識(shí)別方法和裝置。
背景技術(shù)：
：在地震勘探領(lǐng)域中，對(duì)裂縫的識(shí)別和預(yù)測(cè)有著重要的意義。因?yàn)椋趯?shí)際地震勘探中，裂縫既可以充當(dāng)流體的儲(chǔ)集空間，又可以作為流體的運(yùn)移通道，將儲(chǔ)層中分散的儲(chǔ)集空間聯(lián)通起來。所以，有效裂縫的預(yù)測(cè)在實(shí)際勘探中很重要，是決定儲(chǔ)層勘探開發(fā)的關(guān)鍵。具體勘探時(shí)，由于地層上覆地層的壓實(shí)作用，水平裂縫和低角度裂縫幾乎消失，保留下來的一般都是高角度裂縫，而這類裂縫使儲(chǔ)層具有方位各向異性特性，即，地震波的物理性質(zhì)(例如：振幅、速度和頻率等)會(huì)隨著觀測(cè)方位的變換而變化。所以，通過分析裂縫誘導(dǎo)的方位各向異性特性可以得到裂縫的相關(guān)參數(shù)，如裂縫密度、裂縫方位和裂縫流體指示因子，進(jìn)而通過這些數(shù)據(jù)，確定和評(píng)價(jià)儲(chǔ)層內(nèi)油氣的運(yùn)移方向、聚散程度和開發(fā)遠(yuǎn)景等。實(shí)際實(shí)施時(shí)，過去主要采用疊后屬性裂縫檢測(cè)法對(duì)裂縫進(jìn)行識(shí)別。其中，疊后屬性裂縫檢測(cè)法主要是基于疊后三維地震數(shù)據(jù)體，通過研究道間不連續(xù)變化以及構(gòu)造變形來預(yù)測(cè)斷層和構(gòu)造突變區(qū)，從而間接預(yù)測(cè)裂縫發(fā)育帶。但這種方法不能針對(duì)有效裂縫(即開啟并填充流體的裂縫)進(jìn)行預(yù)測(cè)，這也是疊后預(yù)測(cè)方法的局限性。為提高預(yù)測(cè)精度，目前地震工業(yè)大多用疊前方位AVO反演方法進(jìn)行裂縫識(shí)別，它主要基于Rüger反射系數(shù)方程的前兩項(xiàng)求解，不使用第三項(xiàng)，即含sin2θtan2θ項(xiàng)。實(shí)際上，在用Rüger反射系數(shù)方程描述地下某個(gè)成像點(diǎn)的各向異性特征時(shí)，受第三項(xiàng)的影響很大，尤其是當(dāng)入射角相對(duì)較大時(shí)，這種影響更為明顯。并且，隨著地震采集技術(shù)的提高，具有較大偏移距的地震資料越來越多，相應(yīng)地，反射波入射角范圍也較大。當(dāng)入射角范圍較大(如大于30°)時(shí)，Rüger反射系數(shù)方程的第三項(xiàng)，即含sin2θtan2θ項(xiàng)，對(duì)縱波反射系數(shù)的貢獻(xiàn)會(huì)很大，不能忽略。但由于實(shí)施時(shí)未考慮這一項(xiàng)，通過現(xiàn)有技術(shù)進(jìn)行裂縫識(shí)別，往往存在預(yù)測(cè)精度不高且不能對(duì)裂縫流體進(jìn)行定量識(shí)別的技術(shù)問題針對(duì)上述問題，目前尚未提出有效的解決方案。技術(shù)實(shí)現(xiàn)要素：本發(fā)明實(shí)施例提供了一種裂縫識(shí)別方法和裝置，通過求解三項(xiàng)Rüger反射系數(shù)方程，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)，利用上述數(shù)據(jù)，達(dá)到對(duì)裂縫的高精度識(shí)別和對(duì)裂縫內(nèi)裂縫流體定量識(shí)別的目的。本發(fā)明實(shí)施例提供了一種裂縫識(shí)別方法，包括：獲取n個(gè)方位的偏移距道集；根據(jù)層速度將所述n個(gè)方位的偏移距道集轉(zhuǎn)化為n個(gè)方位的角道集，其中，n為大于等于4的正整數(shù)；將所述n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加，得到所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)；根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)；根據(jù)所述裂縫密度和所述各向異性參數(shù)，求解得到裂縫流體指示因子；根據(jù)所述裂縫密度和所述裂縫方位角識(shí)別裂縫，根據(jù)所述裂縫流體指示因子識(shí)別裂縫內(nèi)的裂縫流體類型。在一個(gè)實(shí)施方式中，根據(jù)所述裂縫密度和所述裂縫方位角識(shí)別裂縫，根據(jù)所述裂縫流體指示因子識(shí)別裂縫內(nèi)的裂縫流體類型，包括：根據(jù)所述裂縫密度識(shí)別裂縫的密度；根據(jù)所述裂縫方位角識(shí)別裂縫的方位；根據(jù)所述裂縫流體指示因子定量識(shí)別裂縫內(nèi)的裂縫流體類型。在一個(gè)實(shí)施方式中，將所述n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加，得到所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，包括：將所述n個(gè)方位的角道集中每個(gè)方位的角度范圍等分為多個(gè)分角度范圍；將所述n個(gè)方位的角道集中每個(gè)方位的角道集按照所述分角度范圍進(jìn)行疊加，得到所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)。在一個(gè)實(shí)施方式中，根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)，包括：根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到第一待反演數(shù)據(jù)和第二待反演數(shù)據(jù)；根據(jù)所述第一待反演數(shù)據(jù)，求解得到裂縫方位角和裂縫密度；根據(jù)所述裂縫方位角和所述第二待反演數(shù)據(jù)，求解得到各向異性參數(shù)。在一個(gè)實(shí)施方式中，根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，按照以下方程求解所述第一待反演數(shù)據(jù)和所述第二待反演數(shù)據(jù)：其中，為分角度疊加數(shù)據(jù)，為采集線的方位角，θ為入射角，為第一待反演數(shù)據(jù)，為第二待反演數(shù)據(jù)，A為法向入射的反射系數(shù)。在一個(gè)實(shí)施方式中，根據(jù)所述第一待反演數(shù)據(jù)，按照以下方程求解所述裂縫方位角和所述裂縫密度：其中，為第一待反演數(shù)據(jù)，Biso為各向同性梯度，Bani為裂縫密度，φs為裂縫方位角，為采集線的方位角。在一個(gè)實(shí)施方式中，根據(jù)所述裂縫方位角和所述第二待反演數(shù)據(jù)，按照以下方程求解所述各向異性參數(shù)εν和δν：其中，為第二待反演數(shù)據(jù)，為反射界面兩側(cè)介質(zhì)中縱波速度平均值，Δα為反射界面兩側(cè)介質(zhì)中縱波速度的差值，εν和δν為各向異性參數(shù)，φs為裂縫方位角，為采集線的方位角。在一個(gè)實(shí)施方式中，根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)，包括：根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，通過噪音衰減算法，求解得到所述裂縫密度、所述裂縫方位角和所述各向異性參數(shù)。基于相同的發(fā)明構(gòu)思本發(fā)明實(shí)施例還提供了一種裂縫識(shí)別裝置，包括：獲取模塊，用于獲取n個(gè)方位的偏移距道集；轉(zhuǎn)化模塊，用于根據(jù)層速度將所述n個(gè)方位的偏移距道集轉(zhuǎn)化為n個(gè)方位的角道集，其中，n為大于等于4的正整數(shù)；疊加模塊，用于將所述n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加，得到所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)；第一求解模塊，用于根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)；第二求解模塊，用于根據(jù)所述裂縫密度和所述各向異性參數(shù)，求解得到裂縫流體指示因子；識(shí)別模塊，用于根據(jù)所述裂縫密度、所述裂縫方位角識(shí)別裂縫，根據(jù)所述裂縫流體指示因子識(shí)別裂縫內(nèi)的裂縫流體類型。在一個(gè)實(shí)施方式中，所述疊加模塊，用于將所述n個(gè)方位的角道集中每個(gè)方位的角度范圍按照預(yù)設(shè)規(guī)則等分為多個(gè)分角度范圍；再將所述n個(gè)方位的角道集中每個(gè)方位的角道集按照所述分角度范圍進(jìn)行疊加，得到所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)。在一個(gè)實(shí)施方式中，所述第一求解模塊，包括：求取待反演數(shù)據(jù)單元，用于根據(jù)所述n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到第一待反演數(shù)據(jù)和第二待反演數(shù)據(jù)；求取裂縫方位角和裂縫密度單元，用于根據(jù)所述第一待反演數(shù)據(jù)，求解得到裂縫方位角和裂縫密度；求取各向異性參數(shù)單元，用于根據(jù)所述裂縫方位角和所述第二待反演數(shù)據(jù)，求解得到各向異性參數(shù)。在一個(gè)實(shí)施方式中，所述第一求解模塊，通過噪音衰減算法求解得到所述裂縫密度、所述裂縫方位角和所述各向異性參數(shù)。在一個(gè)實(shí)施方式中，所述識(shí)別模塊，用于根據(jù)所述裂縫密度識(shí)別裂縫的密度；根據(jù)所述裂縫方位角識(shí)別裂縫的方位；根據(jù)所述裂縫流體指示因子定量識(shí)別裂縫流體類型。在本發(fā)明實(shí)施例中，考慮到裂縫的存在使得儲(chǔ)層具有方位各向異性特性，因此通過基于三項(xiàng)Rüger反射系數(shù)方程的疊前分步方位AVO方法，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)，利用這些數(shù)據(jù)便可以實(shí)現(xiàn)對(duì)裂縫的高精度識(shí)別和對(duì)裂縫流體的定量識(shí)別，從而解決了現(xiàn)有技術(shù)中裂縫識(shí)別精度低，不能識(shí)別裂縫流體的問題，從而達(dá)到了有效提高裂縫識(shí)別精度和對(duì)裂縫流體進(jìn)行定量識(shí)別的技術(shù)效果。附圖說明此處所說明的附圖用來提供對(duì)本發(fā)明的進(jìn)一步理解，構(gòu)成本申請(qǐng)的一部分，并不構(gòu)成對(duì)本發(fā)明的限定。在附圖中：圖1是根據(jù)本發(fā)明實(shí)施例的裂縫識(shí)別方法的處理流程圖；圖2是根據(jù)本發(fā)明實(shí)施例的裂縫識(shí)別方法中求解裂縫密度、裂縫方位角和各向異性參數(shù)的處理流程圖；圖3是根據(jù)本發(fā)明實(shí)施例的裂縫識(shí)別方法具體識(shí)別裂縫的過程流程圖；圖4是根據(jù)本發(fā)明實(shí)施例的裂縫識(shí)別裝置的結(jié)構(gòu)示意圖；圖5是根據(jù)本發(fā)明實(shí)施例的裂縫識(shí)別裝置中第一求解模塊的結(jié)構(gòu)示意圖；圖6是根據(jù)本發(fā)明實(shí)施例的25°入射角對(duì)應(yīng)的方位角道集示意圖；圖7是根據(jù)本發(fā)明實(shí)施例的不含噪音的情況下基于三項(xiàng)Rüger反射系數(shù)方程的疊前分步方位AVO方法和基于兩項(xiàng)Rüger反射系數(shù)方程的常規(guī)方位AVO反演方法的結(jié)果對(duì)比圖；圖8是根據(jù)本發(fā)明實(shí)施例的含15％隨機(jī)噪音的情況下基于三項(xiàng)Rüger反射系數(shù)方程的疊前分步方位AVO方法和基于兩項(xiàng)Rüger反射系數(shù)方程的常規(guī)方位AVO反演方法的結(jié)果對(duì)比圖；圖9是根據(jù)本發(fā)明實(shí)施例求解得到的裂縫流體指示因子的理論值、分步方位AVO反演結(jié)果的對(duì)比圖；圖10是根據(jù)本發(fā)明實(shí)施例的過井疊前時(shí)間偏移剖面圖；圖11是根據(jù)本發(fā)明實(shí)施例的目標(biāo)層段的成像測(cè)井解釋結(jié)果和統(tǒng)計(jì)分析的目標(biāo)層段的裂縫的傾角示意圖；圖12是根據(jù)本發(fā)明實(shí)施例的基于常規(guī)方位AVO反演方法得到的裂縫密度剖面圖；圖13是根據(jù)本發(fā)明實(shí)施例的基于三項(xiàng)Rüger反射系數(shù)方程的疊前分步方位AVO方法得到的裂縫密度剖面圖；圖14是根據(jù)本發(fā)明實(shí)施例的裂縫流體指示因子示意圖。具體實(shí)施方式為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下面結(jié)合實(shí)施方式和附圖，對(duì)本發(fā)明做進(jìn)一步詳細(xì)說明。在此，本發(fā)明的示意性實(shí)施方式及其說明用于解釋本發(fā)明，但并不作為對(duì)本發(fā)明的限定。目前大多采用的疊前方位AVO反演方法識(shí)別裂縫之所以存在識(shí)別精度低，且不能對(duì)裂縫流體進(jìn)行定量識(shí)別的技術(shù)問題，其根本原因是由于現(xiàn)有技術(shù)利用的是基于Rüger反射系數(shù)方程的前兩項(xiàng)，沒有使用第三項(xiàng)，即含sin2θtan2θ項(xiàng)。實(shí)際上，在用Rüger反射系數(shù)方程描述地下某個(gè)成像點(diǎn)的各向異性特征時(shí)，受第三項(xiàng)的影響很大，尤其是當(dāng)入射角相對(duì)較大時(shí)，這種影響更為明顯。目前隨著地震采集技術(shù)的提高，具有較大偏移距的地震資料越來越多，相應(yīng)地，反射波入射角范圍也較大。當(dāng)入射角范圍較大(如大于30°)時(shí)，Rüger反射系數(shù)方程中的第三項(xiàng)，即含sin2θtan2θ項(xiàng)，對(duì)縱波反射系數(shù)的貢獻(xiàn)也很會(huì)大，這時(shí)不能忽略。因此，為了獲得較精確的地下裂縫的特征參數(shù)，解決現(xiàn)有技術(shù)中存在的識(shí)別精度低，且不能對(duì)裂縫流體進(jìn)行定量識(shí)別的技術(shù)問題，考慮到關(guān)鍵是要充分利用Rüger反射系數(shù)方程中的第三項(xiàng)，即含sin2θtan2θ項(xiàng)。基于上述思路，本發(fā)明實(shí)施例提供了一種裂縫識(shí)別方法，如圖1所示，可以包括：步驟101：獲取n個(gè)方位的偏移距道集；步驟102：根據(jù)層速度將n個(gè)方位的偏移距道集轉(zhuǎn)化為n個(gè)方位的角道集，其中，n為大于等于4的正整數(shù)；上述實(shí)施例中的層速度，指的是在層狀地層中地震波傳播的速度。它能直接反映地層的巖性，可以用來劃分地層。一般層速度是用地震測(cè)井或聲波測(cè)井測(cè)得的，也可以利用反射記錄計(jì)算得到。在地震勘探中，一般把層速度低于1400米/秒的地層叫做低速層，把高于3500米/秒的地層視為高速層。上述實(shí)施例中的偏移距道集，是地下偏移距共成像道集的簡(jiǎn)稱，是角道集在逆Radon變換下的對(duì)偶量，不同的偏移據(jù)對(duì)應(yīng)不同的角道集。且上述偏移距道集一般由疊前偏移生成。利用偏移距道集與偏移速度模型的相互關(guān)系可以確定速度分析的標(biāo)準(zhǔn)和目標(biāo)函數(shù)。對(duì)應(yīng)的，角道集是根據(jù)偏移距和深度得到樣點(diǎn)的入射角度，把一定角度范圍的地震道集合一起就是角道集，具體過程可能還要考慮地層傾角因素。步驟103：將n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加，得到n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)；具體實(shí)施過程中，在進(jìn)行上述發(fā)明實(shí)施例步驟103：將n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加時(shí)，需要對(duì)每個(gè)方位內(nèi)的角度范圍按照預(yù)設(shè)規(guī)則進(jìn)行劃分，得到多個(gè)分角度范圍，然后再將n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加，得到n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)。即，步驟103可以包括：將n個(gè)方位的角道集中每個(gè)方位的角度范圍等分為多個(gè)分角度范圍；將n個(gè)方位的角道集中每個(gè)方位的角道集按照對(duì)應(yīng)的分角度范圍進(jìn)行疊加，得到n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)。其中，對(duì)每個(gè)方位內(nèi)入射角劃分時(shí)，每個(gè)方位的分角度范圍個(gè)數(shù)t不是唯一的，考慮到具體求解三項(xiàng)Rüger反射系數(shù)方程時(shí)，至少需要3個(gè)以上方程才能求解，這里的t必須為大于等于3的正整數(shù)。具體劃分的預(yù)設(shè)規(guī)則可以包括：將入射角的角度范圍等分為t個(gè)：第一個(gè)分角度范圍為第二個(gè)分角度范圍為第i個(gè)分角度范圍為第t個(gè)分角度范圍為其中，Imax為最大入射角度，Imin為最小入射角度，t為大于等于3的正整數(shù)，i為1到t之間的任意整數(shù)。這里，需要說明的是對(duì)于每個(gè)方位的入射角劃分得到的分角度范圍個(gè)數(shù)t可以是不同的，但必須是符合要求的正整數(shù)。另外，值得注意的是，上述所列舉的預(yù)設(shè)規(guī)則僅是一種示意性描述，是為了更好地說明本發(fā)明，然而，在實(shí)際執(zhí)行的過程中，還可以包括其它衍生預(yù)設(shè)規(guī)則，可以根據(jù)實(shí)際需要進(jìn)行選取，本申請(qǐng)對(duì)此不作限定。在一個(gè)具體實(shí)施方式中，利用層速度將n個(gè)方位的偏移距道集轉(zhuǎn)化為n個(gè)方位的角道集，在本實(shí)施例中，取n＝4，并對(duì)各個(gè)方位的角道集按照角度范圍：3°～14°，14°～25°，25°～36°進(jìn)行劃分，每個(gè)方位得到3個(gè)分角度范圍，對(duì)每個(gè)方位中的3個(gè)分角度范圍的角道集分別進(jìn)行疊加，得到每個(gè)方位的3個(gè)分角度疊加數(shù)據(jù)。然而，值得注意的是，上述所列舉的處理步驟僅是一種示意性描述，是為了更好地說明本發(fā)明，然而，在實(shí)際執(zhí)行的過程中，還可以包括其它衍生步驟，可以根據(jù)實(shí)際需要進(jìn)行選取調(diào)整，本申請(qǐng)對(duì)此不作限定。步驟104：根據(jù)n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)；步驟105：根據(jù)裂縫密度和各向異性參數(shù)，求解得到裂縫流體指示因子；具體地，可以將裂縫密度和各向異性參數(shù)，代入以下的裂縫流體指示因子表達(dá)式，求解得到裂縫流體指示因子：f≈2ϵv-δv-Baniδv-Bani]]>其中，f為裂縫流體指示因子，εν和δν為各向異性參數(shù)，Bani為裂縫密度。上述實(shí)施例中的裂縫流體指示因子是能夠表征地下裂縫內(nèi)的裂縫流體情況的數(shù)據(jù)，具體實(shí)施的，在裂縫飽含液體(油或水)時(shí)，裂縫流體指示因子接近0；在裂縫飽含氣體時(shí)，裂縫流體指示因子遠(yuǎn)遠(yuǎn)大于0。這些數(shù)據(jù)在一定程度上給出了定量預(yù)測(cè)裂縫內(nèi)的裂縫流體類型的標(biāo)準(zhǔn)。在一個(gè)具體的實(shí)施例中，檢測(cè)一口油井的生產(chǎn)數(shù)據(jù)，井累計(jì)產(chǎn)液13400t，其中，產(chǎn)液26.5t/天，產(chǎn)油23.2t/天，這時(shí)，裂縫流體指示因子在井位置表現(xiàn)為低值，這與實(shí)際生產(chǎn)情況吻合，驗(yàn)證了利用本發(fā)明實(shí)施例求解得到的裂縫流體因子可以較好地識(shí)別裂縫流體類型，較準(zhǔn)確地描述實(shí)際情況。步驟106：根據(jù)裂縫密度、裂縫方位角識(shí)別裂縫，根據(jù)裂縫流體指示因子識(shí)別裂縫內(nèi)的裂縫流體的類型。具體實(shí)施時(shí)，通過求解得到的裂縫密度，確定裂縫密度；通過求解得到的裂縫方位角，確定裂縫方位；通過求解得到的裂縫流體指示因子，對(duì)裂縫流體類型進(jìn)行定量識(shí)別。上述發(fā)明實(shí)施例提供的主要方法實(shí)際上就是一種基于三項(xiàng)Rüger反射系數(shù)方程的疊前方位AVO反演方法，通過該方法求解得到裂縫密度、裂縫方位角和各向異性參數(shù)，然后再利用上述裂縫密度、裂縫方位角和各向異性參數(shù)實(shí)現(xiàn)對(duì)裂縫的高精度識(shí)別和對(duì)裂縫內(nèi)裂縫流體類型的定量識(shí)別，解決了現(xiàn)有技術(shù)中裂縫識(shí)別精度低，不能識(shí)別裂縫流體類型的技術(shù)問題，從而達(dá)到了有效提高裂縫識(shí)別精度和對(duì)裂縫流體進(jìn)行定量識(shí)別的技術(shù)效果。具體實(shí)施過程中，考慮到求解裂縫密度、裂縫方位角和各向異性參數(shù)的各個(gè)過程實(shí)際是相互聯(lián)系的，是一個(gè)整體求解的過程，因此上述發(fā)明實(shí)施例中的步驟104根據(jù)n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)的具體求解過程如圖2所示，可以包括：步驟201：根據(jù)n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到第一待反演數(shù)據(jù)和第二待反演數(shù)據(jù)；步驟202：根據(jù)第一待反演數(shù)據(jù)，求解得到裂縫方位角和裂縫密度；步驟203：根據(jù)裂縫方位角和第二待反演數(shù)據(jù)，求解得到各向異性參數(shù)。需要說明的是，上述三個(gè)求解步驟是一個(gè)求解整體，互相有著緊密的聯(lián)系，是一個(gè)不能分開的整體求解過程。結(jié)合上述具體求解步驟201、步驟202和步驟203，本發(fā)明實(shí)施例中一個(gè)具體的裂縫識(shí)別方法處理流程可以如圖3所示包括：獲取n個(gè)方位的共反射點(diǎn)道集，并基于n個(gè)方位的共反射點(diǎn)道集得到n個(gè)方位的部分疊加剖面；然后，通過求解上述三項(xiàng)Rüger反射系數(shù)方程得到n個(gè)方位的第一待反演數(shù)據(jù)和第二待反演數(shù)據(jù)再根據(jù)第一待反演數(shù)據(jù)求解得到裂縫方位角φs和裂縫密度Bani；根據(jù)第二待反演數(shù)據(jù)和裂縫密度Bani，求解得到各向異性參數(shù)εν和δν；根據(jù)裂縫方位角φs和裂縫密度Bani對(duì)裂縫進(jìn)行識(shí)別；根據(jù)裂縫密度Bani、各向異性參數(shù)εν和δν確定用于裂縫流體的定量識(shí)別的裂縫流體指示因子。具體實(shí)施時(shí)，在上述步驟104中，可以通過以下公式求解得到裂縫密度、裂縫方位角和各向異性參數(shù)：其中，α為縱波速度，為反射界面兩側(cè)介質(zhì)中縱波速度平均值，β為橫波速度，為反射界面兩側(cè)介質(zhì)中橫波速度平均值，εv、δv、γ為各向異性參數(shù)，為采集線的方位角，φs為裂縫方位角，θ為入射角，Z縱波入射時(shí)的波阻抗，G為橫波的切向模量，Δα為反射界面兩側(cè)介質(zhì)中縱波速度的差值，ΔZ為反射界面兩側(cè)介質(zhì)中縱波入射波阻抗的差值，ΔG為反射界面兩側(cè)介質(zhì)中橫波的切向模量的差值。因?yàn)樯鲜龉揭还灿腥?xiàng)，即稱為三項(xiàng)Rüger反射系數(shù)方程，具體將得到的分角度疊加數(shù)據(jù)代入上述三項(xiàng)Rüger反射系數(shù)方程，便可以求得裂縫密度、裂縫方位角和各向異性參數(shù)。正如前文所提及的通過分析裂縫誘導(dǎo)的方位各向異性特性就可以獲取裂縫的相關(guān)數(shù)據(jù)，在獲取相關(guān)數(shù)據(jù)后也就相當(dāng)于實(shí)現(xiàn)了裂縫的識(shí)別。進(jìn)一步的，考慮到上述的三項(xiàng)Rüger反射系數(shù)方程在現(xiàn)有技術(shù)使用時(shí)存在求解難度大，計(jì)算耗時(shí)長(zhǎng)，求解實(shí)用性低的問題，針對(duì)上述存在的問題，基于對(duì)上述三項(xiàng)Rüger反射系數(shù)方程的理解，對(duì)其進(jìn)行了變形、推導(dǎo)。在本發(fā)明實(shí)施例中，通過變形、推導(dǎo)，將上述三項(xiàng)Rüger反射系數(shù)方程轉(zhuǎn)換為如下的三個(gè)相關(guān)的線性方程：1)2)3)通過對(duì)這三個(gè)相關(guān)的線性方式的求解就可以得到裂縫識(shí)別時(shí)所需要的裂縫密度、裂縫方位角和各向異性參數(shù)，從而實(shí)現(xiàn)了簡(jiǎn)單、快捷、高效地求解目的。下面對(duì)上述三個(gè)線性方程的推導(dǎo)和應(yīng)用方向進(jìn)行具體描述：方程1：其中，為分角度疊加數(shù)據(jù)，為采集線的方位角，θ為入射角，為第一待反演數(shù)據(jù)，為第二待反演數(shù)據(jù)，A為法向入射的反射系數(shù)。步驟201：根據(jù)n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)求解得到第一待反演數(shù)據(jù)和第二待反演數(shù)據(jù)，即為通過求解變形后的線性方程：得到第一待反演數(shù)據(jù)和第二待反演數(shù)據(jù)。方程2：其中，為第一待反演數(shù)據(jù)，Biso為各向同性梯度，Bani為裂縫密度，φs為裂縫方位角，為采集線的方位角。步驟202根據(jù)第一待反演數(shù)據(jù)，求解得到裂縫方位角和裂縫密度，即為通過求解變形后的線性方程：得到裂縫密度和裂縫方位角。方程3：其中，為第二待反演數(shù)據(jù)，為反射界面兩側(cè)介質(zhì)中縱波速度平均值，Δα為反射界面兩側(cè)介質(zhì)中縱波速度的差值，εν和δν為各向異性參數(shù)，φs為裂縫方位角，為采集線的方位角。步驟203根據(jù)裂縫方位角和第二待反演數(shù)據(jù)，求解得到各向異性參數(shù)，即為通過求解變形后的線性方程：得到各向異性參數(shù)。在具體實(shí)施過程中，考慮到常規(guī)地震資料常常由于存在采集方位范圍窄，并受隨機(jī)噪聲影響大的情況，很容易導(dǎo)致對(duì)上述的方程具體進(jìn)行求解的過程中出現(xiàn)求解穩(wěn)定性差的問題。為了增加方程求解的穩(wěn)定性，減少隨機(jī)噪聲的影響，在本發(fā)明實(shí)施例中，通過噪音衰減算法，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)。下面以通過噪聲衰減方法對(duì)上述3個(gè)線性方程進(jìn)行求解為例進(jìn)行說明，因上述三個(gè)線性方程都是線性的，因此，上述3個(gè)線性方程都可以表示為MX＝R的形式。在實(shí)際計(jì)算的過程中，可以將地震數(shù)據(jù)劃分為多個(gè)采集方位，從而將上述3個(gè)線性方程轉(zhuǎn)換為了超定線性方程組。本發(fā)明實(shí)施例的主要工作之一就是分別求取上述方程組的未知量X。其中。要求解得到合理的X，就要求不僅野外寬(全)方位采集，而且要求地震資料的信噪比要高，然后，實(shí)際資料往往不能滿足這些要求。為此，在本例中采用了噪音衰減算法來解決上述問題，具體的，就是在方程MX＝R基礎(chǔ)上，加入噪音項(xiàng)N，使原來的方程轉(zhuǎn)化為MX＝R+N，然后，再利用奇異值分解方法來解目標(biāo)函數(shù)。通過噪音衰減法求解，實(shí)現(xiàn)了穩(wěn)定求解，并且增強(qiáng)了解的唯一性?；谏鲜龅陌l(fā)明構(gòu)思，本發(fā)明實(shí)施例中還提供了一種裂縫識(shí)別裝置，如下面所述。由于裂縫識(shí)別裝置原理與裂縫識(shí)別方法原理相似，因此涉及到一些相同的技術(shù)和原理在此不再贅述。以下所使用的，術(shù)語“單元”或者“模塊”可以實(shí)現(xiàn)預(yù)定功能的軟件和/或硬件的組合。盡管以下實(shí)施例所描述的裝置較佳地以軟件來實(shí)現(xiàn)，但是硬件，或者軟件和硬件的組合的實(shí)現(xiàn)也是可能并被構(gòu)想的。圖4是本發(fā)明實(shí)施例的一種裂縫識(shí)別裝置的結(jié)構(gòu)框圖。如圖4所示，可以包括：獲取模塊401，轉(zhuǎn)化模塊402，疊加模塊403，第一求解模塊404，第二求解模塊405，識(shí)別模塊406。下面結(jié)合圖4對(duì)該結(jié)構(gòu)進(jìn)行具體說明。獲取模塊401，用于獲取n個(gè)方位的偏移距道集；轉(zhuǎn)化模塊402，用于根據(jù)層速度將n個(gè)方位的偏移距道集轉(zhuǎn)化為n個(gè)方位的角道集，其中，n為大于等于4的正整數(shù)；疊加模塊403，用于將n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加，得到n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)；第一求解模塊404，用于根據(jù)n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到裂縫密度、裂縫方位角和各向異性參數(shù)；第二求解模塊405，用于根據(jù)裂縫密度和各向異性參數(shù)，求解得到裂縫流體指示因子；識(shí)別模塊406，用于根據(jù)裂縫密度、裂縫方位角識(shí)別裂縫，根據(jù)裂縫流體指示因子識(shí)別裂縫內(nèi)的裂縫流體。考慮到具體操作過程中，在將n個(gè)方位的角道集按照預(yù)設(shè)分角度范圍疊加，得到n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)時(shí)，需要將n個(gè)方位每個(gè)方位的角度范圍按照預(yù)設(shè)規(guī)則等分為多個(gè)分角度范圍，因此，上述實(shí)施例中的疊加模塊403，用于將n個(gè)方位的角道集中每個(gè)方位的角度范圍按照預(yù)設(shè)規(guī)則等分為多個(gè)分角度范圍；再將n個(gè)方位的角道集中每個(gè)方位的角道集按照分角度范圍進(jìn)行疊加，得到n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)。考慮到具體實(shí)施時(shí)，第一求解模塊404要處理的是3個(gè)方程，并且對(duì)這3個(gè)方程的求解過程是一個(gè)相互聯(lián)系的整體。為了在一個(gè)整體中，分步求解上述的3個(gè)方程，本發(fā)明實(shí)施例中的第一求解模塊404如圖5所示，可以包括：求取待反演數(shù)據(jù)單元501：用于根據(jù)n個(gè)方位的角道集中每個(gè)方位的分角度疊加數(shù)據(jù)，求解得到第一待反演數(shù)據(jù)和第二待反演數(shù)據(jù)；求取裂縫方位角和裂縫密度單元502：用于根據(jù)第一待反演數(shù)據(jù)，求解得到裂縫方位角和裂縫密度；求取各向異性參數(shù)單元503，用于根據(jù)裂縫方位角和第二待反演數(shù)據(jù)，求解得到各向異性參數(shù)?？紤]到具體求解過程中，求解相關(guān)方程可能會(huì)出現(xiàn)求解不穩(wěn)定的問題，為了提高求解的穩(wěn)定性，上述實(shí)施例中的第一求解模塊404中的求取待反演數(shù)據(jù)單元501，求取裂縫方位角和裂縫密度單元502和求取各向異性參數(shù)單元503具體都是通過噪音衰減算法求解得到裂縫密度、裂縫方位角和各向異性參數(shù)?？紤]到實(shí)際實(shí)施時(shí)，通過第一求解模塊404求解得到了裂縫密度、裂縫方位角和各向異性參數(shù)，通過第二求解模塊405求解得到了裂縫流體指示因子后，還需要對(duì)上述數(shù)據(jù)進(jìn)行處理，從而對(duì)裂縫進(jìn)行綜合識(shí)別，因此，本發(fā)明實(shí)施例中的識(shí)別模塊406，用于根據(jù)裂縫密度識(shí)別裂縫的密度；根據(jù)裂縫方位角識(shí)別裂縫的方位；根據(jù)裂縫流體指示因子定量識(shí)別裂縫流體的類型。從而，實(shí)現(xiàn)了對(duì)裂縫的識(shí)別。根據(jù)本發(fā)明實(shí)施例，將裂縫識(shí)別方法和裝置應(yīng)用到具體的施工中，在一個(gè)具體實(shí)施例中，實(shí)際實(shí)施過程還包含數(shù)值模型的應(yīng)用和實(shí)際資料的應(yīng)用，本發(fā)明實(shí)施例通過數(shù)值模型來驗(yàn)證本發(fā)明實(shí)施例的可行性和有效性，進(jìn)而指導(dǎo)應(yīng)用本發(fā)明實(shí)施例處理實(shí)際資料。如圖6所示，是入射角為25°時(shí)用數(shù)值模型正演模擬方位角道集用于反演的結(jié)果示意圖。其實(shí)施過程可以包括：首先進(jìn)行模型正演，設(shè)計(jì)一雙層介質(zhì)模型，上層為均勻各向同性介質(zhì)層，下層為裂縫介質(zhì)層，模型參數(shù)如表1，并且裂縫中充填氣體。同時(shí)，設(shè)定裂縫縱橫比為0.001，裂縫密度范圍為0.005～0.15，間隔為0.005，即共有30個(gè)CDP點(diǎn)，每個(gè)CDP點(diǎn)對(duì)應(yīng)不同的裂縫密度值。利用各向異性等效介質(zhì)理論求得各個(gè)CDP點(diǎn)的各向異性參數(shù)。設(shè)定入射角為10°、20°、30°和40°，方位角為0°、20°、40°、60°、80°、100°、120°、140°、160°、180°。利用三項(xiàng)Rüger反射系數(shù)方程求解得到反射系數(shù)數(shù)據(jù)，再將30個(gè)CDP點(diǎn)的反射系數(shù)序列分別與40Hz的零相位雷克子波進(jìn)行褶積，得到的合成地震記錄。通過如圖6所示入射角為25°時(shí)的方位道集，可以看出不同方位間反射振幅有明顯差異。表1各向同性/各向異性雙層單界面模型參數(shù)縱波速度(m/s)橫波速度(m/s)密度(kg/m3)各向同性層3670.02000.02400裂縫介質(zhì)層4500.02530.02750具體實(shí)施過程中，當(dāng)不考慮隨機(jī)噪聲時(shí)，本發(fā)明實(shí)施例由于采用基于三項(xiàng)Rüger反射系數(shù)方程分步方位AVO反演方法，相較于基于兩項(xiàng)Rüger反射系數(shù)方程常規(guī)方位AVO反演方法，求解得到的數(shù)據(jù)精度更高，更加接近真實(shí)值。在一個(gè)具體實(shí)施例中，不考慮隨機(jī)噪聲時(shí)，采用基于三項(xiàng)Rüger反射系數(shù)方程分步方位AVO反演方法和基于兩項(xiàng)Rüger反射系數(shù)方程常規(guī)方位AVO反演方法求解結(jié)果比較如圖7所示。其中，圖7是不含噪音的情況下基于三項(xiàng)Rüger反射系數(shù)方程的分步方位AVO反演和基于兩項(xiàng)三項(xiàng)Rüger反射系數(shù)方程的常規(guī)方位AVO反演的結(jié)果對(duì)比圖。圖中分步方位AVO反演結(jié)果中包括裂縫密度Bani、裂縫方位φs，各向異性參數(shù)εv和δv。需要說明的是本發(fā)明實(shí)施例所展示的均是反射界面處的結(jié)果。具體，圖7(a)是通過不同方法求得的裂縫密度Bani對(duì)比結(jié)果，其中，實(shí)線表示理論值，帶星實(shí)線表示常規(guī)反演的結(jié)果，即方位AVO反演的結(jié)果，帶圓圈虛線表示分步方位AVO反演的結(jié)果?？梢钥闯?，在裂縫內(nèi)充填流體的情況下，用常規(guī)方法反演的Bani值偏離理論值很多，尤其是裂縫密度較大的情況，而通過基于三項(xiàng)反射系數(shù)方程的分步方位AVO反演的Bani值非常接近于理論值。圖7(b)是兩種方法反演的裂縫方位φs對(duì)比圖，可以看出在沒有噪音干擾的情況下，兩種反演方法得到的裂縫方位均和理論值吻合。圖7(c)和圖7(d)分別是通過分步方位AVO反演得到的各向異性參數(shù)εv和δv與理論值對(duì)比圖，其中，實(shí)線表示理論值，帶星點(diǎn)的虛線表示分步方位AVO反演的結(jié)果。從圖7(c)和圖7(d)中可以看出，通過本發(fā)明實(shí)施例反演得到的εv和δv都比較接近于理論值。具體實(shí)施過程中，當(dāng)考慮了隨機(jī)噪聲后，本發(fā)明實(shí)施例由于采用基于三項(xiàng)Rüger反射系數(shù)方程分步方位AVO反演方法，相較于基于兩項(xiàng)Rüger反射系數(shù)方程常規(guī)方位AVO反演方法，求解得到的數(shù)據(jù)精度更高，更加接近真實(shí)值。在一個(gè)具體實(shí)施例中，考慮隨機(jī)噪聲時(shí)，采用基于三項(xiàng)Rüger反射系數(shù)方程分步方位AVO反演方法和基于兩項(xiàng)Rüger反射系數(shù)方程常規(guī)方位AVO反演方法求解結(jié)果比較如圖8所示。其中，圖8是含15％隨機(jī)噪音的情況下基于三項(xiàng)Rüger反射系數(shù)方程的分步方位AVO反演和基于兩項(xiàng)Rüger反射系數(shù)方程的常規(guī)方位AVO反演的結(jié)果對(duì)比圖。圖8不僅說明了本發(fā)明實(shí)施例的反演方法比常規(guī)方位AVO反演更精確，更說明了本發(fā)明實(shí)施例的方法有一定的抗噪音能力。其具體求解步驟同不含噪音的情況相同。具體，圖8(a)是通過不同方法求得的裂縫密度Bani對(duì)比結(jié)果，其中，實(shí)線表示理論值，帶星實(shí)線表示常規(guī)方位AVO反演的結(jié)果，帶圓圈虛線表示分步方位AVO反演的結(jié)果。可以看出，所求得的裂縫密度Bani的整體變化特征和無噪音時(shí)候相同，由于有噪音的存在常規(guī)方位AVO反演結(jié)果、分步方位AVO反演結(jié)果均表現(xiàn)不同程度的抖動(dòng)，但是通過本發(fā)明實(shí)施例提出的反演方法得到的Bani更接近于理論值。類似于圖8(a)，圖(b)是含15％隨機(jī)噪音的情況下，通過不同方法求得的裂縫方位角φs對(duì)比結(jié)果圖。圖8(c)和圖8(d)分別是通過分步方位AVO反演得到的各向異性參數(shù)εv和δv與理論值對(duì)比圖，其中，實(shí)線表示理論值，帶星點(diǎn)的虛線表示分步方位AVO反演的結(jié)果?？梢钥闯觯ㄟ^本發(fā)明實(shí)施例的方法反演得到的εv和δv受到噪音的影響都表現(xiàn)不同程度地偏離理論值，但是都在合理范圍內(nèi)浮動(dòng)，也就是說通過本發(fā)明實(shí)施例得到的反演結(jié)果有一定的可信度。在一個(gè)具體實(shí)施例中，將本發(fā)明實(shí)施例分別用于以下四種情況中，通過求解裂縫流體指示因子，確定裂縫內(nèi)的裂縫流體類型：情況1、無噪聲下100％飽含水的情況，情況2、無噪聲下100％包含氣的情況，情況3、含15％噪聲下100％飽含水的情況，情況4、含15％噪聲下100％飽含氣的情況。如圖9所示是裂縫流體指示因子的理論值、分步方位AVO反演結(jié)果的對(duì)比圖。針對(duì)上述四種情況，在其他模型參數(shù)不變的條件下按著上面的步驟求出了裂縫內(nèi)充填水(包括無噪音和有噪音兩種情況)的裂縫密度Bani，各向異性參數(shù)εv和δv，繼而可以求出裂縫流體指示因子。理論研究表明當(dāng)裂縫內(nèi)充填油或水時(shí)，裂縫流體指示因子的值接近0；當(dāng)裂縫內(nèi)充填氣時(shí)裂縫流體指示因子不為0，且其值遠(yuǎn)大于0，據(jù)此可以用裂縫流體指示因子判斷裂縫內(nèi)的裂縫流體類型。通過圖9可以看出，當(dāng)裂縫內(nèi)飽含水時(shí)，裂縫流體指示因子較小接近于0，當(dāng)裂縫充填氣時(shí)，裂縫流體指示因子的值較大，通過裂縫流體指示因子，上述兩種情況可以明顯被區(qū)分開來。雖然受隨機(jī)噪音的影響，裂縫流體指示因子在一定程度上偏離理論值，但是這并不影響利用根據(jù)本發(fā)明實(shí)施例計(jì)算得到的裂縫流體因子區(qū)分裂縫流體類型。通過上述三個(gè)實(shí)施例，根據(jù)數(shù)值模擬表明，相比于基于兩項(xiàng)Rüger反射系數(shù)方程的常規(guī)方位AVO反演方法，基于三項(xiàng)Rüger反射系數(shù)方程的分步方位AVO反演方法得到的裂縫密度和裂縫方位的精度更高；且利用疊前分步方位AVO反演的結(jié)果得到的裂縫流體指示因子確實(shí)能夠很好地識(shí)別裂縫流體類型。在一個(gè)實(shí)施方式中，將本發(fā)明實(shí)施例應(yīng)用到某地區(qū)實(shí)際三維數(shù)據(jù)的一條測(cè)線，根據(jù)具體實(shí)施情況，可以說明本發(fā)明實(shí)施例的實(shí)際應(yīng)用效果。具體可以通過圖10-圖14所示的相關(guān)內(nèi)容，說明本發(fā)明實(shí)施例相較于原有技術(shù)，確實(shí)實(shí)現(xiàn)了對(duì)裂縫的高精度識(shí)別和對(duì)裂縫內(nèi)裂縫流體類型的定量識(shí)別。其中，所述實(shí)施例中圖10表示過井的疊前時(shí)間偏移剖面。圖11是研究井在目標(biāo)層段的成像測(cè)井信息。圖10中橢圓圈的位置對(duì)應(yīng)的井段如圖11(a)所示的深度范圍。從圖11(a)可以看出該井段裂縫非常發(fā)育。圖11(b)表示統(tǒng)計(jì)的該井段內(nèi)發(fā)育的裂縫的傾角，可以看出發(fā)育高陡裂縫，各向異性程度強(qiáng)，測(cè)井評(píng)價(jià)結(jié)果也顯示目標(biāo)層內(nèi)裂縫孔隙度比較高。圖12表示基于兩項(xiàng)Rüger反射系數(shù)方程的常規(guī)方位AVO反演的裂縫密度。圖13表示利用本發(fā)明實(shí)施例求出的過井裂縫密度。從圖12和圖13可以看出，用上述兩種方法得到的裂縫密度基本都和井?dāng)?shù)據(jù)吻合，但是兩者也有一定的差別。通過常規(guī)方位AVO反演的強(qiáng)各向異性值主要分布在井的左側(cè)，而利用本發(fā)明實(shí)施例求出的裂縫密度的高值區(qū)正對(duì)井位置。相比于常規(guī)方法，本發(fā)明實(shí)施例得到的結(jié)果更合理，并與井資料吻合。圖14表示利用本發(fā)明實(shí)施例求出的裂縫流體指示因子。利用本發(fā)明實(shí)施例反演得到的各向異性參數(shù)εv、δv和裂縫密度Bani，求解得到裂縫流體指示因子。理論研究表明當(dāng)裂縫內(nèi)充填油或水時(shí)，裂縫流體指示因子的值較低，而充填氣時(shí)，裂縫流體指示因子值較高。生產(chǎn)數(shù)據(jù)表明該井是口油井，井累計(jì)產(chǎn)液13400t，產(chǎn)液26.5t/天，產(chǎn)油23.2t/天，裂縫流體指示因子在井位置表現(xiàn)低值，與生產(chǎn)情況吻合，這表明利用本發(fā)明實(shí)施例求得的裂縫流體指示因子可以較好地識(shí)別裂縫內(nèi)裂縫流體類型。通過上述實(shí)施例可以看出，將本發(fā)明實(shí)施例應(yīng)用到碳酸鹽巖儲(chǔ)層中，基于三項(xiàng)Rüger反射系數(shù)方程的分步方位AVO反演方法求出的裂縫參數(shù)與成像測(cè)井結(jié)果吻合，且比基于兩項(xiàng)Rüger反射系數(shù)方程的常規(guī)方位AVO反演結(jié)果更精確，得到的裂縫流體指示因子也和井生產(chǎn)情況吻合，充分說明了本發(fā)明實(shí)施例所采用方法的可行性。關(guān)于上述實(shí)施例，必須說明的是本發(fā)明實(shí)施例中裂縫探測(cè)的適用范圍如下：當(dāng)儲(chǔ)層不發(fā)育裂縫或者發(fā)育多組裂縫(裂縫表現(xiàn)不同的走向)時(shí)，各向異性程度弱，建立在三項(xiàng)Rüger反射系數(shù)方程基礎(chǔ)上的分步方位AVO反演方法無法預(yù)測(cè)出裂縫及其內(nèi)的裂縫流體的分布情況；只有對(duì)發(fā)育單組裂縫(裂縫方位呈現(xiàn)某一優(yōu)勢(shì)方向)的儲(chǔ)層，其表現(xiàn)較強(qiáng)的各向異性時(shí)，該方法才準(zhǔn)確能預(yù)測(cè)其中的裂縫和其內(nèi)的裂縫流體類型。從以上的描述中，可以看出，本發(fā)明實(shí)施例實(shí)現(xiàn)了如下技術(shù)效果：通過基于三項(xiàng)Rüger反射系數(shù)方程的疊前分步方位AVO反演方法，實(shí)現(xiàn)了對(duì)裂縫的高精度識(shí)別，同時(shí)，由于通過該方法可以求解得到各向異性參數(shù)，使用所得到的各向異性參數(shù)可以求解得到裂縫流體指示因子，因而還可以通過裂縫流體指示因子實(shí)現(xiàn)對(duì)裂縫流體類型的定量識(shí)別，解決了目前主要采用的裂縫檢測(cè)法中存在的裂縫識(shí)別精度不高，不能對(duì)裂縫流體進(jìn)行定量識(shí)別的技術(shù)問題；通過對(duì)原本所要處理的形式復(fù)雜的三項(xiàng)Rüger反射系數(shù)方程變形轉(zhuǎn)換，將上述方程轉(zhuǎn)換成相關(guān)的3個(gè)線性方程，實(shí)現(xiàn)了對(duì)所需要數(shù)據(jù)分步、簡(jiǎn)單、快捷的求解；通過采用噪音衰減算法進(jìn)行求解，針對(duì)性地解決了線性反演中的病態(tài)問題，克服了具體實(shí)施時(shí)所使用的常規(guī)地震資料中采集方位范圍窄與噪聲大對(duì)求解過程的不利影響，提高了求解過程的穩(wěn)定性，并增強(qiáng)了解的唯一性；另外，由于本發(fā)明實(shí)施例提出的疊前分步方位AVO反演方法是基于三項(xiàng)Rüger反射系數(shù)方程，使之能夠充分利用大炮檢距地震資料，從而不僅提高了對(duì)地震資料的利用率，也提高了裂縫識(shí)別的準(zhǔn)確性。顯然，本領(lǐng)域的技術(shù)人員應(yīng)該明白，上述的本發(fā)明實(shí)施例的各模塊或各步驟可以用通用的計(jì)算裝置來實(shí)現(xiàn)，它們可以集中在單個(gè)的計(jì)算裝置上，或者分布在多個(gè)計(jì)算裝置所組成的網(wǎng)絡(luò)上，可選地，它們可以用計(jì)算裝置可執(zhí)行的程序代碼來實(shí)現(xiàn)，從而，可以將它們存儲(chǔ)在存儲(chǔ)裝置中由計(jì)算裝置來執(zhí)行，并且在某些情況下，可以以不同于此處的順序執(zhí)行所示出或描述的步驟，或者將它們分別制作成各個(gè)集成電路模塊，或者將它們中的多個(gè)模塊或步驟制作成單個(gè)集成電路模塊來實(shí)現(xiàn)。這樣，本發(fā)明實(shí)施例不限制于任何特定的硬件和軟件結(jié)合。以上所述僅為本發(fā)明的優(yōu)選實(shí)施例而已，并不用于限制本發(fā)明，對(duì)于本領(lǐng)域的技術(shù)人員來說，本發(fā)明實(shí)施例可以有各種更改和變化。凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。當(dāng)前第1頁1 2 3