多元約束的古巖溶碳酸鹽巖洞穴儲層分布建模方法
——以塔河油田奧陶系油藏為例

劉鈺銘，侯加根，李永強,2，馬曉強，董越，宋隨宏

1 中國石油大學(北京)地球科學學院，北京 102249 2 中國石化石油勘探開發(fā)研究院，北京 100083

0 引言

古巖溶碳酸鹽巖儲層以大型巖溶洞穴為主要儲集空間，非均質性嚴重，儲層建模難度大[1-9]。為提高建模過程中井間預測的準確性，已有學者對該類特殊儲層的地質建模進行了探索研究。最初，王根久等(2002)[10]、趙敏等(2008)[11]采用地震屬性約束建立了儲層孔隙度三維模型；楊輝延等(2004)[12]、張淑品等(2007)[13]在巖溶帶控制下建立了儲層屬性模型。后來，趙彬等(2011)[14]、劉鈺銘等(2012、2018)[15-16]提出了基于垂向巖溶帶劃分和平面巖溶古地貌分區(qū)的建模思路，加強了地質因素在建模中的控制作用；在此基礎上，侯加根等(2012、2013)[17-18]、馬曉強等(2013)[19]提出了基于儲集體尺度分類的“多類多尺度建?！狈椒?。這些方法大多采用單一因素約束建模，而古巖溶洞穴儲層空間分布復雜且隨機性強，單一因素約束建模難以保證建模結果的準確性；此外，隨著油田生產(chǎn)對地質模型的要求日益提升，也亟需研究更可靠的建模方法。因此，本文擬以塔河油田奧陶系油藏為研究對象，細分古巖溶洞穴儲層類型，建立不同類型古巖溶洞穴儲集體三維分布模型及相應的儲層物性模型，探討針對古巖溶碳酸鹽巖洞穴儲層的建模方法，以期進一步揭示該類儲層非均質特征。

1 研究區(qū)概況

塔河油田位于塔里木盆地北部坳陷區(qū)沙雅隆起南翼的阿克庫勒凸起斜坡帶上，西鄰哈拉哈塘凹陷，東靠草湖凹陷，南接滿加爾坳陷，北部為雅克拉—輪臺斷凸(圖1)。塔河油田奧陶系碳酸鹽巖巖溶屬風化殼巖溶，主要發(fā)育于海西早期，海西早期強烈的構造抬升、褶皺和斷裂活動，造成志留—泥盆系、上奧陶統(tǒng)的普遍剝蝕、缺失，中—下奧陶統(tǒng)頂部也受到部分剝蝕，形成區(qū)域性不整合面(風化殼)[3]。1997年塔河油田奧陶系油藏發(fā)現(xiàn)并投入開發(fā)[9]，根據(jù)油田開發(fā)管理，塔河油田奧陶系古巖溶油藏分為塔河主體區(qū)的二、三、四、六、七、八區(qū)及外圍的十、十一、十二區(qū)、S72井區(qū)等多個開發(fā)區(qū)，每個開發(fā)區(qū)又進一步分為若干個單元。本文選取塔河油田四區(qū)S48單元為研究對象。S48單元面積為13.9 km2，共有鉆井24口，是塔河油田四區(qū)最大的開發(fā)單元，也是塔河油田奧陶系油藏最早的開發(fā)單元。目前該單元已進入產(chǎn)量遞減階段，注水、注氣開發(fā)是其主要開發(fā)方式，開發(fā)方案的調整迫切需要建立精確的三維地質模型。

圖1 塔河油田構造位置圖Fig. 1 Structural map showing the location of Tahe Oil fi eld

2 洞穴儲層類型與識別

2.1 洞穴儲層類型細分

從儲集單元成因和形態(tài)分析出發(fā)，將洞穴儲集體分為暗河型洞穴、落水洞型洞穴、內(nèi)部離散洞穴3種儲層類型。

2.1.1 暗河型洞穴

暗河，又稱地下河，是由地下水匯集，或地表水沿地下巖石裂隙滲入地下，經(jīng)過巖石溶蝕、坍塌以及水的搬運而形成的地下巖溶管道，多發(fā)育于水平潛流帶上部(又稱“飽水帶”)。暗河往往橫向發(fā)育，在平面上可分支分叉，組成管道網(wǎng)絡，垂向上由于地下水位波動，亦可形成多層河道[20]。

2.1.2 落水洞型洞穴

落水洞是地表水沿垂直裂隙進行溶蝕、沖蝕，使裂隙擴大并伴隨部分崩塌作用的產(chǎn)物，也是巖溶水從地表到地下暗河的主要通道。主要發(fā)育于近地表的垂向滲流帶(又稱“包氣帶”)，表面形態(tài)與漏斗相似，整體上可呈裂隙狀、筒狀、錐狀或袋狀。

2.1.3 內(nèi)部離散洞穴

暗河型洞穴與落水洞型洞穴的發(fā)育都需要相對穩(wěn)定的巖溶條件和一定的巖溶作用時間，若局部巖溶作用由于某些條件發(fā)生變化而終止的話，比如潛水面的變化，則這些區(qū)域發(fā)育的洞穴就會停止擴大，而相對孤立地保存下來。它可以是暫時性水流沿裂縫垂直下滲在碳酸鹽巖內(nèi)部形成的溶蝕擴大洞穴，也可能是短時性水平流動形成的地下河的雛形。但考慮到地下資料的局限，而無法將其清楚分開，本文統(tǒng)稱為內(nèi)部離散洞穴。

2.2 洞穴儲層類型單井識別

單井上識別的洞穴儲層數(shù)據(jù)，是建模過程中最直接的數(shù)據(jù)輸入。人們已能從鉆井放空漏失信息，并結合地震剖面“串珠狀”反射等特征在單井上識別洞穴的發(fā)育層段[21-22]，但仍需進一步判別洞穴的類型。本文提出了在單井上進一步判斷洞穴儲層類型的3點依據(jù)。一是洞穴充填物的差異，暗河型洞穴鉆井取心常見砂泥等外部搬運沉積物，其對應自然伽馬曲線發(fā)育明顯高值段(圖2a)，而落水洞型洞穴和內(nèi)部離散洞穴一般無砂泥等搬運沉積物；二是地震響應差異，暗河型洞穴的響應(如反演波阻抗數(shù)據(jù))在剖面上往往橫向發(fā)育，在平面上有一定的延伸規(guī)模，而落水洞型洞穴與內(nèi)部離散洞穴規(guī)模相對小(圖2b，2c，2d)，且落水洞型洞穴沿縱向發(fā)育特征較明顯；三是發(fā)育深度的差異，落水洞型洞穴大多從風化殼(不整合面)開始向下發(fā)育，暗河型洞穴往往距離風化殼有一定距離(塔河四區(qū)為60 m左右)，而內(nèi)部離散洞穴可發(fā)育于各個深度段。

圖2 洞穴儲層類型綜合識別(據(jù)參考文獻[16,19]，有修改)Fig. 2 Comprehensive identi fi cation of cavernous reservoirs (modi fi ed from reference[16, 19])

如T403井，鉆至5488 m時，發(fā)生鉆井液漏失，5488～5554 m段取心明顯見到大量砂泥沉積和垮塌角礫的互層，砂泥沉積段對應自然伽馬明顯高值(圖2a)，而從過該井的波阻抗剖面來看，洞穴橫向展布呈連續(xù)條帶狀(圖2b)，綜合判定其為暗河型洞穴。

3 建模數(shù)據(jù)準備

建模數(shù)據(jù)主要包括井數(shù)據(jù)、地震數(shù)據(jù)和大型斷裂數(shù)據(jù)。井數(shù)據(jù)包括上述24口井的單井洞穴儲層識別結果(類型及發(fā)育層段)以及孔隙度等測井曲線數(shù)據(jù)。地震數(shù)據(jù)為深度域反演波阻抗數(shù)據(jù)。本文“大型斷裂”是相對測井或巖心規(guī)模的斷裂而言，特指通過地震資料能直觀識別的斷裂，即油藏規(guī)模的斷裂。之所以要準備大型斷裂數(shù)據(jù)，是因為多數(shù)巖溶洞穴在成因上與斷裂有關，先期斷裂一定程度上控制巖溶洞穴的發(fā)育[23-24]。大型斷裂數(shù)據(jù)主要通過人工地震解釋和螞蟻體等屬性分析[25]獲取，已有成熟方法，在此不再贅述。

此外，還需要網(wǎng)格化的地層模型作為建模載體。針對S48單元，以奧陶系頂面為地質模型的頂面，頂面以下240 m為底面，設計網(wǎng)格大小為平面25 m×25 m，垂向為1 m，建立了地層網(wǎng)格數(shù)據(jù)體，總網(wǎng)格數(shù)約為580萬。

4 多元約束建模方法

成因控制是儲層地質建模的一貫策略。常用的“等時建?！薄ⅰ?沉積相相控建?！钡炔呗詢H適用于沉積型的常規(guī)碎屑巖儲層。針對“改造型”古巖溶碳酸鹽巖儲層，有學者嘗試引入地質控制因素的干預，提出縱向上按“巖溶帶”分層建模，平面上按照“巖溶古地貌”分區(qū)建模的建模策略[14-15]，具有積極的指導意義，但其內(nèi)涵還需進一步商榷，具體做法也需進一步完善。本文提出了大型斷裂、巖溶帶和地震波阻抗等多元約束的建模方法。

4.1 斷裂約束建模

前人按“巖溶古地貌”分區(qū)建模的依據(jù)是將古地貌視為控制縫洞儲集體發(fā)育的基本因素，不同的古地貌單元具有不同的縫洞儲集體發(fā)育特征[14-15]。嚴格來說，古巖溶作用開始時候的古地貌確實一定程度上控制著巖溶的發(fā)育。但是 “峰叢、洼地”等巖溶地貌往往是巖溶作用的結果，是巖溶作用在近地表的物質體現(xiàn)，即巖溶地貌是古巖溶作用過程的結果，其與地下碳酸鹽巖地層中儲集體的形態(tài)、規(guī)模和發(fā)育程度是否有明確的相關關系還需進一步商榷；此外，有學者嘗試用地層回撥、古構造恢復等方法重構巖溶古地貌，但由于古巖溶作用時間長，地層剝蝕差異程度大，該研究也局限于盆地規(guī)模，對于油藏規(guī)模的巖溶古地貌尚難以獲得較為準確的結果。

在大規(guī)模巖溶發(fā)育的地質背景下，巖溶儲集體發(fā)育的控制因素是地質構造，尤其是斷裂構造[26]，古巖溶洞穴儲集體發(fā)育往往與大型斷裂密切相關[24，26]。S48單元識別出的47個各類洞穴中，38個與大型斷裂直接相交，約占81%；8個在緊鄰斷裂100 m范圍內(nèi)，占17%；只有1個與大型斷裂距離超過100 m，占2%(圖2d)。油氣生產(chǎn)實踐也證明了這一點，塔河油田四區(qū)有產(chǎn)能的生產(chǎn)井中，在斷裂附近300 m以內(nèi)的占90%以上。具體來說，落水洞型洞穴往往發(fā)育于斷裂相交處、拐彎處或斷裂末端，S48單元識別出的19個落水洞型洞穴中，4個處于大型斷裂相交處，8個處于斷裂拐彎處，5個處于斷裂末端(圖2d)；暗河型洞穴上部往往與大型斷裂或落水洞型洞穴相連，S48單元識別出的6個暗河型洞穴累計長度約10 km，共有24條大型斷裂與之相交，有12個落水洞型洞穴與暗河上部相連，即平均約420 m暗河就發(fā)育一條大型斷裂，約830 m發(fā)育一個落水洞型洞穴(圖2d)。

因此，從成因和可操作角度考慮，“斷裂約束建模”比“古地貌分區(qū)建?！备幸饬x?！皵嗔鸭s束建模”的具體做法分為三步：第一步，基于上述洞穴儲層的識別成果(圖2d)，統(tǒng)計S48單元內(nèi)所有洞穴離最近斷裂的距離(簡稱最近距離)，得到洞穴發(fā)育隨最近距離變化的概率密度曲線(圖3)，即最近距離與洞穴發(fā)育的概率關系。第二步，基于大型斷裂空間分布，計算每一個地層網(wǎng)格的最近距離，得到基于三維網(wǎng)格的最近距離數(shù)據(jù)體。第三步，依據(jù)第一步中得到的概率關系將第二步中得到最近距離數(shù)據(jù)體換算為溶洞發(fā)育概率數(shù)據(jù)體，稱為“斷裂約束溶洞發(fā)育概率數(shù)據(jù)體”在接下來的模擬中將此概率作為建模的約束參數(shù)之一。

4.2 巖溶帶約束建模

現(xiàn)代巖溶理論依據(jù)地表水和地下水的運動方式在垂向上劃分巖溶帶，巖溶帶一定程度上控制著洞穴的類型。基于此，提出了“巖溶帶”分層建模方法[14-15]，具體做法是首先在單井上劃分巖溶帶，將巖溶帶視為某種“層”進行井間對比，將地層按照巖溶帶分層，如將不整合面以下約0～60 m視為垂直滲流巖溶帶，60～140 m視為水平潛流巖溶帶；然后在不同的帶內(nèi)分別模擬不同洞穴的分布。然而，巖溶帶的界面是由地下水位決定的，而在漫長的地質歷史時期，地下水位頻繁波動，在油田覆蓋區(qū)由于地下資料的不完備性，巖溶帶劃分主觀性大；另外，巖溶帶不同于地層單元，井間對比也缺乏嚴謹?shù)睦碚撘罁?jù)。

本文對這一方法進行了改進，并不對地層進行嚴格的巖溶帶劃分，而是針對S48單元內(nèi)已經(jīng)證實的洞穴，統(tǒng)計其離奧陶系頂部不整合面的距離(即巖溶發(fā)育深度)，得到不同類型洞穴儲層隨巖溶發(fā)育深度的概率密度曲線(圖4)，據(jù)此可將地層網(wǎng)格的垂向深度數(shù)據(jù)換算出不同類型洞穴儲層的發(fā)育概率，稱為“巖溶帶約束溶洞發(fā)育概率數(shù)據(jù)體”，將此概率作為建模的約束參數(shù)之一。

4.3 地震波阻抗約束建模

井震數(shù)據(jù)整合是現(xiàn)代油藏地質建模研究的重要手段[27]。由于儲層分布的不均勻性，塔河油田開發(fā)井網(wǎng)亦不規(guī)則，且井距較大，大多在400～1000 m。因此，充分整合利用三維地震資料，以降低建模井間預測的不確定性，彌補井數(shù)據(jù)的不足，就尤為重要。

塔河油田奧陶系油藏埋深超過5000 m，其地震資料品質并不高。目前，對洞穴響應較為敏感的是波阻抗反演數(shù)據(jù)。碳酸鹽巖本身致密堅硬，呈高速度、高密度即高波阻抗特征，當內(nèi)部發(fā)育縫洞時，就會引起速度、密度的降低即波阻抗的降低[28]。以塔河油田四區(qū)為例，統(tǒng)計56口井上識別的81個溶洞所對應的波阻抗數(shù)據(jù)，約87%的洞穴對應的絕對波阻抗值在1.0×107～1.6×107kg/m3·m/s。但波阻抗反演數(shù)據(jù)并不能作為反映洞穴的確定性信息，這是因為波阻抗反演數(shù)據(jù)具有地震信息固有的多解性，其對地下情況的解釋不是唯一的[28]，即洞穴與相應波阻抗信息并不是一一對應的，具有類似的波阻抗信息特征并不都是洞穴。

因此，分析不同類型洞穴發(fā)育與波阻抗信息間的相關關系，將波阻抗數(shù)據(jù)體轉化為不同類型洞穴單元發(fā)育概率數(shù)據(jù)體。具體做法是讀取S48單元內(nèi)所有洞穴型儲層對應的波阻抗數(shù)據(jù)，統(tǒng)計得到不同類型洞穴儲層隨波阻抗值變化的發(fā)育概率曲線(圖5)，據(jù)此將地震波阻抗數(shù)據(jù)體轉化為不同類型洞穴儲層的發(fā)育概率數(shù)據(jù)體，稱為“波阻抗約束溶洞發(fā)育概率數(shù)據(jù)體”，在模擬中將此概率作為建模的約束參數(shù)之一。

4.4 多元約束條件整合算法

圖3 洞穴儲層發(fā)育概率與最近斷裂距離關系曲線Fig. 3 Developmental probability curve of cavernous reservoir with the distance to the closest fault

圖4 洞穴儲層發(fā)育概率與距頂不整合面深度關系曲線Fig. 4 Developmental probability curve of cavernous reservoir with the distance to the top unconformity

圖5 洞穴儲層概率與波阻抗關系曲線Fig. 5 Developmental probability curve of cavernous reservoir with the P-impedance

斷裂約束、巖溶帶約束和波阻抗約束都是通過各自的概率體數(shù)據(jù)來實現(xiàn)的，在建模中需將其進一步整合到建模算法中。

目前，用于建立儲層離散變量(如沉積相類型)三維模型的隨機模擬算法主要有序貫指示模擬、多點地質統(tǒng)計學方法和基于目標的建模方法等等?；谀繕说慕７椒ㄒ蟠M儲層屬性具有明確量化的幾何學特征，如河道長、寬、高以及形狀方向等等；多點地質統(tǒng)計學方法是新興的模擬算法，但其要求得到與待模擬目標地質體在空間展布上吻合程度高的訓練圖像。但目前對于地下碳酸鹽巖巖溶儲層的先驗認識還遠遠不足以得到準確的幾何形態(tài)參數(shù)或訓練圖像，導致基于目標的建模方法和多點地質統(tǒng)計學模擬方法較難適用。

序貫指示模擬算法可用于模擬復雜各向異性的地質現(xiàn)象，對先驗認識的要求相對寬松，并可以通過帶有趨勢的克里金或同位協(xié)同指示克里金等模擬算法整合次級約束數(shù)據(jù)。因此，本次建模采用帶有趨勢的序貫指示克里金模擬算法。該模擬方法的關鍵是條件累積概率分布函數(shù)(ccdf)[30]。對于三維空間的每一網(wǎng)格節(jié)點，條件累積概率分布函數(shù)(ccdf)取決于該節(jié)點處各類型的條件概率P( u, k )。該概率滿足：

式中，P( u, k )為變量k(洞穴類型變量)在位置點u的概率，即在位置點u，所有變量類型的概率之和為1。

在整個方法體系中，帶有趨勢的序貫指示克里金模擬算法已經(jīng)非常成熟且廣泛使用，其核心在于如何將前文所述的“波阻抗概率數(shù)據(jù)體”、“ 斷裂約束概率數(shù)據(jù)體”和“巖溶帶約束概率數(shù)據(jù)體”3個次級約束數(shù)據(jù)整合為單一數(shù)據(jù)體進而作為序貫指示克里金模擬算法的趨勢約束。本文采用基于最優(yōu)估計的數(shù)據(jù)融合理論[31]將這3個不同來源的條件概率數(shù)據(jù)體融合為一個最優(yōu)化的概率數(shù)據(jù)體，以實現(xiàn)這三個因素的共同約束。

式中，P( u, k)為變量k在位置點u的概率；P( u, k)波阻抗為波阻抗概率數(shù)據(jù)體中，變量k在位置點u的概率；P( u, k)斷裂為斷控概率數(shù)據(jù)體中，變量k在位置點u的概率；P( u, k)巖溶帶為巖溶帶控概率數(shù)據(jù)體中，變量k在位置點u的概率；F為最優(yōu)估計的數(shù)據(jù)融合算法。

5 建模結果及驗證

通過上述建模數(shù)據(jù)及方法，建立了塔河油田四區(qū)S48單元洞穴儲層分布模型(圖6)，并在此模型的控制下，采用傳統(tǒng)方法，將單井孔隙度數(shù)據(jù)作為條件數(shù)據(jù)，地震反演波阻抗數(shù)據(jù)作為約束數(shù)據(jù)，采用序貫高斯模擬算法，模擬了不同類型洞穴儲集體內(nèi)部的孔隙度參數(shù)分布，建立了儲層孔隙度三維模型。

基于S48單元孔隙度三維模型，結合生產(chǎn)動態(tài)揭示的油水界面信息，計算該單元的地質儲量為1386×104t，將該單元石油地質儲量豐度平面圖與單井累計產(chǎn)油量對比發(fā)現(xiàn)，產(chǎn)能高的區(qū)域石油地質儲量豐度也高(圖7)，間接證實建立的洞穴儲層三維分布模型和孔隙度模型具有較高的可靠性。

圖7 S48單元奧陶系油藏儲量豐度與累計產(chǎn)油量平面對比圖(據(jù)參考文獻[16]，有修改)Fig. 7 Comparison between reserve abundance and cumulative production in S48 unit (modi fi ed from reference[16])

6 結論與展望

建立更加可靠實用的儲層地質模型，進一步表征儲層內(nèi)部非均質性，是油田開發(fā)地質工作的一貫要求。本文提出了斷裂、巖溶帶和波阻抗數(shù)據(jù)聯(lián)合約束建立洞穴型儲層三維分布模型的新方法，建立了塔河油田四區(qū)S48單元儲層模型。但由于洞穴充填物質多樣，充填程度不均[21-22]，該類儲層建模，尤其是不同類型洞穴內(nèi)部屬性分布建模，還需繼續(xù)努力突破。明晰不同類型巖溶洞穴在后期埋藏環(huán)境中經(jīng)歷的垮塌等地質作用過程對儲層質量的影響機制，建立不同類型洞穴儲層內(nèi)部儲集空間類型及儲集性能參數(shù)的空間分布模式，是更準確表征儲層內(nèi)部儲層物性差異的前提。

此外，古巖溶儲層儲集空間的特殊性導致了儲層內(nèi)部油氣富集程度規(guī)律和流體流動機制的復雜性，傳統(tǒng)的含油氣飽和度和滲透率等參數(shù)難以適用，亟需選取合適的物理參數(shù)表征儲層內(nèi)部油氣富集程度和流體流動能力，并揭示這些參數(shù)的在儲層內(nèi)部的分布規(guī)律。

[1] 杜金虎, 周新源, 李啟明, 等. 塔里木盆地碳酸鹽巖大油氣區(qū)特征與主控因素[J]. 石油勘探與開發(fā), 2011, 38(06): 652-661. [DU J H, ZHOU X Y, LI Q M, et al. Characteristics and controlling factors of the large carbonate petroleum province in the Tarim Basin,China[J]. Petroleum Exploration and Development, 2011, 38(06): 652-661.]

[2] 張寶民, 劉靜江. 中國巖溶儲集層分類與特征及相關的理論問題[J]. 石油勘探與開發(fā), 2009, 36(1): 12-29. [ZHANG B M, LIU J J. Classi fi cation and characteristics of karst reservoirs in China and related theories[J]. Petroleum Exploration and Development, 2009,36(1): 12-29.]

[3] 焦方正, 翟曉先. 海相碳酸鹽巖非常規(guī)大油氣田—塔河油田勘探研究與實踐[M]. 北京: 石油工業(yè)出版社, 2008. [JIAO F Z, ZHAI X X. Giant oil fi eld of marine carbonate origin[M]. Beijing: Petroleum Industry Press, 2008.]

[4] WALKDEN G M. Palaeokarstic surfaces in upper Visean (Carboniferous) limestones of the Derbyshire Block, England [J]. Journal of Sedimentary Petrology, 1974, 44(4): 1232-1247.

[5] LOUCKS R G. Paleocave carbonate reservoirs: Origins, burial-depth modi fi cations, spatial complexity, and reservoir implications [J].AAPG bulletin 1999, 83: 1795-1834.

[6] ZHANG K, WANG D R, BRYAN G H. Reservoir characterization of the Odovician oil and gas pools in the Tahe Oil fi eld, Tarim Basin,Northwest China [J]. Petroleum Exploration and Development, 2004, 31(1): 123-126.

[7] LOUCKS R G, MESCHER P K, MCMECHAN G A. Three-dimensional architecture of a coalesced, collapsed-paleocave system in the Lower Ordovician Ellenburger Group, central Texas [J]. AAPG bulletin, 2004, 88(5): 545-564.

[8] SMOSNA R, BRUNER K R, RILEY R A. Paleokarst and reservoir porosity in the Ordovician Beekmantown dolomite of the central Appalachian basin [J]. Carbonates and Evaporites, 2005, 20(1): 50-63.

[9] 張抗. 塔河油田的發(fā)現(xiàn)及其地質意義[J]. 石油與天然氣地質, 1999, 20(2): 24-28. [ZHANG K. The discovery of Tahe Oil fi eld and its geologic implication[J]. Oil and Gas Geology, 1999, 20(2): 24-28.]

[10] 王根久, 王桂宏, 余國義, 等. 塔河碳酸鹽巖油藏地質模型[J]. 石油勘探與開發(fā), 002, 29(1): 109-111. [WANG G J, WANG G H,YU G Y, et al. Geological model of carbonate reservoirs in Tahe Oil fi eld[J]. Petroleum Exploration and Development, 2002, 29(1): 109-111.]

[11] 趙敏, 康志宏, 劉潔. 縫洞型碳酸鹽巖儲集層建模與應用[J]. 新疆石油地質, 2008, 29(3): 318-320. [ZHAO M, KANG Z H, LIU J.Modeling and application of fractured-vuggy carbonate reservoirs[J]. Xinjiang Petroleum Geology, 2008, 29(3): 318-320.]

[12] 楊輝廷, 江同文, 顏其彬, 等. 縫洞型碳酸鹽巖儲層三維地質建模方法初探[J]. 大慶石油地質與開發(fā), 2004, 23(4): 11-12. [YANG H T, JIANG T W, YAN Q B, et al. Establishment of 3D geological models of fracture-cave carbonate reservoir[J]. Petroleum Geology and Oil fi eld Development in Daqing, 2004, 23(4): 11-12.]

[13] 張淑品, 陳福利, 金勇. 塔河油田奧陶系縫洞型碳酸鹽巖儲集層三維地質建模[J]. 石油勘探與開發(fā), 2007, 34(2): 175-180.[ZHANG S P, CHEN F L, JIN Y. 3D modeling of Ordovician fractured-vug carbonate reservoir, Tahe Oil fi eld[J]. Petroleum Exploration and Development, 2007, 34(2): 175-180.]

[14] 趙彬, 侯加根, 劉鈺銘. 塔河油田奧陶系碳酸鹽巖溶洞型儲層三維地質建模與應用[J]. 石油天然氣學報, 2011, 33(5): 12-16. [ZHAO B, HOU J G, LIU Y M. Three dimensional modeling and application of Ordovician carbonate cavity reservoirs in Tahe Oil fi eld[J]. Journal of Oil and Gas Technology, 2011, 33(5): 12-16.]

[15] 劉鈺銘, 侯加根, 胡向陽, 等. 塔河油田古巖溶儲集體三維建模[J]. 中國石油大學學報(自然科學版) , 2012, 36(2): 34-38. [LIU Y M, HOU J G, HU X Y, et al. 3D modeling of paleokarst reservoir in Tahe Oil fi eld[J]. Journal of China University of Petroleum, 2012,36(2): 34-38.]

[16] LIU Y M, HOU J G, LI Y Q, et al. Characterization of architectural elements of Ordovician fractured-cavernous carbonate reservoirs,Tahe Oil fi eld, China[J]. Journal of the Geological Society of India, 2018, 91(3): 315-322.

[17] 侯加根, 馬曉強, 劉鈺銘, 等. 縫洞型碳酸鹽巖儲層多類多尺度建模方法研究—以塔河油田四區(qū)奧陶系油藏為例[J]. 地學前緣,2012, 19(2): 59-66. [HOU J G, MA X Q, LIU Y M, et al. A case study of Ordovician reservoir of 4th block in Tahe Oil fi eld[J]. Earth Science Frontiers, 2012, 19(2): 59-66.]

[18] 侯加根, 馬曉強, 胡向陽, 等. 碳酸鹽巖溶洞型儲集體地質建模的幾個關鍵問題[J]. 高校地質學報, 2013, 19(01): 64-69. [HOU J G, MA X Q, HU X Y, et al. Key issues of 3D geological modeling of paleokarst-cave carbonate reservoir[J]. Geological Journal of China Universities, 2013, 19(01): 64-69.]

[19] 馬曉強, 侯加根, 胡向陽, 等. 論古巖溶洞穴型儲層三維地質建模—以塔河油田四區(qū)奧陶系儲層為例[J]. 地質論評, 2013, 59(02):315-324. [MA X Q, HOU J G, HU X Y, et al. Discussion on 3D geological modeling of palaeokarst cave reservoir: A case study on ordovician reservoir of the 4th block in the Tahe Oil fi eld[J]. Geological Review, 2013, 59(02): 315-324.]

[20] 金強, 田飛. 塔河油田巖溶型碳酸鹽巖縫洞結構研究[J]. 中國石油大學學報(自然科學版) , 2013, 37(5): 15-21. [JIN Q, TIAN F.Investigation of fracture-cave constructions of karsted cabonate reservoirs of Ordovician in Tahe Oil fi eld, Tarim Basin[J]. Journal of China University of Petroleum, 2013, 37(5): 15-21.]

[21] 胡向陽, 權蓮順, 齊得山, 等. 塔河油田縫洞型碳酸鹽巖油藏溶洞充填特征[J]. 特種油氣藏, 2014, 21(01): 18-21. [HU X Y,QUAN L S, QI D S, et al. Features of of cavern fi lling in fractured/vuggy carbonate oil reservoirs, Tahe Oil fi eld[J]. Special oil and gas reservoirs, 2014, 21(01): 18-21.]

[22] 田飛, 金強, 李陽, 等. 塔河油田奧陶系縫洞型儲層小型縫洞及其充填物測井識別[J]. 石油與天然氣地質, 2012, 33(06): 900-908. [TIAN F, JIN Q, LI Y, et al. Identi fi cation of small fracture-vugs and their fi llings through log interpretation in fractured-vuggy Ordovician reservoirs in Tahe oil fi eld[J]. Oil and Gas Geology, 2012, 33(06): 900-908.]

[23] 顧家裕. 塔里木盆地輪南地區(qū)下奧陶統(tǒng)碳酸鹽巖巖溶儲層特征及形成模式[J]. 古地理學報, 1999, 1(1): 54-60. [GU J Y. Char-acteristics of and evolution model of karst reservoirs of lower Ordovician carbonate rocks in Lunnan area of Tarim Basin[J]. Journal of Paleogeology, 1999, 1(1): 54-60.]

[24] 周文, 李秀華, 金文輝, 等. 塔河奧陶系油藏斷裂對古巖溶的控制作用[J]. 巖石學報, 2011, 27(8): 2339-2348. [ZHOU W, LI X H,JIN W H, et al. The control action of fault to paleokarst in view of Ordovician reservoir in Tahe area[J]. Acta Petrologica Sinica, 2011,27(8): 2339-2348.]

[25] 張欣. 螞蟻追蹤在斷層自動解釋中的應用—以平湖油田放鶴亭構造為例[J]. 石油地球物理勘探, 2010, 45(02): 278-281. [ZHANG X. Application of ant tracing algorithm in fault auto-matic interpretation: A case study on Fangheting structure in Pinghu Oil fi eld[J]. Oil Geophysical Prospecting, 2010, 45(02): 278-281.]

[26] 譚周地. 構造對巖溶的控制意義[J]. 長春地質學院學報, 1978, 8(2): 43-47. [TAN Z D. Controlling effect of structure on karst[J].Journal of Jilin University, 1978, 8(2): 43-47.]

[27] 印興耀, 劉永社. 儲層建模中地質統(tǒng)計學整合地震數(shù)據(jù)的方法及研究進展[J]. 石油地球物理勘探, 2002, 37(4): 423-430. [YIN X Y, LIU Y S. Methods and development of integrating seismic data in reservoir model-building[J]. Oil Geophysical Prospecting, 2002,37(4): 423-430.]

[28] 溫志新, 王紅漫, 漆立新, 等. 塔河油田奧陶系縫洞型碳酸鹽巖儲層預測研究[J]. 地學前緣, 2008, 15(1): 94-100. [WEN Z X,WANG H M, QI L X, et al. Research on the prediction of Ordovician carbonatite reservoir of fi ssure-cave type in the Tahe oil fi eld[J].Earth Science Frontiers, 2008, 15(1): 94-100.]

[29] 黃捍東, 張如偉, 趙迪, 等. 塔河奧陶系碳酸鹽巖縫洞預測[J]. 石油地球物理勘探, 2009, 44(2): 213-218. [HUANG H D, ZHANG R W, ZHAO D, et al. Prediction of Ordovician carbonate fracture and cavern in Tahe area[J]. Oil Geophysical Prospecting, 2009, 44(2):213-218.]

[30] 吳勝和, 李宇鵬. 儲層地質建模的現(xiàn)狀與展望[J]. 海相油氣地質, 2007, 12(3): 53-60. WU S H, LI Y P. Reservoir modeling: Current situation and developmental prospect[J]. Maine Oil and Gas Geology, 2007, 12(3): 53-60.]

[31] 王炯琦, 周海銀, 吳翊. 基于最優(yōu)估計的數(shù)據(jù)融合理論[J]. 應用數(shù)學, 2007, 20(02): 392-399. [WANG J Q, ZHOU H Y, WU Y. The theory of data fusion based on state optimal estimation[J]. Methemeica Applicata, 2007, 20(02): 392-399.]