CN111119851B - 一种非对称远探测测井方法 - Google Patents

一种非对称远探测测井方法 Download PDF

Info

Publication number
CN111119851B
CN111119851B CN201811265490.2A CN201811265490A CN111119851B CN 111119851 B CN111119851 B CN 111119851B CN 201811265490 A CN201811265490 A CN 201811265490A CN 111119851 B CN111119851 B CN 111119851B
Authority
CN
China
Prior art keywords
logging method
asymmetric
wave
emitting electrode
plane
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
Application number
CN201811265490.2A
Other languages
English (en)
Other versions
CN111119851A (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.)
Shengli Logging Co Of Sinopec Jingwei Co ltd
China Petrochemical Corp
Sinopec Oilfield Service Corp
Sinopec Shengli Petroleum Engineering Corp
Sinopec Jingwei Co Ltd
Original Assignee
Shengli Logging Co Of Sinopec Jingwei Co ltd
China Petrochemical Corp
Sinopec Oilfield Service Corp
Sinopec Shengli Petroleum Engineering Corp
Sinopec Jingwei Co Ltd
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 Shengli Logging Co Of Sinopec Jingwei Co ltd, China Petrochemical Corp, Sinopec Oilfield Service Corp, Sinopec Shengli Petroleum Engineering Corp, Sinopec Jingwei Co Ltd filed Critical Shengli Logging Co Of Sinopec Jingwei Co ltd
Priority to CN201811265490.2A priority Critical patent/CN111119851B/zh
Publication of CN111119851A publication Critical patent/CN111119851A/zh
Application granted granted Critical
Publication of CN111119851B publication Critical patent/CN111119851B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

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

Abstract

本发明属于声波测井资料处理及解释方法技术领域,尤其涉及一种非对称远探测测井方法。利用该远探测测井方法可以实现对横波的有效探测,最终实现对不同方位的地质异常体反射信号灵敏区分,为井周全方位反射横波成像提供新的思路。一种非对称测量的远探测测井方法,所述非对称远探测测井方法采用对称偶极阵列作为发射组合,并为每一个发射极匹配设置n个接收极,n为大于等于1的自然数;接收极所处平面与发射极所处平面两者平行设置,且任意一个发射极与其相匹配设置的接收极之间存在有不为0的夹角。

Description

一种非对称远探测测井方法
技术领域
该发明属于声波测井资料处理及解释方法技术领域,尤其涉及一种非对称远探测测井方法。
背景技术
声波测井作为一种较为常用的测井方法,其利用的是声波在不同性质岩层中传播时,速度、幅度及频率等声学特性的变化来测定钻井地质剖面的特性,最终判断得到固井的质量。
举例而言,如图1-3所示,根据井孔声学的基本原理,可以确定得到波动方程及其位移解:
(λ+μ)▽(▽·u)+μ▽2u+ρω2u (1)
Figure GDA0001966484330000011
其中,Φ为纵波位移势函数;
Figure GDA0001966484330000012
为z方向的单位矢量;Γ为SV横波位移势函数;χ为SH横波位移势函数。SV横波在竖直平面内偏振,SH横波在水平平面内偏振。纵波、SV横波和SH横波的波长分别满足下式中的一个方程:
Figure GDA0001966484330000013
其中,kp=ω/α和ks=ω/β分别是纵波、横波波数,α、β是纵波、横波速度,各势函数在波数域的解为:
Figure GDA0001966484330000014
其中
Figure GDA0001966484330000015
Figure GDA0001966484330000016
分别为纵波径向波数和横波径向波数。
反射系数An、Bn、Cn、Dn、En、Fn可以由以下井孔流体和地层应力位移连续来确定。
Figure GDA0001966484330000021
Figure GDA0001966484330000022
Figure GDA0001966484330000023
Figure GDA0001966484330000024
其中,τij是孔隙地层受得总应力,井孔内外用上标中的(1)和(2)来区分。将各相关表达式代入,求解反射系数,即可进一步得到井孔中得到波形。声波测井响应的声压时域波形函数为:
Figure GDA0001966484330000025
而进一步根据反射声学基本原理可知:
声反射探测,记录的反射波信号受一些重要的因素控制:
Figure GDA0001966484330000026
其中,RWV代表接收到的反射波;S是声波仪器的***传递函数,其包括信号产生和记录;RD是井眼辐射指向性;RC是接收方位模式;RF是反射体反射系数。它们都能随角频率ω变化,如图1所示。其中S、RD、RC这些参数与井孔声学密切相关,反射声学决定参数RF。
除了这些因素,反射波主要受沿着波程的传播损耗控制。传播损耗由两个部分组成:几何扩散1/D,其中D是传播到反射体又返回来的总距离;波程的幅度衰减exp(-ωT/2Q),其中T是沿着波程D的旅行时,Q是品质因子。方程(10)表明,如果传播到远处反射体并反射回来的T在记录时间内,并且地层衰减1/D没有把反射波振幅减弱到噪声水平,在适当的辐射和接受情况下,则记录的数据能够对反射体成像。
①发射部分
方程(10)的理论模型中另一个重要的因素是声源辐射特性RD。虽然单极声源的辐射特性已经详细研究过(Meredith,1990),但关于偶极S波反射探测的辐射特性需要更进一步分析。
充满流体井眼中偶极源的辐射是具有方向的。偶极子声波测井常被用于测量地层S波速度以及S波方位各向异性。偶极子声源或接收器***的一个非常有用的特性是方向性:激发的或接收的声波振幅依赖于声波质点运动的方向(偏振方向)与声源或接收器方向之间的角度φ。假定直井情况下,偶极源激发的SV和SH波的方位决定于(Schmitt,1988):
Figure GDA0001966484330000031
其中φ是方位角,θ与垂井的交角;uφ和uθ分别是SH和SV波位移。附图2了SH和SV波偏振方向的例子。
偶极源的最大优势是激发的SH和SV波分别对cosφ和sinφ的方位敏感,这就提供了用偶极源确定反射体方位的基础。
②方位接收
在理论模型中一个重要的因素是井眼对反射波的响应,也就是方程(10)的RC。虽然这种影响已经由Schoenberg(1986)和Peng(1993)详细研究过,但它对反射横波探测的重要性最近才被意识到。图3展示了快地层中井径为0.2m充液井眼的SH和SV波接收模式。在入射平面内SH波响应在所有入射角显示宽模式,这表明井眼可接收对来自各方向的SH波。相比之下,当入射角接近90°,SV波模式接近0°,这意味着不能被探测到垂直入射井眼的SV波。SH和SV波的这种响应特征在下面被用来S波成像分析和现场资料解释。
因此发明人发现,现有远探测测井过程中仍然存在有如下缺陷:在利用横波反射信号对不同方位的地质异常体进行成像时(使用如图4所示对称结构设置的偶极子和接收偶极子时),技术人员无法解决偶极反射横波远探测左右分辨不清问题。
发明内容
本发明提供了一种非对称远探测测井方法,利用该远探测测井方法可以实现对横波的有效探测,最终实现对不同方位的地质异常体反射信号灵敏区分,为井周全方位反射横波成像提供新的思路。
为解决上述技术问题,本发明采用了如下技术方案:
一种非对称测量的远探测测井方法,所述非对称远探测测井方法采用对称偶极阵列作为发射组合,并为每一个发射极匹配设置n个接收极,n为大于等于1的自然数;
接收极所处平面与发射极所处平面两者平行设置,且任意一个发射极与其相匹配设置的接收极之间存在有不为0的夹角。
较为优选的,所述非对称远探测测井方法的应用场景为水平井或大斜度井。
本发明提供了一种非对称远探测测井方法,其中该方法对称偶极阵列作为发射组合,每一个发射极匹配设置有n个接收极,任意一个发射极与其相匹配设置的任意接收极之间存在有不为0的夹角。该非对称远探测测井方法通过采用交错设置发射极与接收极的技术手段,使得横波探测时可产生非对称的偏心效果,从而能够更为灵敏的获取不同方位地质异常体的反射信号。
附图说明
图1为现有技术一种利用多分量偶极测井仪进行横波成像的示意图;
图2为现有技术偶极S波反射探测x向偶极源远场指向性图;
图3为快地层SH和SV入射波的井眼接收方位图;
图4为现有技术发射偶极子和接收偶极子结构对称设置示意图;
图5为本发明提供实施例一的示意图;
图6为实施例一模型;
图7为实施例一发射器状态;
图8为实施例一接收器状态;
图9为实施例二发射器状态;
图10为实施例二研究应力状态点相对位置;
图11为实施例一全波波形;
图12为实施例一反射波波形放大;
图13为实施例二地层对应位置波形。
具体实施方式
本发明提供了一种非对称远探测测井方法,利用该远探测测井方法可以实现对横波的有效探测,最终实现对不同方位的地质异常体反射信号灵敏区分,为井周全方位反射横波成像提供新的思路。
具体的,该非对称远探测测井方法采用对称偶极阵列作为发射组合,并为每一个发射极匹配设置n个接收极,n为大于等于1的自然数;接收极所处平面与发射极所处平面两者平行设置,且任意一个发射极与其相匹配设置的接收极之间存在有不为0的夹角。尤其优选的,该非对称远探测测井方法的应用场景为水平井或大斜度井。
值得注意的是,在实际使用时该非对称远探测测井方法有多种不同应用方式,以下述具体实施例为例进行示例说明。
应用场景一:对称极发射,多分量非对称极接收
参考图5,以对称偶极阵列作为发射组合,并为每一组发射极匹配设置一组接收极为例。具体的,将X1、X2作为一组发射极,将Y1、Y2作为另一组发射极端;与此相对应的,(A1、A4)、(A3、A2)、(A3、A1)、(A2、A4)作为接收极组合,且接收极所处平面(X1、X2、Y1、Y2)与发射极所处平面(A1、A2、A3、A4)两者平行设置;定义仪器中轴点为O点,那么∠A1OX1=∠A3OY1=∠A2OX2=∠A4OY2=θ;其中,θ为非0角度。具体的,当X1、X2为发射极时,(A1、A4)和(A3、A2)用以接收X1、X2的平行分量,(A3、A1)和(A2、A4)用以接收X1、X2的交叉分量;当Y1、Y2为发射极,(A3、A1)和(A2、A4)用以接收X1、X2的平行分量,(A1、A4)和(A3、A2)用以接收X1、X2的交叉分量。由于发射极、接收极是交错设置的,在横波探测时会产生非对称偏心波形,因而可以更为灵敏的反应不同方位地质异常体的反射信号到时和幅度信息,从而进一步计算到时和方位能量,区别区分地质异常体的具体方位。
这里以图6模型为例,图6模型所处地层为硬地层,井外有一裂缝,偶极发射声源产生横波信息如图7所示,采用n2、n4、n6和n8号接收器接收井孔流体的压力信息(见图8),(n2、n8)与(n4、n6)可以形成偏心的偶极接收器(如图11所示);最前面部分是直达波,0.0018s附近开始出现反射波信息。因为接收器距离裂缝的远近各有不同,到时会有差异,n28到时会早于n46(如图12所示)。因而,可以通过到时先后区分地质异常***于仪器的方位,削除成像方位的多解性。对于仪器尺寸相对波长不能忽略时,反射波信号受仪器阻挡作用也较大,正对反射体反射信号幅度强,远离反射体反射信号幅度弱,也可以通过反射信号幅度来判断地质异常体方位。在地质异常体方位未知时,可以通过两两正交分量进行合成,形成相应方位的接收波形。
应用场景二:偏心极发射,正交分量极接收
同样参考图5,采用偏心发射极作为发射组合,并为每一组发射极匹配设置一组接收极。具体的,将(Y1、X1)的组合、(X2、Y1)的组合、(Y2、X2)的组合、(X1、Y2)的组合分别作为发射极端;与此相对应的,(A2、A1)和(A3、A2)分别为接收极组合,并且接收极所处平面(X1、X2、Y1、Y2)与发射极所处平面(A1、A2、A3、A4)两者平行设置,定义仪器中轴点为O点,有∠A1OX1=∠A3OY1=∠A2OX2=∠A4OY2=θ,其中,θ为非0角度。具体的,当X1、X2作为发射极时,(A1、A4)和(A3、A2)用以接收X1、X2的平行分量,(A3、A1)和(A2、A4)用以接收X1、X2的交叉分量;当Y1、Y2作为发射极时,(A3、A1)和(A2、A4)用以接收X1、X2的平行分量,(A1、A4)和(A3、A2)用以接收X1、X2的交叉分量。同样的,由于发射极、接收极是交错设置的,在横波发射时会产生非对称偏心波形,因而可以更为灵敏的反应不同方位地质异常体的反射信号到时和幅度信息,从而进一步计算到时和方位能量,区别区分地质异常体的具体方位。
这里以图9模型为例,图9模型所处地层为硬地层,发射声源为偏心偶极声源(见图9)。为研究仪器发射能量进入地层情况,在井外4m处布置n1-n8切应变探测点;n1和n5所对应波形如图13所示,图中n1到偏心发射器较近,其到时明显要早于n5;而从反射体到接收,因为对称设置,反射信号传递过程两侧变化一致;而由于偏心发射时靠近异常体一侧传播距离小,波形较远离异常体一侧传播距离小,从而可以利用反射信号到时区分地质异常体方位。同时当随钻等仪器尺寸与波场相比不可忽略时,近发射端到达地层相同位置信号强,远发射端到达地层相同位置信号相对弱。因而,也可以通过对比到时和振幅进行异常体方位的确定。另外,若异常体走向与发射方向不一致,可以通过对应方位正交分量进行合成。
需要解释说明的是,限于篇幅有限,发射极、接收极交错设置的方式以及具体夹角情况无法一一列举,本领域技术人员可根据实际情况按需设置,在此不做赘述。
本发明提供了一种非对称远探测测井方法,其中该方法对称偶极阵列作为发射组合,每一个发射极匹配设置有n个接收极,任意一个发射极与其相匹配设置的任意接收极之间存在有不为0的夹角。该非对称远探测测井方法通过采用交错设置发射极与接收极的技术手段,使得横波探测时可产生非对称的偏心效果,从而能够更为灵敏的获取不同方位地质异常体的反射信号。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (2)

1.一种非对称远探测测井方法,其特征在于,所述非对称远探测测井方法采用对称偶极阵列作为发射组合,并为每一个发射极匹配设置n个接收极,n为大于等于1的自然数;
所述发射组合包括两组具有正交关系的对称偶极阵列,接收极所处平面与发射极所处平面两者平行设置,且任意一个发射极与其相匹配设置的接收极之间存在有不为0的夹角,且所有接收极处于同一平面,每个发射极均与所有接收极交错设置,其中,
通过对称极发射且多分量非对称极接收或者偏心极发射且正交分量极接收的方式,使得在横波探测时产生非对称偏心波形,以反应不同方位地质异常体的反射信号到时和幅度信息,从而通过获得计算到时和方位能量来区别区分地质异常体的具体方位。
2.根据权利要求1所述的一种非对称远探测测井方法,其特征在于,所述非对称远探测测井方法的应用场景为水平井或大斜度井。
CN201811265490.2A 2018-10-29 2018-10-29 一种非对称远探测测井方法 Active CN111119851B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811265490.2A CN111119851B (zh) 2018-10-29 2018-10-29 一种非对称远探测测井方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811265490.2A CN111119851B (zh) 2018-10-29 2018-10-29 一种非对称远探测测井方法

Publications (2)

Publication Number Publication Date
CN111119851A CN111119851A (zh) 2020-05-08
CN111119851B true CN111119851B (zh) 2023-03-14

Family

ID=70483914

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811265490.2A Active CN111119851B (zh) 2018-10-29 2018-10-29 一种非对称远探测测井方法

Country Status (1)

Country Link
CN (1) CN111119851B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA1295725C (en) * 1989-05-18 1992-02-11 Sen-Tsuen Chen Method and system for measuring azimuthal anisotropy effects using acoustic multipole transducers
CN1690729A (zh) * 2000-11-13 2005-11-02 贝克休斯公司 声学测井装置及钻井环箍上使用声学测井装置的方法
CN101029565A (zh) * 2006-11-30 2007-09-05 中国石油天然气集团公司 裸眼井方位声波测井方法
CN102508299A (zh) * 2011-09-29 2012-06-20 中国石油大学(华东) 在井中发射并接收偶极横波的远探测方法
CN102720484A (zh) * 2012-05-18 2012-10-10 中国海洋石油总公司 一种随钻声波测井仪器及测井方法
CN103603656A (zh) * 2013-08-21 2014-02-26 中国石油大学(北京) 一种基于相控圆弧阵的声波测井方位接收方法及装置
CN104481526A (zh) * 2014-12-08 2015-04-01 中国科学院声学研究所 一种双交叉偶极子测井方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA1295725C (en) * 1989-05-18 1992-02-11 Sen-Tsuen Chen Method and system for measuring azimuthal anisotropy effects using acoustic multipole transducers
CN1690729A (zh) * 2000-11-13 2005-11-02 贝克休斯公司 声学测井装置及钻井环箍上使用声学测井装置的方法
CN101029565A (zh) * 2006-11-30 2007-09-05 中国石油天然气集团公司 裸眼井方位声波测井方法
CN102508299A (zh) * 2011-09-29 2012-06-20 中国石油大学(华东) 在井中发射并接收偶极横波的远探测方法
CN102720484A (zh) * 2012-05-18 2012-10-10 中国海洋石油总公司 一种随钻声波测井仪器及测井方法
CN103603656A (zh) * 2013-08-21 2014-02-26 中国石油大学(北京) 一种基于相控圆弧阵的声波测井方位接收方法及装置
CN104481526A (zh) * 2014-12-08 2015-04-01 中国科学院声学研究所 一种双交叉偶极子测井方法

Also Published As

Publication number Publication date
CN111119851A (zh) 2020-05-08

Similar Documents

Publication Publication Date Title
US7392135B2 (en) Adaptive equalization of downhole acoustic receivers
US6714480B2 (en) Determination of anisotropic moduli of earth formations
US7623412B2 (en) Anisotropy measurement while drilling
US8755248B2 (en) Unipole and bipole acoustic logging while drilling tools
US20070153628A1 (en) Enhanced noise cancellation in VSP type measurements
US8797825B2 (en) Method and apparatus for measuring formation anisotropy while drilling
US10782432B2 (en) Monopole acoustic logging while drilling instrument used together with bottom hole assembly, method for measuring shear wave velocity of slow formations
CN202170793U (zh) 一种随钻声波测井装置和发射换能器
US4713968A (en) Method and apparatus for measuring the mechanical anisotropy of a material
US6188961B1 (en) Acoustic logging apparatus and method
US9494704B2 (en) Estimating formation stresses using sonic data
US11808908B2 (en) Real-time reconfiguration of phased array operation
EP2304474A2 (en) Method and system for generating a beam of acoustic energy from a borehole, and applications thereof
US10768329B2 (en) Seismic sensing systems and processes for using same
US10605944B2 (en) Formation acoustic property measurement with beam-angled transducer array
CN101029565B (zh) 裸眼井方位声波测井方法
US9030911B2 (en) Method and system for delineating a second wellbore from a first wellbore
Che et al. An experimental study on azimuthal reception characteristics of acoustic well-logging transducers based on phased-arc arrays
CN103558637B (zh) 基于三分量传感器的远探测方法
US6176344B1 (en) Method and system for determining the azimuth position and distance of a reflecting subsurface formation
CN111119851B (zh) 一种非对称远探测测井方法
Walker et al. Logging services: Towards the 3D measurement of formation properties in high-resolution with a continuous depth of investigation
Che et al. Azimuthal cement evaluation with an acoustic phased-arc array transmitter: numerical simulations and field tests
CN114779346A (zh) 一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置
Boonen et al. A Dual-Frequency Lwd Sonic Tool Expands Existing Unipolar Transmitter Technology To Supply Shear Wave Data In Soft Formations.

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
TA01 Transfer of patent application right

Effective date of registration: 20200630

Address after: 100027 Beijing, Chaoyangmen, North Street, No. 22, No.

Applicant after: SINOPEC Group

Applicant after: SINOPEC OILFIELD SERVICE Corp.

Applicant after: SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

Applicant after: WELL LOGGING COMPANY, SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

Address before: 100101 Beijing, Beichen West Road, No. 8 Beichen world center, block A, No. 703, No.

Applicant before: SINOPEC OILFIELD SERVICE Corp.

Applicant before: SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

Applicant before: WELL LOGGING COMPANY, SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20220218

Address after: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen

Applicant after: SINOPEC Group

Applicant after: SINOPEC OILFIELD SERVICE Corp.

Applicant after: SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

Applicant after: Sinopec Jingwei Co.,Ltd.

Applicant after: Shengli logging company of Sinopec Jingwei Co.,Ltd.

Address before: 100027 Chaoyangmen North Street, Chaoyang District, Chaoyang District, Beijing

Applicant before: SINOPEC Group

Applicant before: SINOPEC OILFIELD SERVICE Corp.

Applicant before: SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

Applicant before: WELL LOGGING COMPANY, SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant