CN110764139B - 一种各向异性纵横波走时高次幂线性组合定位方法 - Google Patents
一种各向异性纵横波走时高次幂线性组合定位方法 Download PDFInfo
- Publication number
- CN110764139B CN110764139B CN201810840579.0A CN201810840579A CN110764139B CN 110764139 B CN110764139 B CN 110764139B CN 201810840579 A CN201810840579 A CN 201810840579A CN 110764139 B CN110764139 B CN 110764139B
- Authority
- CN
- China
- Prior art keywords
- event
- opj
- shoot
- positioning
- longitudinal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/129—Source location
Abstract
本发明提供一种各向异性纵横波走时高次幂线性组合定位方法。该方法利用每个纵横波走时的定位特点,将纵波走时定位方程、横波走时定位方程、纵横波走时差值定位方程三者按照定位贡献的高低进行高次幂线性组合,建立目标方程;然后利用目标方程,开展射孔定位分析,获得射孔各向异性参数;然后利用射孔各向异性参数,对微地震事件做初步定位分析,获得微地震事件初始定位结果;最后,再利用目标方程,在微地震事件初始定位结果的基础上做进一步的层析定位分析,通过不断地修正各向异性参数,反演出定位精度更高的微地震事件定位结果。
Description
技术领域
本发明属于井中微地震监测信号处理技术领域,尤其涉及一种各向异性纵横波走时高次幂线性组合定位方法。
背景技术
微地震压裂监测技术是非常规致密砂岩气、页岩气藏储层油气田开发中一项关键技术之一,根据反演定位出的震源信息,可以获得裂缝属性(主应力走向、裂缝宽度、密度等),用来评价压裂效果,分析裂缝诱发规律,优化布井等。因此,在微地震信号处理中,最终目的是震源定位,亦称为微地震信号处理最核心技术。
井中微地震监测是微地震观测方式之一,其特点是利用井下三分量检波器来接收微地震全波场信号。相对于地面微地震监测,井中接收到的数据信噪比较高、微地震事件个数与类型较丰富。但是,与地面成百上千个检波器网状监测不同,由于井中微地震检波器个数有限(一般12至32级三分量井中检波器),通常井中微地震检波器放于垂直井段且检波器之间的间距一般为10米。也就是说,井中微地震监测微地震事件视角范围非常小,利用常规定位方法容易出现不稳定、精度不高的微地震定位结果。为了解决这一技术问题,需开发定位精度更高的新的定位方法。
目前,井中微地震定位方法主要有两种思路:一是基于P波、S波事件旅行时正演,代表算法有网络搜索法、模拟退火法、geiger法等,其优点是容易实现,其缺点是由于初至相位信号弱导致P波、S波事件旅行时难以准确拾取,影响定位结果;二是基于波动方程褶积,代表算法有干涉法、逆时偏移法、被动源成像法,其优点是不需要拾取事件初至,其缺点是对资料信噪比、速度模型要求高、检波器个数要求较多数量,计算成本高。
以上方法中,基于走时定位方法简单、实用且基于各向同性均匀介质假设,常用于实际井中微地震事件定位处理。但是随着非常规致密砂岩气、页岩气藏储层压裂微地震开发,地层具有各向异性非均质性,微地震纵横波走时及传播路径不同于以往各向同性,使用现有的定位方法不能满足其定位精度的要求。
发明内容
为了解决上述技术问题,本发明提出了一种新的各向异性介质纵横波走时高次幂线性组合定位方法,考虑地层各向异性对纵横波走时影响,由此消除或减少微地震事件定位的不稳定性、不准确性,以满足煤层气、页岩气等非常规微地震监测需要。该方法主要包括以下步骤:
S1:将各向异性纵波走时定位方程OPJP、各向异性横波走时定位方程OPJS和各向异性纵横波走时时差定位方程OPJPS按照定位贡献的高低进行高次幂线性组合,建立目标方程OPJ;
S2:拾取射孔的纵横波走时Tshoot,P、Tshoot,S作为观测值,结合已知的射孔空间位置和声波测井纵横波速度,利用所述目标方程OPJ对射孔开展定位分析,反演出射孔各向异性参数Δ={ε,δ},其中ε、δ为描述射孔各向异性的两个参数;
S3:拾取微地震事件纵横波走时Tevent,P、Tevent,S作为观测值,将反演出的射孔各向异性参数Δ={ε,δ}作为微地震事件初始各向异性参数,结合已知的声波测井纵横波速度,利用所述目标方程OPJ对微地震事件开展定位分析,获得微地震事件初始定位结果RTevent={Levent,Zevent},其中Levent、Zevent为描述微地震事件空间位置的径向坐标、深度坐标;
S4:以微地震事件初始定位结果为空间约束,以微地震事件初始各向异性参数为参数约束,在扰动范围内,利用所述目标方程OPJ对微地震事件做进一步的层析定位处理,获得定位精度更高的微地震事件最终定位结果。
根据本发明的实施例,上述步骤S4获得的微地震事件最终定位结果需要满足以下条件:根据步骤S4获得的微地震事件最终定位结果确定的微地震事件纵横波走时与步骤S3拾取的微地震事件纵横波走时之间的误差达到最小。
根据本发明的实施例,上述目标方程OPJ可以是:
OPJ=OPJ4 PS+OPJ3 PS+OPJ3 P+OPJ2 P+OPJ2 S+OPJS。
根据本发明的实施例,上述步骤S2具体包括以下步骤:
将拾取的射孔的纵横波走时Tshoot,P、Tshoot,S、已知的射孔空间位置和声波测井纵横波速度代入所述目标方程OPJ,获得射孔目标方程OPJshoot:
OPJshoot=OPJ4 shoot,PS+OPJ3 shoot,PS+OPJ3 shoot,P+OPJ2 shoot,P+OPJ2 shoot,S+OPJshoot,S
OPJshoot,P=|Tshoot,P-Tshoot,Pi|
OPJshoot,S=|Tshoot,S-Tshoot,Si|
OPJshoot,PS=|(Tshoot,P-Tshoot,S)-(Tshoot,Pi-Tshoot,Si)|
其中,OPJshoot,P为射孔各向异性纵波走时定位方程、OPJshoot,S为射孔各向异性横波走时定位方程、OPJshoot,PS为射孔各向异性纵横波时差定位方程,Tshoot,Pi、Tshoot,Si为反演的射孔纵横波走时;
反演射孔各向异性,利用以下偏导方程求解射孔各向异性参数Δ={ε,δ}:
根据本发明的实施例,求解出的射孔各向异性参数Δ={ε,δ}需要满足以下条件:求解出的射孔各向异性参数Δ={ε,δ}能够使得所述射孔目标方程OPJshoot的值达到最小。
根据本发明的实施例,上述步骤S3具体包括以下步骤:
将反演出的射孔各向异性参数作为微地震事件初始各向异性参数;
将拾取的微地震事件纵横波走时Tevent,P、Tevent,S,所述微地震事件初始各向异性参数和已知的声波测井纵横波速度代入所述目标方程OPJ,获得微地震事件目标方程OPJevent,
OPJevent=OPJ4 event,PS+OPJ3 event,PS+OPJ3 event,P+OPJ2 event,P+OPJ2 event,S+OPJevent,S
OPJevent,P=|Tevent,P-Tevent,Pi|
OPJevent,S=|Tevent,S-Tevent,Si|
OPJevent,PS=|(Tevent,P-Tevent,S)-(Tevent,Pi-Tevent,Si)|
其中,OPJevent,P为微地震各向异性纵波走时定位方程、OPJevent,S为微地震各向异性横波走时定位方程、OPJevent,PS为微地震各向异性纵横波时差定位方程,Tevent,Pi、Tevent,Si为反演的微地震事件纵横波走时;
反演微地震事件空间位置,利用以下偏导方程求解微地震事件初始定位结果RTevent={Levent,Zevent}:
根据本发明的实施例,求解出的微地震事件初始定位结果RTevent={Levent,Zevent}需要满足以下条件:求解出的微地震事件初始定位结果RTevent={Levent,Zevent}能够使得所述微地震事件目标方程OPJevent的值达到最小。
根据本发明的实施例,优选地利用最小二乘法或网格搜索法求解各偏导方程。
根据本发明的实施例,上述步骤S4具体包括以下步骤:
以微地震事件初始定位结果为空间位置中心,以微地震事件初始各向异性参数为参数中心,在扰动范围内,利用所述微地震事件目标方程OPJevent,对表征微地震事件空间位置的径向坐标、深度坐标以及微地震事件各向异性参数进行反演,利用以下偏导方程求解微地震事件最终定位结果:
根据本发明的实施例,优选地利用最小二乘法或网格搜索法求解步骤S4中的偏导方程,具体包括以下步骤:
在扰动范围内,针对每个微地震事件各向异性参数可能值,反演出与之对应的表征微地震事件空间位置的径向坐标、深度坐标,并据此推演出相应的微地震事件纵横波走时,计算推演出的微地震事件纵横波走时与拾取的微地震事件纵横波走时之间的误差;
寻找所有可能误差中的最小值,与所述最小值对应的径向坐标、深度坐标以及各向异性参数即为微地震事件最终定位结果以及相应的修正后的微地震事件各向异性参数。
与现有技术相比,本发明具有如下优点或有益效果:
本发明提供的一种基于各向异性介质的井中微地震定位方法,是一种井中微地震定位精度更高的、效果更好的微地震处理方法。首先,相比各向同性,本发明提出的定位方法考虑到地层各向异性,反演的纵横波走时更接近实际观测值。其次,本发明提出的定位方法利用P波走时、S波走时、PS波时差高次幂线性组合,与目前常见纵横波组合比较,能够充分利用纵横波定位敏感性,同时通过层析再定位,大幅度地提高了计算效率和定位精度。本发明提供的定位方法的简单易行、误差可控,能够为井中微地震定位处理提供可靠的技术保障。
本发明的其他优点、目标,和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书,权利要求书,以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本申请的技术方案或现有技术的进一步理解,并且构成说明书的一部分。其中,表达本申请实施例的附图与本申请的实施例一起用于解释本申请的技术方案,但并不构成对本申请技术方案的限制。
图1是本发明实施例一所采用的定位方法的工作原理图;
图2是本发明实施例二的井中微地震观测系统侧视图:射孔、检波器、事件;
图3是本发明实施例二的井中微地震模型事件空间坐标、各向异性参数的列表;
图4是本发明实施例二的射孔各向异性参数ε、δ的真实值和反演值的列表;
图5是基于图4反演出的射孔各向异性参数所获得的微地震事件初始定位结果的示意图;
图6是图5反演出的微地震事件初始定位结果与微地震事件真实空间位置的误差统计;
图7是基于图5反演出的微地震事件初始定位结果通过做进一步层析分析所获得的微地震事件最终定位结果的示意图;
图8是图7反演出的微地震事件最终定位结果与微地震事件真实空间位置的误差统计。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征,在不相冲突前提下可以相互结合,所形成的技术方案均在本发明的保护范围之内。
实施例一
本发明充分利用纵横波走时联合定位,提供了一种基于各向异性介质纵横波走时高次幂线性组合定位方法。图1示出了该定位方法的工作原理图。该方法大致可以划分为以下四个步骤:
S1,建立纵横波走时高次幂线性组合定位目标方程OPJ(以下简称目标方程OPJ),即:将各向异性纵波走时定位方程OPJP、各向异性横波走时定位方程OPJS、各向异性纵横波走时时差定位方程OPJPS三者根据定位贡献的高低进行高次幂复杂线性组合,建立目标方程OPJ;
S2,利用目标方程OPJ,通过射孔定位分析,获得射孔各向异性参数ε、δ;
S3,利用目标方程OPJ,采用S2获得的射孔各向异性参数ε、δ作为微地震事件初始各向异性参数,开展微地震事件初步定位分析,获得微地震事件初始定位结果RTevent={Levent,Zevent};
S4,利用目标方程OPJ,开展微地震事件层析定位分析,即:以微地震事件初始定位结果作为空间约束、以微地震事件初始各向异性参数为参数约束,设置相应的扰动范围,根据目标方程OPJ,反演微地震事件空间位置,期间不断修正各向异性参数ε、δ,直至反演出的纵横波走时无限地接近拾取的纵横波走时(即无限地接近观测的真实值),由此获得定位精度更高的微地震事件定位结果RT* event={L* event,Z* event}。
下面具体介绍上述各步骤的实施过程。
首先,对各向异性纵波走时定位方程OPJP、各向异性横波走时定位方程OPJS、各向异性纵横波走时定位方程OPJPS三者进行高次幂线性组合,建立本发明目标方程OPJ:
OPJ=OPJ4 PS+OPJ3 PS+OPJ3 P+OPJ2 P+OPJ2 S+OPJS (1)
OPJP=|TP-TPi| (2)
OPJS=|TS-TSi| (3)
OPJPS=|(TP-TS)-(TPi-TSi)| (4)
其中,TP、TS分别为拾取的射孔或微地震事件纵横波走时,TPi、TSi分别为反演出的射孔或微地震事件纵横波走时。在此,目标方程OPJ是与射孔或微地震事件空间位置(径向坐标L、深度坐标Z)、各向异性参数(ε、δ)、纵横波速度(VP0、VS0)相关的函数。换言之,本发明目标方程OPJ有6个自变量。
然后,输入射孔信息,利用目标方程OPJ(式(1)~(4)),开展射孔各向异性参数反演。具体地,将拾取的射孔纵横波走时Tshoot,P、Tshoot,S作为实际观测值,将已知的声波测井纵横速度和已知的射孔空间位置,代入上述目标方程OPJ(式(1)~(4))中:
OPJshoot=OPJ4 shoot,PS+OPJ3 shoot,PS+OPJ3 shoot,P+OPJ2 shoot,P+OPJ2 shoot,S+OPJshoot,S (5)
OPJshoot,P=|Tshoot,P-Tshoot,Pi| (6)
OPJshoot,S=|Tshoot,S-Tshoot,Si| (7)
OPJshoot,PS=|(Tshoot,P-Tshoot,S)-(Tshoot,Pi-Tshoot,Si)| (8)
此时式(5)仅仅是各向异性参数ε、δ的函数,通过以下偏导方程实现对射孔各向异性的反演:
在此,可以优选地利用最小二乘法或网格搜索法求解式(9),寻找一组各向异性参数Δ={ε,δ},使得反演出的射孔纵横波走时无限接近实际观测值,即无限接近拾取的射孔纵横波走时。
其次,再次利用目标方程OPJ(式(1)~(4)),输入射孔各向异性参数,开展微地震事件定位处理,以此获得微地震事件初始定位结果RTevent。具体地,将拾取的微地震事件纵横波走时Tevent,P、Tevent,S作为实际观测值,结合已知的声波测井纵横速度,输入上一步反演出的射孔各向异性参数Δ={ε,δ},代入上述目标方程OPJ(式(1)~(4)):
OPJeventt=OPJ4 event,PS+OPJ3 event,PS+OPJ3 event,P+OPJ2 event,P+OPJ2 event,S+OPJevent,S(10)
OPJevent,P=|Tevent,P-Tevent,Pi| (11)
OPJevent,S=|Tevent,S-Tevent,Si| (12)
OPJevent,PS=|(Tevent,P-Tevent,S)-(Tevent,Pi-Tevent,Si)| (13)
此时式(13)仅仅是微地震事件径向坐标L、深度坐标Z的函数,通过以下偏导方程实现对微地震事件空间位置的反演:
同样地,在此可以优选地利用最小二乘法或网格搜索法求解式(14),当反演出的微地震事件纵横波走时与观测值(即拾取的微地震事件纵横波走时)之间的误差最小时,输出对应的微地震事件空间位置RTevent={Levent,Zevent},即为微地震事件初始定位结果。
最后,再次利用本发明目标方程OPJ(式(1)~(4)),对微地震事件做进一步的层析定位分析,以此获得定位精度更高的微地震事件空间位置和相应的修正后的各向异性参数。
具体地,由于微地震事件所在位置各向异性参数与射孔所在位置存在细微差异,这种差异会导致利用射孔各向异性参数进行微地震事件定位的结果可能会存在较大误差。为了减小此误差,本发明在此提出利用层析思路,进一步进行精确定位。求解本发明微地震事件具体表达式OPJevent。
与上一步骤不同的是,在此除了反演微地震事件空间位置的径向坐标L、深度坐标Z之外,同时还反演微地震事件各向异性参数ε、δ。换言之,此时式(13)不仅是微地震事件径向坐标L、深度坐标Z的函数,而且还是微地震事件各向异性参数ε、δ的函数。因此,在此需要通过以下四个偏导方式实现层析反演:
同样地,可以优选地利用最小二乘法或网格搜索法来求解上述四个偏导方程。具体过程为:
以微地震事件初始定位位置RTevent={Levent,Zevent}为中心,根据精度要求,建立微地震事件空间网格[Levent-ΔL,Levent+ΔL]、[Zevent-ΔZ,Zevent+ΔZ];
同样地,以射孔各向异性参数Δ={ε,δ}为中心,建立微地震事件各向异性参数网格[ε-Δε,ε+Δε]、[δ-Δδ,δ+Δδ];
针对微地震事件,每个各向异性参数范围内所有可能值,反演对应的微地震事件空间位置的径向坐标、深度坐标,然而射线追踪正演出相应的纵横波初至走时,带入式(10),计算正演出的纵横波走时与拾取的纵横波走时之间的误差;寻找所有可能误差中的最小值,与所述最小值对应的微地震事件空间位置以及各向异性参数,即为最终的微地震事件层析定位结果RT* event={L* event,Z* event}以及相应的修正后的各向异性参数Δ*={ε*,δ*}。
实施例二
在本实施例中,利用现有的井中微地震模型数据来验证利用本发明的方法获得的微地震事件定位结果的准确性。
本次井中微地震观测例子为14级井下检波器监测微地震信号、1个已知射孔信号与21个微地震事件信号,其观测系统的径向、深度坐标几何关系如图2所示。图3所示的列表是井中微地震模型事件不同震源点位置的空间坐标以及与之对应的各向异性参数ε、δ。在实施本发明的定位方法之前,已经利用已有高精度射线追踪算法,根据图2观测方式与图3的各向异性参数表正演出每个震源点达到检波器的纵横波走时,作为已知观测值来开展反演输入走时真实值。需要说明的是,在此不考虑横向变化速度影响,即不考虑VTI介质,另外声波测井提供纵横波垂向速度作为已知输入。
首先,根据本发明式(1)~(4),开展基于已知射孔位置的纵横波高次幂线性组合定位分析(如式(5)~(8)),反演出射孔各向异性参数。在具体操作时,各向异性参数ε的扰动范围选取[0,0.7]且精度要求是0.003,各向异性参数δ的扰动范围选取[-0.1,0.1]且精度要求是0.001。利用式(9)进行求解。图4的列表示出了射孔各向异性参数ε、δ的真实值和反演值。由图4可以看出,反演值与真实值之间误差非常小,这对进行下一步的微地震事件定位处理十分有利。
然后,将上述射孔定位反演出的各向异性参数作为微地震事件初始各向异性参数输入,以微地震事件纵横波的真实走时为观测对象,直接对微地震事件进行纵横波走时高次幂线性组合定位处理(如式(10)~(13))。通过求解式(14),寻找微地震事件的空间位置,使得对应的纵横波走时最接近真实的观测值,将此时空间位置作为微地震事件初始定位结果。
从图5可以看出,基于射孔定位分析出的各向异性,利用本发明的目标方程获得的微地震事件定位结果,存在一定的误差(图6显示了反演值与真实值的对比,其中dL_deta为径向误差,dZ_deta为深度误差)。这说明,不同位置各向异性差异性会引起一定程度的定位误差,尤其是在径向上。
最后,为了进一步提高定位精度,在上一步获得的微地震事件初始定位结果的基础上,再次利用本发明的目标方程(10)~(13),进行进一步的层析定位,求解式(15)。即,以微地震事件初始定位结果为中心,给定径向+/-25米的扰动范围、深度+/-10米的扰动范围且精度为1米,给定各向异性参数ε+/-0.3的扰动范围、各向异性参数δ+/-0.05的扰动范围且精度为0.001,同时反演微地震事件空间位置与各向异性参数,直至纵横波走时误差进一步减小且无限接近观测真实值,从而实现微地震事件最终精确定位(如图7所示)。
图8为上述层析再定位结果的误差统计结果。从图8可以看出,径向、深度的定位误差显著地进一步减小。这充分验证了利用本发明的方法可以实现各向异性介质微地震定位处理,并且获得更高精度的微地震定位结果。
在此,应该理解的是,尽管上文已经描述了所公开系统和方法的示例性实施例,但是它们仅用于示例性的目而并非加以限制。本公开并非是详尽的,并且不限制所公开的精确形式。例如,本发明的高次幂线性组合目标方程OPJ也可以是以下带有系数的线性组合:
OPJ=αOPJ4 PS+βOPJ3 PS+χOPJ3 P+δOPJ2 P+εOPJ2 S+φOPJS。
其中,α,β,χ,δ,ε,φ为常数系数。
因此,在不脱离本公开广度或范围的情况下,可根据上述教导进行修改及变型,或者从本公开的实践中获知修改及变型。
Claims (9)
1.一种各向异性纵横波走时高次幂线性组合定位方法,包括以下步骤:
S1:将各向异性纵波走时定位方程OPJP、各向异性横波走时定位方程OPJS和各向异性纵横波走时时差定位方程OPJPS按照定位贡献的高低进行高次幂线性组合,建立目标方程OPJ;
S2:拾取射孔的纵横波走时Tshoot,P、Tshoot,S,结合已知的射孔空间位置和声波测井纵横波速度,利用所述目标方程OPJ对射孔开展定位分析,反演出射孔各向异性参数Δ={ε,δ},其中ε、δ为描述射孔各向异性的两个参数;
S3:拾取微地震事件纵横波走时Tevent,P、Tevent,S,将反演出的射孔各向异性参数Δ={ε,δ}作为微地震事件初始各向异性参数,结合已知的声波测井纵横波速度,利用所述目标方程OPJ对微地震事件开展定位分析,获得微地震事件初始定位结果RTevent={Levent,Zevent},其中Levent、Zevent为描述微地震事件空间位置的径向坐标、深度坐标;
S4:以微地震事件初始定位结果为空间约束,以微地震事件初始各向异性参数为参数约束,在扰动范围内,利用所述目标方程OPJ对微地震事件做进一步的层析定位处理,获得定位精度更高的微地震事件最终定位结果;
其中,所述目标方程OPJ为:
OPJ=OPJ4 PS+OPJ3 PS+OPJ3 P+OPJ2 P+OPJ2 S+OPJS
OPJP=|TP-TPi|
OPJS=|TS-TSi|
OPJPS=|(TP-TS)-(TPi-TSi)|
其中,TP、TS分别为拾取的射孔或微地震事件纵横波走时,TPi、TSi分别为反演出的射孔或微地震事件纵横波走时。
2.根据权利要求1所述的定位方法,其特征在于,所述步骤S4获得的微地震事件最终定位结果满足以下条件:根据步骤S4获得的微地震事件最终定位结果推演出的微地震事件纵横波走时与步骤S3拾取的微地震事件纵横波走时之间的误差达到最小。
3.根据权利要求1所述的定位方法,其特征在于,其特征在于,所述步骤S2具体包括以下步骤:
将拾取的射孔的纵横波走时Tshoot,P、Tshoot,S、已知的射孔空间位置和声波测井纵横波速度代入所述目标方程OPJ,获得射孔目标方程OPJshoot:
OPJshoot=OPJ4 shoot,PS+OPJ3 shoot,PS+OPJ3 shoot,P+OPJ2 shoot,P+OPJ2 shoot,S+OPJshoot,S
OPJshoot,P=|Tshoot,P-Tshoot,Pi|
OPJshoot,S=|Tshoot,S-Tshoot,Si|
OPJshoot,PS=|(Tshoot,P-Tshoot,S)-(Tshoot,Pi-Tshoot,Si)|
其中,OPJshoot,P为射孔各向异性纵波走时定位方程,OPJshoot,S为射孔各向异性横波走时定位方程,OPJshoot,PS为射孔各向异性纵横波时差定位方程,Tshoot,Pi、Tshoot,Si为反演的射孔纵横波走时;
反演射孔各向异性,利用以下偏导方程求解射孔各向异性参数Δ={ε,δ}:
4.根据权利要求3所述的定位方法,其特征在于,反演出的射孔各向异性参数Δ={ε,δ}满足以下条件:所述射孔各向异性参数Δ={ε,δ}能够使得所述射孔目标方程OPJshoot的值达到最小。
5.根据权利要求1所述的定位方法,其特征在于,所述步骤S3具体包括以下步骤:
将反演出的射孔各向异性参数作为微地震事件初始各向异性参数;
将拾取的微地震事件纵横波走时Tevent,P、Tevent,S,所述微地震事件初始各向异性参数和已知的声波测井纵横波速度代入所述目标方程OPJ,获得微地震事件目标方程OPJevent,
OPJevent=OPJ4 event,PS+OPJ3 event,PS+OPJ3 event,P+OPJ2 event,P+OPJ2 event,S+OPJevent,S
OPJevent,P=|Tevent,P-Tevent,Pi|
OPJevent,S=|Tevent,S-Tevent,Si|
OPJevent,PS=|(Tevent,P-Tevent,S)-(Tevent,Pi-Tevent,Si)|
其中,OPJevent,P为微地震各向异性纵波走时定位方程,OPJevent,S为微地震各向异性横波走时定位方程,OPJevent,PS为微地震各向异性纵横波时差定位方程,Tevent,Pi、Tevent,Si为反演出的微地震事件纵横波走时;
反演微地震事件空间位置,利用以下偏导方程求解微地震事件初始定位结果RTevent={Levent,Zevent}:
6.根据权利要求5所述的定位方法,其特征在于,获得的微地震事件初始定位结果RTevent={Levent,Zevent}满足以下条件:所述微地震事件初始定位结果RTevent={Levent,Zevent}能够使得所述微地震事件目标方程OPJevent的值达到最小。
7.根据权利要求3或5所述的定位方法,其特征在于,利用最小二乘法或网格搜索法求解各偏导方程。
9.根据权利要求8所述的定位方法,其特征在于,利用最小二乘法或网格搜索法求解步骤S4中的偏导方程,具体包括以下步骤:
在扰动范围内,针对每个微地震事件各向异性参数可能值,反演出与之对应的表征微地震事件空间位置的径向坐标、深度坐标,并据此推演出相应的微地震事件纵横波走时,计算推演出的微地震事件纵横波走时与拾取的微地震事件纵横波走时之间的误差;
寻找所有可能误差中的最小值,与所述最小值对应的径向坐标、深度坐标以及各向异性参数即为微地震事件最终定位结果以及相应的修正后的微地震事件各向异性参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810840579.0A CN110764139B (zh) | 2018-07-27 | 2018-07-27 | 一种各向异性纵横波走时高次幂线性组合定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810840579.0A CN110764139B (zh) | 2018-07-27 | 2018-07-27 | 一种各向异性纵横波走时高次幂线性组合定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110764139A CN110764139A (zh) | 2020-02-07 |
CN110764139B true CN110764139B (zh) | 2021-05-25 |
Family
ID=69327968
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810840579.0A Active CN110764139B (zh) | 2018-07-27 | 2018-07-27 | 一种各向异性纵横波走时高次幂线性组合定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110764139B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2014136780A (ru) * | 2014-09-11 | 2016-03-27 | Открытое акционерное общество "Акустический институт имени академика Н.Н. Андреева" | Способ сейсмического мониторинга процесса освоения месторождения углеводородов на акваториях |
CN105510880A (zh) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | 一种基于双差法的微地震震源定位方法 |
CN105549068A (zh) * | 2015-12-09 | 2016-05-04 | 中国科学院地质与地球物理研究所 | 一种三维各向异性微地震干涉逆时定位方法及系统 |
CN106353792A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种适用于水力压裂微震震源定位的方法 |
CN108254780A (zh) * | 2018-01-22 | 2018-07-06 | 河海大学 | 一种微地震定位及各向异性速度结构层析成像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8494777B2 (en) * | 2008-04-09 | 2013-07-23 | Schlumberger Technology Corporation | Continuous microseismic mapping for real-time 3D event detection and location |
-
2018
- 2018-07-27 CN CN201810840579.0A patent/CN110764139B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2014136780A (ru) * | 2014-09-11 | 2016-03-27 | Открытое акционерное общество "Акустический институт имени академика Н.Н. Андреева" | Способ сейсмического мониторинга процесса освоения месторождения углеводородов на акваториях |
CN105510880A (zh) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | 一种基于双差法的微地震震源定位方法 |
CN106353792A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种适用于水力压裂微震震源定位的方法 |
CN105549068A (zh) * | 2015-12-09 | 2016-05-04 | 中国科学院地质与地球物理研究所 | 一种三维各向异性微地震干涉逆时定位方法及系统 |
CN108254780A (zh) * | 2018-01-22 | 2018-07-06 | 河海大学 | 一种微地震定位及各向异性速度结构层析成像方法 |
Non-Patent Citations (7)
Title |
---|
Joint microseismic location and anisotropic tomography using differential arrival times and differential backazimuths;Junlun Li et al.;《Geophysical Journal International》;20131231;第1917-1931页 * |
Locating Bakken microseismic events with simultaneous anisotropic tomography and extended double-difference method;Junlun Li et al.;《SEG Houston 2013 Annual Meeting》;20131231;第2073-2078页 * |
Microseismic event location using an inverse method of joint P–S phase arrival difference and P-wave arrival difference in a borehole system;Wen Zhou et al.;《Journal of Geophysics and Engineering》;20150310;第220-226页 * |
Microseismic joint location and anisotropic velocity inversion for hydraulic fracturing in a tight Bakken reservoir;Junlun Li et al.;《GEOPHYSICS》;20141031;第79卷(第5期);第C111-C122页 * |
各向异性介质中的微震监测和声搜索定位方法;黄麟淇等;《东北大学学报(自然科学版)》;20150531;第36卷;第238-243页 * |
地层各向异性对微地震反演定位的影响;徐奔奔;《中国地球科学联合学术年会 2015》;20151231;第1528-1529页 * |
水平旋转与极性判断结合的井中微地震资料震相识别方法;余波等;《石油物探》;20131130;第52卷(第6期);第586-第595页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110764139A (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110133715B (zh) | 一种基于初至时差和波形叠加的微地震震源定位方法 | |
CA2542577C (en) | Radial profiling of slowness: methods and apparatus for near-wellbore alteration estimation | |
US10768323B2 (en) | Methods and systems for seismic data analysis using a tilted transversely isotropic (TTI) model | |
AU2014407527B2 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
CN102943669B (zh) | 适用于水平井轨迹设计的高精度深度域随机模拟反演方法 | |
CN110531410A (zh) | 一种基于直达波场的最小二乘逆时偏移梯度预条件方法 | |
CN103852785A (zh) | 地层各向异性的评价方法 | |
CN111665556B (zh) | 地层声波传播速度模型构建方法 | |
CN110515122B (zh) | 正演网格搜索定位及微地震信号识别方法及装置 | |
CN111208564A (zh) | 一种深度域层位标定方法及装置 | |
CN109085642B (zh) | 一种各向异性介质微地震事件定位方法 | |
CN105607119A (zh) | 近地表模型构建方法与静校正量求取方法 | |
CN102877828A (zh) | 一种三维多井联合井地ct成像方法 | |
CN110764136B (zh) | 各向异性纵横波走时线性组合与非线性组合联合定位方法 | |
CN110764140B (zh) | 基于射孔双差各向异性纵横波非线性联合定位方法 | |
CN110764139B (zh) | 一种各向异性纵横波走时高次幂线性组合定位方法 | |
CN110764138B (zh) | 一种各向异性纵横波走时非线性组合定位方法 | |
CN109212594B (zh) | 一种各向异性介质纵横波联合定位方法 | |
CN110764137B (zh) | 基于射孔混合时差各向异性纵横波非线性联合定位方法 | |
CN110764148B (zh) | 一种各向异性矢量波场井地联合定位方法 | |
CN109655892B (zh) | 各向异性纵横波走时三元多项式组合定位方法及系统 | |
CN113267830A (zh) | 基于非结构网格下二维重力梯度与地震数据联合反演方法 | |
CN110888158B (zh) | 一种基于rtm约束的全波形反演方法 | |
CN109655919B (zh) | 各向异性纵横波走时点乘定位方法及系统 | |
CN109212593B (zh) | 一种基于多射孔双差的纵横波联合定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |