CN108732625A - 一种井周岩石非均匀各向异性的识别方法及系统 - Google Patents

一种井周岩石非均匀各向异性的识别方法及系统 Download PDF

Info

Publication number
CN108732625A
CN108732625A CN201810488255.5A CN201810488255A CN108732625A CN 108732625 A CN108732625 A CN 108732625A CN 201810488255 A CN201810488255 A CN 201810488255A CN 108732625 A CN108732625 A CN 108732625A
Authority
CN
China
Prior art keywords
inhomogeneous
module
dipole
yanshi
anisotropics
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.)
Granted
Application number
CN201810488255.5A
Other languages
English (en)
Other versions
CN108732625B (zh
Inventor
黑创
唐晓明
苏远大
庄春喜
罗明璋
李修权
张磊
陈伟国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangtze University
Original Assignee
Yangtze University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yangtze University filed Critical Yangtze University
Priority to CN201810488255.5A priority Critical patent/CN108732625B/zh
Publication of CN108732625A publication Critical patent/CN108732625A/zh
Application granted granted Critical
Publication of CN108732625B publication Critical patent/CN108732625B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Acoustics & Sound (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种井周岩石非均匀各向异性的识别方法及系统，所述方法通过在深度区间内进行陈列声波测井，获取深度区间内正交偶极四分量陈列波形数据，运用声波信号的散射衰减品质因子识别岩石均匀各向异性，岩石非均匀各向异性的结果可用于指导油气开发，为油气的勘探、开采和生产提供重要的信息。

Description

一种井周岩石非均匀各向异性的识别方法及系统
技术领域
本发明涉及油气和矿产资源勘探技术领域，特别涉及一种井周岩石非均匀各向异性的识别方法及系统。
背景技术
地下岩石普遍存在非均匀性。由于岩石圈中岩石的成分、密度、相态和缝洞中流体的组成、温度和压力等不均匀性，实际地层介质弹性参数，如弹性模量、速度和泊松比等，在空间分布上往往是非均匀的。岩石非均匀的各向异性识别对于油气和矿产资源的精细勘探至关重要。
测井尺度的岩石非均匀各向异性可以通过不同方位的时差差异、层析成像等方法进行识别，这些方法利用的主要是声波测井数据中的直达波信息，只能评价井壁附近的岩石非均匀各向异性。目前对于远离井壁的岩石非均匀各向异性的识别方法研究较少。事实上，当远离井壁的岩石存在非均匀各向异性时，不同方位会产生不同强度的散射波，可以利用不同方向声波测井数据中散射波的强弱识别岩石非均匀各向异性。而在声波测井常规数据处理中，散射波和尾波往往被视为噪声而加以压制。现代声波测井仪器已广泛应用于油气勘探，采集了大量正交偶极的阵列声波测井数据。如果能从声波测井数据正确分析和提取岩石的非均匀各向异性信息，将会提供很有价值的地层信息，并为声波测井的处理和解释开辟新的应用前景。
发明内容
针对现有技术的不足，本发明提供一种井周岩石非均匀各向异性的识别方法及系统，计算不同方位偶极声波散射衰减品质因子识别岩石非均匀各向异性，为计算声波测井频率下井周岩石非均匀各向异性识别提供一种实用可行的方法和技术。
为实现上述目的，本发明提供一种井周岩石非均匀各向异性的识别方法，所述方法包括以下步骤：
S1、在深度区间内进行陈列声波测井，获取深度区间内正交偶极四分量陈列波形数据；
S2、根据步骤S1获取深度区间多个深度位置的正交偶极四分量声波测井数据，并记录每一个深度位置处的仪器方位；
S3、将步骤S2中获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
S4、将方位为的偶极声波测井数据进行带通滤波，消除测井随机噪声，得到滤波后的偶极声波测井数据v(t)；
S5、利用希尔伯特变换计算v(t)的能量包络A(t)；
S6、创建目标函数ΔE，所述目标函数ΔE定义为理论的能量包络AC(t)和实测的能量包络A(t)的函数：
ΔE(t,QI,QS)＝|AC(t)-A(t)|2
其中QS为散射衰减品质因子，QI为固有衰减品质因子；
S7、根据步骤S6求解目标函数，得到方位为的全局极小值和进而得到不同方位散射衰减品质因子的集合根据集合中元素的差异判断岩石是否具有非均匀各向异性；
S8、重复步骤S3～步骤S7，获得深度区间内岩石非均匀各向异性结果。
在上述技术方案中，步骤S2中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
在上述技术方案中，步骤S4中所述方位取值区间为0度～360度。
在上述技术方案中，步骤S6中所述理论的能量包络AC(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，ID是直达波振幅AD在时间窗(t1,t2)的积分，表示为
在上述技术方案中，步骤S7所述集合中元素的差异判断分为以下步骤：
S71、若集合中任意两个元素数值相等，则岩石为非均匀各向同性；
S72、若集合中存在两个数值不相等的元素，则岩石为非均匀各向异性；
本发明还提供一种井周岩石非均匀各向异性的识别系统，其特征在于，所述系统包括采集模块、记录模块、转换模块、滤波模块、计算模块、创建模块、判断模块和循环模块；
采集模块，用于在深度区间内进行陈列声波测井，采集深度区间内正交偶极四分量陈列波形数据；
记录模块，用于在深度区间多个深度位置的正交偶极四分量声波测井数据的获取点，记录深度位置处的仪器方位；
转换模块，用于将获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
滤波模块，用于将方位为的偶极声波测井数据进行带通滤波，消除测井随机噪声，得到滤波后的偶极声波测井数据v(t)；
计算模块，用于运用希尔伯特变换计算v(t)的能量包络A(t)；
创建模块，创建目标函数ΔE，所述目标函数ΔE定义为理论的能量包络AC(t)和实测的能量包络A(t)的函数：
ΔE(t,QI,QS)＝|AC(t)-A(t)|2
其中QS为散射衰减品质因子，QI为固有衰减品质因子；
判断模块，用于求解目标函数，得到方位为的全局极小值和进而得到不同方位散射衰减品质因子的集合根据集合中元素的差异判断岩石是否具有非均匀各向异性；
循环模块，用于循环重复转换模块～判断模块，获得深度区间内岩石非均匀各向异性结果。
在上述技术方案中，记录模块中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
在上述技术方案中，滤波模块中所述方位取值区间为0度～360度。
在上述技术方案中，创建模块中所述理论的能量包络AC(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，ID是直达波振幅AD在时间窗(t1,t2)的积分，表示为
在上述技术方案中，判断模块中所述集合中元素的差异判断分为以下模块：
相同模块，若集合中任意两个元素数值相等，则岩石为非均匀各向同性；
相异模块，若集合中存在两个数值不相等的元素，则岩石为非均匀各向异性。
本发明一种井周岩石非均匀各向异性的识别方法及系统，具有以下有益效果：本发明采用不同方位散射衰减品质因子的差异来识别岩石非均匀各向异性，与现有技术方案只能判断井壁附近岩石非均匀各向异性相比，本发明可判断远离井壁岩石的非均匀各向异性。岩石非均匀各向异性的结果可用于指导油气开发如井位优化、酸化压裂等，为油气的勘探、开采和生产提供重要的信息。
附图说明
图1为本发明一种井周岩石非均匀各向异性的识别方法流程图；
图2为本发明实施例计算得到的某一深度位置不同方位散射衰减品质因子示意图；
图3为本发明一种井周岩石非均匀各向异性的识别系统模块图；
具体实施方式
下面结合附图对本发明作进一步详细描述，本发明提供一种井周岩石非均匀各向异性的识别方法，所述方法包括以下步骤，如图1所示：
S1、在深度区间内进行陈列声波测井，获取深度区间内正交偶极四分量陈列波形数据；
所述正交偶极四分量阵列数据是用多极子测井仪器的偶极模式采集的数据，正交偶极是相互垂直的两组换能器X和Y，与X1和Y1发射对应的接收X2和Y2，采集的数据为四组数据分别为：X1发X2收、X1发Y2收、Y1发X2收、Y1发Y2收。
S2、根据步骤S1获取深度区间多个深度位置的正交偶极四分量声波测井数据，并记录每一个深度位置处的仪器方位；
其中，仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
S3、将步骤S2中获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
其中，通过四组数据组合计算得到不同方位的偶极声波测井数据即为横波数据。
S4、将方位为的偶极声波测井数据进行带通滤波，消除测井随机噪声，得到滤波后的偶极声波测井数据v(t)；
其中，所述方位取值区间为0度～360度，为方位的横波数据。
S5、利用希尔伯特变换计算v(t)的能量包络A(t)；
S6、创建目标函数ΔE，所述目标函数ΔE定义为理论的能量包络AC(t)和实测的能量包络A(t)的函数：
ΔE(t,QI,QS)＝|AC(t)-A(t)|2
其中QS为散射衰减品质因子，QI为固有衰减品质因子；
其中，声波信号的散射衰减品质因子是描述岩石非均匀性的重要参数，不同方位声波测井信号的散射衰减品质因子的差异预示着岩石的非均匀各向异性。
具体的，所述理论的能量包络AC(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，ID是直达波振幅AD在时间窗(t1,t2)的积分，表示为
S7、根据步骤S6求解目标函数，得到方位为的全局极小值和进而得到不同方位散射衰减品质因子的集合根据集合中元素的差异判断岩石是否具有非均匀各向异性；
具体的，所述集合中元素的差异判断分为以下步骤：
S71、若集合中任意两个元素数值相等，则岩石为非均匀各向同性；
S72、若集合中存在两个数值不相等的元素，则岩石为非均匀各向异性；
S8、重复步骤S3～步骤S7，获得深度区间内岩石非均匀各向异性结果。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明；
如图2所示利用本发明的处理步骤计算得到的某一深度位置不同方位(0度～360度)的散射衰减品质因子，从图中可看出在45度方位的散射衰减品质因子最大(约为435)，150度方位对应的散射衰减品质因子最小(约为145)，利用步骤S7的判断条件，该深度点岩石具有非均匀各向异性。
如图3所示本发明还提供一种井周岩石非均匀各向异性的识别系统，所述系统包括采集模块、记录模块、转换模块、滤波模块、计算模块、创建模块、判断模块和循环模块；
采集模块，用于在深度区间内进行陈列声波测井，采集深度区间内正交偶极四分量陈列波形数据；
记录模块，用于在深度区间多个深度位置的正交偶极四分量声波测井数据的获取点，记录深度位置处的仪器方位；
其中，记录模块中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
转换模块，用于将获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
滤波模块，用于将方位为的偶极声波测井数据进行带通滤波，消除测井随机噪声，得到滤波后的偶极声波测井数据v(t)；
其中，所述方位取值区间为0度～360度。
计算模块，用于运用希尔伯特变换计算v(t)的能量包络A(t)；
创建模块，创建目标函数ΔE，所述目标函数ΔE定义为理论的能量包络AC(t)和实测的能量包络A(t)的函数：
ΔE(t,QI,QS)＝|AC(t)-A(t)|2
其中QS为散射衰减品质因子，QI为固有衰减品质因子；
具体的，所述理论的能量包络AC(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，ID是直达波振幅AD在时间窗(t1,t2)的积分，表示为
判断模块，用于求解目标函数，得到方位为的全局极小值和进而得到不同方位散射衰减品质因子的集合根据集合中元素的差异判断岩石是否具有非均匀各向异性；
具体的，所述集合中元素的差异判断分为以下模块：
相同模块，若集合中任意两个元素数值相等，则岩石为非均匀各向同性；
相异模块，若集合中存在两个数值不相等的元素，则岩石为非均匀各向异性。
循环模块，用于循环重复转换模块～判断模块，获得深度区间内岩石非均匀各向异性结果。
说明书中未阐述的部分均为现有技术或公知常识。本实施方式仅用于说明该发明，而不用于限制本发明的范围，本领域技术人员对于本发明所做的等价置换等修改均认为是落入该发明权利要求书所保护范围内。

Claims (10)

1.一种井周岩石非均匀各向异性的识别方法，其特征在于，所述方法包括以下步骤:
S1、在深度区间内进行陈列声波测井，获取深度区间内正交偶极四分量陈列波形数据；
S2、根据步骤S1获取深度区间多个深度位置的正交偶极四分量声波测井数据，并记录每一个深度位置处的仪器方位；
S3、将步骤S2中获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
S4、将方位为的偶极声波测井数据进行带通滤波，消除测井随机噪声，得到滤波后的偶极声波测井数据v(t)；
S5、利用希尔伯特变换计算v(t)的能量包络A(t)；
S6、创建目标函数ΔE，所述目标函数ΔE定义为理论的能量包络AC(t)和实测的能量包络A(t)的函数：
ΔE(t,QI,QS)＝|AC(t)-A(t)|2
其中QS为散射衰减品质因子，QI为固有衰减品质因子；
S7、根据步骤S6求解目标函数，得到方位为的全局极小值和进而得到不同方位散射衰减品质因子的集合根据集合中元素的差异判断岩石是否具有非均匀各向异性；
S8、重复步骤S3～步骤S7，获得深度区间内岩石非均匀各向异性结果。
2.根据权利要求1所述一种井周岩石非均匀各向异性的识别方法，其特征在于，步骤S2中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
3.根据权利要求1所述一种井周岩石非均匀各向异性的识别方法，其特征在于，步骤S4中所述方位取值区间为0度～360度。
4.根据权利要求1所述一种井周岩石非均匀各向异性的识别方法，其特征在于，步骤S6中所述理论的能量包络AC(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，ID是直达波振幅AD在时间窗(t1,t2)的积分，表示为
5.根据权利要求1所述一种井周岩石非均匀各向异性的识别方法，其特征在于，步骤S7所述集合中元素的差异判断分为以下步骤：
S71、若集合中任意两个元素数值相等，则岩石为非均匀各向同性；
S72、若集合中存在两个数值不相等的元素，则岩石为非均匀各向异性。
6.一种井周岩石非均匀各向异性的识别系统，其特征在于，所述系统包括采集模块、记录模块、转换模块、滤波模块、计算模块、创建模块、判断模块和循环模块；
采集模块，用于在深度区间内进行陈列声波测井，采集深度区间内正交偶极四分量陈列波形数据；
记录模块，用于在深度区间多个深度位置的正交偶极四分量声波测井数据的获取点，记录深度位置处的仪器方位；
转换模块，用于将获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
滤波模块，用于将方位为的偶极声波测井数据进行带通滤波，消除测井随机噪声，得到滤波后的偶极声波测井数据v(t)；
计算模块，用于运用希尔伯特变换计算v(t)的能量包络A(t)；
创建模块，创建目标函数ΔE，所述目标函数ΔE定义为理论的能量包络AC(t)和实测的能量包络A(t)的函数：
ΔE(t,QI,QS)＝|AC(t)-A(t)|2
其中QS为散射衰减品质因子，QI为固有衰减品质因子；
判断模块，用于求解目标函数，得到方位为的全局极小值和进而得到不同方位散射衰减品质因子的集合根据集合中元素的差异判断岩石是否具有非均匀各向异性；
循环模块，用于循环重复转换模块～判断模块，获得深度区间内岩石非均匀各向异性结果。
7.根据权利要求6所述一种井周岩石非均匀各向异性的识别系统，其特征在于，记录模块中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
8.根据权利要求6所述一种井周岩石非均匀各向异性的识别系统，其特征在于，滤波模块中所述方位取值区间为0度～360度。
9.根据权利要求6所述一种井周岩石非均匀各向异性的识别系统，其特征在于，创建模块中所述理论的能量包络AC(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，ID是直达波振幅AD在时间窗(t1,t2)的积分，表示为
10.根据权利要求6所述一种井周岩石非均匀各向异性的识别系统，其特征在于，判断模块中所述集合中元素的差异判断分为以下模块：
相同模块，若集合中任意两个元素数值相等，则岩石为非均匀各向同性；
相异模块，若集合中存在两个数值不相等的元素，则岩石为非均匀各向异性。
CN201810488255.5A 2018-05-21 2018-05-21 一种井周岩石非均匀各向异性的识别方法及系统 Active CN108732625B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810488255.5A CN108732625B (zh) 2018-05-21 2018-05-21 一种井周岩石非均匀各向异性的识别方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810488255.5A CN108732625B (zh) 2018-05-21 2018-05-21 一种井周岩石非均匀各向异性的识别方法及系统

Publications (2)

Publication Number Publication Date
CN108732625A true CN108732625A (zh) 2018-11-02
CN108732625B CN108732625B (zh) 2020-02-07

Family

ID=63938777

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810488255.5A Active CN108732625B (zh) 2018-05-21 2018-05-21 一种井周岩石非均匀各向异性的识别方法及系统

Country Status (1)

Country Link
CN (1) CN108732625B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112859163A (zh) * 2021-01-07 2021-05-28 中国石油大学(北京) 利用散射波确定压裂区地层品质因子变化的方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104007462A (zh) * 2014-04-16 2014-08-27 孙赞东 基于衰减各向异性的裂缝预测方法
WO2016032353A1 (ru) * 2014-08-26 2016-03-03 Закрытое Акционерное Общество "Технологии Обратных Задач" Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам
CN107678061A (zh) * 2017-10-24 2018-02-09 中国地质大学（北京） 基于横波分裂振幅差属性的流体检测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104007462A (zh) * 2014-04-16 2014-08-27 孙赞东 基于衰减各向异性的裂缝预测方法
WO2016032353A1 (ru) * 2014-08-26 2016-03-03 Закрытое Акционерное Общество "Технологии Обратных Задач" Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам
CN107678061A (zh) * 2017-10-24 2018-02-09 中国地质大学（北京） 基于横波分裂振幅差属性的流体检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XIAO-MING TANG EL: ""Elastic wave scattering to characterize heterogeneities in the borehole environment"", 《GEOPHYSICAL JOURNAL INTERNATIONAL》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112859163A (zh) * 2021-01-07 2021-05-28 中国石油大学(北京) 利用散射波确定压裂区地层品质因子变化的方法和装置

Also Published As

Publication number Publication date
CN108732625B (zh) 2020-02-07

Similar Documents

Publication Publication Date Title
CN104656157B (zh) 一种识别页岩气甜点区的方法及装置
CN108107478B (zh) 大地电磁同步探测与实时反演方法及系统
BRPI0909445B1 (pt) Método para obtenção de anisotropia acústica e formação de imagem por meio de amostragem azimutal de alta resolução
CN102866417A (zh) 一种地下溶洞地震跨孔ct探测及层析成像装置及方法
CN108957548B (zh) 一种多波多分量联合观测地震页岩气富集区预测方法
CN109407144B (zh) 一种基于多波的单孔孤石立体探测方法
CN107479098B (zh) 一种水力压裂过程中同井微地震监测工艺
CN106908837B (zh) 一种压裂裂缝形态和裂缝高度确定方法
CN202837558U (zh) 一种地下溶洞地震跨孔ct探测及层析成像装置
CN105735971B (zh) 一种基于弹性波的钻孔深度检测系统及其检测方法
CN109001839A (zh) 一种城市车载综合地球物理探测系统及工作方法
CN104950327B (zh) 确定地面微地震观测系统的检波器的位置的方法
AU2020101105A4 (en) A method for grounded short-offset transient electromagnetic geophysical surveying
CN105093314B (zh) 一种测定微地震震源的方法
CN102720484B (zh) 一种随钻声波测井仪器及测井方法
CN106324682B (zh) 一种应用于冻土层地区的表层结构调查方法
EP4089446A1 (en) A method for seismic frequency resonance exploration technology
CN102508310B (zh) 一种煤田火区上部地层孔隙度分布的探测方法
CN101100940A (zh) 一种阵列化声信号检测系统及其工程应用
CN108732625A (zh) 一种井周岩石非均匀各向异性的识别方法及系统
RU2650084C2 (ru) Способ мониторингового контроля физического состояния геологической среды
CN106646668B (zh) 一种雷达测井标准井模型的建立方法
CN203630367U (zh) 一种地面探测地质的系统
CN110927632A (zh) 一种频率域水平x方向磁场分量观测及资料处理方法
CN105929445B (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