基于高頻衰減梯度的碳酸鹽巖溶洞儲集體規模識別方法———以塔河油田巖溶古河道溶洞儲集體為例

發布時間：2020-01-14所屬分類：工程師職稱論文瀏覽：1次

摘 要： 摘要：塔河奧陶系碳酸鹽巖油藏儲集空間以溶洞儲集體為主。溶洞儲集體地震反射特征受溶洞規模、溶洞充填物類型及充填程度的共同影響，僅依據地震反射特征難以準確刻畫溶洞規模。從地震、地質及鉆井資料出發，建立了不同規模、不同充填物的溶洞模型，利用非均

　　摘要：塔河奧陶系碳酸鹽巖油藏儲集空間以溶洞儲集體為主。溶洞儲集體地震反射特征受溶洞規模、溶洞充填物類型及充填程度的共同影響，僅依據地震反射特征難以準確刻畫溶洞規模。從地震、地質及鉆井資料出發，建立了不同規模、不同充填物的溶洞模型，利用非均勻介質波動方程正演模擬技術，對比分析了不同規模、不同充填物的溶洞模型對地震反射振幅及頻率響應的影響，提出了基于高頻衰減梯度屬性的溶洞儲集體規模識別技術。以塔河油田巖溶古河道溶洞儲集體為例，利用高頻衰減梯度屬性識別巖溶古河道規模，結果表明，高頻衰減梯度不受溶洞儲集體充填物的物性差異的影響，直觀反映了溶洞發育的規模。基于高頻衰減梯度屬性的巖溶古河道溶洞儲集體規模識別結果與已鉆井結果高度吻合，證明了該方法的有效性。

　　關鍵詞：溶洞儲集體;正演模擬;高頻衰減梯度;巖溶古河道;儲層預測;塔河油田

　　與巖溶作用相關的古洞穴碳酸鹽巖縫洞系統是一種重要的油氣儲集體[1]，塔河油田奧陶系碳酸鹽巖油藏是一種與風化殼古巖溶相關的縫洞型油氣藏[2]，溶洞儲集體是該油藏的主要儲集體，具有儲集空間大，鉆遇概率高的特點。塔河油田已進入全面開發階段，如何提高溶洞儲集體地震預測與油藏描述精度是這一階段亟待解決的問題。碳酸鹽巖溶洞儲集體的外部形狀不規則，內部結構及縫洞連通性復雜，溶洞的發育規模不等以及溶洞內填充物類型多樣等導致了溶洞儲集體地震反射特征復雜，根據地震振幅屬性對溶洞發育規模及內部結構的解釋具有多解性。因此，針對溶洞儲集體的地震預測技術必須建立在正確認識溶洞儲集體地震反射特征的基礎上，利用地震屬性與溶洞儲集體儲層參數之間的關系，實現對溶洞儲集體結構特征的精細描述。

　　正演模擬是研究溶洞儲集體地震響應特征的重要手段。董良國等[3]認為，隨著溶洞直徑的增大，溶洞反射波振幅呈現先增強后減弱再增強最后趨于穩定的變化規律。姚姚[4]基于隨機介質模型及非均質彈性波波動方程，研究了不同寬度和高度溶洞儲集體的地震反射波振幅特征。曲壽利等[5]基于正演模擬定量分析了不同尺度溶洞儲集體的地震反射特征。學者們還采用振幅聚類、不連續性、頻率及疊后波阻抗等地震屬性刻畫了溶洞的發育情況[6-9]。王立華等[10]設計了溶洞物理模型，進行了正演模擬，分析了不同大小溶洞模型的地震響應特征，并在屬性分析時指出利用振幅極值及頻率屬性可較好地識別溶洞尺度，但對于溶洞的形態及充填物等因素未做進一步研究。胡向陽等[11]綜合利用塔河油田鉆井和測井等資料，揭示了溶洞充填物的主要類型及充填模式，并提出了不同充填物的識別方法。

　　頻率屬性是廣泛應用于溶洞規模判斷及充填物識別的屬性，如何利用頻率屬性判斷溶洞規模以及識別充填物一直是我們目前研究的重點與難點。地震剖面的振幅異常綜合反映了地下溶洞儲集體的規模及充填特征，因頻率異常易受振幅異常的影響，基于頻率屬性的溶洞儲集體識別存在多解性或誤區。為了更好地利用頻率屬性判識溶洞儲集體規模及充填情況，有必要基于正演模擬開展不同規模和充填特征的溶洞儲集體振幅及頻率響應特征分析。

　　1正演模擬分析

　　1.1靶區溶洞規模及充填特征統計分析

　　我們依據鉆井的放空、漏失情況，巖心分析、成像測井分析，密度、聲波、中子、伽馬等電測曲線異常特征判斷溶洞是否充填并識別填充物。溶洞測井響應特征通常表現為低密度、低速度、高中子異常，機械沉積充填(砂泥巖)為高伽馬異常。由以上溶洞發育判別標識，分析塔河油田主體工區鉆遇溶洞井溶洞規模及充填特征。圖1顯示了塔河油田主體工區鉆遇溶洞井溶洞厚度統計結果，可以看出超過半數的溶洞厚度小于20m，圖2顯示了塔河油田主體工區鉆遇溶洞井溶洞充填情況統計結果，可以看出溶洞以未充填和半充填為主。本文正演模擬時主要考慮溶洞的充填物類型、充填程度及發育規模3個因素。

　　巖石物理參數是溶洞模型的數據基礎，明確不同充填特征的溶洞模型巖石物理參數是準確建立不同充填特征溶洞模型的前提。圖3及表1均為塔河油田主體工區不同充填特征的溶洞模型縱波速度，溶洞未充填時速度較低，砂泥質充填溶洞時根據充填物類型及充填程度不同，縱波速度變化較大，為2700～4100m/s�？傮w而言，泥質充填溶洞時的縱波速度小于砂質充填溶洞時的縱波速度，即隨著泥質含量的增加，縱波速度逐漸降低。

　　1.2不同規模、不同充填特征的溶洞模型反射特征及頻率響應特征分析

　　常規溶洞儲集體地震波場正演模擬是基于等效介質建立的地質模型，等效介質是對地下介質的一種粗略近似，通常給定溶洞儲集體儲層的平均等效速度，該方法在描述沉積相對穩定的大套地層時有一定的適用性，但對于非均質性較強的碳酸鹽巖儲層，利用該方法建立的等效溶洞介質模型與實際情況存在較大差異，不能真實反映地下介質的非均勻性。本文采用測井約束下的非均質建模方法[12]，在測井數據的約束下，基于地質認識和測井解釋成果，將溶洞儲集體的速度分為兩部分，即背景速度和空間隨機擾動，將二者相加，可構建出縱、橫向具有一定非均質特性的溶洞儲集體模型。采用空間網格與時間步長均可任意變化的高階有限差分聲波正演模擬方法，進行不同規模和不同充填特征的溶洞儲集體地震波場特征研究。

　　1.2.1規模相同，充填特征不同的溶洞儲集體頻率響應特征分析

　　模型組1如圖4所示，該組模型中6個溶洞的橫向寬度均為200m，縱向厚度均為10m，6個溶洞內部充填物巖石物理參數如表2所示。

　　圖5a為信噪比為5時，對模型組1中6個溶洞模型正演模擬得到的Kirchhoff疊前深度偏移剖面，圖5b為根據6個溶洞模型正演模擬的結果拾取的均方根振幅曲線。在溶洞規模相同的情況下，溶洞儲集體內部充填速度越低，其對應“串珠”能量越強。

　　如圖6所示，采用短時傅里葉變換分別獲取6個溶洞模型的單頻能量，可以看出單頻能量的異常與“串珠”能量正相關，即“串珠”能量強，單頻異常能量就強。將頻率為50Hz時的6個溶洞模型橫向比較可知，溶洞①的高頻能量要強于溶洞⑥的高頻能量。所以僅根據高頻能量強度來確定溶洞規模易造成錯誤的地質解釋成果，這種假象是由“串珠”能量差異引起的。

　　1.2.2規模及充填特征均不相同，振幅能量相近的溶洞頻率響應特征分析

　　從模型組1可以看出，溶洞充填物縱波速度越低，地震反射越強。溶洞規模越小，地震反射振幅越弱，故同時調整模型組2中6個溶洞模型的規模和充填物縱波速度使它們振幅一致，并分析其頻率響應特征。模型組2如圖7所示，該組模型為6個寬度為200m，厚度分別為5，10，15，20，25，40m的溶洞，充填速度如表2所示。

　　相關期刊推薦：《石油物探》創刊于1962年，最初是一份原地質礦產部石油地質海洋地質局系統內部交流的學術刊物。所刊登的內容主要為國內地震、測井、電法、磁法和重力等地球物理勘探新理論、新技術和新成果，以地震勘探為主。欄目設有科技論文、勘探開發史例、經驗交流、綜合評述以及學術討論等。

　　圖8a為信噪比為5時，對模型組2中6個溶洞模型正演模擬得到的偏移剖面，圖8b為根據模型組2中6個溶洞模型正演模擬的結果拾取的均方根振幅曲線。隨著溶洞規模的增大，地震反射振幅增大，同時充填物縱波速度的增加又會導致振幅的減弱。在該組模型中，隨著溶洞規模的增大，充填物縱波速度也在增加，這使得6個溶洞所產生的“串珠”能量幾乎相等，溶洞規模和充填物縱波速度這兩個因素對振幅的影響相互抵消，排除了由振幅不同引起頻率響應特征差異的可能。

　　圖9為模型組2中6個振幅相當的溶洞模型的頻率響應特征，可以看出，發育規模不同時，高頻異常對薄溶洞響應特征更為明顯，尤其當頻率大于50Hz時，相較于厚溶洞，薄溶洞頻率異常更明顯;厚溶洞頻率響應特征通常表現為低頻異常，高頻異常不明顯。由此可見，對于“串珠”能量相當的溶洞，可將頻率異常作為溶洞規模判識的依據之一。

　　1.2.3充填特征相同、規模不同的溶洞頻率響應特征分析

　　為了進一步分析頻率屬性與溶洞規模的關系，設計一組溶洞模型(模型組3)，使其充填物速度相同，模型橫向寬度不變，縱向厚度逐漸增加，如圖10所示。該模型組包括6個寬度為200m，縱向厚度分別為5，10，15，20，25，40m的溶洞。溶洞充填物縱波速度范圍為3201～3399m/s，橫波速度范圍1261～1339m/s，密度為1.759～1.906g/cm3。

　　圖11a為信噪比為5時，對模型組3中6個溶洞模型正演模擬得到的偏移剖面，圖11b為根據模型組3中6個溶洞模型正演模擬的結果拾取的均方根振幅曲線。隨著溶洞規模的增大，其振幅也相應地增大，在調諧厚度內，地震反射振幅隨厚度的增加而增加。圖12為模型組3中6個溶洞模型的頻率響應特征，較厚的溶洞主要為低頻異常，較薄的溶洞無高頻異常，這是由能量差異引起的頻率異常的不確定性。

　　1.3高頻衰減梯度屬性

　　對地震道進行譜分析[13-14]，獲得各單頻能量;采用最小二乘方法，對各單頻能量{yf1，yf2，…，yfn｝隨頻率{f1，f2，…，fn｝的變化規律進行線性擬合或對數擬合。單頻能量y隨頻率f的變化規律滿足線性擬合公式y=a×f+b，其中，系數a為單頻能量隨頻率變化的梯度屬性，系數b為擬合關系的截距。

　　分析模型組3中不同規模溶洞的頻率響應特征，可以看出不同規模的溶洞，其頻率變化的梯度特征不同。規模相對較小的溶洞，其高頻衰減梯度絕對值小，而隨著溶洞規模的增加其高頻衰減梯度絕對值越大。模型組1和模型組2的頻率響應特征也遵循該規律。圖13為基于模型組2和模型組3計算得到的高頻衰減梯度隨溶洞厚度變化的情況，可以看出，高頻衰減梯度不受溶洞充填物的影響，而與溶洞規模直接相關。

　　2應用效果

　　研究區巖溶古河道為發育在地層T47頂面以下(一間房組或鷹山組)的巖溶水道，巖溶古河道是塔河油田溶洞儲集體的重要組成部分。優選地震屬性，結合巖溶古河道地震反射特征，可精細刻畫巖溶古河道平面分布特征，利用高頻衰減梯度屬性可進一步確定巖溶古河道的發育規模。如圖14所示，高頻衰減梯度高絕對值處(粉色)對應河道發育規模較大，高頻衰減梯度低絕對值處(藍色)對應河道發育規模較小，圖14a中箭頭處TH12138井鉆遇的溶洞較厚，測井解釋厚度達到50m，其高頻衰減梯度絕對值高，而圖14b中箭頭處TH10112井鉆遇的溶洞較薄，測井解釋厚度為14m，其高頻衰減梯度絕對值低。

　　將工區21口鉆遇巖溶古河道井的高頻衰減梯度與測井解釋得到的溶洞厚度進行交會分析，發現高頻衰減梯度絕對值與河道厚度呈正相關，其中TH12164和TH12112井離擬合曲線較遠(圖15)。原因可能是這兩口井處于巖溶斜坡高陡部位，水力梯度較大，下切深度明顯;河道的縱向厚度大，橫向寬度小，綜合尺度小，從而造成高頻衰減梯度絕對值大。

　　圖16疊合顯示了預測得到的巖溶古河道發育情況與斷裂體系，工區主要發育近東西向及北北東向展布的兩個斷裂，與巖溶古河道的展布特征匹配關系較好，說明巖溶古河道發育與斷裂密不可分，分支河規模較小，河流曲折處、分支河與主體河道交匯處斷裂發育規模較大，并沿河流方向規模逐漸增大。

　　3結論

　　本文針對碳酸鹽巖儲層非均質性強的特點，摒棄基于層狀介質或等效介質建立地震地質模型的傳統方法，利用測井約束下的高精度非均質建模技術，構建了不同尺度和不同充填特征的碳酸鹽巖溶洞儲集體模型，利用波動方程正演模擬技術，研究了不同規模和不同充填物的溶洞儲集體地震反射振幅和頻率屬性，提出了基于高頻衰減梯度的碳酸鹽巖溶洞儲集體規模識別方法。正演模擬結果表明，高頻衰減梯度不受溶洞充填物影響，直觀反映了溶洞規模。基于高頻衰減梯度屬性的塔河油田巖溶古河道溶洞儲集體規模識別表明，隨著溶洞規模的增加，高頻衰減梯度絕對值隨之增大，溶洞儲集體規模識別結果與已鉆井結果高度吻合。