CN107045145B - 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法 - Google Patents

地震层序控制下的叠前振幅随偏移距变化检测缝洞方法 Download PDF

Info

Publication number
CN107045145B
CN107045145B CN201611254793.5A CN201611254793A CN107045145B CN 107045145 B CN107045145 B CN 107045145B CN 201611254793 A CN201611254793 A CN 201611254793A CN 107045145 B CN107045145 B CN 107045145B
Authority
CN
China
Prior art keywords
seismic
fracture
avo
fracture hole
attribute
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
CN201611254793.5A
Other languages
English (en)
Other versions
CN107045145A (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
BGP Inc
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 BGP Inc filed Critical BGP Inc
Priority to CN201611254793.5A priority Critical patent/CN107045145B/zh
Publication of CN107045145A publication Critical patent/CN107045145A/zh
Application granted granted Critical
Publication of CN107045145B publication Critical patent/CN107045145B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging

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

本发明提供了一种地震层序控制下的叠前振幅随偏移距变化检测缝洞方法。所述方法通过地震资料解释,建立研究区主要目的层的地震层序格架;利用测井资料建立关键井AVO正演模型,分析缝洞储层AVO响应特征,提取共反射点道集的截距、梯度属性;对截距、梯度属性沿地震层序框架进行方差体提取,得目的层段AVO属性方差体;根据缝洞储层标定结果,划定截距、梯度识别缝洞的方差门槛值,选择单个属性或两个属性交会识别门槛值以上缝洞发育体系,建立缝洞与AVO属性方差体关系。本发明能够充分利用AVO响应和地震等时地层格架的横向分辨能力,不仅可快速识别缝洞发育有利区,且可间接识别油气富集区,并且具有快速、有效检测特点。

Description

地震层序控制下的叠前振幅随偏移距变化检测缝洞方法
技术领域
本发明涉及地球物理勘探技术中的储层预测领域,具体来讲提供一种地震层序控制下的叠前振幅随偏移距变化(叠前AVO)检测缝洞方法,本发明的地震缝洞检测方法能够利用三维地震叠前AVO道集提取的AVO属性,在地震层序控制下生成梯度、截距属性方差体,进而判别缝洞发育情况。
背景技术
目前,叠前AVO裂缝检测是基于垂直裂缝发育带的存在会引起地下介质的方位各向异性,地下介质的方位各向异性要引起地震波特征的方位各向异性,可以通过地震资料提取地震波特征(速度、振幅、衰减)的方位各向异性来达到检测垂直裂缝发育带的目的,主要有AVA技术和及FVA技术。
AVA是指振幅随方位的变化,主要利用反射振幅随入射方位的变化特征来检测裂缝发育带,通过求取某点各个方位的振幅信息把它放到座标中,如果各个方位的振幅信息都相同,表现为圆形,则该点是各向同性的,不存在高角度裂缝;如果各个方位的振幅信息不相同,呈现椭圆特征,则该点是各向异性的,且椭圆的扁率代表裂缝密度大小。FVA是指纵波频率属性随入射方位的变化,是利用三维纵波的频率属性随入射方位的变化特征来检测裂缝发育带。根据裂缝储层的地震散射理论研究,地震衰减和裂缝密度场的空间变化有关,沿裂缝走向传播时衰减慢,而垂直裂缝方向衰减快,因此依据裂缝引起的地震波高频衰减的方位各向异性所导致的地震反射频率的方位各向异性特征来检测裂缝。
叠前AVO(Amplitude variation with offset,振幅随偏移距变化)裂缝预测技术是利用多个不同方位角的叠加数据体,提取各方位角叠加数据的振幅和频率属性,分析其在同一位置各方位角的频率属性变化,根据井点裂缝标定结果,判别地层缝洞发育的叠前裂缝检测技术。通常包括如下步骤:
步骤一:偏移距的选择,根据面元小、覆盖次数高且分布均匀、各方位偏移距分布均匀的原则按实际地震资料的偏移距-方位角分布选取适当的偏移距范围数据;
步骤二:方位划分,根据断裂主要方位和地震数据覆盖次数均匀的原则划分数据叠加方位,并QC(质量控制)各方位叠加数据体在裂缝发育位置是否具有相应振幅、频率异常响应;
步骤三:地震波场正演模拟,采用地震波场数值模拟技术,模拟裂缝的方向、密度和所含流体变化产生的地震响应特征,确定可用于解决研究区地下裂缝的地震属性。
步骤四:属性优选,提取了目的层能量、频率相关等多种属性,在剖面和平面上观察裂缝发育段的属性响应,选择敏感属性半定量识别裂缝密度、强度和方位信息。
目前的叠前AVO裂缝检测在利用叠前分方位角数据进行裂缝识别时,要求三维地震资料是宽方位采集,及各方位角地震数据的入射角和覆盖次数是均匀分布的,但受采集成本等多方面因数限制,一般地震资料的采集长宽比大于0.67就被称为宽方位资料,但具体在叠加方位角数据时,为了保证各方位叠加数据的入射角和覆盖次数均匀分布,减少地震资料的不均匀带来的影响,一般都会舍去部分小入射角和大入射角数据,尤其是长宽比在0.7左右的地震资料,舍去的大入射角数据更多一些,如果根据测井建立的波动方程正演模型指示的裂缝响应是在中、大入射角更明显时,根据该方法预测的裂缝效果就会被大打折扣,影响对缝洞系统的真实表征。
发明内容
本发明的目的在于解决现有技术存在的上述不足中的至少一项。
本发明提供了一种地震层序控制下的叠前振幅随偏移距变化检测缝洞方法。所述检测缝洞方法包括以下步骤:通过地震资料解释,建立研究区主要目的层的地震层序格架;利用测井资料建立关键井AVO正演模型,分析缝洞储层的AVO响应特征,提取共反射点道集的截距、梯度属性;对截距、梯度属性沿地震层序框架进行方差体提取,得到目的层段的AVO属性方差体;根据关键井缝洞储层标定结果,划定截距、梯度识别缝洞的方差门槛值,选择单个属性或两个属性交会识别门槛值以上缝洞发育体系,建立缝洞与AVO属性方差体关系。
在本发明的一个示例性实施例中,所述地震资料解释的步骤还可包括根据反射终止关系,识别地震层序。例如,反射终止关系可包括削蚀、顶超、下超和/或上超。
在本发明的一个示例性实施例中,所述测井资料可包括纵波测井数据和/或横波测井数据。
在本发明的一个示例性实施例中,所述检测缝洞方法可包括在所述建立缝洞与AVO属性方差体关系的步骤之后,在平面上和剖面上预测缝洞发育情况,绘制目的层段缝洞发育带预测平面图。
在本发明的一个示例性实施例中,所述建立研究区主要目的层的地震层序格架的步骤和所述建立关键井AVO正演模型的步骤无先后顺序之分。
与现有技术相比,本发明的有益效果包括:能够充分利用AVO响应和地震等时地层格架的横向分辨能力,不仅可以快速识别缝洞发育有利区,而且由于油气充注的缝洞体会进一步增大其AVO响应,也可以间接识别油气富集区,与生产结合紧密,具有快速、有效检测特点;理论体系完整可靠,流程清楚,使缝洞检测过程更加直接;能够适用三维地震资料,具有一定钻井数目区块(最好保证有缝洞发育有差异的模式井2口或以上)的缝洞预测;易于实施,具有极大的推广和实用价值。
附图说明
图1示出了本发明的一个示例性实施例的地震层序控制下的叠前AVO检测缝洞方法的流程图;
图2示出了本发明的一个示例性实施例的缝洞发育井AVO正演模型和井旁CRP道集的响应图;
图3示出了本发明的一个示例性实施例的地震资料精细解释剖面图;
图4示出了本发明的一个示例性实施例的连井AVO属性剖面图;
图5示出了本发明的一个示例性实施例的AVO梯度属性方差体连井剖面图。
具体实施方式
在下文中,将结合示例性实施例来详细说明本发明的地震层序控制下的叠前振幅随偏移距变化(叠前AVO)检测缝洞方法。
总体来讲,本发明的地震层序控制下的叠前AVO检测缝洞方法,能够在精细地震解释的基础上,建立地震层序格架,利用叠前共反射点道集提取梯度、截距AVO属性,在地震层序格架下提取沿层梯度、截距方差体。本发明能够充分利用各入射角信息,判别缝洞发育情况,而且可间接反映含油气的缝洞体储层发育情况,其值高低是缝洞发育的指示指标,根据缝洞标定门槛值,超过该门槛值表示缝洞发育,值越大指示其缝洞密度越高,油气饱和度越大。
在本发明的一个示例性实施例中,地震层序控制下的叠前AVO检测缝洞方法可以包括如下步骤:
(1)以地震资料为基础,根据削蚀、顶超、下超和上超等反射终止关系,识别地震层序,建立研究区目的层段的地震层序格架。
(2)利用纵、横波测井数据建立AVO正演模型,分析缝洞储层的AVO响应。提取梯度、截距等AVO属性体。
(3)地震层序格架下,对梯度、截距开展目的层段方差体分析,提取沿层梯度、截距沿层方差属性,依据测井缝洞标定结果,选取其中一种或两种交会划定方差体识别门槛值,根据方差体值大小识别缝洞发育区。
下面,将结合附图,并以GST工区为例详细描述地震层序控制下的叠前AVO检测缝洞技术。
如图1所示,在本示例性实施例中,地震层序控制下的叠前AVO检测缝洞方法可通过以下步骤实现:
1、通过精细地震资料解释,根据削蚀、顶超、下超和上超等反射终止关系,识别地震层序,建立研究区主要目的层的地震层序格架。
2、利用测井资料建立关键井AVO正演模型,分析缝洞储层的AVO响应特征;提取共反射点道集的截距、梯度属性。
3、对截距、梯度属性沿地震层序框架进行方差体提取,得到目的层段的AVO属性方差体。
4、根据关键井缝洞储层标定结果,划定截距、梯度识别缝洞的方差门槛值,选择单个属性或两个属性交会识别门槛值以上缝洞发育体系,建立缝洞与AVO属性方差体关系。
5、在平面上和剖面上预测缝洞发育情况,绘制目的层段缝洞发育带预测平面图。
图2示出了利用缝洞发育的井的电测资料制作的AVO正演模型和井旁共反射点道集的响应图。如图2所示,在缝洞发育段顶部模型与近旁地震道截距均正值,都出现了振幅随偏移距增大的AVO响应,具有Ⅰ类AVO异常现象,说明AVO响应较好地反映了缝洞储层发育情况。其中,AVO正演模型是指振幅随偏移距变化模型,井旁CRP道集指共反射点道集。
如图3所示,根据层位标定结果,在三维地震资料上精细解释地震反射层位,不仅解释出地质界面(例如,图3中A所示点线界面),还依据超覆和削蚀现象解释识别内部层序界面(例如,图3中B所示界面)。
图4是根据叠前共反射点道集提取的AVO截距、梯度属性剖面,缝洞发育井高截距、高梯度明显,值域横向变化大,延伸范围小,与模型正演结果对应。
图5是目的层提取的AVO梯度属性方差体连井剖面,图中方差越大表示变化率越高,井轨迹上的方块(例如,井Gaoshi8(即GS8)、Gaoshi8(即,GS9)、Gaoshi10(即,GS10)等上的方块)表示储层发育段。剖面上大于0.65以上的高方差异常值与井缝洞储层发育段对应,方差体纵横向变化大的,表明AVO属性体变化较大,缝洞体系非均质性强;方差值低于0.4的区域基本呈现较高异常较发育,横向有层状分布趋势,但相互间连通性差;方差值高于0.6的区域基本呈现高值异常横向不连续,纵向上目的层上端方差值高于0.6的区域基本呈现高值异常较下段发育。其中,在目的层测试获90万方/日以上的高产工业井GS6、GS7、GS9井附近的诸多区域呈现高值异常发育。
此外,通过在研究区应用本发明的地震层序控制下的叠前AVO检测缝洞方法,预测了震旦系灯四段、灯二段缝洞发育有利区,预测结果表明缝洞储层发育的高产井和工业气井多位于AVO属性体方差高值区,位于属性体方差低值区的井多为储层不发育井的微气井,灯四段和灯二段的缝洞检测符合分别为87%、84%,具有较高的预测精度,可对缝洞发育有利区快速圈定,具有良好的应用前景。
综上所述,本发明提出地震层序控制下的叠前AVO检测缝洞技术,利用AVO模型正演分析含油气缝洞储层的AVO响应,提取AVO的梯度、截距属性体,在精细解释基础上,建立地震层序框架,对AVO属性进行方差分析处理,得到目的层段梯度和截距的方差体,利用正演模型缝洞体标定,结合测井解释成果、测试成果,划定缝洞体方差体识别门槛,提取门槛值以上方差体平均值,在平面上快速划分缝洞储层发育区,或直接在数据体上进行缝洞体异常地质体追踪。
本发明充分利用了AVO响应和地震等时地层格架的横向分辨能力,不仅可以快速识别缝洞发育有利区,而且由于油气充注的缝洞体会进一步增大其AVO响应,也可以间接识别油气富集区,与生产结合紧密,具有快速、有效检测特点。此外,本发明的特点还包括:理论体系完整可靠,流程清楚,使缝洞检测过程更加直接;能够适用三维地震资料,具有一定钻井数目区块(最好保证有缝洞发育有差异的模式井2口或以上)的缝洞预测;易于实施,具有极大的推广和实用价值。
尽管上面已经结合示例性实施例及附图描述了本发明,但是本领域普通技术人员应该清楚,在不脱离权利要求的精神和范围的情况下,可以对上述实施例进行各种修改。

Claims (6)

1.一种地震层序控制下的叠前振幅随偏移距变化检测缝洞方法,其特征在于,所述检测缝洞方法包括以下步骤:
通过地震资料解释,建立研究区主要目的层的地震层序格架;
利用测井资料建立关键井AVO正演模型,分析缝洞储层的AVO响应特征,提取共反射点道集的截距、梯度属性;
对截距、梯度属性沿地震层序框架进行方差体提取,得到目的层段的AVO属性方差体;
根据关键井缝洞储层标定结果,划定截距、梯度识别缝洞的方差门槛值,选择单个属性或两个属性交会识别门槛值以上缝洞发育体系,建立缝洞与AVO属性方差体关系。
2.根据权利要求1所述的地震层序控制下的叠前振幅随偏移距变化检测缝洞方法,其特征在于,所述地震资料解释的步骤还包括根据反射终止关系,识别地震层序。
3.根据权利要求2所述的地震层序控制下的叠前振幅随偏移距变化检测缝洞方法,其特征在于,所述反射终止关系包括削蚀、顶超、下超和/或上超。
4.根据权利要求1所述的地震层序控制下的叠前振幅随偏移距变化检测缝洞方法,其特征在于,所述测井资料包括纵波测井数据和/或横波测井数据。
5.根据权利要求1所述的地震层序控制下的叠前振幅随偏移距变化检测缝洞方法,其特征在于,所述检测缝洞方法包括在所述建立缝洞与AVO属性方差体关系的步骤之后,在平面上和剖面上预测缝洞发育情况,绘制目的层段缝洞发育带预测平面图。
6.根据权利要求1所述的地震层序控制下的叠前振幅随偏移距变化检测缝洞方法,其特征在于,所述建立研究区主要目的层的地震层序格架的步骤和所述建立关键井AVO正演模型的步骤无先后顺序之分。
CN201611254793.5A 2016-12-30 2016-12-30 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法 Active CN107045145B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611254793.5A CN107045145B (zh) 2016-12-30 2016-12-30 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611254793.5A CN107045145B (zh) 2016-12-30 2016-12-30 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法

Publications (2)

Publication Number Publication Date
CN107045145A CN107045145A (zh) 2017-08-15
CN107045145B true CN107045145B (zh) 2019-02-15

Family

ID=59543388

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611254793.5A Active CN107045145B (zh) 2016-12-30 2016-12-30 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法

Country Status (1)

Country Link
CN (1) CN107045145B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109991661B (zh) * 2019-04-08 2020-03-06 成都理工大学 油气检测方法及装置
CN110174701A (zh) * 2019-06-05 2019-08-27 中国煤炭地质总局地球物理勘探研究院 一种火成岩侵入范围的预测方法及装置
CN112198549A (zh) * 2019-07-08 2021-01-08 中国石油天然气集团有限公司 一种基于地震正演模板的叠前裂缝确定方法及系统
CN113325468B (zh) * 2020-02-28 2023-06-30 中国石油天然气集团有限公司 储层缝洞分布范围预测方法及装置
CN111665562B (zh) * 2020-07-20 2022-03-01 西南石油大学 一种高精度地震层序划分方法
CN112363221A (zh) * 2020-10-27 2021-02-12 中国石油天然气集团有限公司 Walkaway VSP测线的布设方法及装置
CN112462421A (zh) * 2020-10-30 2021-03-09 中国石油天然气集团有限公司 储层信息预测方法、装置、电子设备及存储介质
CN114966851B (zh) * 2022-05-13 2023-05-05 广州海洋地质调查局 一种储层预测方法、装置及存储介质

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5999486A (en) * 1998-07-23 1999-12-07 Colorado School Of Mines Method for fracture detection using multicomponent seismic data
CN101551466B (zh) * 2008-04-03 2011-04-20 中国石油天然气集团公司 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法
CN102109611B (zh) * 2009-12-23 2012-08-15 中国石油天然气集团公司 一种利用低频地震属性预测油藏优质储层的方法
CN104375178B (zh) * 2014-11-05 2017-10-17 中国石油天然气股份有限公司 一种碳酸盐岩缝洞储层预测方法和装置
CN104818982B (zh) * 2015-04-20 2017-12-05 中国石油天然气股份有限公司 一种缝洞型储层连通性定量分析方法及装置
CN104820239B (zh) * 2015-05-13 2018-07-06 中国石油大学(华东) 一种方位叠前地震属性解耦提取方法
CN105182424B (zh) * 2015-08-03 2017-09-01 中国石油天然气股份有限公司 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置
CN106019370B (zh) * 2016-05-06 2018-01-02 中国石油天然气股份有限公司 岩溶层断珠体油藏的井点确定方法及装置

Also Published As

Publication number Publication date
CN107045145A (zh) 2017-08-15

Similar Documents

Publication Publication Date Title
CN107045145B (zh) 地震层序控制下的叠前振幅随偏移距变化检测缝洞方法
CN104635269B (zh) 一种基于叠前方位角道集预测火成岩裂缝型储层的方法
WO2018010628A1 (zh) 一种基于大面积致密储层地震岩石物理反演方法
CN105653815B (zh) 一种基于岩石物理模型理论的定量解释油藏流体分布方法
CN104657523A (zh) 砂砾岩综合地质建模方法
CN104678434B (zh) 一种预测储层裂缝发育参数的方法
CN106951660A (zh) 一种海相碎屑岩水平井储层测井解释方法及装置
CN104155701B (zh) 一种利用叠前地震资料及井信息的多尺度裂缝预测方法
CN104636980B (zh) 针对河道砂油藏类型油气汇集条件的地球物理表征方法
US20140343859A1 (en) Time-lapse 4d scattering for imaging hydraulically induced fractures
CN106556869A (zh) 一种拉张环境下定量刻画断层走滑位移的方法
CN102053259A (zh) 一种油气藏识别方法
CN104950327A (zh) 确定地面微地震观测系统的检波器的位置的方法
CN109655905A (zh) 基于全方位角数据叠前avaz反演裂缝预测方法及系统
CN103869366B (zh) 一种确定裂隙裂缝走向的方法及装置
CN106033127B (zh) 基于横波速度变化率的地应力方位地震预测方法
CN106842299B (zh) 一种基于地震属性的裂缝定量化预测的方法
CN106468171B (zh) 多信息耦合的递进式砂砾岩岩相判别方法
CN104375166B (zh) 确定地震勘探激发深度的方法
US9244183B2 (en) Fracture characterization from refraction travel time data
CN106501851B (zh) 一种地震属性优化方法及装置
CN104820237B (zh) 计算地层的方差体的方法
CN114135269B (zh) 一种致密砂岩油层识别方法及装置
CN113093275B (zh) 曲流河及曲流河三角洲油田钻井成功率提高方法及装置
Li et al. Fault and fracture prediction of tight gas reservoir based on seismic likelihood attribute

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
CB02 Change of applicant information

Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant after: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

CB02 Change of applicant information
TA01 Transfer of patent application right

Effective date of registration: 20180402

Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Applicant after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

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

Effective date of registration: 20201113

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Patentee before: BGP Inc., China National Petroleum Corp.

TR01 Transfer of patent right