微地震散射波结构成像
发布时间:2020-09-11 17:27
水力压裂是一种有效的非常规油气开发工程技术,通过在目标储层产生裂缝以改善低渗透性储层,可以更有效地从储层中采收油气。在水力压裂过程中,高压流体通过压裂井注入到储层中,地下应力状态发生改变,进而诱发微地震事件,微地震事件则可能被安置在地表或监测井中的地震台站记录下来。对检测到的微地震事件进行定位,并且根据其分布可以用来评估诱发裂缝的发育情况和储层改造体积。通常认为,检测到的事件越完整,对微地震事件的定位越准确,则越有利于准确评估地下裂缝网络的发育情况以及储层改造体积。但是由于微震信号能量弱,很难完备地检测到微震信号,而且在某些情况下裂缝的产生不一定会伴随微地震的产生,因此只根据微地震事件的分布来描述水力压裂造成的裂缝区域存在偏差。对于压裂产生的裂缝,它们对入射地震波来说是一个强散射体。基于此,本论文提出了一种表征压裂产生裂缝分布的新方法,即利用微地震事件波形中的散射波进行逆时偏移,从而直接对裂缝区域进行成像。对于多级水力压裂系统,前面压裂段诱发产生的裂缝,可以充当后续压裂产生微地震波场的强散射体,因此后段压裂诱发的微震事件的波形则可能包含前段裂缝引起的散射波。因此,对某段压裂的微地震波形进行逆时偏移成像,通过对波形中强散射波进行成像,则可以描述前段压裂诱发的裂缝以及天然存在的裂缝。基于一个井下微地震监测系统,发展了基于声波方程的井下微地震逆时偏移算法并基于正演模型进行了算法的验证。微地震逆时偏移的成像结果不仅会受到微地震监测系统的影响,也可能受到周围地下结构的成像干扰,从而无法较好的照明目标裂缝区域。因此,本论文提出基于染色算法的针对目标成像的微地震逆时偏移算法。通过对目标成像区域进行染色标记,可以创建一个仅与被标记的地下结构及其周围结构有关的虚波场,从而获得一个仅与被标记结构及其周围结构有关的成像结果,实现改善目标压裂段的成像照明情况,提高目标区域的成像效果。上述发展的微地震偏移成像算法基于声波波动方程,需要对台站接收到的波场进行反传和对震源产生的波场进行正传,通过P波或者S波的波形互相关为成像条件对地下结构进行成像,因此在某种程度上受震源函数和辐射花样的影响。为了消除震源对成像结果的影响,本论文进一步发展了基于弹性波动方程的微地震逆时成像算法,只对台站记录的波形沿着时间反方向传播,并对分离的P、S波场实施成像条件就可以实现。这种基于逆时成像的算法不需要知道微震的位置、震源函数和震源机制,可以确定微震的位置和结构。为了提高成像的分辨率,本论文提出分组成像条件来提高成像的能量聚焦。与自相关/互相关成像条件相比,分组成像条件能够较好的压制空间采样不足造成的成像噪声,也能更好的适应稀疏的台站分布。与几何平均成像条件相比,分组成像条件对台站数据的要求更低,不容易产生不稳定的成像结果,更适合实际地震数据处理。论文基于合成和实际数据,进行了算法的测试,发现基于分组成像条件的微地震成像算法,一方面可以对微震进行定位,另一方面也可以对结构进行成像。
【学位单位】:中国科学技术大学
【学位级别】:博士
【学位年份】:2018
【中图分类】:TE357.1
【部分图文】:
逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑图2.1二维井中微地震监测的观测系统与P波速度模型。P波背景速度为4200m/s,逡逑在X邋=邋1000邋m处有断裂带,断裂带设置为比背景速度低15%的低速异常。监测井位于X邋=逡逑445米处,井中安置36级检波器,在深度上自Z=975邋m分布到Z=2025邋m,检波器间距为逡逑30m,由白色三角形表示。8个微地震事件分别位于(850邋m,邋2120邋m),邋(850邋m,邋1975邋m),邋(850逡逑m,邋2260邋m),邋(862邋m,邋2048邋m),邋(862邋m,邋2170邋m),邋(830邋m,邋2107邋m),邋(843邋m,邋2041邋nr)和(842邋m,逡逑2235邋m),用白色星号表示。逡逑(a)逦(b)逡逑。.25逡逑30-5邋-f'逡逑I逡逑0逦0.3逦0.6逦0.9逦1.2逦1.5逦0逦5逦10逦15逦20逦25逦30逦35逡逑X(km)逦Receiver邋number逡逑图2.2位于(850m,邋1975m)的地震事件的P波波场快照与台站记录。(a)通过有限差逡逑分法正演P波波场在传播了邋t邋=邋0.1786s时的波场快照。(b)井中36级检波器的记录,时逡逑间采样间隔为0.595ms。模型设置如图2.1所示。逡逑13逡逑
X(km)逦Receiver邋number逡逑图2.2位于(850m,邋1975m)的地震事件的P波波场快照与台站记录。(a)通过有限差逡逑分法正演P波波场在传播了邋t邋=邋0.1786s时的波场快照。(b)井中36级检波器的记录,时逡逑间采样间隔为0.595ms。模型设置如图2.1所示。逡逑13逡逑
逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑图2.4叠后成像结果。将图2.3所示的8个事件的逆时偏移成像结果归一化并且叠加逡逑得到叠后成像结果。逡逑(a)逡逑0逦—r—,邋ii1逡逑0.6邋10-6逡逑0.9逦04逡逑:::B邋_邋n邋::逡逑0逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑(b)逡逑0邋卜二_*''—ni1逡逑o.3逦I邋°-8逡逑。6邋■。6逡逑!;:|^H::逡逑£卜逡逑壚邋T邋*逦I-0-4逡逑z-1?邋i"邋W-1^.6逡逑:;HIKhHI:8逡逑0逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑图2.5传感器分布在不同深度时逆时偏移成像结果的对比。(a)当36级传感器分布逡逑从Z=1575邋m到2625邋m时的8个事件的叠后成像结果;(b)当36级传感器分布从Z=775逡逑m到1825邋m时的8个事件的叠后成像结果。逡逑16逡逑
本文编号:2816967
【学位单位】:中国科学技术大学
【学位级别】:博士
【学位年份】:2018
【中图分类】:TE357.1
【部分图文】:
逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑图2.1二维井中微地震监测的观测系统与P波速度模型。P波背景速度为4200m/s,逡逑在X邋=邋1000邋m处有断裂带,断裂带设置为比背景速度低15%的低速异常。监测井位于X邋=逡逑445米处,井中安置36级检波器,在深度上自Z=975邋m分布到Z=2025邋m,检波器间距为逡逑30m,由白色三角形表示。8个微地震事件分别位于(850邋m,邋2120邋m),邋(850邋m,邋1975邋m),邋(850逡逑m,邋2260邋m),邋(862邋m,邋2048邋m),邋(862邋m,邋2170邋m),邋(830邋m,邋2107邋m),邋(843邋m,邋2041邋nr)和(842邋m,逡逑2235邋m),用白色星号表示。逡逑(a)逦(b)逡逑。.25逡逑30-5邋-f'逡逑I逡逑0逦0.3逦0.6逦0.9逦1.2逦1.5逦0逦5逦10逦15逦20逦25逦30逦35逡逑X(km)逦Receiver邋number逡逑图2.2位于(850m,邋1975m)的地震事件的P波波场快照与台站记录。(a)通过有限差逡逑分法正演P波波场在传播了邋t邋=邋0.1786s时的波场快照。(b)井中36级检波器的记录,时逡逑间采样间隔为0.595ms。模型设置如图2.1所示。逡逑13逡逑
X(km)逦Receiver邋number逡逑图2.2位于(850m,邋1975m)的地震事件的P波波场快照与台站记录。(a)通过有限差逡逑分法正演P波波场在传播了邋t邋=邋0.1786s时的波场快照。(b)井中36级检波器的记录,时逡逑间采样间隔为0.595ms。模型设置如图2.1所示。逡逑13逡逑
逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑图2.4叠后成像结果。将图2.3所示的8个事件的逆时偏移成像结果归一化并且叠加逡逑得到叠后成像结果。逡逑(a)逡逑0逦—r—,邋ii1逡逑0.6邋10-6逡逑0.9逦04逡逑:::B邋_邋n邋::逡逑0逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑(b)逡逑0邋卜二_*''—ni1逡逑o.3逦I邋°-8逡逑。6邋■。6逡逑!;:|^H::逡逑£卜逡逑壚邋T邋*逦I-0-4逡逑z-1?邋i"邋W-1^.6逡逑:;HIKhHI:8逡逑0逦0.3逦0.6逦0.9逦1.2逦1.5逡逑X(km)逡逑图2.5传感器分布在不同深度时逆时偏移成像结果的对比。(a)当36级传感器分布逡逑从Z=1575邋m到2625邋m时的8个事件的叠后成像结果;(b)当36级传感器分布从Z=775逡逑m到1825邋m时的8个事件的叠后成像结果。逡逑16逡逑
【参考文献】
相关期刊论文 前2条
1 孙文博;孙赞东;朱兴卉;;相对保幅的角度域VSP逆时偏移(英文)[J];Applied Geophysics;2011年02期
2 马德堂,朱光明;弹性波波场P波和S波分解的数值模拟[J];石油地球物理勘探;2003年05期
本文编号:2816967
本文链接:https://www.wllwen.com/kejilunwen/shiyounenyuanlunwen/2816967.html