基于微地震连续裂缝网络模型的SRV研究.pdf
《基于微地震连续裂缝网络模型的SRV研究.pdf》由会员分享,可在线阅读,更多相关《基于微地震连续裂缝网络模型的SRV研究.pdf(8页珍藏版)》请在咨信网上搜索。
1、 第 47 卷第 4 期物 探 与 化 探Vol.47,No.4 2023 年 8 月GEOPHYSICAL&GEOCHEMICAL EXPLORATION Aug.,2023doi:10.11720/wtyht.2023.1139李秋辰,陈冬,许文豪,等.基于微地震连续裂缝网络模型的 SRV 研究J.物探与化探,2023,47(4):1048-1055.http:/doi.org/10.11720/wtyht.2023.1139Li Q C,Chen D,Xu W H,et al.Determining stimulated reservoir volume based on the mic
2、roseismic continuous fracture network modelJ.Geophysical and Geochemical Exploration,2023,47(4):1048-1055.http:/doi.org/10.11720/wtyht.2023.1139基于微地震连续裂缝网络模型的 SRV 研究李秋辰1,陈冬2,许文豪1,易善鑫2,谢兴隆1,关俊朋2,崔芳姿1(1.中国地质调查局 水文地质环境地质调查中心,河北 保定 071051;2.江苏省地质调查研究院,江苏 南京 210008)摘 要:在干热岩的开发和改造阶段,需要对储层进行水力压裂改造,储层改造体积(s
3、timulated reservoir volume,SRV)是评价压裂效果的主要标准。微地震监测作为水力压裂的技术环节之一,可以对储层改造体积进行有效估算。本文探讨了基于微地震连续裂缝网络模型的 SRV 计算方法,以指导干热岩的水力压裂工作。首先,基于微震事件时空分布特征与多项震源参数进行连续裂缝网络建模;其次,提取裂缝网格长度,并选取适当的裂缝高度与裂缝宽度以计算储层改造体积;最后,选取某双井型干热岩微地震监测数据对该方法进行实际应用。实例应用结果表明,该方法可有效估算储层改造体积,为干热岩后续开发与改造提供依据。关键词:干热岩;微地震监测;连续裂缝网络;储层改造体积中图分类号:P631.
4、4 文献标识码:A 文章编号:1000-8918(2023)04-1048-08收稿日期:2022-06-11;修回日期:2023-04-07基金项目:江苏省碳达峰碳中和科技创新专项资金(重大科技示范 BE2022859)第一作者:李秋辰(1988-),男,硕士学位,工程师。Email:874033010 通讯作者:许文豪(1993-),男,硕士学位,工程师。Email:941734230 0 引言随着人类生存发展观念的逐渐转变,绿色低碳已经成为中国的能源战略之一。在此时代背景下,干热岩的勘探开发日渐受到重视1。微地震监测技术可以对干热岩储层改造过程中引起的声发射事件进行识别解析2,从而提取岩
5、石力学参数、进行储层应力分析3与刻画裂缝的空间展布信息并对其成像4-5,进而对储层改造体积(stimulated reser-voir volume,SRV)进行估算。鉴于上述特点,微地震监测已经成为干热岩开采过程中的关键配套技术6-7。干热岩的储层体积改造主要是通过泵入大量低黏度压裂液,形成由天然裂缝和水力裂缝组成的复杂网状裂缝结构8-9。SRV 的大小直接影响干热岩产能的高低10,因此前人对其进行了大量研究。2006 年,Mayerhofer 等11在研究页岩中的微地震事件与裂缝之间的关系时,首次提出了“油藏改造体积”的概念,2010 年验证了页岩气增产效果随着储层改造体积变大而愈发明显的
6、观点7。Warpinskif等12提出一种利用箱状体体积来近似估计 SRV 的方法,Baidurja 等13对该方法进行改进,利用凹包体体积来计算储层改造体积。上述两种基于微震事件包裹体模型的 SRV 计算方法虽然方便快捷,但往往低估了地下裂缝的复杂性,人们逐渐偏向选择基于裂缝网络建模的方法对储层改造体积进行评估10,14。Dershowitz 等15首次引出离散水力裂缝网络(discrete fracture network,DFN)建模的观点,Xu等16-17在此基础上加入流体与裂缝及裂缝与裂缝之间的相互作用,从而对 SRV 进行估算。Meyer等18-19通过模拟地下裂缝的扩展情况和支撑
7、剂在裂缝网络中的运移方式,利用 DFN 模型计算出储层改造体积。Alexandre 等20通过线状 DFN 建模,模拟出由射孔点向储层深处生长的裂缝网络。Seifol-lahi 等21基于微地震事件产生的裂缝面,给出了面状 DFN 模型的研究方法。上述几种 DFN 模型虽然考虑了裂缝的复杂性,4 期李秋辰等:基于微地震连续裂缝网络模型的 SRV 研究但裂缝参数往往是随机分布的,不符合裂缝的真实情况,且反演效率很低。本文在线状 DFN 模型的基础上,基于微震事件的时间空间分布特征进行裂缝网络模型重构,并在此基础上加入微震事件的震源机制解与地应力场分析结果进行约束,形成连续裂缝网络(continu
8、ous fracture network,CFN),对 SRV进行合理估算。最后将该方法运用到某干热岩双井微震监测中,取得良好的实际应用效果。1 算法流程1.1 CFN 模型建立流程在对微震事件进行基本处理(时间空间定位、矩震级计算、震源机制解及地应力反演)之后,可以获取事件的空间坐标、发震时刻、矩震级、震源机制解与局部地应力场等信息,对于所有给定的已经 包 含 上 述 信 息 的 微 震 事 件 点e x,y,z,t,m,s,d,r(),设置事件点集合 M。其中,x、y、z 为震源的空间坐标,t 为发震时刻,m 为矩震级,s、d、r 分别为事件破裂面的走向、倾角和滑移角。对 M 进行连续裂缝
9、网络建模由以下步骤组成10:1)按 t 由小到大的顺序,对 M 中的事件点进行排列。2)定义连续裂缝网络 N 的种子点。当监测到射孔事件时,定义种子点为射孔事件点,否则定义发震时刻最早的事件点为种子点。3)定义微震事件点 e x,y,z,t,m,s,d,r()和连续裂缝网络 N 之间的连接准则 d(e,N):微震事件点以最短路径与网络进行连接,事件可以与事件进行连接,也可以与路径进行连接(图 1)。种子点的距离定义为 0。图 1“事件点网络”连接准则Fig.1“Event-network”connection principle 4)对连接准则 d(e,N)进行约束,具体包括:时间约束、震源机
10、制约束、最长扩展距离约束与最大间距约束。其中,时间约束为在给定时间范围内进行路径连接;震源机制约束为使用震源机制反演的事件破裂面走向、倾角与滑移角对其约束;最长扩展距离约束为给定一个最长扩展距离,大于该距离后路径不再扩展;最大间距约束为给定一个最大间距(参考矩震级设定),大于该间距后事件不再进行路径连接。5)根据约束后的连接准则 d(e,N),确定事件点集 M 中的第 i 个微震事件点 ei至连续裂缝网络 N的连接点 c,生成 ei与 c 之间的路径 di。6)从事件点集 M 中移除 ei,将 ei与 di加入连续裂缝网络 N 后,在该时刻重定义 N。7)增量从 i 至 i+1,循环 3)6)
11、步骤,直至事件点集 M 中没有事件点。至此,完成连续裂缝网络建模。1.2 目标函数设置目标函数一为10,22:f1w()=mi=1d2i,j,(1)j=acgmin d2ei,Fj(),(2)di,j=Axx+Ayy+Azz+M()A2x+A2y+A2z,e(Fj)iFj;+,其他。(3)式中:w 表示等待优化的裂缝变量;m 表示事件点的数量;i、j 表示第 i 个点和第 j 条裂缝;di,j表示微地震事件 ei到其对应的裂缝椭圆面 Fj的垂直投影距离;A 为 m 3 的矩阵;M 为事件点集。该目标函数主要通过所有事件点到裂缝面的最短距离之和,来约束裂缝的数量、中心位置、倾向与倾角,从而确定整
12、个裂缝的空间位置。设置目标函数二为:f2w()=nj=1+ajbj1+mj,(4)式中:n 表示裂缝数量;aj、bj分别表示第 j 条裂缝的长短半轴;为正常数;mj表示与第 j 条裂缝相关联的微地震事件数量,该目标函数通过使建模过程中接受更小的裂缝面积和更多的相关联微地震事件点,来约束裂缝的数量、旋转角与裂缝椭圆面的长短半轴,从而确定裂缝的大小和方位。设置目标函数三为:f3w()=nj=1mji=1(j)i-j2,(5)9401物 探 与 化 探47 卷 式中:n 表示裂缝数量;mj表示与第 j 条裂缝相关联的微地震事件数量;(j)i表示第 j 条裂缝上的第 i 个事件点的倾向与倾角;j是第
13、j 条裂缝的倾向与倾角。该目标函数约束了裂缝上事件点与该条裂缝的产状信息的偏差,从而确保生成的裂缝的产状信息与该条裂缝上事件点对应的信息一致。最终目标函数为:f w()=1f1w()+2f2w()+3f3w(),(6)式中:1、2与 3分别表示相应子目标函数的权重因子。1.3 SRV 计算建立连续裂缝网络模型之后,可统计裂缝分支指数 与 裂 缝 网 络 长 度 等 参 数。裂 缝 分 支 指 数(branch index,BI)表示裂缝网络的复杂程度,将射孔点延伸出去的裂缝视为主缝,其分支指数为 1;直接从主缝分支出去的裂缝,其分支指数为 2;后续分支裂缝的分支指数依次增加 1。裂缝网络长度表
14、示裂缝网络中所有裂缝的累积长度,即所有单条裂缝长度的加和。其算法流程见图 2。图 2 基于连续裂缝网络模型的 SRV 算法流程Fig.2 Flow chart of SRV algorithm based oncontinuous fracture network model 在本文计算 SRV 时,认为震源半径为裂缝高度,震源半径可由震源机制解获得;认为裂缝影响半径为裂缝宽度,裂缝影响半径根据岩石模拟实验获得。结合 CFN 的统计结果,在对计算空间网格化之后,可获得储层改造体积如下所示23:SRV=klkrkwk,(7)式中:l 为第 k 个网格的裂缝长度;r 为第 k 个网格的裂缝高度;w
15、 为第 k 个网格的裂缝宽度;k 为裂缝的分布范围。2 实际应用分析2.1 事件定位本实例从某干热岩工区水力压裂过程的 4 621个双井微震监测事件当中,选取了矩震级 Mw-2且具有高信噪比的 4 180 个微震事件,采用本文所述方法对该事件集进行研究。首先结合工区地质资料与射孔校正速度模型(图 3),对事件集进行定位,并给出微震事件点的时间空间分布(图 4),图 4 中带色圆球表示微震事件的时空演化过程。空间分布显示:微震事件主要分布于水平 x 向-500350 m、水平 y 向-400500 m、地下 z 向-4 810-4 650 m 之间。事件形态为近椭圆形,主轴为近 NW 向,短轴为
16、近 SE 向。定位结果的精度对连续裂缝网络建模的准确性有很大影响,该压裂段使用双井监测,相对于单井监测,定位目标函数更加收敛,且能有效解决多解性问题。选取某一微震事件,其定位目标函数等值线如图 5 所示,图 5a 与图 5b 分别表示单井定位与双井定位目标函数等值线。图 3 射孔校正速度模型Fig.3 Model of perforation correction velocity0501 4 期李秋辰等:基于微地震连续裂缝网络模型的 SRV 研究axy 向平面;bxz 向剖面axy plan;bxz section图 4 微震事件时间空间分布Fig.4 Spatio-temporal dis
17、tribution of microseismic eventsa单井定位目标函数等值线;b双井定位目标函数等值线acontour map of positioning objective function of single well;bcontour map of positioning objective function of double well图 5 定位目标函数等值线Fig.5 Contour map of positioning objective function1501物 探 与 化 探47 卷 2.2 震源机制解、矩震级与地应力分析为获取事件的更多震源参数信息,对实例中
18、的全部事件的震源机制、矩震级与地应力场进行反演(图 6)。反演结果显示:基于震源机制的裂缝方位为近 NW 向,裂缝倾角以高角度为主;垂向应力为该区段最大主应力,最大水平主应力为近 NW 向,最小水平主应力为近 SE 向;矩震级集中分布在 0 级以下。a震源机制解玫瑰图;b应力主轴置信区间图;c矩震级统计直方图afocal mechanism solution rose diagram;bconfidence of principal stress axis;cstatistical histogram of moment magnitude图 6 微震事件震源参数反演结果Fig.6 Inver
19、sion results of source parameters of microseismic events2.3 连续裂缝网络(CFN)建模在获取事件集的空间坐标、发震时刻、矩震级、震源机制解及地应力分析等信息后,采用本文所述方法,以射孔事件为种子点、最大扩展间距为 700 m、最长间距为 75 m,同时加以震源参数约束进行该压裂段的连续裂缝网络建模,最终给出连续裂缝网络模型的时间空间分布(图 7)。模型显示:裂缝水平方向上呈近椭圆形分布,长轴为近 NW 向,短轴为近 SE 向;垂直方向上集中分布在-4 800-4 650 m 之间;形成时间以射孔点为原点、呈均匀向远端扩散的规律。上述特
20、征与事件时间空间特征、震源机制解与地应力分析的结果吻合。裂缝网络长度统计结果为 336.6103 m,裂缝分支指数为 5。2.4 储层改造体积(SRV)计算在该压裂段 CFN 模型的基础上,依据震源机制解中的震源半径参数,裂缝高度取值为 5 m;参考相关岩石模拟实验的经验结果,裂缝影响半径取值为3 m。利用式(1)对该压裂段的储层改造体积进行计算,并给出储层改造体积示意(图 8)。经计算,储层改造体积为 504.9104 m3,证明该压裂段水力压裂效果较好,可部署下一阶段工作。3 结论本文探讨了基于连续网络裂缝模型计算 SRV的方法:利用微震事件的时间空间参数进行线状DFN 模型重建,与此同时
21、加入微震事件的多项震源参数进行约束,形成连续裂缝网络模型。参考模型统计结果设定裂缝网格长度,参考震源半径设定缝高,参考裂缝影响半径设定缝宽,对 SRV 进行合理估算,进而对后续井位的布设、水力压裂及储层改造的效果评价提供依据。利用本文所述方法对某双井型干热岩微震事件进行实例分析,结果表明:基于微震事件的时空分布,依据“事件点网络”的连接准则与事件属性参数建立的连续裂缝网络,不仅可以对水力裂缝网络进行解释,在本文中同时也是储层改造体积计算的基础。基于连续裂缝网络模型,根据震源半径与2501 4 期李秋辰等:基于微地震连续裂缝网络模型的 SRV 研究图 7 连续裂缝网络(CFN)模型Fig.7 C
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 地震 连续 裂缝 网络 模型 SRV 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。