縫洞型儲層聲波遠探測測井響應(yīng)模擬及應(yīng)用

摘要： 塔河油田奧陶系碳酸鹽巖縫洞型儲層縫洞類型多樣、結(jié)構(gòu)復(fù)雜。聲波遠探測技術(shù)能夠探測到井外數(shù)十米范圍內(nèi)的地質(zhì)異常體，在這種縫洞型油藏評價中具有廣闊的應(yīng)用前景。本文首先采用完全匹配層（perfectly matched layer， PML）吸收邊界條件的交錯網(wǎng)格有限差分方法模擬塔河油田縫洞型儲層的聲波遠探測測井響應(yīng)；然后根據(jù)縫洞空間分布建立過井及井旁裂縫不同長度、開度、角度的裂縫分布模型以及井旁洞穴、孔洞尺寸、數(shù)量、遠近變化的縫洞型儲層模型；再用有限差分方法模擬波場，分析裂縫、洞穴、溶蝕孔洞的波場特征，用偏移成像處理方法得到遠探測聲波波場裂縫、洞穴、溶蝕孔洞的成像結(jié)果，建立相應(yīng)的聲波遠探測成像解釋模型；最后將實鉆井聲波遠探測測井資料與常規(guī)測井、電成像測井和地震剖面特征及試油成果進行對比。結(jié)果表明，應(yīng)用此方法建立的聲波遠探測解釋模型對碳酸鹽巖縫洞型儲層的判識效果良好，適用于塔河油田碳酸鹽巖縫洞型儲層的識別與評價。

關(guān)鍵詞：塔河油田；縫洞型儲層；聲波遠探測；數(shù)值模擬；解釋模型；實例應(yīng)用

doi：10.13278/j.cnki.jjuese.20230195 中圖分類號：P631.8 文獻標(biāo)志碼：A

收稿日期： 2023-08-09

作者簡介： 藺學(xué)旻（1968—），男，高級工程師，從事測井綜合研究，E-mail：lin3d@189.cn

基金項目： 中國石化科技攻關(guān)項目（P20039）

Supported by Science and Technology Tackle Project of SINOPEC （P20039）

Simulation and Application of Acoustic Remote Detection Logging Response in Fracture-Cavity Reservoirs： Taking Fracture-Cavity" Carbonate Reservoirs in Tahe Oilfield as an Example

Lin Xuemin， Xiao Honglin

Research Institute of Exploration and Development， Sinopec Northwest Oilfield Company， Urumqi 830011， China

Abstract： Fracture-cavity reservoirs of Ordovician carbonate rock in Tahe oilfield have various types and complex structures. Acoustic remote detection technology can detect geological anomalies in the range of tens of meters outside the well， and has broad application prospects in the evaluation of fracture-cavity reservoirs. Firstly， the staggered grid finite difference method with perfectly matched layer （PML） absorption boundary conditions is used to simulate the acoustic remote detection logging response of fracture-cavity reservoirs in Tahe oilfield. Then， through the simplified modeling of fracture-cavity spatial distribution， the fracture distribution model with different lengths， openings and angles of" fractures passing through and beside the well and the fracture-cavity reservoir model with the changes of cave， hole size， number and distance beside the well are established. Then the wave field is simulated by finite difference， the wave field characteristics of cracks， caves and dissolution holes are analyzed， and the imaging results of cracks， caves and dissolution holes of acoustic remote detection" wave field are obtained by offset imaging processing method， and the corresponding imaging interpretation model of acoustic remote detection is established. At last， the acoustic remote detection logging data of real drilling are compared with conventional logging， electric imaging logging， seismic profile characteristics， and oil test results. The results show that the acoustic remote detection interpretation model established by this method is effective for the identification and evaluation of carbonate fracture-cavity reservoirs.

Key words： Tahe oil field; fracture-cavity reservoirs; acoustic remote detection; numerical simulation; interpretation model; example application

0 前言

塔河油田位于塔里木盆地北部沙雅隆起阿克庫勒凸起南部，是我國投入開發(fā)的儲量最大、產(chǎn)量最高的海相碳酸鹽巖油藏[1-4]，已累計產(chǎn)原油突破1×108 t[5]，其奧陶系碳酸鹽巖先后經(jīng)歷多期構(gòu)造抬升和海平面變化，在加里東晚期—海西期經(jīng)歷了長期風(fēng)化巖溶作用[6]，形成了復(fù)雜多樣的巖溶縫-洞體系，包括大型溶洞、溶蝕孔洞、溶蝕孔隙、構(gòu)造溶蝕裂縫等，其中大型洞、孔是最主要的儲集空間；不同尺度的斷裂和裂縫是主要的連通通道[7]。近年來塔河油田碳酸鹽巖油藏的精細開發(fā)需要對縫洞型儲層進行高精度刻畫，常規(guī)測井由于受徑向探測深度的限制，無法描述遠離井筒縫洞型儲層的結(jié)構(gòu)，因此需要獲得對井周數(shù)十米范圍內(nèi)響應(yīng)信息的描述，搭建測井與地震識別縫洞的橋梁。

近年發(fā)展起來的利用充液井中偶極子聲源進行的遠探測測井[8-14]，從不定向發(fā)展到定向，探測距離也發(fā)展到可以對井外近百米內(nèi)的地質(zhì)異常體進行探測[15]，填補了常規(guī)測井與地震之間的空白，在塔河油田碳酸鹽巖縫洞型油藏評價中具有廣闊的應(yīng)用前景。但縫洞型儲層復(fù)雜多變，聲波遠探測成像特征與儲層發(fā)育類型及發(fā)育程度的對應(yīng)關(guān)系尚不明確，制約了常規(guī)成像測井近井筒儲層評價與聲波遠探測遠井筒縫洞識別及地震資料區(qū)域性縫洞型儲層刻畫的結(jié)合。

本文利用有限差分模擬方法[16]，通過開展縫洞型儲層遠探測特征的正演模擬，建立了過井及井旁裂縫不同長度、開度、傾角的成像特征和井旁洞穴、孔洞尺寸、數(shù)量、遠近變化的成像特征，形成了塔河油田縫洞型儲層的聲波遠探測成像解釋模型，并通過在塔河油田奧陶系實鉆井聲波遠探測測井的實例，將聲波遠探測成像結(jié)果與電成像測井、常規(guī)測井及地震資料縫洞識別結(jié)果進行對比，為建立塔河油田成像測井、常規(guī)測井、聲波遠探測井到地震的全尺度縫洞型儲層評價體系提供依據(jù)。

1 孔隙彈性介質(zhì)陣列聲波測井響應(yīng)有限差分模擬方法

1.1 Biot孔隙彈性介質(zhì)理論及基本方程

利用Biot孔隙彈性介質(zhì)模型近似縫洞型儲層。孔隙彈性介質(zhì)的彈性波傳播理論由Biot[17-18]提出，是一種描述多孔介質(zhì)在受到機械作用時孔隙流體流動行為的數(shù)學(xué)模型，考慮了孔隙流體和固體基質(zhì)之間相對運動引起的能量耗散。它將多孔介質(zhì)的變形描述為固體基質(zhì)的彈性變形和飽和孔隙流體流動引起的變形兩個分量的總和，并考慮了固相和液相之間的耦合及流體壓縮性和流體飽和度變化的影響。圖1給出了Biot孔隙彈性介質(zhì)模型示意圖。

為了處理問題方便， Biot提出了五個基本假設(shè)：1）Biot理論中的“孔隙”為有效孔隙，即液相是連續(xù)的，而密閉孔隙被認為是固相的一部分；2）孔隙在空間中隨機分布，并且是各向同性的，即孔隙占據(jù)面積與固相占據(jù)面積的比率與橫截面的方向無關(guān)；3）孔隙的尺度為宏觀尺度；4）固體骨架被認為具有壓縮性和剪切剛性，即為彈性固體基質(zhì)；5）介質(zhì)單位體積的形變是完全可逆的。基于這些假設(shè)，Biot[17-18]根據(jù)固相與液相的質(zhì)量和動量守恒方程，通過流體流動和力學(xué)邊界條件以及描述多孔介質(zhì)流體流動特性的本構(gòu)關(guān)系的耦合，導(dǎo)出了Biot孔隙彈性介質(zhì)波動方程，用以近似模擬孔隙彈性介質(zhì)中彈性波的傳播特征與傳播規(guī)律。

圖2a給出了軸向柱對稱的井筒流體與地層模型[19-20]。設(shè)極坐標(biāo)下Biot固-液雙相模型中固相質(zhì)點位移矢量為u=（ur，uθ，uz）T（ur、uθ、uz分別為位移矢量u

的徑向、軸向和垂向分量），液相質(zhì)點位移矢量為U=（Ur，Uθ，Uz）T（Ur、Uθ、Uz分別為位移矢量U的徑向、軸向和垂向分量）。根據(jù)廣義胡克定律，考慮平行于波傳播方向流體的流動， Biot理論固-液雙相橫向各向異性孔隙彈性介質(zhì)中的應(yīng)力應(yīng)變和孔隙壓力關(guān)系可以表示為如下本構(gòu)方程[21-22]：

式中：Δr、Δz分別為離散r方向和z方向的步長；vPmax為介質(zhì)中的最大縱波速度；vS為橫波速度；fc為激發(fā)源主頻。

為提高計算效率和精度，在差分方程求解時，還選用了非分裂PML（NPML）吸收邊界條件。圖3給出了不同吸收邊界條件計算精度的對比，可以看出，與無吸收邊界條件相比，NPML吸收邊界條件和阻尼吸收邊界條件均使邊界波場得到了較大的衰減，證實了NPML和阻尼吸收邊界條件的有效性。圖3b給出了圖3a中紅色橢圓內(nèi)NPML和阻尼吸收邊界條件波場放大的對比，可以看出，阻尼吸收邊界條件的波場幅度更大，表明NPML吸收邊界條件對波場的吸收效果好于阻尼吸收邊界條件。

2 縫洞型儲層偶極聲場響應(yīng)

根據(jù)塔河油田實鉆井資料井眼及泥漿數(shù)據(jù)，本次數(shù)值模擬采用井筒直徑為0.2 m，井孔泥漿密度為1.1 g/cm3，泥漿縱波速度為1 500 m/s。從取心實驗數(shù)據(jù)和常規(guī)測井資料得到碳酸鹽巖基質(zhì)地層的縱波速度為 5 800 m/s，橫波速度為3 600 m/s，密度為 2.69 g/cm3。依據(jù)塔河油田數(shù)百口測井常規(guī)測井中的密度數(shù)據(jù)、偶極子聲波中的縱橫波資料，在數(shù)值模擬中對塔河油田縫洞型儲層的縱波速度取值為4 500 m/s，橫波速度為2 750 m/s，密度為2.55 g/cm3。

2.1 單個過井裂縫延伸長度、開度和傾角變化時的波場響應(yīng)

塔河油田不同尺度的斷裂和沿裂縫溶蝕擴大的大型裂縫是重要的油氣儲集空間和通道。根據(jù)塔河油田數(shù)百口井的常規(guī)測井、電阻率掃描成像測井資料可知，這些斷裂和溶蝕擴大裂縫開度為數(shù)厘米到數(shù)十厘米，本次模擬開度選定為20 cm，設(shè)定傾角為60°，走向為0°，背景介質(zhì)為灰?guī)r，模擬分析穿井筒的裂縫長度對偶極反射聲場的影響。全裂縫長度為1、5、10、15 m。圖4中每一個裂縫長度左側(cè)為模擬反射波的變密度圖，右側(cè)為由模擬反射波數(shù)據(jù)處理得到的成像結(jié)果的變密度圖（下同，不再贅述）。從圖4中可以看出，對于過井傾斜裂縫而言，隨著裂縫長度的增加，成像結(jié)果越來越清晰。由此可知，裂縫的長度越長，在偶極遠探測測井中越容易被探測到。

假定裂縫長度為10 m，傾角為60°，走向為0°，背景介質(zhì)為灰?guī)r，模擬分析裂縫開度對偶極反射聲場的影響。根據(jù)塔河油田常規(guī)及電成像測井反映的沿裂縫溶蝕形成的大裂縫和不同尺度斷裂的開度，本次模擬開度選為5、10、20、30 cm。圖5是每一個裂縫開度變化的模擬成像圖。從圖5中可以看出，裂縫的開度越大，成像結(jié)果越清晰，但沒有本質(zhì)差別。

假定裂縫長度為10 m，開度為20 cm，走向為0°，背景介質(zhì)為灰?guī)r，模擬分析裂縫傾角對偶極反射聲場的影響。令裂縫傾角為20°～80°，每隔20°模擬一次。裂縫離井距離（裂縫中點到井軸之間的距離）為10 m。圖6為過井筒裂縫傾角變化模擬成像圖。從圖6中可以看出，隨著裂縫傾角的增大，裂縫的成像長度越來越大，也越來越清晰。這是因為：裂縫與井平行時輻射到地層中的橫波能量最多，井外反射體最容易成像；隨著裂縫傾角變小，聲源輻射到地層中的橫波的能量逐漸減小，裂縫的成像長度和清晰度也逐漸變差。

2.2 多個井旁裂縫延伸長度、開度和傾角變化時的波場響應(yīng)

塔河油田碳酸鹽巖油藏的儲層中除了不同尺度斷裂和沿裂縫溶蝕形成的大裂縫之外，發(fā)育更多的是尺寸較小的裂縫帶，研究裂縫帶的遠探測聲場特征至關(guān)重要。

根據(jù)塔河油田裂縫開度分布及聲波遠探測的分辨率，假定裂縫開度為5 cm，傾角為60°，走向為0°，裂縫長度分別為0.5、1.5、5.0、10.0 m，背景介質(zhì)為灰?guī)r，模擬分析井旁裂縫組裂縫長度對偶極反射聲場的影響，結(jié)果如圖7所示。從圖7中可以看出，由于生成裂縫組的位置有隨機性，裂縫組之間的相互干涉不同，成像有一定的變化；但總體來看，裂縫長度越長越易成像，對于小尺寸的裂縫組，最終成像時呈現(xiàn)一種團簇狀。

根據(jù)塔河油田電成像顯示的裂縫開度范圍并考慮到陣列聲波儀器分辨率，設(shè)定裂縫開度為1、2、5、10 cm，傾角為60°，走向為0°，裂縫長度為0.8 m，背景介質(zhì)為灰?guī)r，模擬分析井旁裂縫組裂縫開度對偶極反射聲場的影響，結(jié)果如圖8所示。從圖8中可以看出，由于生成的裂縫組的位置有隨機性，裂縫組之間的相互干涉成像反映不同，總體來看，裂縫開度越大反射波能量越強，成像越明顯，成像結(jié)果反映的是裂縫組的整體形態(tài)。

假定裂縫長度為0.8 m，開度為5 cm，走向為0°，傾角為0°、30°、60°、90°，背景介質(zhì)為灰?guī)r，模擬分析井旁裂縫組裂縫傾角對偶極反射聲場的影響，結(jié)果如圖9所示。從圖9中可以看出，隨著平行裂縫組傾角的增加，裂縫組的整體形態(tài)更為清晰，由于裂縫組的尺度較小，最終成像時呈現(xiàn)一種團簇狀。這是由于裂縫組的反射波之間相互疊加干涉，不能夠完全還原模型圖，但能夠看出裂縫組的整體區(qū)域輪廓；且由于被前面的裂縫組遮擋，后面的裂縫組的反射能量較弱。

2.3 井旁單個洞穴直徑、中心距井軸徑向距離變化時的波場響應(yīng)

在塔河油田碳酸鹽巖縫洞型油藏中，洞穴是極為重要的一種儲層類型，它儲存油氣的貢獻度超過了50%。為了進一步了解洞穴產(chǎn)生的遠探測聲場特征，需要從洞穴尺寸、離井距離方面出發(fā)，模擬分析洞穴遠探測聲場隨上述因素的變化規(guī)律。

根據(jù)塔河油田常規(guī)及成像測井對264個井層洞穴進行統(tǒng)計而確定的洞穴尺寸范圍，本次模擬假定洞穴離井距離為10 m，洞穴直徑為10、50、100、300 cm，背景介質(zhì)為灰?guī)r，模擬分析洞穴尺寸對偶極反射聲場的影響，結(jié)果如圖10所示。從圖10中可以看出：當(dāng)單個洞穴的直徑較小時，成像結(jié)果呈斑點狀；當(dāng)單個洞穴的直徑較大時，可呈現(xiàn)洞穴的輪廓。

假定洞穴直徑為3 m，離井距離為10、20、40、60 m，背景介質(zhì)為灰?guī)r，模擬分析洞穴離井筒距離變化對偶極反射聲場的影響，結(jié)果如圖11所示。從圖11中可以看出，洞穴離井距離越遠，洞穴的輪廓成像越模糊。這是由于隨著距離的增加，反射波的能量越來越弱。長距離的反射體成像與儀器輻射出去的能量有直接關(guān)系。

2.4 井旁多個孔洞直徑、中心距井軸徑向距離變化時的波場響應(yīng)

塔河油田碳酸鹽巖儲層中的溶蝕孔洞往往不是單個存在，而是多個孔洞發(fā)育成群；所以需要從孔洞的尺度、數(shù)目方面出發(fā)，詳細模擬分析它們對孔洞群遠探測聲場的影響。

根據(jù)塔河油田巖心和電成像資料反映的孔洞尺寸，設(shè)定孔洞直徑為1、2、5、8 cm，孔洞數(shù)目為50個，背景介質(zhì)為灰?guī)r，模擬分析孔洞尺寸對偶極反射聲場的影響，結(jié)果如圖12所示。從圖12中可以看出，孔洞的聲波遠探測成像為“小串珠”的形態(tài)，隨著孔洞尺度的增加，反射波的能量變大，“小串珠”現(xiàn)象越來越明顯。

假定孔洞直徑為5 cm，孔洞數(shù)目為50、100、200、300個，背景介質(zhì)為灰?guī)r，模擬分析孔洞數(shù)目變化對偶極反射聲場的影響，結(jié)果如圖13所示。從圖13中可以看出，隨著孔洞數(shù)目的增加，反射波的能量增大，成像結(jié)果呈現(xiàn)團簇狀，可分辨孔洞組的區(qū)域，但孔洞組內(nèi)孔洞的形態(tài)無法識別。

3 儲層聲場數(shù)值模擬結(jié)果特征分析及實測資料解釋應(yīng)用

3.1 遠探測聲場數(shù)值模擬特征模式

從前文裂縫、洞穴、孔洞模型的模擬結(jié)果分析可以得到如下規(guī)律性認識。

1）對于裂縫型儲層，基本呈現(xiàn)為連續(xù)線性反射。不同的裂縫長度、開度、傾角、離井距離等具有不同的成像結(jié)果，裂縫長度大、開度大、傾角高、離井距離近，成像結(jié)果清晰。

2）對于洞穴型儲層：當(dāng)單個洞穴的直徑較小時，成像結(jié)果呈斑點狀；當(dāng)單個洞穴的直徑較大時，可呈現(xiàn)洞穴的輪廓。洞穴直徑、洞穴離井距離對成像結(jié)果均有影響。

3）孔洞型儲層分別成像為小串珠、團簇狀雜亂斑點。孔洞尺寸、數(shù)目對成像結(jié)果均有影響。

據(jù)此，可以總結(jié)出儲層遠探測聲場的解釋模式及對應(yīng)的識別圖版，以便于解釋塔河油田碳酸鹽巖縫洞型儲層的聲波遠探測反射聲場顯示特征，推測儲層裂縫與溶蝕孔洞的發(fā)育特征。表1給出了由裂縫、洞穴及孔洞模型的聲場模擬響應(yīng)總結(jié)得到的碳酸鹽巖縫洞型儲層遠探測特征解釋模式。表1中圖版類似于傳統(tǒng)的地球物理解釋圖版，可直接據(jù)以推斷對應(yīng)的縫洞分布特征。

3.2 聲場數(shù)值模擬圖版在塔河油田的應(yīng)用

在上述聲波遠探測響應(yīng)特征數(shù)值模擬結(jié)果及解釋模式的基礎(chǔ)上，結(jié)合塔河油田實際測井資料進行了聲波遠探測處理解釋，并結(jié)合常規(guī)綜合測井、電成像測井及地震資料進行對比。

3.2.1 TH124XX井

圖14所示給出了TH124XX井碳酸鹽巖井段常規(guī)測井、電成像測井、聲波遠探測測井及地震綜合圖。該井6 415～6 457 m常規(guī)測井顯示雙側(cè)向電阻率有所降低，聲波時差和中子孔隙度同步增大，表現(xiàn)為孔洞、裂縫型儲層特征（圖14a），在電成像上顯示發(fā)育數(shù)條高導(dǎo)縫，微裂縫為主，溶蝕較發(fā)育，孔隙頻譜顯示孔隙較發(fā)育，總孔隙度為3.7%～4.2%（圖14b），聲波遠探測成像（圖14c）顯示井外65 m范圍內(nèi)見斑點狀反射異常，對照前述解釋模型，此井段為溶蝕孔洞發(fā)育，常規(guī)測井、電成像的儲層識別與聲波遠探測測井一致。6 440～6 495 m電阻率中高（圖14a），儲層不發(fā)育，在電成像上顯示為個別細微裂縫（圖14b），但聲波遠探測成像顯示井外40～60 m范圍有兩組明顯帶狀反射異常（圖14c），對照前述解釋模型為井外發(fā)育裂縫帶，裂縫延伸較遠、規(guī)模較大，45°切片最清晰，裂縫走向沿北東南西向展布（圖14c）。6 537～6 630 m在電成像顯示裂縫不發(fā)育，僅發(fā)育斑點狀溶蝕孔洞，孔隙頻譜顯示儲層非均質(zhì)性較強，總孔隙度為3.0%～3.3%（圖14b），聲波遠探測成像顯示井外60 m范圍見多組近似垂直及斑點狀反射異常（圖14c），根據(jù)前述解釋模型，此段為發(fā)育的溶蝕孔洞帶。6 643～6 657 m電成像顯示高導(dǎo)縫特征，裂縫為中角度，上段裂縫寬度小，主要為微裂縫，下段裂縫寬度大、縫面粗糙，被溶蝕作用改造明顯（圖14b），聲波遠探測成像顯示井外約60 m近似垂直的帶狀反射異常（圖14c），根據(jù)前述解釋模型，此段為裂縫發(fā)育帶，45°切片最清晰，推測裂縫走向沿北東向展布（圖14c）。地震資料（圖14d、e）顯示“串珠”規(guī)模小，相干弱，與常規(guī)測井、電成像測井和聲波遠探測測井解釋本井主要為裂縫、孔洞型儲層一致。本井全井酸壓后12 mm油嘴電潛泵生產(chǎn)，油壓1.59 MPa，日產(chǎn)油24 m3。

3.2.2 TH122XX井

圖15給出了TH122XX井常規(guī)測井、聲波遠探測測井及地震綜合圖。6 424.5～6 425.5 m電阻率急劇降低，聲波時差稍有增大，為單組系裂縫發(fā)育（圖15a）。6 425.5～6 454.5 m電阻率逐漸降低，聲波時差局部增大，為洞穴頂部連續(xù)的裂縫發(fā)育帶（圖15a）；此段聲波遠探測成像顯示為與數(shù)值模擬一致的穿井壁的井旁裂縫帶特征（圖15 b）。6 458.5～6 463.0 m雙側(cè)向電阻率急劇降低，井徑擴大超過14 in，聲波時差急劇增大，表現(xiàn)為放空洞穴的曲線特征（圖15a），聲波遠探測成像顯示為徑向深度為20～40 m、高度為20余m的微弧狀反射特征，為洞穴型儲層的成像特征（圖15b），與常規(guī)測井反映一致。從6 485 m到井底電阻率略有降低，聲波時差微幅增大，為孔洞型儲層特征（圖15a）；聲波遠探測成像顯示為從井筒周圍延伸到徑向深度40 m范圍

內(nèi)的團簇狀雜亂斑，為溶蝕孔洞發(fā)育帶（圖15b）。圖15 c、d地震資料顯示“串珠狀”明顯，儲層發(fā)育，與常規(guī)測井及聲波遠探測解釋的洞穴相吻合。本井全井段酸壓，5 mm油嘴開井，油壓8.2 MPa，日產(chǎn)油45 m3。

3.2.3 TH1021XX井

圖16給出了TH1021XX井常規(guī)測井、聲波遠探測測井及地震綜合圖。常規(guī)測井6 025～6 085 m顯示儲層不發(fā)育，且以孔洞型儲層為主（圖16 a）；聲波遠探測成像顯示距離井筒20～40 m見強反射近弧狀成像特征信號（圖16 b）。由于常規(guī)測井探測深度有限，無法反映遠離井筒的這個井旁儲集體，根據(jù)模擬成像特征分析為井旁洞穴，近弧狀成像為洞穴邊界。偶極聲波遠探測成像在東西向反射最強烈，成像最清晰（圖16 b）。地震資料顯示本井井旁有東西向斷裂發(fā)育（圖16 c、d），聲波遠探測成像顯示的是此斷裂破碎并溶蝕形成的溶蝕洞穴邊界。本井全井段酸壓，4.5 mm油嘴，油壓17 Mpa，日產(chǎn)油78 m3。

4 結(jié)論

1）本文根據(jù)陣列聲波測井響應(yīng)有限差分模擬的原理和方法，對塔河油田碳酸鹽巖縫洞型儲層的聲波遠探測聲場進行了數(shù)值模擬，得到了過井及井旁裂縫不同長度、開度、角度的成像特征和井旁洞穴、孔洞尺寸、數(shù)量、遠近變化的成像特征，建立了塔河油田縫洞型儲層聲波遠探測的基礎(chǔ)解釋模板。

2）塔河油田三口井實測陣列聲波資料的聲波遠探測測井與常規(guī)測井、電成像測井和地震剖面特征的對比證明，聲波遠探測測井對碳酸鹽巖縫洞型儲層的判識具有較好效果，其結(jié)果與電成像測井、常規(guī)測井和地震屬性剖面結(jié)果對應(yīng)性較好。

3）[JP1]通過聲波遠探測數(shù)值模擬建立的塔河油田縫洞型儲層聲波反射聲場特征、解釋模式與相應(yīng)特征圖像，為識別與評價塔河油田縫洞型油藏遠離井筒的儲層發(fā)育類型與發(fā)育程度及延展范圍奠定了基礎(chǔ)。通過微尺度成像測井、小尺度常規(guī)測井和中尺度聲波遠探測測井與大尺度地震資料的有機結(jié)合，未來有望建立起塔河油田從成像測井、常規(guī)測井、聲波遠探測測井到地震資料的全尺度縫洞型儲層評價體系。

參考文獻（References）：

[1] 康玉柱，康志宏. 塔河大油田勘探實踐與技術(shù)創(chuàng)新[M]. 北京： 中國石化出版社，2009.

Kang Yuzhu，Kang Zhihong. Exploration Practice in Tahe Huge Oilfield and Technical Innovation[M]. Beijing： China Petrochemical Press， 2009.

[2] 李陽，范智慧. 塔河奧陶系碳酸鹽巖油藏縫洞系統(tǒng)發(fā)育模式與分布規(guī)律[J]. 石油學(xué)報，2011，32（1）：101-106.

Li Yang，F(xiàn)an Zhihui. Developmental Pattern and Distribution Rule of the Fracture-Cavity System of Ordovician Carbonate Reservoirs in the Tahe Oilfield[J]. Acta Petrolei Sinica，2011，32（1）：101-106.

[3] 翟曉先，云露.塔里木盆地塔河大型油田地質(zhì)特征及勘探思路回顧[J]. 石油與天然氣地質(zhì)，2008，29（5）：565-573.

Zhai Xiaoxian，Yun Lu. Geology of Giant Tahe Oilfield and a Review of Exploration Thinking in the Tarim Basin[J]. Oil ＆ Gas Geology，2008，29（ 5） ： 565-573.

[4] 竇之林. 碳酸鹽巖縫洞型油藏描述與儲量計算[J]. 石油實驗地質(zhì)，2014，36（ 1） ： 9-15.

Dou Zhilin. Description and Reserves Calculation of Fractured-Vuggy Carbonate Reservoirs[J]. Petroleum Geology ＆ Experiment，2014，36（1）：9-15.

[5] 郭春濤，史江濤，劉亮，等. 塔里木盆地塔河地區(qū)中下奧陶統(tǒng)沉積特征及其演化模式[J]. 吉林大學(xué)學(xué)報（地球科學(xué)版），2024，54（1）：68-82.

Guo Chuntao，Shi Jiangtao，Liu Liang， et al. Sedimentary Characteristics and Evolution Model of Middle and Lower Ordovician in Tahe Area，Tarim Basinn[J]. Journal of Jilin University （Earth Science Edition），2024，54（1）：68-82.

[6] 何治亮，彭守濤，張濤. 塔里木盆地塔河地區(qū)奧陶系儲層形成的控制因素與復(fù)合-聯(lián)合成因機制[J]. 石油與天然氣地質(zhì)，2010，31（6）：743-752.

He Zhiliang， Peng Shoutao， Zhang Tao. Controlling Factors of Ordovician Reservoir Formation in Tahe Area， Tarim Basin and Its Composite and Joint Genetic Mechanism[J]. Oil and Gas Geology， 2010，31（6）：743-752.

[7] 竇之林. 塔河油田碳酸鹽巖縫洞型油藏開發(fā)技術(shù)[M]. 北京：石油工業(yè)出版社，2012.

Dou Zhilin. Development Technology of Carbonate Fracture Cave Type Reservoir in Tahe Oilfield[M]. Beijing： Petroleum Industry Press，2012.

[8] Tang X M， Glassman H， Patterson D. Single-Well Acoustic Imaging in Anisotropic Formations[J]. Geophysics， 2008， 73（4）： D11-D16.

[9] Tang X M， Patterson D. Single-Well S-Wave Imaging Using Multi-Component Dipole Acoustic Log Data[J]. Geophysics， 2009， 74（6）： 211-223.

[10] Wei Z T， Tang X M. Numerical Simulation of Radiation， Reflection， and Reception of Elastic Waves from a Borehole Dipole Source[J]. Geophysics， 2012， 77（6）： D253-D261.

[11] Tang X M， Cao J J， Wei Z T. Shear-Wave Radiation， Reception， and Reciprocity of a Borehole Dipole Source： With Application to Modeling of Shear-Wave Reflection Survey[J]. Geophysics， 2014， 79（2）： T43-T50.

[12] 唐曉明， 魏周拓. 利用井中偶極聲源遠場輻射特性的遠探測測井[J]. 地球物理學(xué)報， 2012， 55（8）： 2798-2807.

Tang Xiaoming， Wei Zhoutuo. Remote Detection Logging Using Far-Field Radiation Characteristics of Dipole Sound Source in Well[J]. Chinese Journal of Geophysics， 2012， 55（8）： 2798-2807.

[13] 魏周拓. 反射聲波測井?dāng)?shù)值與物理模擬研究[D]. 青島： 中國石油大學(xué) （華東）， 2011.

Wei Zhoutuo. Numerical and Physical Simulation of Reflective Acoustic Logging[D]. Qingdao： China University of Petroleum （East China）， 2011.

[14] 李丹.反射聲波測井的成像研究[D]. 北京：中國石油大學(xué)（北京），2021.

Li Dan. Study on Imaging for Acoustic Reflection Imaging Logging[D]. Beijing： China University of Petroleum， Beijing，2021.

[15] 董經(jīng)利，許孝凱，張晉言，等. 聲波遠探測技術(shù)概述及發(fā)展[J]. 地球物理學(xué)進展，2020， 35（2）：566-572.

Dong Jingli， Xu Xiaokai， Zhang Jinyan， et al. Overview and Development of Acoustic Remote Detection Technology[J]. Advances in Geophysics， 2020， 35（2）： 566-572.

[16] 李航，孫宇航，李佳慧，等. 基于波動方程的地震波數(shù)值模擬研究綜述[J/OL]. 吉林大學(xué)學(xué)報（地球科學(xué)版），[2024-07-29]. https：//doi.org/10.13278/j.cnki.jjuese. 20230308.

Li Hang， Sun Yuhang， Li Jiahui，et al. Comprehensive Review of Numerical Simulation of Seismic Waves Based on the Wave Equation[J/OL]. Journal of Jilin University （Earth Science Edition）， [2024-07-29]. https：//doi.org/10. 13278/j.cnki.jjuese.20230308.

[17] Biot M A. Theory of Propagation of Elastic Waves in a Fluid-Saturated Porous Solid： I： Low-Frequency Range[J]. The Journal of the Acoustical Society of America， 1956， 28（2）： 168-178.

[18] Biot M A. Theory of Propagation of Elastic Waves in a Fluid-Saturated Porous Solid： II： Higher Frequency Range[J]. The Journal of the Acoustical Society of America， 1956， 28（2）： 179-191.

[19] 唐曉明，鄭傳漢. 定量測井聲學(xué)[M]. 北京：石油工業(yè)出版社，2004.

Tang Xiaoming， Zheng Chuanhan. Quantitative Log Acoustics[M]. Beijing： Petroleum Industry Press， 2004.

[20] 曹景記， 唐曉明， 魏周拓.偶極聲源在裸眼井及套管井外的橫波輻射特征[J]. 地球物理學(xué)報， 2014， 57（5）： 1683-1692.

Cao Jingji， Tang Xiaoming， Wei Zhoutuo. S-Wave Radiation Characteristics of Dipole Sound Source Outside Open Hole and Cased-Hole[J]. Chinese Journal of Geophysics， 2014， 57（5）： 1683-1692.

[21] 楊頂輝，滕吉文. 各向異性介質(zhì)中三分量地震記錄的FCT有限差分模擬[J]. 石油地球物理勘探， 1997， 32（2）： 181-190.

Yang Dinghui， Teng Jiwen. FCT Finite-Difference Simulation of Three-Component Seismic Records in Anisotropic Media[J]. Petroleum Geophysical Exploration， 1997， 32（2）： 181-190.

[22] Ursula Iturrarán-Viveros， Molero Miguel. Simulation of Sonic Waves Along a Borehole in a Heterogeneous Formation： Accelerating 2.5-D Finite Differences Using [Py]OpenCL[J]. Computers amp; Geosciences， 2013， 56： 161-169.

[23] Pan G A， Abubakar A， Habashy T M. An Effective Perfectly Matched Layer Design for Acoustic Fourth-Order Frequency-Domain Finite-Difference Scheme[J]. Geophysical Journal International， 2012， 188（1）： 211-222.

[24] Hiroshi Takenaka， Tanaka Hiroki， Okamoto Taro， et al. Quasi-Cylindrical 2.5D Wave Modeling for Large-Scale Seismic Surveys[J]. Geophysical Research Letters， 2003， 30（21）： 1.

[25] 徐冰，李文博，王易安，等. 不規(guī)則井眼的斯通利波反射數(shù)值模擬[J]. 科學(xué)技術(shù)與工程， 2012， 12（9）： 2029-2032.

Xu Bing， Li Wenbo， Wang Yian， et al. Numerical Simulation of Stoneleigh Wave Reflection in Irregular Borehole[J]. Science， Technology and Engineering， 2012， 12（9）： 2029-2032.

[26] Wang T， Tang X. Finite-Difference Modeling of Elastic Wave Propagation： A Nonsplitting Perfectly Matched Layer Approach[J]. Geophysics， 2003， 68（5）： 1749-1755.

[27] Wang X， Hornby B， Dodds K. Dipole Sonic Response in Deviated Borehole Penetrating an Anisotropic Formation[C]//72nd Annual Meeting. Salt Lake City： Society of Exploration Geophysicists， 2002： 360-363.

[28] Kurkjian A L， Chang S K. Acoustic Multipole Sources in Fluid-Filled Boreholes[J]. Geophysics， 1986， 51（1）： 148-163.