-
伊拉克 BU 油田一直采用天然能量开发，随着开发的深入，天然能量衰竭严重，采出程度仅为 6% 左右[1]，迫切需要注水开发，以确保地层压力平衡。为制定合理可行的注水方案，需对研究区不同沉积期次储层展布范围及演化规律进行更加详细地研究，确定优势储层发育带。为此，以BU油田Mishirif 组MB2 1 段为重点研究层位，在前人相关认识的基础上[1-2]，从 MB2 1 段沉积控制因素及沉积演化规律入手，综合应用岩心、地震、测井及分析化验等资料，选取沉积过程控制因素的定量参数，对各沉积期次的储层沉积过程进行数值模拟，获得不同岩性储层在时间和空间上的演化特征，以期为 BU 油田注水方案编制提供地质依据，也为类似厚层碳酸盐岩储层勘探开发提供指导和借鉴。
-
1 区域地质概况
-
BU油田属于伊拉克米桑（Missan）油田群，位于伊拉克东南部米桑省，毗邻伊朗边界，西北距巴格达约350 km，南距巴士拉约175 km。区域构造位于扎格罗斯造山带和阿拉伯台地东部边缘的过渡带上，属于扎格罗斯构造低角度褶皱带，BU 油田位于西背斜构造带，分为南、北2个高点，断层不发育[2]。 Mishirif组MB2 1 段是油田重要储油层之一，地质储量占整个油田群的 70% 以上[1]；沉积类型为碳酸盐岩台地，储层厚度大，平均厚度为80 m；依据生物碎屑含量将其岩性划分为颗粒灰岩、泥粒灰岩、粒泥灰岩、含粒泥晶灰岩和泥晶灰岩，根据岩性和电性特征，垂向上由上到下可以划分为Ⅰ—Ⅵ共6个小层； 沉积能量整体呈增高的趋势，但在Ⅴ小层和Ⅲ小层发育2个短时高能沉积带，因此，形成了高低能量交替的沉积环境。处于高能沉积环境的Ⅰ，Ⅲ和Ⅴ小层以颗粒灰岩和泥粒灰岩为主，其中Ⅰ小层储层厚度最大，平均为15 m，是油田主力产层；处于低能沉积环境的Ⅱ，Ⅳ和Ⅵ小层以泥晶灰岩、含粒泥晶灰岩和粒泥灰岩为主，其中Ⅵ小层储层厚度最大，平均为 25 m。储集空间类型包括原生粒间孔、铸模孔、次生粒间溶孔、基质微孔及少量溶蚀孔洞和微裂缝。受沉积和成岩作用影响，储层非均质性较强，孔隙度和渗透率变化范围较大；其中孔隙度为 1.4%～29.7%，平均为 17.2%，渗透率为 0.2～3 212 mD，平均为 16.2 mD；Ⅰ小层为高孔高渗透储层，其余均为中孔中渗透储层。
-
2 沉积过程数值模拟
-
2.1 技术流程
-
沉积过程数值模拟技术已在中外油田取得了一些应用成果，证明其是可行的[3-6]。但对于碳酸盐岩储层的沉积过程模拟实例较少，本次研究依据BU 油田的地震、测井、岩心和分析化验等资料，分析可容空间和碳酸盐岩产率对沉积过程数值模拟的控制作用。其中可容空间用古水深、地层沉降和海平面升降表征；碳酸盐岩产率用产率随沉积时间、深度和波浪的变化表征。由于难以获得准确的关键参数，因此，以单井地层厚度和岩性拟合程度为标准，通过优化参数，进行多次模拟实验，直至得出较为合理的沉积数值模型，分析碳酸盐岩台地沉积演化特征，具体技术流程如图1所示。
-
图1 沉积过程数值模拟技术流程
-
Fig.1 Technical flow chart of numerical simulation results of reservoir sedimentary process of
-
2.2 沉积数值模型建立
-
依据上述技术流程，开展 BU 油田碳酸盐岩储层沉积过程数值模拟。建立的沉积数值模型长度为36.6 km，宽度为14.4 km，面积为527.04 km2，平面网格间距为 0.2 km。模拟目的层为 MB2 1 段，储层平均厚度为80 m，沉积时间为距今96～90 Ma，总沉积时长为6 Ma，时间步长为0.1 Ma，岩性设置为颗粒灰岩、泥粒灰岩、粒泥灰岩、含粒泥晶灰岩和泥晶灰岩 5种。
-
2.3 沉积过程控制参数分析
-
可容空间变化、碳酸盐岩产率是影响碳酸盐岩沉积的主控因素[7-9]，因此这2个因素是沉积过程数值模拟的关键参数。基于沉积背景、古生物和岩心分析化验资料，通过井震结合定量化表征关键参数的初始值。
-
2.3.1 可容空间变化
-
可容空间是指沉积物表面与基准面之间可供潜在沉积物充填的空间，由早期未被充填遗留下来的古地形和后期形成的沉积地层组成，并随地质年代发生变化。早期遗留下来的古地形为初始可容空间，由基底初始水深决定，古地貌越低，基底初始水体越深，初始可容空间越大；反之，古地貌越高，基底初始水体越浅，初始可容空间越小。沉积过程中可容空间的变化是形成地层沉积的基础，受控于基底沉降与基准面的变化。基底沉降可由现今地层厚度与水深变化来反推，基准面的变化用海平面升降曲线来表征，二者的综合即为可容空间的总变化。
-
由研究区古生物和沉积环境资料可知，研究区碳酸盐岩储层沉积时期的水深为 0～50 m。以地震解释层位和井点分层为硬数据，建立研究区三维构造模型，对该模型进行褶皱恢复和去压实作用，恢复研究区古地貌，确定碳酸盐岩形成的初始地形。综合古生物和古地形分析，获得初始水深平面分布图（图2），确定初始可容空间。沉积过程中水深的变化幅度远小于地层厚度，因此，水深变化对可容空间的影响可以忽略。用 MB2 1 段的地层厚度与沉积时间的比值来表征基底沉降速率，平均沉降速率为14.3 m/Ma，最大沉降速率为22.7 m/Ma，最小沉降速率为7.3 m/Ma。由岩性旋回变化规律和古生物分布特征可知，底部Ⅵ小层以泥晶灰岩和含粒泥晶灰岩为主，顶部Ⅰ小层以颗粒灰岩为主，因此，从距今 96～90 Ma，海平面整体是一个由深变浅的过程，在中间有2期短暂的旋回变化，岩性以颗粒灰岩、泥粒灰岩和粒泥灰岩为主。由BU-3井的单井岩性解释结果和沉积厚度变化可得BU油田MB2 1 段相对海平面变化曲线。
-
图2 BU油田初始水深平面分布
-
Fig.2 Map of initial water depth of BU Oilfield
-
2.3.2 碳酸盐岩产率
-
碳酸盐岩在孤立台地上生长，其产率受沉积时间、水深、波浪改造作用共同影响。不同的沉积环境形成不同的沉积产物。通过薄片鉴定识别出岩性，标定测井曲线，并与自然伽马、中子、密度曲线相结合，采用神经网络算法，进行单井岩性划分[10]。通过统计实钻井点各小层不同岩性含量来表征沉积时期每种岩性的产率。按沉积时间单元统计各井点岩性含量与沉积时间的比值来定义各岩性随沉积时间变化的产率（表1）。同时，从单井沉积旋回划分结果来看，水体越浅，越容易沉积颗粒灰岩等高能粗粒物质；水体越深，越容易沉积泥晶灰岩等低能细粒物质，因此，可以定义岩性随水深的变化，即从水深为 0～50 m，每 10 m 范围内，由深到浅发育的沉积岩性依次为颗粒灰岩、泥粒灰岩、粒泥灰岩、含粒泥晶灰岩和泥晶灰岩。波浪改造作用包括波浪深度、角度及频率，不同的波浪改造作用对碳酸盐岩建造的影响不同，由岩性划分结果可见，研究区台地四周均有颗粒灰岩发育，说明在每个方向均受到波浪能量的改造，因此，波浪角度分别定义为 45°，135°，225°，315°，波浪改造作用深度范围定义为 5～20 m，碳酸盐岩台地呈 ES—WN 向延伸，其产率较高。
-
2.4 模拟结果分析
-
采用法国石油研究院研发的 Dionisos软件开展沉积过程数值模拟。沉积过程是复杂的，获取的资料具有不完备性，模拟结果必然会存在一定的不确定性，因此，需要对模拟结果进行优选，重点检查模拟结果与实际地层厚度及岩性之间的关系。检查各沉积单元单井模拟地层厚度与实钻地层厚度的误差，当与局部井点误差较大时，对模拟地层沉降进行微调，直到与实钻井地层厚度误差控制在 10% 以内（表2），当进行大幅度调整时，需检查地层对比结果的等时性。当模拟结果与实钻井岩性有较大差异时，调整不同水深岩性发育及其产率，同时，考虑波浪能量的方向与大小，波浪能量作用较强的部位，颗粒灰岩发育，反之泥晶灰岩发育。
-
从模拟结果（图3）来看，Ⅰ，Ⅲ和Ⅴ小层主要发育颗粒灰岩和泥粒灰岩，Ⅱ，Ⅳ和Ⅵ小层主要发育含粒泥晶灰岩和泥晶灰岩，与实钻井岩性基本一致。用模拟结果定义沉积相，BU油田高部位发育生物礁，周缘发育碎屑滩，其中碎屑滩分为滩核、滩缘和开阔海。模拟结果还可以表征储层的非均质性，不同沉积时期岩性具有差异。对模拟结果做等时切片（图4），可以看出不同沉积时期沉积物的建造过程。从下向上依次为：第 1 期（距今 96～94 Ma），沉积时期水体较深时，生物礁不发育，以开阔海为主；第 2期（距今 94～93.8 Ma），水体短时间下降，礁体发育，在高部位形成生物礁，周缘发育小范围滩体；第 3 期（距今 93.8～93.1 Ma），水体上升，生物礁被淹没，仅在局部高点发育生物礁，大面积发育开阔海；第 4 期（距今 93.1～92.5 Ma），水体下降明显，大面积发育生物礁和滩体，第 5 期（距今 92.5～91.2 Ma），水体缓慢上升，礁体分布范围缩小，伴有周缘滩体和开阔海大面积发育；第 6 期（距今 91.2～90 Ma），水体大幅下降，生物礁大面积发育，开阔海不发育，是最好的储层形成时期。
-
图3 BU油田模拟结果与实钻井岩性划分对比
-
Fig.3 Comparison of simulated lithology with drilled data of BU Oilfield
-
图4 BU油田储层沉积过程数值模拟结果
-
Fig.4 Numerical simulation results of reservoir sedimentary process of BU Oilfield
-
3 模拟结果应用
-
三维地质建模是利用随机模拟算法对储层分布进行预测，已在中外众多油田获得良好的应用效果，其关键是等时约束建模和成因控制建模[11-16]。因此，以沉积过程数值模拟结果作为地质建模的约束条件，是一种新的建模思路。将模拟结果按沉积时间做切片，输出颗粒灰岩概率分布体，对每个时间切片所反映的颗粒灰岩概率分布在三维地质模型中进行赋值，并建立颗粒灰岩三维地质模型，以该模型为约束条件，以井点解释的颗粒灰岩含量数据为基础数据建立岩相概率分布模型（图5）。在此基础上，采用相控建模的策略，运用序贯高斯随机模拟算法，进行孔隙度和渗透率模型的建立，这种建模方法既遵循地质统计学原理，又符合沉积储层分布规律。根据该研究成果，BU油田注水开发方案采用底注顶采的模式，由于底部相对致密的泥晶灰岩发育高角度微裂缝，在纵向上沟通储层，当底部注水时，水沿着微裂缝向上流动，可将顶部颗粒灰岩中的油驱出，该油田已完成4口井方案实施，恢复了地层压力，缓解了局部压力亏空现象，在剩余油富集区，日产油量能力增幅近 50%，达 0.38×104 m3 / d，BU油田整体日产油能力提高近 19%，达 0.94×104 m3 /d。
-
图5 BU油田颗粒灰岩含量岩相概率分布模型剖面
-
Fig.5 Lithofacies probability distribution model profile of grain limestone content in BU Oilfield
-
4 结论
-
本文提出了一种基于沉积过程的碳酸盐岩储层模拟方法，用该方法获取的模拟结果作为约束条件，与序贯高斯随机模拟算法相结合，进行三维地质模型的建立，可以得到既符合沉积规律又遵循地质统计学特征的岩相模型。对碳酸盐岩台地沉积影响因素进行分析，认为可容纳空间变化、碳酸盐岩产率等环境参数对储层分布范围的影响较大。沉积过程数值模拟可以为储层分布提供预测模型，在少井及地震资料分辨率受限的情况下，可以提供定量地质模型。
-
参考文献
-
[1] 张义楷，王志松，史长林，等.伊拉克米桑油田碳酸盐岩储层成岩作用［J］.科学技术与工程，2016，16（5）：45-53． ZHANG Yikai，WANG Zhisong，SHI Changlin，et al.Carbonate reservoir characteristics and diagenesis in Missan oilfield，Iraq ［J］.Science Technology and Engineering，2016，16（5）：45-53．
-
[2] 张义楷，康安，闵小刚，等.伊拉克米桑油田群Mishrif组MB2 1段碳酸盐岩储层特征及成因［J］.石油实验地质，2016，38（3）：360-365． ZHANG Yikai，KANG An，MIN Xiaogang，et al.Characteristics and genesis of carbonate reservoirs in the Mishrif MB2 1 member in the Missan oil fields，Iraq［J］.Petroleum Geology & Experi⁃ ment，2016，38（3）：360-365.
-
[3] 何文军，刘敏珠，吴俊军，等.准噶尔盆地阜东斜坡阜19井区三叠系韭菜园子组沉积正演模拟［J］.油气地质与采收率，2018，25（6）：7-15.HE Wenjun，LIU Minzhu，WU Junjun，et al.Forward modeling of sedimentation in the Triassic Jiucaiyuanzi Formation in Well Fu19 area of the Fudong slope，Junggar Basin［J］.Petroleum Geol⁃ ogy and Recovery Efficiency，2018，25（6）：7-15.
-
[4] 王鹏飞，叶小明，霍春亮，等.BZ油田古近系储层沉积过程数值模拟［J］.断块油气田，2017，24（5）：604-607.WANG Pengfei，YE Xiaoming，HUO Chunliang，et al.Numerical simulation of sedimentary process for Paleogene reservoir in BZ oilfield［J］.Fault-Block Oil & Gas Field，2017，24（5）：604-607.
-
[5] 秦国省，邹存友，赵亮，等.基于沉积过程的储层表征研究进展及发展趋势［J］.科学技术与工程，2017，17（30）：110-120.QIN Guosheng，ZOU Cunyou，ZHAO Liang，et al.Research prog⁃ ress and development tendency of reservoir characterization based on sedimentary process［J］.Science Technology and Engi⁃ neering，2017，17（30）：110-120.
-
[6] 张文彪，段太忠，刘彦锋，等.综合沉积正演与多点地质统计模拟碳酸盐岩台地——以巴西Jupiter油田为例［J］.石油学报，2017，38（8）：925-934.ZHANG Wenbiao，DUAN Taizhong，LIU Yanfeng，et al.Integrat⁃ ed sedimentary forward modeling and multipoint geostatistics in carbonate platform simulation：a case study of Jupiter oilfield in Brazil［J］.Acta Petrolei Sinica，2017，38（8）：925-934.
-
[7] 张娟.沾化凹陷长堤地区沙一段下部生物灰岩发育模式及分布特征［J］.油气地质与采收率，2018，25（3）：36-42.ZHANG Juan.Development mode and distribution characteristics of biolimestone at the lower Es 1 in Changdi area of Zhanhua Sag ［J］.Petroleum Geology and Recovery Efficiency，2018，25（3）：36-42.
-
[8] 徐康，刘光祥，胡文瑄，等.塔中地区上奥陶统层序格架内礁滩体展布及沉积相［J］.大庆石油地质与开发，2018，37（2）：26-33.XU Kang，LIU Guangxiang，HU Wenxuan，et al.Distribution and sedimentary facies of the reef flats in Upper Ordovician sequence framework of Tazhong area［J］.Petroleum Geology & Oilfield De⁃ velopment in Daqing，2018，37（2）：26-33.
-
[9] 付育璞，郑强，庞雯，等.叙利亚O油田Shiranish组碳酸盐岩储层主控因素分析［J］.石油实验地质，2017，39（3）：355-361.FU Yupu，ZHENG Qiang，PANG Wen，et al.Carbonate reservoir controls in the Shiranish Formation of O oil field，Syria［J］.Petro⁃ leum Geology & Experiment，2017，39（3）：355-361.
-
[10] 王振洲，张春雷，高世臣.利用决策树方法识别复杂碳酸盐岩岩性——以苏里格气田苏东 41-33 区块为例［J］.油气地质与采收率，2017，24（6）：25-33.WANG Zhenzhou，ZHANG Chunlei，GAO Shichen.Lithology identification of complex carbonate rocks based on decision tree method：An example from Block Sudong41-33 in Sulige Gas Field ［J］.Petroleum Geology and Recovery Efficiency，2017，24（6）：25-33.
-
[11] 赵习，刘波，郭荣涛，等.储层表征技术及应用进展［J］.石油实验地质，2017，39（2）：287-294.ZHAO Xi，LIU Bo，GUO Rongtao，et al.Reservoir characterization and its application to development［J］.Petroleum Geology & Ex⁃ periment，2017，39（2）：287-294.
-
[12] 程媛，张冲，陈雨龙，等.基于压汞资料的碳酸盐岩储层渗透率预测模型——以扎纳若尔油田 KT-Ⅰ和 KT-Ⅱ含油层系灰岩储层为例［J］.油气地质与采收率，2017，24（3）：10-17.CHENG Yuan，ZHANG Chong，CHEN Yulong，et al.Permeability prediction model of carbonate reservoir based on mercury injec⁃ tion data：A case study of oil-bearing limestone reservoir in KT-Ⅰand KT-Ⅱmembers of Zahnanor oilfield［J］.Petroleum Geolo⁃ gy and Recovery Efficiency，2017，24（3）：10-17.
-
[13] 李滨.大庆长垣萨中中部断层区三维地质建模［J］.大庆石油地质与开发，2017，36（4）：68-72.LI Bin.3D geological modeling of Middle Sazhong Fault Block in Daqing Placanticline［J］.Petroleum Geology & Oilfield Develop⁃ ment in Daqing，2017，36（4）：68-72.
-
[14] 王鸣川，段太忠，杜秀娟，等.沉积相耦合岩石物理类型的孔隙型碳酸盐岩油藏建模方法［J］.石油实验地质，2018，40（2）：253-259.WANG Mingchuan，DUAN Taizhong，DU Xiujuan，et al.Geologi⁃ cal modeling method based on sedimentary facies coupled rock type for porous carbonate reservoirs［J］.Petroleum Geology & Ex⁃ periment，2018，40（2）：253-259.
-
[15] 高振南，霍春亮，罗成栋，等.双模迭代技术在裂缝性碳酸盐岩油藏中的应用［J］.特种油气藏，2018，25（3）：108-112.GAO Zhennan，HUO Chunliang，LUO Chengdong，et al.Applicia⁃ tion of dual-mode iterative technology in fractured carbonate res⁃ ervoir［J］.Special Oil & Gas Reservoirs，2018，25（3）：108-112.
-
[16] 刘江丽，徐美茹，吕学菊，等.地震约束鲕粒滩型碳酸盐岩储层三维地质建模［J］.中国石油勘探，2017，22（3）：63-70.LIU Jiangli，XU Meiru，LÜ Xueju，et al.3D geologic modeling of oolitic-beach carbonate reservoirs based on seismic constraint［J］.China Petroleum Exploration，2017，22（3）：63-70.
-
摘要
综合利用岩心、地震、测井及分析化验等资料，分析BU油田碳酸盐岩沉积环境的控制因素，开展储层沉积过程数值模拟研究，定量表征不同沉积时期储层展布范围及沉积演化规律。研究结果表明，碳酸盐岩储层沉积过程数值模拟的可容空间变化、碳酸盐岩产率等环境参数对储层分布范围的影响较大。BU油田Mishirif组MB2 1 段纵向上分为Ⅰ—Ⅵ共6个小层，高能沉积带与低能沉积带交替出现，其中高能沉积带岩性主要为颗粒灰岩，在顶部Ⅰ小层较发育，而低能沉积带岩性主要为泥晶灰岩，并有高角度微裂缝发育，在Ⅵ小层底部发育。该研究成果为BU油田底注顶采的注水开发方案提供地质依据，已完成4口井方案实施，恢复了地层压力，缓解了局部压力亏空现象，整体日产油能力提高近19%，达0.94×104 m3 /d。
Abstract
Based on the core，seismic，logging and experiment data，the control factors of sedimentary environment were an- alyzed，and the numerical simulation of reservoir sedimentary process was carried out to quantitatively characterize the res- ervoir distribution and evolution at different stages. The study results show that the initial parameters of numerical simula- tion of carbonate reservoir sedimentary process have a great influence on the reservoir distribution，and main parameters in- clude accommodating space change，carbonate productivity，etc. In the vertical direction of BU Oilfield，the MB2 1 Member of Mishirif Formation can be divided into 6 layers includingⅠ-Ⅵ. High-energy sedimentary belt and low-energy sedimen- tary belt appear alternately. The lithology of high-energy sedimentary belt is mainly grain limestone on the top of LayerⅠ， while the lithology of low-energy sedimentary belt is mainly micrite limestone with micro-fracture on the bottom of Layer Ⅵ. The research results provide the geologic basis for the water injection development scheme of bottom injection and top production in BU Oilfield. Four wells have been drilled to build up the formation pressure and to alleviate the local pressure deficit，and the daily production rate of the Oilfield has increased by nearly 19%，and increases to 0.94×104 m3 /d.
