CN107037484A - 一种提取偶极反射横波的方法 - Google Patents

一种提取偶极反射横波的方法 Download PDF

Info

Publication number
CN107037484A
CN107037484A CN201610079487.6A CN201610079487A CN107037484A CN 107037484 A CN107037484 A CN 107037484A CN 201610079487 A CN201610079487 A CN 201610079487A CN 107037484 A CN107037484 A CN 107037484A
Authority
CN
China
Prior art keywords
wave
slowness
reflection
dipole
stratum
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
CN201610079487.6A
Other languages
English (en)
Other versions
CN107037484B (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
Logging Co Of Triumph Petroleum Engineering Co Ltd Of China Petrochemical Industry
Sinopec Oilfield Service Corp
Sinopec Shengli Petroleum Engineering Corp
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 Logging Co Of Triumph Petroleum Engineering Co Ltd Of China Petrochemical Industry, Sinopec Oilfield Service Corp, Sinopec Shengli Petroleum Engineering Corp filed Critical Logging Co Of Triumph Petroleum Engineering Co Ltd Of China Petrochemical Industry
Priority to CN201610079487.6A priority Critical patent/CN107037484B/zh
Publication of CN107037484A publication Critical patent/CN107037484A/zh
Application granted granted Critical
Publication of CN107037484B publication Critical patent/CN107037484B/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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/21—Frequency-domain filtering, e.g. band pass

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种提取偶极反射横波的方法，涉及应用地球物理测井及勘探技术领域，该方法基于频率以及慢度差异，去除了直达弯曲波、漏能p波、反射斯通利波等噪声波，保留了有用的偶极反射横波。一种提取偶极反射横波的方法，包括：对偶极阵列声波进行低通滤波；对偶极阵列声波进行组合，得到共发射器组合以及共接受器组合，使用频率波数域滤波方法将共发射器组合中不同慢度波场进行开窗滤波，形成下行反射横波；使用频率波数域滤波方法将共接收器组合中不同慢度波场进行开窗滤波，形成上行反射横波；对偶极阵列声波进行组合，得到共源距组合，使用频率波数域滤波方法将共源距组合中不同慢度波场进行开窗滤波，形成上行以及下行反射横波。

Description

一种提取偶极反射横波的方法
技术领域
本发明涉及应用地球物理测井及勘探技术领域，尤其涉及一种提取偶极反射横波的方法。
背景技术
本领域技术人员在偶极反射横波测井成像过程中发现，在接收到的偶极反射横波波形中，还存在有直达弯曲波、漏能p波、斯通利波等噪声波形，且上述噪声波形振幅较高。而现有技术未能对上述噪声波作有效的处理，使得反射横波其信噪比较低，最终导致反射成像结果中具有较强的干扰。因此，设计一种可以消除上述噪声波的方法就显得十分必要。
发明内容
本发明提出了一种提取偶极反射横波的方法，该方法基于频率以及慢度差异，去除了直达弯曲波、漏能p波、反射斯通利波等噪声波，保留了有用的偶极反射横波。
为解决上述技术问题，本发明的实施例采用如下技术方案：
一种提取偶极反射横波的方法，包括：
步骤(1)：对偶极阵列声波进行低通滤波，所述低通滤波的截止频率为4.7kHz；
步骤(2)：判断地层与井轴是否平行；若地层与井轴平行，则继续执行步骤(3)；若地层与井轴不平行，则继续执行步骤(4)；
步骤(3)：对偶极阵列声波进行组合，得到共发射器组合以及共接受器组合；
使用频率波数域滤波方法将共发射器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成下行反射横波；
使用频率波数域滤波方法将共接收器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成上行反射横波；
步骤(4)：对偶极阵列声波进行组合，得到共源距组合；
使用频率波数域滤波方法将共源距组合中不同慢度波场进行开窗滤波，保留慢度范围为-2s～0的波场，形成上行反射横波；保留慢度范围为0～2s的波场，形成下行反射横波；
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
优选的，当地层与井轴平行且偶极反射横波检测仪器位在界面上方时，偶极反射横波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
优选的，当地层与井轴平行且当地层与井轴平行且仪器位于界面下方时，偶极反射横波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
优选的，当地层与井轴不平行且当地层与井轴平行且仪器位于界面上方时，偶极反射横波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
优选的，当地层与井轴不平行且当地层与井轴平行且仪器位于界面下方时，偶极反射横波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
本发明提供了一种提取偶极反射横波的方法，该方法通过频率域滤波以及速度域滤波方式，消除了对成像干扰较大的漏能p波、直达波以及反射波(尤其是指斯通利波)，使得反射横波的信噪比大幅提高。
附图说明
图1为本发明提取偶极反射横波方法的流程图；
图2为进行低通滤波处理后形成的波形图；
图3为偶极反射横波发射仪器与界面的示意图(偶极反射横波发射仪器位于界面上方)；
图4为偶极反射横波发射仪器与界面的示意图(偶极反射横波发射仪器位于界面下方)；
图5为根据共发射器组合形成的下行反射横波的波形图；
图6为偶极反射横波接收仪器与界面的示意图(偶极反射横波接收仪器位于界面上方)；
图7为根据共接收器组合形成的上行反射横波的波形图；
图8为偶极反射横波检测仪器与界面的示意图(偶极反射横波检测仪器按照共源距进行组合后，位于界面上方)。
具体实施方式
本发明提出了一种提取偶极反射横波的方法，该方法基于频率以及慢度差异，去除了直达弯曲波、漏能p波、反射斯通利波等噪声波，保留了有用的偶极反射横波。
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合下述附图对本发明实施例做详细描述。
本发明提供一种提取偶极反射横波的方法，如图1所示，该方法包括：
步骤(1)：对偶极阵列声波进行低通滤波，低通滤波的截止频率为4.7kHz。
在提取偶极反射横波的过程中，首先对偶极阵列声波进行低通滤波，该低通滤波的截止频率为4.7kHz。通过该低通滤波，消除噪声漏能p波影响；尤其对于软底层来说，绝大多说漏能p波可被消除。其中，图2为进行低通滤波处理后形成的波形图。
步骤(2)：判断地层与井轴是否平行；
其中，当地层与井轴平行时，对偶极阵列声波进行组合，分别得到共发射器组合和共接收器组合，并进一步根据共发射器组合和共接收器组合形成上行反射横波以及下行反射横波(即继续执行步骤3)；若地层与井轴不平行时，则对偶极阵列声波进行组合，得到共源距组合，并根据共源距组合形成上行反射横波以及下行反射横波(即继续执行步骤4)。
步骤(3)：对偶极阵列声波进行组合，得到共发射器组合以及共接受器组合；
使用频率波数域滤波方法将共发射器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成下行反射横波；
使用频率波数域滤波方法将共接收器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成上行反射横波。
具体的，以图3或图4为例，对共发射器组合进行开窗滤波处理。首先假设偶极反射横波发射仪器位于界面的上方，将界面与井轴夹角定义为α，发射器与界面距离定义为Z，发射器与第一个接收器距离定义为z，地层速度定义为v，慢度定义为s。如图3所示，其中发射器以三角符号表示，接收器以圆形符号表示。
由基本的三角关系可以很容易得出反射波的传播距离为且两接收器接收的反射波传播路程小于接收器间距dz，即接收到波形的到达时间满足：其中v为在地层中传播的横波速度。
相邻两接收器接收到波形视横波速度为：
其中dz为接收器的间距，当dz→0时，上式(2)为精确视速度，表达式为：
即视横波慢度为：
而当偶极反射横波发射仪器移动至界面以下时，如图4所示，(2)式和(3)式中反射波由上行波变为下行波，公式相应变化如下：
其中dz为接收器的间距，当dz→0时，上式为精确视速度，表达式为：
视横波慢度为：
值得注意的是，其中，上述公式中α其取值范围为0≤α≤90°，也就是说满足：
对于而言，当α＝0时，其取值最小，当α＝90°时，其取得最大值。
因此，s01′或s02′满足：
另一方面，反射斯通利波-sst其慢度范围满足：-sst<-sf；而Sf为井孔流体声速，其满足：-sf<-s。因此通过公式(9)可以确定，通过开窗滤波手段，对一定慢度范围内的共发射器组合(具体的，即保留慢度为内的共发射器组合信息)进行滤波，可将反射斯通利波噪声消除，从而达到提取下行反射横波信息的目的。其中，形成的下行反射横波的波形图如图5所示。
同样道理，对共接收器组合进行开窗滤波处理，如图6所示，其中当偶极反射横波接收仪器在界面的上方时，将其视横波速度表示为s01″；当偶极反射横波接收仪器在界面的下方时，将其视横波速度表示为s02″。
与上述共发射器组合视横波慢度推导过程相类似，最终计算可得s01″或s02″满足：
同样的，由于斯通利波sst其慢度范围满足：sf<sst；Sf为井孔流体声速，其满足：s<sf。因此，通过开窗滤波手段，对一定慢度范围内的共接收器组合(具体的，即保留慢度为内的共接收器组合信息)进行滤波，可将斯通利波噪声消除，从而达到提取上行反射横波信息的目的。其中，形成的上行反射横波的波形图如图7所示。
至此，在地层与井轴平行的情况下，通过共发射器组合以及共接受器组合，通过开窗滤波处理，分别形成了下行反射横波以及上行反射横波，且消除了斯通利波噪声的影响。
步骤(4)：对偶极阵列声波进行组合，得到共源距组合；
使用频率波数域滤波方法将共源距组合中不同慢度波场进行开窗滤波，保留慢度范围为-2s～0的波场，形成上行反射横波；保留慢度范围为0～2s的波场，形成下行反射横波；
具体的，步骤(3)给出的是地层与井轴平行时，形成上行反射横波以及下行反射横波的过程；而步骤(4)则给出的是地层与井轴不平行时形成上行反射横波以及下行反射横波的过程。
如图8所示，将偶极反射横波检测仪器按照共源距进行组合，并将不同深度的相同源距波形组成道集。以共源距组合作为参考，在仪器位于界面上方时，偶极反射横波的视慢度为：
通过计算可以确定的是，在未知地层信息下，s03′满足：0≤s03′≤2s；此时直达波慢度均为0，上行反射斯通利波视慢度为sst′＝2sst>2s。因此，通过开窗滤波手段，对一定慢度范围内的共源距组合(具体的，即保留慢度为0～2s内的共源距组合信息)进行滤波，可将斯通利波以及直达波等噪声消除，从而达到提取上行反射横波信息的目的。
与公式(11)相类似的是，当仪器位于界面下方时，偶极反射横波慢度满足公式：
通过计算可以确定的是，在未知地层信息下，s03″满足：-2s≤s03″≤0；此时直达波慢度均为0，下行反射斯通利波视慢度为sst″＝-2sst<-2s。因此，通过开窗滤波手段，对一定慢度范围内的共源距组合(具体的，即保留慢度为-2s～0内的共源距组合信息)进行滤波，可将斯通利波以及直达波等噪声消除，从而达到提取下行反射横波信息的目的。
至此，在地层与井轴不平行的情况下，通过共源距组合，通过开窗滤波处理，分别形成了下行反射横波以及上行反射横波，且消除了斯通利波噪声的影响。
本发明提供了一种提取偶极反射横波的方法，该方法通过频率域滤波以及速度域滤波方式，消除了对成像干扰较大的漏能p波、直达波以及反射波(尤其是指斯通利波)，使得反射横波的信噪比大幅提高。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (5)

1.一种提取偶极反射横波的方法，其特征在于，包括：
步骤(1)：对偶极阵列声波进行低通滤波，所述低通滤波的截止频率为4.7kHz；
步骤(2)：判断地层与井轴是否平行；若地层与井轴平行，则继续执行步骤(3)；若地层与井轴不平行，则继续执行步骤(4)；
步骤(3)：对偶极阵列声波进行组合，得到共发射器组合以及共接受器组合；
使用频率波数域滤波方法将共发射器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成下行反射横波；
使用频率波数域滤波方法将共接收器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成上行反射横波；
步骤(4)：对偶极阵列声波进行组合，得到共源距组合；
使用频率波数域滤波方法将共源距组合中不同慢度波场进行开窗滤波，保留慢度范围为-2s～0的波场，形成上行反射横波；保留慢度范围为0～2s的波场，形成下行反射横波；
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
2.根据权利要求1所述的一种提取偶极反射横波的方法，其特征在于，
当地层与井轴平行且偶极反射横波检测仪器位在界面上方时，偶极反射横 波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
3.根据权利要求1所述的一种提取偶极反射横波的方法，其特征在于，
当地层与井轴平行且当地层与井轴平行且仪器位于界面下方时，偶极反射 横波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
4.根据权利要求1所述的一种提取偶极反射横波的方法，其特征在于，
当地层与井轴不平行且当地层与井轴平行且仪器位于界面上方时，偶极反 射横波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
5.根据权利要求1所述的一种提取偶极反射横波的方法，其特征在于，
当地层与井轴不平行且当地层与井轴平行且仪器位于界面下方时，偶极反 射横波慢度满足公式：
其中，s表示慢度，z表示发射器与第一个接收器之间的距离，Z表示发射器与界面之间的距离。
CN201610079487.6A 2016-02-04 2016-02-04 一种提取偶极反射横波的方法 Active CN107037484B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610079487.6A CN107037484B (zh) 2016-02-04 2016-02-04 一种提取偶极反射横波的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610079487.6A CN107037484B (zh) 2016-02-04 2016-02-04 一种提取偶极反射横波的方法

Publications (2)

Publication Number Publication Date
CN107037484A true CN107037484A (zh) 2017-08-11
CN107037484B CN107037484B (zh) 2020-02-21

Family

ID=59532359

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610079487.6A Active CN107037484B (zh) 2016-02-04 2016-02-04 一种提取偶极反射横波的方法

Country Status (1)

Country Link
CN (1) CN107037484B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100208552A1 (en) * 2009-02-13 2010-08-19 Camwell Paul L Acoustic telemetry stacked-ring wave delay isolator system and method
CN102012524A (zh) * 2010-09-29 2011-04-13 中国海洋石油总公司 一种海上三维地震观测系统羽状漂移定量评估方法
CN102830433A (zh) * 2012-08-16 2012-12-19 中国石油化工股份有限公司 基于频率域用偶极横波测井资料计算各向异性方位角的方法
CN103437756A (zh) * 2013-09-12 2013-12-11 北京环鼎科技有限责任公司 一种偶极子声波测井仪
CN103439739A (zh) * 2013-04-08 2013-12-11 中国石油集团东方地球物理勘探有限责任公司 地球物理勘探用纵横波匹配方法及匹配装置
CN103726836A (zh) * 2012-10-12 2014-04-16 中国石油集团长城钻探工程有限公司 基于声波测井资料提取模式波慢度的方法
CN204877453U (zh) * 2015-07-02 2015-12-16 中石化石油工程技术服务有限公司 一种全井眼扫描偶极声波阵列声系

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100208552A1 (en) * 2009-02-13 2010-08-19 Camwell Paul L Acoustic telemetry stacked-ring wave delay isolator system and method
CN102012524A (zh) * 2010-09-29 2011-04-13 中国海洋石油总公司 一种海上三维地震观测系统羽状漂移定量评估方法
CN102830433A (zh) * 2012-08-16 2012-12-19 中国石油化工股份有限公司 基于频率域用偶极横波测井资料计算各向异性方位角的方法
CN103726836A (zh) * 2012-10-12 2014-04-16 中国石油集团长城钻探工程有限公司 基于声波测井资料提取模式波慢度的方法
CN103439739A (zh) * 2013-04-08 2013-12-11 中国石油集团东方地球物理勘探有限责任公司 地球物理勘探用纵横波匹配方法及匹配装置
CN103437756A (zh) * 2013-09-12 2013-12-11 北京环鼎科技有限责任公司 一种偶极子声波测井仪
CN204877453U (zh) * 2015-07-02 2015-12-16 中石化石油工程技术服务有限公司 一种全井眼扫描偶极声波阵列声系

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
庄春喜 等: "偶极横波远探测测井数据处理及应用", 《测井技术》 *

Also Published As

Publication number Publication date
CN107037484B (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
Tang Imaging near-borehole structure using directional acoustic-wave measurement
CN109143374B (zh) 一种井周散射体成像方法及系统
CN110529087B (zh) 一种评价地层水力压裂效果的方法及装置
CN104216008A (zh) 一种井中压裂微地震事件识别方法
CN106526678B (zh) 一种反射声波测井的波场分离方法及装置
CN103376469B (zh) 一种基于超声成像测井的裂缝定量评价方法
RU2015134848A (ru) Удаление псевдоэхосигналов из акустических волн
CN104265277B (zh) 一种利用管波与地层声波干涉原理提取地层声速的方法
CN113156515A (zh) 一种声波远探测成像降噪处理方法和装置
CN106950600B (zh) 一种近地表散射面波的去除方法
CN103913770A (zh) 基于vsp资料对地震数据进行处理的方法
CN109164492B (zh) 一种提取套管井地层声波速度的方法
Liu et al. Detection of karst cavity beneath cast-in-place pile using the instantaneous phase difference of two receiver recordings
Wang et al. Generalized collar waves in acoustic logging while drilling
CN105301658A (zh) 利用微侧井记录提取近地表反褶积算子的方法
CN104237944B (zh) 一种适用于交错网格有限差分的全吸收pml方法
CN107037484A (zh) 一种提取偶极反射横波的方法
CN107037483B (zh) 一种基于时间慢度域的偶极反射横波三维成像方法
Tao et al. Processing of array sonic logging data with multi-scale STC technique
NO20210922A1 (en) Waveform processing utilizing an amplitude adaptive data mask
CN105158799B (zh) 基于相位准则的物理探测回波信号识别与提取方法
CN112764109B (zh) 偶极横波反射波的分离提取方法及装置
CN105301653A (zh) 一种多因素联合识别的钻机噪音压制方法
CN112817042B (zh) 沉积储层中油气识别方法及装置
CN104459772A (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
TR01 Transfer of patent right

Effective date of registration: 20200629

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

Co-patentee after: SINOPEC OILFIELD SERVICE Corp.

Patentee after: SINOPEC Group

Co-patentee after: SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

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

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

Co-patentee before: SINOPEC SHENGLI PETROLEUM ENGINEERING Co.,Ltd.

Patentee before: SINOPEC OILFIELD SERVICE Corp.

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

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220221

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

Patentee after: SINOPEC Group

Patentee after: SINOPEC OILFIELD SERVICE Corp.

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

Patentee after: Sinopec Jingwei Co.,Ltd.

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

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

Patentee before: SINOPEC Group

Patentee before: SINOPEC OILFIELD SERVICE Corp.

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

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

TR01 Transfer of patent right