365体育官网入口

  • (010)84657853 84658665 84657900
  • mt@aolipumps.org
  • 北京市朝阳区芍药居35号中煤信息大厦(100029)

365体育官网入口

★ 生孩子开发建设 ★

地质统计学反演在淮南地区致密砂岩气储层的地质建模研究

刘 俊1,崔 凡(fan)2,王 瑞3,刘 波4,陈柏平2

(1.淮河生物质能控股企业集团电话煤业分集团厂家,湖南省淮南市,232000;2.中矿山机械师范大学本科(济南),济南市海淀区,100083;3.中是由化工环保股不多集团厂家是由探矿研发深入分析院,济南市海淀区,100083;4.等等工业师范大学本科,等等省邯郸市,056038)

摘(zhai) 要 针对(dui)致密(mi)(mi)砂(sha)(sha)岩(yan)(yan)气(qi)储(chu)层(ceng)(ceng)(ceng)较(jiao)薄(bo)(bo)(bo)、砂(sha)(sha)泥互(hu)层(ceng)(ceng)(ceng)地(di)(di)(di)质(zhi)(zhi)现象(xiang)普(pu)遍,常规(gui)的(de)(de)地(di)(di)(di)震反演方(fang)法无法精细分辨薄(bo)(bo)(bo)储(chu)层(ceng)(ceng)(ceng)空(kong)间展布(bu)情(qing)况,利用基(ji)于马尔科夫链-蒙特卡洛算法的(de)(de)地(di)(di)(di)质(zhi)(zhi)统计学(xue)反演方(fang)法在薄(bo)(bo)(bo)储(chu)层(ceng)(ceng)(ceng)预(yu)测(ce)上(shang)(shang)的(de)(de)优势,预(yu)测(ce)了淮(huai)南(nan)地(di)(di)(di)区(qu)煤(mei)田的(de)(de)煤(mei)系地(di)(di)(di)层(ceng)(ceng)(ceng)致密(mi)(mi)砂(sha)(sha)岩(yan)(yan)气(qi)储(chu)层(ceng)(ceng)(ceng)岩(yan)(yan)性,并建立三维地(di)(di)(di)质(zhi)(zhi)模型。研究结果(guo)发(fa)现,研究区(qu)煤(mei)系地(di)(di)(di)层(ceng)(ceng)(ceng)内砂(sha)(sha)岩(yan)(yan)储(chu)层(ceng)(ceng)(ceng)大多较(jiao)薄(bo)(bo)(bo),且具有(you)砂(sha)(sha)泥互(hu)层(ceng)(ceng)(ceng)的(de)(de)特点,通(tong)过与实际测(ce)井资料对(dui)比,利用地(di)(di)(di)质(zhi)(zhi)统计学(xue)反演方(fang)法预(yu)测(ce)致密(mi)(mi)砂(sha)(sha)岩(yan)(yan)储(chu)层(ceng)(ceng)(ceng)厚度的(de)(de)相对(dui)误(wu)差(cha)(cha)小于17%,对(dui)于5 m以上(shang)(shang)的(de)(de)砂(sha)(sha)岩(yan)(yan)储(chu)层(ceng)(ceng)(ceng),相对(dui)误(wu)差(cha)(cha)控制在5%以内。反演结果(guo)可为预(yu)测(ce)致密(mi)(mi)砂(sha)(sha)岩(yan)(yan)气(qi)藏提(ti)供精细的(de)(de)岩(yan)(yan)层(ceng)(ceng)(ceng)分布(bu)信息,为致密(mi)(mi)砂(sha)(sha)岩(yan)(yan)气(qi)储(chu)层(ceng)(ceng)(ceng)“甜点”提(ti)供参考。

重(zhong)要词 地质(zhi)统计学反演(yan);致密(mi)砂岩气;地质(zhi)建模

0 引言

致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)作为一种非常规(gui)天然气(qi),是一种清洁的化石能源。20世纪80年(nian)代初,美国(guo)(guo)的致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)勘探(tan)开发(fa)(fa)技术获(huo)得重大(da)突破,之后(hou)开始快速发(fa)(fa)展起来。目前,致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)已在全球各国(guo)(guo)得到(dao)开发(fa)(fa)和(he)研究。我国(guo)(guo)盆地发(fa)(fa)育具有多(duo)期(qi)多(duo)类型的特点,致(zhi)(zhi)密(mi)储层(ceng)广泛发(fa)(fa)育,致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)储量较(jiao)丰富。近(jin)年(nian)来,我国(guo)(guo)致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)产量日益增长(zhang),2017年(nian)产量约1×108 m3,2019年(nian)产量约6×108 m3 [1]。我国(guo)(guo)致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)勘探(tan)和(he)开发(fa)(fa)尚处(chu)于(yu)发(fa)(fa)展时期(qi),如(ru)何(he)将地球物理探(tan)测手段与地质资料(liao)有机结合,准确预测致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)储层(ceng)位置及厚(hou)度对致(zhi)(zhi)密(mi)砂(sha)(sha)(sha)岩(yan)(yan)(yan)气(qi)开采和(he)利(li)用具有重要的意义。

区别的的堆积的环境享有区别的水趋势特征英文,所确立的砂体在岩相主成部分、板材厚度、内壁非均质性及其砂岩残渣碎屑成分表主成部分、泥质水分含量、粒状磨料粒度、分选性状等几角度各具独特,会造成区别的堆积的环境所确立的砂体享有区别的原有融入系数度和融入率。虽说成岩效应对原有融入系数度的更新改造不强,虽然成岩效应是在的堆积效应的基础性努力上进行的,因为前兆的成岩效应也遭遇的堆积的环境的影响力到,进而进这一步影响力到成岩效应的类别、屈服强度,对砂岩的融入系数演变起必须的操作效应。

沉积(ji)环境是致密(mi)(mi)砂(sha)岩气(qi)储(chu)(chu)(chu)层(ceng)(ceng)形成(cheng)的(de)(de)重(zhong)要因素(su)之一,不同的(de)(de)沉积(ji)环境影响了砂(sha)体(ti)的(de)(de)厚(hou)度、形态及空间展(zhan)布,了解致密(mi)(mi)砂(sha)岩气(qi)储(chu)(chu)(chu)层(ceng)(ceng)所处地质环境的(de)(de)岩性、储(chu)(chu)(chu)层(ceng)(ceng)厚(hou)度等(deng)信息能为致密(mi)(mi)砂(sha)岩气(qi)的(de)(de)预(yu)测提供有价值的(de)(de)参考。致密(mi)(mi)砂(sha)岩气(qi)储(chu)(chu)(chu)层(ceng)(ceng)一般具有纵向砂(sha)体(ti)厚(hou)度薄(bo)、砂(sha)泥互(hu)层(ceng)(ceng),横向砂(sha)体(ti)厚(hou)度变化快的(de)(de)特点[2]。目前,常用(yong)的(de)(de)约束(shu)系(xi)数(shu)脉冲(chong)反(fan)演(yan)三维地震(zhen)反(fan)演(yan)技术(shu)由(you)于(yu)地震(zhen)数(shu)据分(fen)辨率(lv)限(xian)制,无法对致密(mi)(mi)砂(sha)岩气(qi)储(chu)(chu)(chu)层(ceng)(ceng)进行精细刻(ke)画(hua)。利用(yong)曲(qu)线(xian)(xian)重(zhong)构技术(shu)虽然能进一步提高(gao)约束(shu)系(xi)数(shu)脉冲(chong)反(fan)演(yan)的(de)(de)精度,但仍旧无法突破约束(shu)系(xi)数(shu)脉冲(chong)反(fan)演(yan)的(de)(de)理论(lun)极限(xian),对超(chao)薄(bo)储(chu)(chu)(chu)层(ceng)(ceng)预(yu)测效果不理想,且曲(qu)线(xian)(xian)重(zhong)构技术(shu)目前仍存在争议(yi)[3]。

地质(zhi)统计(ji)学反演(yan)方法(fa)融合了测井纵向分辨率高和(he)三维地震横向分辨率高的(de)(de)(de)优(you)点(dian),目前已被广(guang)泛应用到油气(qi)(qi)(qi)勘探领(ling)域的(de)(de)(de)储(chu)层(ceng)预测[4-7]。针(zhen)对(dui)淮南致(zhi)密(mi)砂(sha)岩(yan)气(qi)(qi)(qi)储(chu)层(ceng)特点(dian),笔者(zhe)采用地质(zhi)统计(ji)学反演(yan)方法(fa)对(dui)致(zhi)密(mi)砂(sha)岩(yan)气(qi)(qi)(qi)储(chu)层(ceng)所(suo)在地层(ceng)进(jin)行反演(yan),并构建(jian)地质(zhi)模型(xing),一方面可为致(zhi)密(mi)砂(sha)岩(yan)气(qi)(qi)(qi)储(chu)层(ceng)预测提供(gong)精(jing)细的(de)(de)(de)地质(zhi)背景信息,进(jin)一步(bu)保证(zheng)致(zhi)密(mi)砂(sha)岩(yan)气(qi)(qi)(qi)“甜点(dian)”预测的(de)(de)(de)准确性(xing),并促进(jin)致(zhi)密(mi)砂(sha)岩(yan)气(qi)(qi)(qi)资源的(de)(de)(de)高效开发和(he)利用,另一方面对(dui)预测其他非常规(gui)天然气(qi)(qi)(qi)并建(jian)立精(jing)确地质(zhi)模型(xing)提供(gong)借鉴。

1 研究区地质和原始数据概况

研(yan)究区(qu)(qu)位于安徽省凤台县张(zhang)集(ji)煤矿某采(cai)区(qu)(qu)内,经钻探揭露研(yan)究区(qu)(qu)内地层(ceng)(ceng)(ceng)(ceng)有(you)寒(han)武系(xi)(xi)、奥(ao)陶(tao)系(xi)(xi)、石炭系(xi)(xi)、二叠系(xi)(xi)、新近系(xi)(xi)、第四系(xi)(xi)地层(ceng)(ceng)(ceng)(ceng)。致密砂岩(yan)(yan)储层(ceng)(ceng)(ceng)(ceng)一(yi)般(ban)位于二叠系(xi)(xi)的(de)山西组和石河子组内的(de)煤系(xi)(xi)地层(ceng)(ceng)(ceng)(ceng)中(zhong)。煤系(xi)(xi)地层(ceng)(ceng)(ceng)(ceng)主(zhu)(zhu)要有(you)砂岩(yan)(yan)、泥岩(yan)(yan),以及1、6、8、11-2、13-1号共5个主(zhu)(zhu)要可(ke)采(cai)煤层(ceng)(ceng)(ceng)(ceng)。煤层(ceng)(ceng)(ceng)(ceng)厚度为(wei)0~10.74 m,工作面瓦斯绝(jue)对(dui)涌出量为(wei)1.300 8~9.680 0 m3/min。

还包括由于叠后立体立体强震总计数剧数剧资料、测井总计数剧数剧资料(还包括岩性和纵波特性阻抗匹配)开设的岩土工程勘察总计学反演。强震总计数剧数剧资料在反演过程中 中起还包括的作用,反演以后报告需用忠于自己强震总计数剧数剧资料,测井数剧资料还包括为反演提供了反演叁数的倡导知识基础总计数剧数剧资料和反演帮助先决条件,以后使反演以后报告在横排可提供强震总计数剧数剧资料的展布大趋势,在垂向可提供测井相关信息的表现形式。立体立体强震叠后总计数剧数剧资料剖面和几率段领域如1如图甲如图所示,叠后强震总计数剧数剧资料的几率网络带宽为18~78 Hz,主频约为45 Hz。钻井队岩性和纵波特性阻抗匹配总计如2如图甲如图所示,从图例可发现,煤系的地层中砂岩储层较薄,最薄砂岩储层厚薄不能1 m,且发生砂泥互层现状,这也增重了强震反演的難度。反演常用反演子波为主导频50 Hz的雷克子波。

2 地质统计学反演方法原理

地质统计(ji)学反(fan)(fan)演(yan)(yan)方法是(shi)基于贝叶斯(si)定理的(de)(de)反(fan)(fan)演(yan)(yan)框架,通过已知三(san)维地震、测井数据(ju)等先验信息统计(ji)分析并获(huo)取(qu)概率(lv)(lv)密度(du)函(han)(han)数、变差函(han)(han)数等反(fan)(fan)演(yan)(yan)参(can)数来反(fan)(fan)演(yan)(yan)在已知地质、测井、地震等信息的(de)(de)条件下储层的(de)(de)岩(yan)(yan)性(xing)或物(wu)性(xing)参(can)数,由于反(fan)(fan)演(yan)(yan)过程中涉及岩(yan)(yan)性(xing)、纵(zong)波阻抗等多维度(du)参(can)数,计(ji)算十分困难,目前,较(jiao)为有效的(de)(de)方法是(shi)采用马(ma)尔科夫链-蒙(meng)特卡(ka)洛算法构建马(ma)尔科夫链进行随(sui)机采样,反(fan)(fan)演(yan)(yan)过程中根(gen)据(ju)概率(lv)(lv)密度(du)函(han)(han)数、变差函(han)(han)数等反(fan)(fan)演(yan)(yan)参(can)数进行模拟计(ji)算,最终可获(huo)得岩(yan)(yan)性(xing)、纵(zong)波阻抗等各种高分辨率(lv)(lv)数据(ju)[8-10]。

图1 3D地震灾害所造成的叠后数据显示剖面和频宽範圍
图2 钻孔岩性和纵波阻抗匹配总计

贝叶斯反演(yan)框架是一种利用概(gai)(gai)率(lv)统(tong)计(ji)原(yuan)理对(dui)多种信息(xi)进(jin)行(xing)整合(he)并判别的(de)方法,能(neng)够公(gong)平地(di)将(jiang)多种不确定(ding)信息(xi)源进(jin)行(xing)统(tong)一分(fen)析(xi)判断,常被用来计(ji)算在(zai)已知(zhi)先验(yan)信息(xi)概(gai)(gai)率(lv)(E)和假设条件概(gai)(gai)率(lv)(H)的(de)前提下,某一事件(X)发生的(de)后验(yan)概(gai)(gai)率(lv)[11-13],其公(gong)式(shi)如下:

(1)

式(shi)中:P(X|H,E)——待求解的目(mu)标后验概率,即最终的纵波阻(zu)抗、岩性参数反演模型;

P(X|H)——在已(yi)知(zhi)假设条件(jian)(H)下事(shi)件(jian)X的先验概率分布(bu),即通过测井(jing)数据获取的岩(yan)性纵波阻抗(kang)、岩(yan)性等信息(xi)得到的预测性岩(yan)性或(huo)弹(dan)性参数的先验信息(xi);

P(E|X)——事件X与已知验证(E)之间的(de)相似(si)性,即(ji)依(yi)赖三维地震数(shu)据(ju)和测井资料的(de)似(si)然概率密度函(han)数(shu)[14];

P(H|E)——归一化参量,表示满(man)足(zu)真(zhen)实的地震(zhen)、测(ce)井等信(xin)息和通过反(fan)演参数(shu)分析得到的预测(ce)信(xin)息,可(ke)近(jin)似看(kan)作全概(gai)率事件(jian)(P(H|E)~1)。

通(tong)过(guo)地(di)质统计学(xue)反(fan)演参数(shu)(shu)求得各储层的(de)纵波(bo)阻抗zp,然后将其转化为反(fan)射(she)系数(shu)(shu)Rp,再与子波(bo)w进(jin)行褶积计算,生成合成记录(lu)sy,其与实际地(di)震(zhen)记录(lu)s之间的(de)差异可以(yi)用(yong)噪(zao)音水平来衡量。反(fan)射(she)系数(shu)(shu)和纵波(bo)阻抗的(de)计算公(gong)式如下:

(2)

式中:ρ1和(he)ρ2——上下储层的密度,kg/m3;

v1和v2——上下储层的速度,m/s;

zp——纵波阻(zu)抗(kang),(kg/m3)·(m/s);

ρ——储层(ceng)的密度,kg/m3;

v——储层的(de)速(su)度,m/s。

褶积测算流程有以下:

sy=w(t)×Rp(t)

(3)

式中:t——时间(jian),s。

这些整换算进程能否觉得为计算方式(4):

P(zp|v,s)∝P(zp|v)P(s|sy(zp))

(4)

式(4)表(biao)示后(hou)验概(gai)(gai)率(lv)密(mi)度(du)(du)函(han)数(shu)为已(yi)知(zhi)变差函(han)数(shu)v、合成(cheng)记录(lu)sy以及地震数(shu)据s的(de)(de)前提条件下,求(qiu)解(jie)纵(zong)波阻(zu)抗zp的(de)(de)概(gai)(gai)率(lv)。这可(ke)以看作求(qiu)解(jie)所有(you)输入(ru)概(gai)(gai)率(lv)密(mi)度(du)(du)函(han)数(shu)(纵(zong)波阻(zu)抗、岩(yan)性)的(de)(de)交集(ji)部(bu)分,所以地质统计学反(fan)演(yan)是(shi)一个(ge)真正意义(yi)上的(de)(de)同时(shi)反(fan)演(yan)算(suan)法(fa)(fa),同时(shi)考(kao)虑所有(you)的(de)(de)已(yi)知(zhi)信息源的(de)(de)影响。求(qiu)解(jie)出各(ge)个(ge)输入(ru)概(gai)(gai)率(lv)密(mi)度(du)(du)函(han)数(shu)的(de)(de)交集(ji)空间,即后(hou)验概(gai)(gai)率(lv)分布(bu)后(hou),由于未知(zhi)参数(shu)的(de)(de)后(hou)验概(gai)(gai)率(lv)复杂、维(wei)度(du)(du)较高,需(xu)要利用马尔科夫链-蒙特卡洛算(suan)法(fa)(fa)从中获取统计意义(yi)上的(de)(de)公平样本(ben)集(ji)[15-16]。

3 致密砂岩储层地质统计学反应实际应用

在事实用中,质地核算表学反演时候主耍包扩的数据显示表格报告分析开始筹备 、反演性能拟合曲线、质地核算表学反演3个方法。里面,的数据显示表格报告分析开始筹备 必备的数据显示表格报告分析包扩测井的数据显示表格报告分析(岩性、纵波阻抗匹配)和3d强震的数据显示表格报告分析。质地核算表学反演是源于贝叶斯反演架构设计,采取马尔科夫链-蒙特卡洛百度算法做出抽样检验、绘制和模以,并最终能够获取如果想要的储层性能或状态体的时候。 下面一开始争对反演关键,对设计区的叠后数据源多角度层位说明,并认定高密度砂岩储层所处的煤系岩层的反演位置,以作打造一开始的反演地层构架,如图甲如图所示3如图所示。
图3 反演地层知识体系剖面

地(di)(di)质统计(ji)(ji)学反演过程中(zhong),选取合适的(de)地(di)(di)质统计(ji)(ji)学反演参(can)数(shu)(shu)尤为(wei)(wei)重(zhong)要(yao),该参(can)数(shu)(shu)主要(yao)为(wei)(wei)概(gai)(gai)率(lv)密(mi)度函(han)数(shu)(shu)和变(bian)差(cha)函(han)数(shu)(shu)。概(gai)(gai)率(lv)密(mi)度函(han)数(shu)(shu)表示(shi)储层参(can)数(shu)(shu)或(huo)属(shu)性(xing)在空(kong)间上的(de)分布(bu)规律[15]。本次应(ying)用中(zhong),利用多元统计(ji)(ji)对(dui)井(jing)数(shu)(shu)据的(de)纵波(bo)阻抗(kang)和岩(yan)(yan)(yan)性(xing)进行统计(ji)(ji),获(huo)取纵波(bo)阻抗(kang)、岩(yan)(yan)(yan)性(xing)直方(fang)图(tu),并采用高斯函(han)数(shu)(shu)对(dui)纵波(bo)阻抗(kang)、岩(yan)(yan)(yan)性(xing)的(de)分布(bu)进行拟(ni)合,用拟(ni)合后(hou)的(de)均(jun)值、方(fang)差(cha)表示(shi)纵波(bo)阻抗(kang)、岩(yan)(yan)(yan)性(xing)的(de)分布(bu),泥(ni)岩(yan)(yan)(yan)拟(ni)合后(hou)的(de)概(gai)(gai)率(lv)密(mi)度函(han)数(shu)(shu)均(jun)值为(wei)(wei)9.41×106(kg/m3)·(m/s),方(fang)差(cha)为(wei)(wei)8.2×105;砂岩(yan)(yan)(yan)拟(ni)合后(hou)的(de)概(gai)(gai)率(lv)密(mi)度函(han)数(shu)(shu)均(jun)值为(wei)(wei)10.85×106(kg/m3)·(m/s),方(fang)差(cha)为(wei)(wei)5.1×105,如图(tu)4所(suo)示(shi)。

图4 煤系土层砂岩和泥岩储层的纵波阻抗匹配概率计算公式密度计算方程

变差(cha)(cha)函(han)(han)数(shu)(shu)是一种距(ju)离函(han)(han)数(shu)(shu),反映(ying)了地质统计学反演建模中(zhong)储层(ceng)参(can)数(shu)(shu)或属性(xing)在(zai)空间距(ju)离上的(de)连(lian)续性(xing),包括横向(xiang)和(he)(he)垂(chui)向(xiang)变差(cha)(cha)函(han)(han)数(shu)(shu)。变差(cha)(cha)函(han)(han)数(shu)(shu)具有一定的(de)主观性(xing),可(ke)(ke)采用(yong)指数(shu)(shu)函(han)(han)数(shu)(shu)或高(gao)斯(si)函(han)(han)数(shu)(shu)进行拟(ni)合得到。垂(chui)向(xiang)变差(cha)(cha)函(han)(han)数(shu)(shu)可(ke)(ke)根据井(jing)资(zi)料采样的(de)两点数(shu)(shu)据间的(de)协(xie)方(fang)差(cha)(cha)拟(ni)合得到,并在(zai)反演过程(cheng)中(zhong)对变差(cha)(cha)函(han)(han)数(shu)(shu)进行不断修改,最终的(de)参(can)数(shu)(shu)如图(tu)5所(suo)示,泥岩(yan)和(he)(he)砂(sha)岩(yan)的(de)变程(cheng)分(fen)别为25 ms和(he)(he)30 ms,公(gong)式(shi)(5)为两点数(shu)(shu)据间的(de)协(xie)方(fang)差(cha)(cha)和(he)(he)变程(cheng)的(de)计算公(gong)式(shi)[11]:

(5)

式中(zhong):i、j——上下采(cai)样(yang)点(dian);

σ——采(cai)样点所在储(chu)层(ceng)对应的概率密(mi)度函数的方差;

Δt——上下采样点的间隔时间;

R——变程。

图5 煤系岩层砂岩储层的下降函数值
水平垂直降低变量可从自然灾害均方根功能提升,其平均值为800 m。 效果,能够质地构造测算学反演的办法领取科研区煤系质地构造的岩性数据库体并创造出一个二维质地构造立体建模方法,图6和图7为二维质地构造立体建模方法的岩性立体建模方法和剖面图,表1为反演效果与盲井(未陆续参与反演的测井)各岩性厚薄测算效果(以紧挨煤层顶板这几天的砂岩和泥岩完成测算),整体研究分析能知质地构造测算学反演效果的导致精度较高,砂岩厚薄相较数据误差率为0.33%~16.67%,泥岩厚薄相较数据误差率为2.50%~14.29%,反演的各储层岩性与测井的资料基本的相符。
图6 顺利通过岩土工程勘察统计表格学反演步骤树立的3d煤系土层岩性整治
从反演毕竟的三维立体岩性体和剖面需要知道,要求研究分析区的煤系地壳构造基本有砂岩、泥岩和煤3种岩性,砂岩储层占基本个部分。在-6和13--6煤层之前基本有11套间隔砂岩储层,抗衡在很多不间隔砂岩储层,砂岩它的厚度面积为0~30 m。13--6和12-2号煤层之前砂岩储层较薄且基本与泥岩轮换出显,是非均质砂岩气藏的通常功能。12-2号和1-6煤层、-6和-6煤层之前都具备设计五套比很厚的泥岩,可对下端砂岩储层起到了“开放式”的功效,为非均质砂岩气具备了不错的储气地壳因素。-6和1-6煤层之前砂岩层比很厚,顶、底的泥岩层较薄,阻碍于储气,含气或许性较低。综上所述可以知道,非均质砂岩气最有弊于的存贮区域性为-6和-6煤层之前,该区域煤层比很厚、埋深较深,有有益于煤层气的导致,-6和-6煤层具备设计五套比很厚的泥岩还是有有益于导致非均质砂岩气存贮密闭式前景。
表1 地层汇总数学反演与测井事实上各岩性壁厚汇总数 煤层偏号砂岩体积尺寸实际存在/m反演/m相比较于不确定度/%泥岩体积尺寸实际存在/m反演/m相比较于不确定度/%13-5号8.28.42.438.58.83.5311-2号2.83.110.711.41.614.29八号2.22.513.602.02.514.2810号1.21.416.678.08.22.505号30.030.10.333.13.23.22
图7 水文地质计算学反演岩性剖面及安全验证

4 结论与建议

(1)地质学勘察汇总学反演方案就能对非均质砂岩气储层开展能够预计,其预计岩性和料厚报告单与测井实际情况文件不一性较高,而言料厚大的储层(>5 m),预计相对应的测量测量误差稳定制在5%时间内,料厚薄的储层(1~5 m),其精确性也较高,预计相对应的测量测量误差在17%时间内,证明地质学勘察汇总学反演方案在淮南的地方煤田的非均质砂岩气储层预计存在能行性。 (2)针对探究区煤系的地层内非均质砂岩储放的反演然而,预测出5号和-6煤层之前为非均质砂岩气储层“糕点”区域划分的将会性不大,意见重心在所在区域域划分内参与非均质砂岩气抽采试验台。 (3)鉴于少了抽采井数值,后期的进十步进十步探索,用含气砂岩和可含气砂岩开始岩相控住,有机会进十步加快地质学统计分析学反演方式方法在高密度砂岩气储层的予测正确性。

基(ji)准文(wen)(wen)献(xian)综述:

[1] 谢玉洪,蔡东升,孙晗森.国海油如此规气勘察搭建合一化探索性与预期效果[J].国石化勘察,2020,25(2):27-32. [2] 寇琳琳,孙现瑶,胡耀,等.低密度砂岩气藏调查构思及地质勘察措施浅谈[J].江西热,2019,46(10):148-149. [3] 蒋振会,齐颖.斜率抽象化反演在愈来愈规油藏储层竖向分析中的技术应用[J].海洋生物地层学术前沿,2016,32(7):64-70. [4] 张国华,张会星.地壳学统计数据学反演与常规性反演在储层予测中的采用视觉效果相比较讲解[J].我国的焦炭地壳学,2015,27(4):57-62. [5] 张国华.地质勘察统计表格学反演还有其在储层预測中的app[D].成都:国内 大海学校,2015. [6] 王伟方.水文地质统计数据学反演在储层預測中的软件应用[J].中州传统的煤炭,2016(8):137-140,144. [7] 陈顺顺,姜子振,刘媛媛.岩土工程勘察数据汇总学反表演功夫术在红岗北扶余油藏中的用[J].关中探铁路工人程,2017,29(3):51-53,57. [8] 郝立亚.依托于蒙特卡洛虚拟的贝叶斯个数价格波动类别及app研究分析[D].威海:湖北大学考研,2011. [9] 郑进京,朱慧明.研究背景MCMC的方式的贝叶斯AR(p)模型工具解析[J].统计学与策略,2005(20):6-8. [10] ARILD B, ODD K, RAGNAR H, et al.Bayesian lithology and fluid prediction from seismic prestack data[J].Society of Exploration Geophysicists, 2008, 73(3):13-21. [11] 王朋岩,李耀华,赵荣.叠后MCMC法岩性反演贝叶斯科研[J].宇宙物理性学发展,2015,30(4):1918-1925. [12] 宫秀军.贝叶斯学会认识论名词解释应用学科调查[D].北京市:中国现代学科院学科调查生院(计算的技能学科调查所), 2002. [13] 袁晶.贝叶斯形式在数据选原因中的采用[D].广州:山西本科大学,2013. [14] 郑卖瓜.贝叶斯金钱日期字段预侧型号及技术应用调查[D].广东药科高校:湖南省高校,2005. [15] 孙月成.为Bayesian-MCMC百度算法的地层核算学反演以及其在油藏模以中的用途[J].宇宙工具学进况,2018,33(2):724-729. [16] 朱嵩,毛根海,刘国华,等.提高效率的MCMC做法非常采用[J].水利水电学报,2009,40(8):1019-1023.

Research on geological modeling of tight sandstone gas reservoir in Huainan area by geostatistics inversion

LIU Jun1, CUI Fan2, Wang Rui3, LIU Bo4, CHEN Boping2

(1.Coal Industry Branch Company of Huaihe Energy (Group) Co., Ltd., Huainan, Anhui 232000,China;2.China University of Mining and Technology-Beijing, Haidian, Beijing 100083, China;3.China Petrochemical Corporation Petroleum Exploration and Development Research Institute, Haidian, Beijing 100083, China;4.Hebei University of Engineering, Handan, Hebei 056038, China)

Abstract In view of the fact that the tight sandstone gas reservoir was thin, the geological phenomenon of sand-mud interbedding was common, and the conventional seismic inversion method cuold not accurately distinguish the spatial distribution of the thin reservoir, using the advantage of the geostatistical inversion method based on Markov chain-Monte Carlo algorithm in thin reservoir prediction, the lithology of the tight sandstone gas reservoir of coal measures in Huainan area was predicted, and three dimensional geological model was established.The results showed that most of the sandstone reservoirs in the coal measure strata in the study area were thin, and have the characteristics of sand-mud interbedding.By comparing with the actual logging data, the relative error of the predicted thickness of tight sandstone gas reservoir by the geostatistical inversion method was less than 17%, and for the sandstone reservoirs more than 5 m, the relative error was less than 5%.The inversion results could provide fine strata distribution information for the prediction of tight sandstone gas reservoir, and provide a reference for the "sweet spot" of tight sandstone gas reservoir.

Key words geostatistics inversion; tight sandstone gas; geological modeling

中图类别(bie)号(hao) TD712

文献资料标志logo码 A

手机端二维码扫描查看

使用格试:刘俊,崔凡,王(wang)瑞,等.地质统计(ji)学反演在淮南地区致密砂岩气(qi)储(chu)层的(de)地质建模研(yan)究[J].中国煤炭,2021,47(5):7-12.doi:10.19880/j.cnki.ccm.2021.05.002

LIU Jun, CUI Fan, Wang Rui, et al.Research on geological modeling of tight sandstone gas reservoir in Huainan area by geostatistics inversion [J].China Coal,2021,47(5):7-12.doi:10.19880/j.cnki.ccm.2021.05.002

私募基(ji)金创业项目(mu):淮南矿(kuang)业集团科(ke)研项目(mu)(HNKYJTJS(2020)16号;HNKYJTJS(2018)181号)

我们(men)简绍:刘俊(1985-),汉族(zu),安(an)徽淮南人,本(ben)科学(xue)历,高级工程师(shi),主要从事矿(kuang)井地质方向研究。E-mail:zgkdlj@126.com

(损失小编 郭东芝)