一種預測儲層裂縫發(fā)育參數的方法
【技術領域】
[0001] 本發(fā)明屬于石油地球物理勘探領域，具體涉及一種預測儲層裂縫發(fā)育參數的方 法，屬于地震資料解釋范疇，是一種利用各個方位角地震數據的目的層的雙程走時值來計 算儲層的裂縫發(fā)育參數的方法。
【背景技術】
[0002] 中國南方探區(qū)的碳酸鹽巖裂縫型儲層相對發(fā)育，具體表現(xiàn)為在四川盆地內，鉆遇 了大量的裂縫型氣藏，并獲得工業(yè)氣流。這些裂縫型碳酸鹽巖儲層基本上是受到構造應力 的作用下形成，具有非常強的非均質性，平面上呈不均勻分布，勘探開發(fā)難度很大，如何識 別出裂縫型儲層是地球物理工作的重點之一。所以，對該類型儲層的預測及技術探索具有 很重要的意義。
[0003] 根據研究發(fā)現(xiàn)地下目的層在發(fā)育裂縫的狀態(tài)下，地震P波振幅、雙程走時值通過 該裂縫體時具有隨方位角變化的特點。例如，XiangyangLi在1999年《Geophysics》上 發(fā)表的〈〈FracturedetectionusingazimuthalvariationofP-wavemoveoutfrom orthogonalseismicsurveylines》。研究資料表明：①裂縫體頂層的反射P波雙程走時 不隨方位角變化，振幅隨方位角有微弱的變化，對各向異性反映不明顯。②裂縫體底層的反 射P波表現(xiàn)出很強的方位各向異性，P波反射振幅及雙程走時與測線和裂縫方位有關，測線 與裂縫平行時振幅最強、雙程走時最短；隨著測線與裂縫夾角的增大，振幅逐漸減弱、雙程 走時逐漸變長；至測線與裂縫方向垂直時，振幅最弱，雙程走時最長。由此可見，P波在通過 裂縫體時表現(xiàn)出很強的方位各向異性特征，裂縫對P波響應的影響程度取決于裂縫方位與 測線之間的夾角及地下一定開啟裂縫系統(tǒng)的規(guī)模，該現(xiàn)象尤其是對底層的走時影響較大。
[0004] 現(xiàn)有的常規(guī)裂縫型儲層預測的技術方法，主要是與疊后及疊前有關的相關技術。 如相干技術、曲率、P波各向異性分析技術、道集時差計算等來實現(xiàn)預測。一些專利如發(fā)明 《裂縫預測方法和裝置》（【申請?zhí)枴?01010205983.4)的專利公開了一種裂縫預測方法和裝 置，主要是以目的層為拾取中心，利用拾取時窗獲取每個地震道的反射振幅；利用獲取的方 位角和反射振幅進行橢圓擬合，以確定裂縫方向和裂縫密度；發(fā)明《一種儲層裂縫的確定方 法》（申請?zhí)?：201210242061.X)的專利公開了一種儲層裂縫的確定方法，包括：根據實際地 震資料類型對實際地震資料進行多種方法裂縫預測分析，獲取儲層的裂縫屬性及分布；利 用實際地質參數進行正演及FMI資料，比對儲層裂縫分布，獲取匹配度及吻合度；根據匹配 度及吻合度從多種裂縫預測分析方法中優(yōu)選選取一種最佳裂縫分析方法。根據優(yōu)選出的裂 縫分析方法，解釋實際地震資料，確定儲層裂縫。但是，這些常規(guī)的技術方法實現(xiàn)起來難度 相當大，操作起來易受各種不利因素影響，主要存在如下幾個問題：
[0005] (1)利用井資料進行裂縫預測的約束，在遠離井的區(qū)域預測效果往往不好。
[0006] (2)利用振幅、頻率等屬性進行橢圓擬合計算裂縫發(fā)育方向及強度，預測結果往往 差異很大，難于對這些技術方法進行取舍及影響成果的判斷。
[0007] (3)對微型裂縫發(fā)育參數的評價，常規(guī)技術方法往往難于準確描述。

【發(fā)明內容】

[0008] 本發(fā)明的目的在于解決上述現(xiàn)有技術中存在的難題，提供一種預測儲層裂縫發(fā)育 參數的方法，采用相對簡單的技術流程來解決對裂縫型儲層的識別，更為有效地對該類型 儲層進行預測。
[0009] 本發(fā)明是通過以下技術方案實現(xiàn)的：
[0010] 一種預測儲層裂縫發(fā)育參數的方法，包括：
[0011] (1)輸入地震數據；
[0012] (2)利用不同方位角的地震數據提取目的層雙程走時值，并進行歸一化處理，得到 經歸一化處理后的各個目的層的雙程走時值；
[0013] (3)對設定網格上各點分別生成地震玫瑰圖和綜合裂縫發(fā)育強度平面圖，（一方 面是經橢圓擬合計算形成每個網格點上的地震玫瑰圖，另一方面是利用網格點上所計算的 離散系數值，再經插值計算后得到綜合裂縫發(fā)育強度平面圖）將此兩圖疊合后進行顯示， 獲得裂縫參數平面圖。
[0014] 所述步驟（1)中輸入的地震數據為經動校正后形成的⑶P道集數據。
[0015] 所述步驟⑵包括以下步驟：
[0016](21)、根據地震資料及實際情況，對處理后的⑶P道集數據進行分方位角、入射角 范圍限定的疊加及偏移處理，得到各個分方位角偏移數據體；
[0017](22)、根據地震資料的實際情況及目的層的反射特征，對目的層使用相同的解釋 參數對各個分方位角偏移數據體進行拾取解釋，得到各個方位角數據體的目的層位的雙程 走時值，并對拾取得到的各個方位角數據體的目的層位的雙程走時值進行網格化插值和平 滑濾波處理，得到各個方位角數據體的處理后的目的層雙程走時數據。
[0018] 所述步驟（22)中的網格化插值和平滑濾波處理具體如下：
[0019] 通過利用測井資料及各個方位角偏移數據體對井上的裂縫型儲層進行準確標定， 并進行目的層的選取。目的層選取原則上是針對裂縫型儲層段的底部進行選取。并確定好 目的層的反射特征，如波峰、波谷或零相位等這些容易進行層位追蹤的反射特征。其中，解 釋后的層位網格化插值、平滑濾波等處理是將解釋網格（通常解釋網格大于1線Xl道）的 層位的雙程走時值根據網格化插值算法內插成1線Xl道，并對網格化計算后的層位雙程走 時數據設定一定的濾波參數進行平滑濾波，使其消除某些走時值突變點，更接近實際的目 的層的解釋。網格化插值算法可有最小平方法、加權平均法、（泛）克里金法等，使用這些 計算方法應根據解釋層位具體情況而定。
[0020](23)、對所述各個方位角數據體的處理后的目的層雙程走時數據進行歸一化處理 得到經歸一化處理后的各個目的層的雙程走時值。具體如下：
[0021] 根據各個方位角數據體的目的層雙程走時的最大值及最小值的實際情況，統(tǒng)一將 各個方位角走時數值歸一化校正到相對較小的n~m值域上。歸一化處理亦即是利用加、 減、乘、除等組合運算，使得雙程走時值均有一定程度上的變小并使各點上的最大雙程走時 值與最小雙程走時值的倍數相對增大，以利于進行橢圓擬合計算。
[0022] 所述步驟（3)包括以下步驟：
[0023] (31)、根據實際情況設置設定網格的網格參數；所述網格參數包括網格間距和網 格數；
[0024](32)、提取設定網格點上所述經歸一化處理后的各個目的層的雙程走時值，并對 該點進行離散系數計算；將該點上的離散系數作為該點綜合裂縫發(fā)育強度參數；
[0025] (33)、計算設定網格點上的裂縫走向及該走向的裂縫相對強度，得到地震玫瑰 圖；
[0026] (34)、利用設定網格點上的離散系數的值經插值計算后生成綜合裂縫發(fā)育強度平 面圖；
[0027] (35)、將所述綜合裂縫發(fā)育強度平面圖和所述地震玫瑰圖疊合在一張平面圖上進 行顯示，從而獲得裂縫參數平面圖。
[0028] 所述步驟（33)中是采用三點橢圓擬合來進行計算的，具體如下：
[0029]