CN107436449A - 一种提取水平层状介质各向异性参数的方法 - Google Patents
一种提取水平层状介质各向异性参数的方法 Download PDFInfo
- Publication number
- CN107436449A CN107436449A CN201610362068.3A CN201610362068A CN107436449A CN 107436449 A CN107436449 A CN 107436449A CN 201610362068 A CN201610362068 A CN 201610362068A CN 107436449 A CN107436449 A CN 107436449A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mfrac
- msup
- mover
- msub
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000013598 vector Substances 0.000 claims abstract description 23
- 238000010276 construction Methods 0.000 claims abstract description 3
- 238000012986 modification Methods 0.000 claims description 10
- 230000004048 modification Effects 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 7
- 239000002245 particle Substances 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 15
- 238000000605 extraction Methods 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- 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. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于微地震监测数据处理技术,公开了一种提取水平层状介质各向异性参数的方法,该方法包括:利用旅行时残差、射线路径长度向量、群慢度与各层中各向异性参数构造超定方程;求解超定方程获得各向异性参数。本发明只利用井中微地震自身的信息提取水平层状介质各向异性参数,特别适合没有三维地震和正交偶极子声波测井的地区。
Description
技术领域
本发明涉及微地震监测数据处理技术,是一种用井中微地震自身的信息提取水平层状介质各向异性参数的计算方法。
背景技术
微地震压裂监测技术现在已成为致密储层油气田开发中一项常用的监测技术,其最主要的应用是水力压裂裂缝分布的监测,在优化压裂方案和开发井网部署方面起到重要作用。微地震监测有地面监测和井中监测两种方式。地面监测就是在监测目标(如压裂井)周围地面布置若干接收点或测线进行微地震监测。井中监测是在监测目标区域周围临近的井中布置接收检波器排列进行微地震监测。
微地震监测的首要目标是对微地震事件进行定位,而定位精度取决于模型的准确。井中微地震一般采用水平层状介质模型。由于观测井离压裂井较近(水平距离约为200米~600米),并且受压裂对地层破坏的影响,井中微地震更易受各向异性因素影响。有些地区受条件限制(无地面三维地震资料、无正交偶极子声波测井)不能提供各向异性参数。
因此,本领域亟需一种提取水平层状介质各向异性参数的方法。
发明内容
为了解决现有技术中存在的问题,本发明提出了一种利用井中微地震信息提取水平层状介质各向异性参数的方法。该方法只利用井中微地震自身的信息提取水平层状介质各向异性参数,特别适合没有三维地震和正交偶极子声波测井的地区。
根据本发明的一个方面,提供一种提取水平层状介质各向异性参数的方法,该方法包括:利用旅行时残差、射线路径长度向量、群慢度与各层中各向异性参数构造超定方程;求解超定方程获得各向异性参数。
进一步地,该方法在求解超定方程获得各向异性参数后获得修正量;判断是否满足模型的修正量充分小,不满足条件则更新模型,在初始模型各个参数加上所获得的修正量,然后再重新构造方程和求解;满足条件则输出模型。
进一步地,所述超定方程为:
式中为旅行时残差,Δτi为发震时刻误差;为射线路径长度向量;为P波或SV波的群慢度,即 为或矢量;分别为垂向速度的修正向量、ε的修正向量和δ的修正向量。
进一步地,在水平层状介质中P波和SV波的精确相速度公式为:
式中P表示纵波,SV表示质点振动的垂直平面内的横波分量,α0为纵波垂向速度,β0为横波垂向速度。
进一步地,在水平层状介质中P波和SV波的群速度计算公式为:
θg=θ+Δθ (5)
式中,θg为群角。V为P波和SV波的相速度。
进一步地,通过给定井中检波器获得到达时利用震源位置、发震时刻和初始模型来计算理论到达时两者旅行时残差为:
进一步地,通过LU分解、SVD分解、共轭梯度算法或LSQR方法求解超定方程。
本发明方法只利用井中微地震自身的信息提取水平层状介质各向异性参数,特别适合没有三维地震和正交偶极子声波测井的地区。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式中,相同的参考标号通常代表相同部件。
图1显示了相速度和群速度之间的关系。
图2显示了井中微地震观测的示例图。
图3显示了根据本发明实施例的利用井中微地震信息提取水平层状介质各向异性参数的方法流程图。
附图仅用于示例性说明,不能理解为对本专利的限制;为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的;相同或相似的标号对应相同或相似的部件;附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
下面参照图1和图2描述本发明方法的原理。首先,利用Thomsen参数表示相速度和群速度。水平层状介质各向异性参数一般使用传统的Thomsen参数化方式。P波和SV波的精确相速度公式为:
式中P表示纵波(正号),SV表示质点振动的垂直平面内的横波分量(负号)。α0为纵波垂向速度,β0为横波垂向速度。θ为入射角,P波和SV波都用一个符号表示,下标表明了何种波。
在水平层状介质中P波和SV波的群速度计算公式为:
θg=θ+Δθ (5)
式中,θg为群角。V为P波和SV波的相速度。
接下来,求取各向异性参数。
在井中观测微地震时,已知的信息包括微地震事件的初至信息、射孔或导爆索的具体位置和记录到的初至信息。利用这些信息可以计算出Thomsen参数 (是一种各向异性参数表示法,指利用ε,δ参数描述各项异性)。
微地震事件i(i=1,…,n)到接收点k(k=1,…,m)的旅行时可以用射线理论表示为:
其中,u为群慢度场。对于微地震监测来讲,实际观测的旅行时为其中τi为发震时刻。
给定井中检波器得到的到达时(P波、SV波),利用震源(射孔、导爆索、微地震事件)位置,发震时刻和初始模型来计算理论到达时初始模型中的信息有P波和SV波的垂向速度(一般可利用声波测井建立)和ε,δ,如果没有地面地震信息或VSP信息,初始模型信息一般为零)。
两者旅行时的残差为
方程(8)为非线性方程,利用泰勒展开来线性化该方程,则残差可以用Thomsen参数的扰动表示为:
式中Δτi为发震时刻误差;为射线路径长度向量,可根据射线追踪获得; 为P波或SV波的群慢度(即);为或矢量;由方程(3)求出;分别为垂向速度的修正向量、ε的修正向量和δ的修正向量。向量元素的个数等于层的个数。
求解方程(9)组成的方程即可获得所需的各向异性参数(实际就是要求解修正向量),计算方法可选LU分解、SVD分解、共轭梯度算法、 LSQR方法。
下面以LSQR方法为例,描述求解过程。
利用方程(9)可构建稀疏线性方程组Ax=b。
LSQR方法具有计算量小的优点,并且能很容易地利用矩阵的稀疏性简化计算,因此适合求解大型稀疏矩阵。
LSQR原理如下:
(1)初始化
(2)双对角化:
βi+1ui+1=Avi-αiui
αi+1vi+1=ATui+1-βi+1vi
(3)正交变换:
(4)解矢量:
wi+1=vi+1-wiθi+1/ρi
上述迭代公式中,αi,βi为对角化矩阵A的标量序列,ui,vi为矢量序列,λ为阻尼因子,其余为中间变量。利用LSQR方法求解方程(9)是本领域技术人员容易实现的,再次不再赘述。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
可以利用声波测井建立水平层状模型,模型信息可包括P波和SV波的垂向速度,各层ε,δ初始化为零;接下来利用观测井三分量检波器记录的数据,拾取微地震事件初至,例如纵波初至和横波初至;然后构建超定方程,利用数值计 算方法求解即可得到各向异性参数。
参照图3,描述根据本发明一个实施例的方法流程。首先,输入数据,所述数据可来自观测井三分量检波器记录的数据,例如三分量数据的旅行时;接下来通过射线追踪获得射线路径长度向量;构造超定方程(9),例如方程(9)里面包括了射线路径的长度,以及各层中各向异性参数与群慢度的关系,换言之,方程(9)就是旅行时残差与各层参数之间的关系;通过LSQR方法求解方程(9)获得各向异性参数,即求解修正向量等;接下来判断是否满足模型的修正量充分小,该过程可根据本领域技术人员的经验判断;不满足条件则更新模型,例如在初始模型各个参数加上这些修正量(加法运算),然后再重新构造方程和求解;满足条件则输出模型,即输出每层中P波和SV波的垂向速度和ε,δ。
下面结合具体实施方式对本发明作进一步的说明。其中,附图仅用于示例性说明,表示的仅是示意图,而非实物图,不能理解为对本专利的限制;为了更好地说明本发明的实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (7)
1.一种提取水平层状介质各向异性参数的方法,其特征在于,该方法包括:
利用旅行时残差、射线路径长度向量、群慢度与各层中各向异性参数构造超定方程;
求解超定方程获得各向异性参数。
2.根据权利要求1所述的方法,其特征在于,该方法进一步包括:
求解超定方程获得各向异性参数后获得修正量;
判断是否满足模型的修正量充分小,不满足条件则更新模型,在初始模型各个参数加上所获得的修正量,然后再重新构造方程和求解;
满足条件则输出模型。
3.根据权利要求1所述的方法,其特征在于,所述超定方程为:
<mrow>
<msubsup>
<mi>r</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msup>
<mi>&Delta;&tau;</mi>
<mi>i</mi>
</msup>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>&epsiv;</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mover>
<mi>&epsiv;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mover>
<mi>&delta;</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mover>
<mi>&delta;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中为旅行时残差,Δτi为发震时刻误差;为射线路径长度向量;为P波或SV波的群慢度,即 为或矢量;分别为垂向速度的修正向量、ε的修正向量和δ的修正向量。
4.根据权利要求3所述的方法,其特征在于,在水平层状介质中P波和SV波的精确相速度公式为:
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>P</mi>
<mo>,</mo>
<mi>S</mi>
<mi>V</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>&epsiv;sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>-</mo>
<mfrac>
<mi>f</mi>
<mn>2</mn>
</mfrac>
<mo>&PlusMinus;</mo>
<mfrac>
<mi>f</mi>
<mn>2</mn>
</mfrac>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>f</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>&alpha;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&beta;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
<mi>D</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msup>
<mi>&epsiv;sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
<mi>f</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>-</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
<mi>f</mi>
</mfrac>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中P表示纵波,SV表示质点振动的垂直平面内的横波分量,α0为纵波垂向速度,β0为横波垂向速度。
5.根据权利要求4所述的方法,其特征在于,在水平层状介质中P波和SV波的群速度计算公式为:
<mrow>
<msubsup>
<mi>V</mi>
<mi>g</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mi>dv</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>d</mi>
<mi>&theta;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mi>&Delta;</mi>
<mi>&theta;</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mfrac>
<mrow>
<msup>
<mi>dV</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>d</mi>
<mi>&theta;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
θg=θ+Δθ (5)
<mrow>
<mfrac>
<mrow>
<msup>
<mi>dV</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>d</mi>
<mi>&theta;</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msubsup>
<mi>&alpha;</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>&lsqb;</mo>
<mi>&epsiv;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>&PlusMinus;</mo>
<mfrac>
<mrow>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>+</mo>
<mn>2</mn>
<msup>
<mi>&epsiv;sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>-</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>-</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>4</mn>
<mi>&theta;</mi>
</mrow>
<mrow>
<mi>f</mi>
<mi>D</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,θg为群角。V为P波和SV波的相速度。
6.根据权利要求3所述的方法,其特征在于,通过给定井中检波器获得到达时利用震源位置、发震时刻和初始模型来计算理论到达时两者旅行时残差为:
<mrow>
<msubsup>
<mi>r</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>t</mi>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>t</mi>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
7.根据权利要求1所述的方法,其特征在于,通过LU分解、SVD分解、共轭梯度算法或LSQR方法求解超定方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362068.3A CN107436449A (zh) | 2016-05-26 | 2016-05-26 | 一种提取水平层状介质各向异性参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362068.3A CN107436449A (zh) | 2016-05-26 | 2016-05-26 | 一种提取水平层状介质各向异性参数的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107436449A true CN107436449A (zh) | 2017-12-05 |
Family
ID=60454464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610362068.3A Pending CN107436449A (zh) | 2016-05-26 | 2016-05-26 | 一种提取水平层状介质各向异性参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107436449A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108387937A (zh) * | 2018-01-26 | 2018-08-10 | 中国科学院声学研究所 | 一种交叉偶极子各向异性快速反演方法及系统 |
CN109025979A (zh) * | 2018-04-24 | 2018-12-18 | 长江大学 | 一种正演模拟蚀变地层随钻多极子声波频散曲线的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014164970A2 (en) * | 2013-03-12 | 2014-10-09 | Ion Geophysical Corporation | Methods and systems for locating seismic events |
CN104407378A (zh) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种各向异性参数反演方法及装置 |
CN104422955A (zh) * | 2013-08-22 | 2015-03-18 | 中国石油化工股份有限公司 | 一种利用旅行时变化量进行各向异性参数提取的方法 |
CN105319590A (zh) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种基于hti介质的各向异性单参数反演方法 |
CN105510880A (zh) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | 一种基于双差法的微地震震源定位方法 |
-
2016
- 2016-05-26 CN CN201610362068.3A patent/CN107436449A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014164970A2 (en) * | 2013-03-12 | 2014-10-09 | Ion Geophysical Corporation | Methods and systems for locating seismic events |
CN104422955A (zh) * | 2013-08-22 | 2015-03-18 | 中国石油化工股份有限公司 | 一种利用旅行时变化量进行各向异性参数提取的方法 |
CN105319590A (zh) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种基于hti介质的各向异性单参数反演方法 |
CN105510880A (zh) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | 一种基于双差法的微地震震源定位方法 |
CN104407378A (zh) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种各向异性参数反演方法及装置 |
Non-Patent Citations (3)
Title |
---|
吴国忱: "《各向异性介质地震波传播与成像》", 31 July 2006, 中国石油大学出版社 * |
李河昭 等: "起伏地表下VTI介质射线追踪与多参数联合反演", 《SPG/SEG北京2016国际地球物理会议电子文集》 * |
郭立鹏 等: "复杂各向异性介质初至波射线追踪", 《石油物探》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108387937A (zh) * | 2018-01-26 | 2018-08-10 | 中国科学院声学研究所 | 一种交叉偶极子各向异性快速反演方法及系统 |
CN108387937B (zh) * | 2018-01-26 | 2019-09-03 | 中国科学院声学研究所 | 一种交叉偶极子各向异性快速反演方法及系统 |
CN109025979A (zh) * | 2018-04-24 | 2018-12-18 | 长江大学 | 一种正演模拟蚀变地层随钻多极子声波频散曲线的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101893720B (zh) | 一种地震波的矢量波场分离与合成的方法和系统 | |
CN105093274B (zh) | 一种水力压裂裂缝震源机制的反演方法及系统 | |
CA2935604C (en) | Methods and systems of determining a fault plane of a microseismic event | |
CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
CN105842735A (zh) | 具有复杂速度分布的区域岩体微震震源定位方法 | |
CN108181653B (zh) | 针对vti介质逆时偏移方法、设备及介质 | |
CN104101902B (zh) | 地震属性聚类方法及装置 | |
CN103675916B (zh) | 一种三分量检波器埋置方向高精度校正的方法 | |
CN110261902B (zh) | 一种基于多谱能量合成的地下浅层震源定位方法 | |
CN104237937B (zh) | 叠前地震反演方法及其系统 | |
CN102053261A (zh) | 一种地震数据处理方法 | |
CN104730574A (zh) | 构建近地表结构模型的方法 | |
CN104199090A (zh) | 一种地面监测微地震定位的速度模型构建与求解方法 | |
CN108873063A (zh) | 一种微地震矩张量反演的方法及装置 | |
CN106199697A (zh) | 模拟微地震的弹性波正演方法 | |
CN106249295A (zh) | 一种井中微地震p、s波联合快速定位方法及系统 | |
CN108717201A (zh) | 一种隧道围岩微震源定位方法 | |
CN106249297A (zh) | 基于信号估计的水力压裂微地震震源定位方法及系统 | |
CN105242313B (zh) | 一种弹性波逆时偏移极性反转的校正方法及系统 | |
CA2935465A1 (en) | Methods and systems of generating a velocity model | |
CN107436449A (zh) | 一种提取水平层状介质各向异性参数的方法 | |
CN106289371A (zh) | 水下爆破安全测定装置与方法 | |
CN104166159A (zh) | 四维微地震监测的裂缝形态处理方法和系统 | |
CN103412325A (zh) | 一种地震纵、横波波场分离与去噪的方法 | |
CN104597489A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171205 |