CN103837893A - 基于vsp资料的声波测井曲线校正方法 - Google Patents
基于vsp资料的声波测井曲线校正方法 Download PDFInfo
- Publication number
- CN103837893A CN103837893A CN201410085309.5A CN201410085309A CN103837893A CN 103837893 A CN103837893 A CN 103837893A CN 201410085309 A CN201410085309 A CN 201410085309A CN 103837893 A CN103837893 A CN 103837893A
- Authority
- CN
- China
- Prior art keywords
- acoustic logging
- vsp data
- time
- acoustic
- dark
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000012937 correction Methods 0.000 title claims abstract description 10
- 239000002131 composite material Substances 0.000 description 13
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000001915 proofreading effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 102000000584 Calmodulin Human genes 0.000 description 1
- 108010041952 Calmodulin Proteins 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于VSP资料的声波测井曲线校正方法，所述方法包括：将声波测井曲线按VSP资料标定的初至的时深关系校正到同一基准面上；对目的层段中的每一层利用VSP资料对声波测井曲线的时深关系进行校正，计算VSP资料的时深关系点与声波测井曲线的时深关系点的时间差；计算VSP资料和声波测井曲线每一层的时深关系点的时间差的平均值；将声波测井曲线整体拉伸所述时间差的平均值；调整声波速度值，当声波测井曲线的平均层速度与VSP资料的平均层速度一致时，得到校正后的声波测井曲线。本发明利用VSP资料准确的时-深关系，并结合测井分层和区域地质特征，对声波测井曲线进行校正，使校正后的声波合成记录与地面地震剖面波组特征更匹配。
Description
技术领域
本发明属于石油天然气地震勘探调查领域，主要应用于石油地震勘探的室内地震资料处理解释，更具体地讲，涉及一种基于VSP资料的声波测井曲线校正方法。
背景技术
地震勘探技术的应用从研究地质构造发展到研究地层岩性，随着寻找隐蔽性复杂油藏技术的发展，对岩性解释的精度要求越来越高。这时，地震地质层位的精细标定就显得尤为重要。
目前地球物理工程师常利用声波测井资料对地震地质层位进行标定，而声波测井资料由于受地层岩性、井孔、泥浆渗入地层和对泥岩的蚀变等诸多因素的影响，会使得部分井的井段或全井的声波测井资料有不同程度的失真；另一方面声波测井是以深度来采样的，为了与时间域的地面地震记录相对比，要先求得反射系数，做成合成记录，合成记录与地震记录是否匹配很大程度上与声波测井资料密切相关。分析合成记录与井旁地震记录，当确定了主要目的层与地震剖面相对应的相位，可以发现波组相位相对应，而合成记录的层时间厚度小于地震记录的层时间厚度。此外，人为读数的误差也会影响到声波测井资料的读数。
目前，利用声波测井资料对地震地质层位进行标定时，需将合成记录与地面地震记录进行匹配，现有的对声波测井资料进行校正的方法往往采用手工拉伸方法，其结果是薄层标定精度较低，且存在人为误操作风险。
发明内容
本发明的目的在于利用VSP资料的准确的时深关系对声波测井曲线进行校正，有效解决声波测井曲线的失真问题，并使得校正后的声波合成记录能够精细准确地标定地震地质层位。
本发明的一方面提供一种基于VSP资料的声波测井曲线校正方法，所述方法包括：将声波测井曲线按VSP资料标定的初至的时深关系校正到同一基准面上；从声波测井曲线中选取目的层段，对所述目的层段中的每一层执行以下步骤：利用VSP资料对声波测井曲线的时深关系进行校正，将实际地质的界面深度作为声波测井曲线时深关系的校正点，计算VSP资料的时深关系点与声波测井曲线的时深关系点的时间差；计算VSP资料和声波测井曲线每一层的时深关系点的时间差的平均值；将声波测井曲线整体拉伸所述时间差的平均值；根据VSP资料的平均层速度和区域地层平均速度，调整声波速度值，当声波测井曲线的平均层速度与VSP资料的平均层速度一致时，得到校正后的声波测井曲线。
此外，通过将VSP资料的每一层与声波测井曲线的每一层进行一一对应，来利用VSP资料对声波测井曲线的时深关系进行校正。
此外，利用时间漂移法将声波测井曲线整体拉伸所述时间差的平均值。
根据本发明的对声波测井曲线进行校正的方法，能够方便、快捷的利用VSP资料获得的时深关系对声波测井曲线进行校正，使校正后的声波合成记录与地面地震剖面波组特征更匹配。
附图说明
通过下面结合附图进行的详细描述，本发明的上述和其它目的、特点和优点将会变得更加清楚，其中：
图1是示出根据本发明的实施例的对声波测井曲线进行校正的方法的流程图。
具体实施方式
现在对本发明实施例进行详细的描述，其示例表示在附图中，其中，相同的标号始终表示相同部件。下面通过参照附图对实施例进行描述以解释本发明。
为了解决声波测井曲线存在失真，声波合成记录与地面地震记录匹配不好，不能准确标定地震地质层位的问题，本发明提出了一种基于VSP资料的声波测井曲线校正方法，获得的校正后的声波测井曲线能更好地反映地层速度特征，声波合成记录与地面地震剖面波组特征更匹配，提高了声波合成记录品质。
图1是示出根据本发明的实施例的对声波测井曲线进行校正的方法的流程图。
参照图1，在步骤101中，将声波测井曲线按VSP资料标定的初至的时深关系校正到同一基准面上。
具体地讲，由于声波测井曲线与VSP资料的深度均为测量深度，但是声波测井通常不是从井口开始测量，因此，声波测井曲线缺少井口到声波测井起测段的信息，所以可利用VSP资料标定的初至的时深关系对声波测井曲线进行校正，使声波测井曲线的起始点的时深关系与VSP资料标定的初至的时深关系统一到同一基准面。
在步骤102中，从声波测井曲线中选取目的层段，对所述目的层段中的每一层执行以下步骤：利用VSP资料对声波测井曲线的时深关系进行校正，将实际地质的界面深度作为声波测井曲线时深关系的校正点，计算VSP资料的时深关系点与声波测井曲线的时深关系点的时间差。
例如，利用VSP资料对声波测井曲线的时深关系进行校正的步骤可为：将VSP资料的每一层与声波测井曲线的每一层进行一一对应，完成对声波测井曲线与VSP资料标定地震地质层位。例如，选取的目的层段中可能包含多层，将VSP资料的第一层与声波测井曲线的第一层进行对应，VSP资料的第二层与声波测井曲线的第二层进行对应，以此类推，使声波测井曲线每一层与VSP资料的每一层逐层对应，以完成对声波测井曲线的时深关系的校正。
可选地，由于实际地质的地质分层比较可靠，因此利用实际地质的界面深度作为声波测井曲线的时深关系的校正点进行校正，由于声波测井曲线校正后的时深关系仍与实际地质的时深关系存在细微的偏差，因此需对声波测井曲线的时间进行进一步的校正。
具体地讲，利用VSP资料对声波测井曲线的时深关系进行校正时，在所述选取的目的层段中，VSP资料的第一层的时深关系点与声波测井曲线的第一层的时深关系点之间具有一个第一时间差，VSP资料的第二层的时深关系点与声波测井曲线的第二层的时深关系点之间具有一个第二时间差，以此类推。这里，VSP资料与声波测井曲线的每一层的时深关系点的时间差不同，有细微差异。这里，由于计算VSP资料的时深关系点与声波测井曲线的时深关系点的时间差的方法为本领域的公知常识，本发明对此部分的内容不再详细说明。
在步骤103中，计算VSP资料和声波测井曲线每一层的时深关系点的时间差的平均值。
例如，由于VSP资料与声波测井曲线的每一层的时间差不同，因此在对声波测井曲线进行拉伸时，如果每一层的时间差之间的差异较小，则在对声波测井曲线进行拉伸时，可计算出所述目的层段的VSP资料与声波测井曲线的每一层的时间差的平均值，将声波测井曲线整体拉伸所述时间差的平均值，以完成对声波测井曲线的时间的校正。如果其中某一层的时间差与其他层的时间差的差异较大，则需对步骤102进行检查（例如，检查地震地质层位标定的是否不准确）。
在步骤104中，将声波测井曲线整体拉伸所述时间差的平均值。
例如，可利用时间漂移法将声波测井曲线整体拉伸所述时间差的平均值。这里，由于时间漂移法为本领域的公知常识，本发明对此部分的内容不再详细描述。
在步骤105中，根据VSP资料的平均层速度和区域地层平均速度，调整声波速度值，当声波测井曲线的平均层速度与VSP资料的平均层速度一致时，得到校正后的声波测井曲线。
具体地讲，在声波测井曲线的时深关系改变的同时，声波速度值也会发生变化，因此，可利用VSP资料的平均层速度，并结合区域地质特征（例如，区域地层平均速度），调整声波速度值。由于可利用声波速度值计算得到声波测井曲线的平均层速度，因此当调整声波速度值使声波测井曲线的平均层速度与VSP资料的平均层速度一致时，校正结束，得到校正后的声波测井曲线。
这里，由于计算VSP资料的平均层速度以及利用声波速度值计算声波测井曲线的平均层速度的方法为本领域的公知常识，本发明对此部分的内容不再详细说明。
本发明利用VSP资料准确的时-深关系，结合测井分层和区域地质特征，很好地解决了对声波测井曲线的校正问题。
本发明提供一种基于VSP资料的声波测井曲线校正方法，该方法能够方便、快捷的利用VSP获得的时-深关系对声波测井曲线进行校正，使校正后的声波合成记录与地面地震地层剖面波组特征更匹配。解决了目前声波合成记录与地面地震记录标定地震地质层位效果不好，人工拉伸声波测井曲线，声波合成记录不能精细准确地标定地震地质层位的问题，在地震资料处理和地震资料解释方面具有广阔的应用前景。
尽管已经参照其示例性实施例具体显示和描述了本发明，但是本领域的技术人员应该理解，在不脱离权利要求所限定的本发明的精神和范围的情况下，可以对其进行形式和细节上的各种改变。
Claims (3)
1.一种基于VSP资料的声波测井曲线校正方法，所述方法包括：
（1）将声波测井曲线按VSP资料标定的初至的时深关系校正到同一基准面上；
（2）从声波测井曲线中选取目的层段，对所述目的层段中的每一层执行以下步骤：利用VSP资料对声波测井曲线的时深关系进行校正，将实际地质的界面深度作为声波测井曲线时深关系的校正点，计算VSP资料的时深关系点与声波测井曲线的时深关系点的时间差；
（3）计算VSP资料和声波测井曲线每一层的时深关系点的时间差的平均值；
（4）将声波测井曲线整体拉伸所述时间差的平均值；
（5）根据VSP资料的平均层速度和区域地层平均速度，调整声波速度值，当声波测井曲线的平均层速度与VSP资料的平均层速度一致时，得到校正后的声波测井曲线。
2.根据权利要求1所述的方法，其中，通过将VSP资料的每一层与声波测井曲线的每一层进行一一对应，来利用VSP资料对声波测井曲线的时深关系进行校正。
3.根据权利要求1所述的方法，其中，利用时间漂移法将声波测井曲线整体拉伸所述时间差的平均值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410085309.5A CN103837893B (zh) | 2014-03-10 | 2014-03-10 | 基于vsp资料的声波测井曲线校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410085309.5A CN103837893B (zh) | 2014-03-10 | 2014-03-10 | 基于vsp资料的声波测井曲线校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103837893A true CN103837893A (zh) | 2014-06-04 |
| CN103837893B CN103837893B (zh) | 2017-06-30 |
Family
ID=50801592
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410085309.5A Active CN103837893B (zh) | 2014-03-10 | 2014-03-10 | 基于vsp资料的声波测井曲线校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103837893B (zh) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104133249A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
| CN104156624A (zh) * | 2014-08-27 | 2014-11-19 | 中国电子科技集团公司第二十二研究所 | 测井设备的时深数据的处理方法和装置 |
| CN104656134A (zh) * | 2015-03-11 | 2015-05-27 | 中国电建集团北京勘测设计研究院有限公司 | 一种地震波初至读时与时距曲线动态关联的方法和系统 |
| CN104678440A (zh) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | 井约束二维地震变速度场非线性误差校正方法 |
| CN105445802A (zh) * | 2014-09-04 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于时频谱约束的高精度时深校正方法 |
| CN106054254A (zh) * | 2016-08-01 | 2016-10-26 | 中国石油天然气集团公司 | 一种测井曲线校正方法及装置 |
| CN106154329A (zh) * | 2015-04-16 | 2016-11-23 | 中国石油化工股份有限公司 | 利用vsp速度校正声波均方根速度的方法及其系统 |
| CN106324680A (zh) * | 2016-08-18 | 2017-01-11 | 中国石油天然气集团公司 | 地层破裂压力预测方法 |
| CN106569280A (zh) * | 2016-10-28 | 2017-04-19 | 中国石油化工股份有限公司 | 一种提高断阶构造地震解释精度的标定方法 |
| CN112649876A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 建立地震偏移速度模型的方法和装置 |
| CN112987098A (zh) * | 2021-04-02 | 2021-06-18 | 中油奥博(成都)科技有限公司 | 基于vsp的微地震速度模型优化方法 |
| CN113514876A (zh) * | 2021-06-18 | 2021-10-19 | 中国石油化工股份有限公司 | 一种基于地层速度校正的设计井深度计算方法 |
| WO2023184055A1 (en) * | 2022-03-28 | 2023-10-05 | Saudi Arabian Oil Company | A method and system for determining first breaks of sonic waveforms |
| CN117077068A (zh) * | 2023-10-18 | 2023-11-17 | 中国科学院地质与地球物理研究所 | 深地油气精准导航随钻声波测井数据实时标定方法与系统 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107765309B (zh) * | 2017-10-09 | 2019-03-29 | 中国石油天然气集团有限公司 | 一种基于垂直地震数据的钻前深度预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129084A (zh) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
| CN103018778A (zh) * | 2012-12-12 | 2013-04-03 | 中国石油天然气股份有限公司 | 校正声波测井曲线进行储层预测的方法及设备 |
-
2014
- 2014-03-10 CN CN201410085309.5A patent/CN103837893B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129084A (zh) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
| CN103018778A (zh) * | 2012-12-12 | 2013-04-03 | 中国石油天然气股份有限公司 | 校正声波测井曲线进行储层预测的方法及设备 |
Non-Patent Citations (7)
| Title |
|---|
| 云美厚 等: "地震与声波测井匹配校正方法研究", 《石油地球物理勘探》, vol. 37, no. 6, 31 December 2002 (2002-12-31), pages 601 - 605 * |
| 侯伯刚 等: "地震反演中测井数据的预处理", 《物探与化探》, vol. 33, no. 3, 30 June 2009 (2009-06-30), pages 331 - 336 * |
| 公维民: "声波测井曲线深度校正方法", 《石油仪器》 * |
| 张卫红: "VSP和声波测井资料在岩性解释中的联合应用", 《石油物探》, vol. 41, no. 2, 30 June 2002 (2002-06-30) * |
| 曾照荣 等: "声波测井与VSP测井平均速度的对比分析", 《新疆石油地质》 * |
| 朱广生 等: "塔里木盆地巨型速度场的建立", 《新疆石油地质》, vol. 23, no. 6, 31 December 2002 (2002-12-31), pages 507 - 510 * |
| 高志凌 等: "用VSP时深关系对声波测井资料进行校正", 《中国地球物理年会2007》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104133249B (zh) * | 2014-07-31 | 2016-08-31 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
| CN104133249A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
| CN104156624B (zh) * | 2014-08-27 | 2017-08-01 | 中国电子科技集团公司第二十二研究所 | 测井设备的时深数据的处理方法和装置 |
| CN104156624A (zh) * | 2014-08-27 | 2014-11-19 | 中国电子科技集团公司第二十二研究所 | 测井设备的时深数据的处理方法和装置 |
| CN105445802A (zh) * | 2014-09-04 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于时频谱约束的高精度时深校正方法 |
| CN104678440A (zh) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | 井约束二维地震变速度场非线性误差校正方法 |
| CN104656134A (zh) * | 2015-03-11 | 2015-05-27 | 中国电建集团北京勘测设计研究院有限公司 | 一种地震波初至读时与时距曲线动态关联的方法和系统 |
| CN106154329A (zh) * | 2015-04-16 | 2016-11-23 | 中国石油化工股份有限公司 | 利用vsp速度校正声波均方根速度的方法及其系统 |
| CN106054254A (zh) * | 2016-08-01 | 2016-10-26 | 中国石油天然气集团公司 | 一种测井曲线校正方法及装置 |
| CN106324680A (zh) * | 2016-08-18 | 2017-01-11 | 中国石油天然气集团公司 | 地层破裂压力预测方法 |
| CN106569280A (zh) * | 2016-10-28 | 2017-04-19 | 中国石油化工股份有限公司 | 一种提高断阶构造地震解释精度的标定方法 |
| CN112649876A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 建立地震偏移速度模型的方法和装置 |
| CN112987098A (zh) * | 2021-04-02 | 2021-06-18 | 中油奥博(成都)科技有限公司 | 基于vsp的微地震速度模型优化方法 |
| CN113514876A (zh) * | 2021-06-18 | 2021-10-19 | 中国石油化工股份有限公司 | 一种基于地层速度校正的设计井深度计算方法 |
| CN113514876B (zh) * | 2021-06-18 | 2024-01-23 | 中国石油化工股份有限公司 | 一种基于地层速度校正的设计井深度计算方法 |
| WO2023184055A1 (en) * | 2022-03-28 | 2023-10-05 | Saudi Arabian Oil Company | A method and system for determining first breaks of sonic waveforms |
| CN117077068A (zh) * | 2023-10-18 | 2023-11-17 | 中国科学院地质与地球物理研究所 | 深地油气精准导航随钻声波测井数据实时标定方法与系统 |
| CN117077068B (zh) * | 2023-10-18 | 2024-03-08 | 中国科学院地质与地球物理研究所 | 深地油气精准导航随钻声波测井数据实时标定方法与系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103837893B (zh) | 2017-06-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103837893A (zh) | 基于vsp资料的声波测井曲线校正方法 | |
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| Treadgold et al. | Eagle Ford shale prospecting with 3D seismic data within a tectonic and depositional system framework | |
| CN108680951A (zh) | 一种基于地震信息判断煤层气富集沉积控制作用的方法 | |
| CN106842313A (zh) | 基于方位叠前地震数据的各向异性参数反演方法 | |
| CN108663713B (zh) | 一种建立深度域构造模型的方法 | |
| CN103592680B (zh) | 一种基于正反演的测井数据和深度域地震剖面合成方法 | |
| CN105093306A (zh) | 一种地球物理勘探中储层自动解释与厚度求取方法 | |
| CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
| CN106526671B (zh) | 一种碳酸盐岩储层地震波响应特征分析方法 | |
| CN111208564A (zh) | 一种深度域层位标定方法及装置 | |
| Jenkinson et al. | Joint PP-PS angle-stack analysis and AVA inversion in Grane Field, offshore Norway | |
| Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
| Altenbernd et al. | Geophysical evidence for the extent of crustal types and the type of margin along a profile in the northeastern Baffin Bay | |
| CN107238859B (zh) | 计算虚拟井合成地震记录的方法及系统 | |
| CN116500679A (zh) | 一种井震标定方法、装置和相关设备 | |
| Jackson et al. | Seismic refraction data in the Gulf of Saint Lawrence: implications for the lower-crustal blocks | |
| Luo et al. | Well-to-seismic calibration in depth domain using dynamic depth warping techniques | |
| CN114578452B (zh) | 一种定量计算地下反射系数振幅占比的方法 | |
| Luo et al. | A well-to-seismic calibration method for seismic data in depth domain | |
| Liu et al. | Estimation of Thomsen’s anisotropic parameters from walkaway VSP and applications | |
| Stein et al. | A successful statics methodology for land data | |
| O’Brien et al. | Multicomponent VSP imaging of tight-gas sands | |
| Wu et al. | Characteristic reflection wavefield based wave-equation traveltime inversion | |
| Constance et al. | Integration of Multicomponet Surface Seismic, Multicomponent VSP and Microseismic in Reservoir characterization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180208 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp. 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 Patentee before: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200918 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: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee before: BGP Inc., China National Petroleum Corp. |