CN102353990B - Method for manufacturing seismic slice with overthrust fault - Google Patents
Method for manufacturing seismic slice with overthrust fault Download PDFInfo
- Publication number
- CN102353990B CN102353990B CN 201110145744 CN201110145744A CN102353990B CN 102353990 B CN102353990 B CN 102353990B CN 201110145744 CN201110145744 CN 201110145744 CN 201110145744 A CN201110145744 A CN 201110145744A CN 102353990 B CN102353990 B CN 102353990B
- Authority
- CN
- China
- Prior art keywords
- fault
- layer position
- tomography
- horizon
- section
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 16
- 239000010410 layer Substances 0.000 claims description 72
- 238000003325 tomography Methods 0.000 claims description 40
- 239000011229 interlayer Substances 0.000 claims description 10
- 238000007789 sealing Methods 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 9
- 238000011161 development Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention provides a method for manufacturing a seismic slice with an overthrust fault, which comprises the following steps of: carrying out explanation on a geologic horizon (with an overthrust fault) and a fault; connecting the geologic horizon with the fault so as to form a new horizon according to the explanation results of the geologic horizon and the fault; carrying out up-and-down horizon connection on the obtained new horizon at the position of the fault so as to form a closed fault block; and according to the formed fault block and the connected horizon, executing an operation ofslicing.
Description
Technical field
The invention belongs to the oil-gas exploration and development field, more particularly, relate to a kind of method that contains the earthquake microsection manufacture of overthrust fault.
Background technology
Earthquake section analytical technology is important analytical technology in the 3-d seismic exploration, has important effect for the identification in reservoir space distribution zone.Earthquake section is to reconfigure in the horizontal with original earthquake data or through the geological data of certain special processing, forms new data, and new data are representation of data Changing Pattern in the horizontal preferably.Generally, earthquake section is divided into three major types, and first kind is section when waiting, and namely along the identical seismic reflection time geological data is extracted, and forms new data; Second kind is horizon slice, comprises along the geologic horizon of a certain seismic interpretation geological data being extracted with geologic horizon along this seismic interpretation identical time of floating up or down carrying out data pick-up, forms new data; The third is the interlayer stratum section that developed recently gets up, and namely carries out data pick-up along time thickness or depth thickness between the geologic horizon hierarchical of certain two seismic interpretation, forms new data.Second kind and the third section analytical approach all can relate to the tomography of seismic interpretation, traditional microsection manufacture method is calculated under trap-down or unfaulted model, when the stratum contains trap-up, separately twice calculating of upper lower burrs, cast out trap-up information, make and can't finish the slice of data that contains trap-up.Yet, development along with the seismic prospecting exploitation, carried out the transition to the complex structure exploration stage now, containing the contrary tectonic structure of covering fracture has become geological phenomenon very general in the exploratory development, the needs that conventional microsection manufacture and analytical approach can't satisfy research and production.
Summary of the invention
To contain overthrust fault upper lower burrs earthquake section and can't once calculate and finish in order to solve, the problem of footwall information dropout, carry out at the tectonic structure that contains overthrust fault and carry out the earthquake microsection manufacture, the reservoir prediction work that section is analyzed, the invention provides a kind of method that contains the earthquake microsection manufacture of overthrust fault, this method can be made the earthquake section simultaneously to the geologic model upper lower burrs that contains trap-up, comprehensive utilization upper lower burrs data, increased the measurable area of earthquake reservoir prediction, saved seismic interpretation worker's working time, and made and predict the outcome more comprehensively, more reasonable.
In order to address the above problem, a kind of method that contains the earthquake microsection manufacture of overthrust fault is provided, this method can comprise: geologic horizon and the tomography that contains overthrust fault made an explanation; The geologic horizon of the explanation results of geologic horizon and tomography is connected with tomography forms new layer position; The fault block that the levels position is connected to form sealing is carried out at position of fault in the new layer position that forms; Carry out section according to the fault block that forms with the layer position that is connected.
Can adopt to contain that the contrary seismic interpretation method of covering fracture is finished geologic horizon and tomography makes an explanation, contain the overthrust fault zone, same geologic horizon was repeated to pick up and show in the different time.
Be connected with tomography at the geologic horizon with the explanation results of geologic horizon and tomography and form in the step of new layer position, can be with the explanation results of geologic horizon and tomography according to a certain end of main profile direction from survey line, each section of the layer position of explaining coupled together with the tomography that intersects with its, forms one new layer.
Carrying out the levels position in the new layer position that will form at position of fault is connected to form in the step of fault block of sealing, when tomography penetrates adjacent when two-layer up and down, be that the basis connects with these two adjacent layer position with the line of rent, it is the fault block on border one by one with the tomography that the layer position that connects will form.
The method of carrying out section comprises: section when waiting, horizon slice and interlayer section.
In horizon slice, adding, subtract a certain fixed value along the main profile direction with newly-generated layer position or this layer position according to the layer position of sealing is time or degree of depth reference mark, extract corresponding geological data, form the data volume of a new layer position type, this layer bit data body shows that with planar fashion the section that is this geological data shows.
In the interlayer section, between certain two layer position, calculate the Along ent of two interlayers according to the number of slices that will make, in containing the trap-up zone, after calculating Along ent respectively and calculate all Along ents according to upper lower burrs, extract the data volume that corresponding geological data forms one by one new layer position type according to Along ent, these layers bit data body shows that with planar fashion the section that is this geological data shows.
Description of drawings
By the description of carrying out below in conjunction with accompanying drawing, above-mentioned and other purposes of the present invention and characteristics will become apparent, wherein:
Fig. 1 illustrates according to the geologic horizon of overthrust fault and fault interpretation result's the diagram of containing of the present invention;
Fig. 2 is the diagram that connects according to layer of the present invention position;
Fig. 3 is the diagram according to tomography piecemeal of the present invention;
Fig. 4 is the diagram of the example of the tomography piecemeal of making according to the present invention;
Fig. 5 is the diagram of the stratum section cutting position of the acquisition according to the present invention;
Fig. 6 and Fig. 7 illustrate the stratum slice plane figure of the upper lower burrs that obtains according to the present invention respectively;
Fig. 8 is the process flow diagram that illustrates according to the method for the earthquake microsection manufacture that contains overthrust fault of exemplary embodiment of the present invention.
Embodiment
The present invention is described below with reference to accompanying drawings.
In order more clearly to describe technical scheme of the present invention, shown in the accompanying drawings by adopt software "
Seismic interpretation slip-stick artist " carry out situation of the present invention, thus exemplary embodiment of the present invention is clearly described.
In accompanying drawing of the present invention, the DKH work area is the name in the earthquake work area of adopting among the figure, and 452 surveys line are survey line numberings in this work area, and CDP is a kind of numbering of geological data, and the cdp table in the literary composition shows the numbering of the data point on 452 surveys line.
The invention provides a kind of method that contains the earthquake microsection manufacture of overthrust fault.This method can comprise the steps: the geologic horizon that contains overthrust fault and tomography are made an explanation; The geologic horizon of the explanation results of geologic horizon and tomography is connected with tomography forms new layer position; The new layer position that forms is connected to form the fault block of sealing in row levels position, position of fault limit; Carry out section according to the fault block that forms with the layer position that is connected.
Fig. 8 is the process flow diagram that illustrates according to the method for the earthquake microsection manufacture that contains overthrust fault of exemplary embodiment of the present invention.
With reference to Fig. 8, at step S1, geologic horizon and the tomography that contains overthrust fault made an explanation.According to the present invention, adopt to contain against the seismic interpretation method of covering fracture and finish geologic horizon explanation and fault interpretation.Namely containing the overthrust fault zone, this interpretation procedure repeats to pick up and show in the different time to same geologic horizon., explain that good layer position and tomography are continuous in the horizontal, jump phenomena can not appear.
Geologic horizon is explained and to be referred to will demarcate the geochron horizon calibration to geological data by well-log information, on geological data, laterally follow the trail of according to Seismic reflection character then, complete layer position identification on whole geological data, this process is called geologic horizon and explains.
Fault interpretation refers to the Seismic reflection character by geological data, determines the position of each seismic fault, spatially tomography is made up then, finishes plane and the space spread of tomography and sketches, and this process is called fault interpretation.
Fig. 1 illustrates according to the geologic horizon of overthrust fault and fault interpretation result's the diagram of containing of the present invention.Fig. 1 finishes geologic horizon with the three-dimensional work area of DKH and explains and fault interpretation that the result who obtains such as Fig. 1 are a survey line explanation results wherein.Wherein transverse axis is CDP number, and the longitudinal axis is that (unit: ms), T1f4, T1f3ct, T1p11, T1p21 are layer positions to time shaft among Fig. 1, and fault1 and fault2 are tomographies.
Subsequently, at step S2, the geologic horizon among the step S1 and fault interpretation result's geologic horizon is connected with tomography forms new layer position.More particularly, with the explanation results among the step S1 according to a certain end of main profile direction from survey line, each section of the layer position of explaining coupled together with the tomography that intersects with its, forms a new layer position, new layer position will replace the geologic horizon of explanation originally to finish follow-up work.Wherein, survey line is a unit in charge of construction in the earthquake data acquisition work progress, constructs for two dimension, be exactly the data that once construction is gathered, show as a two-dimentional line, be called survey line, the data of a survey line are exactly a two-dimension earthquake data section, in two-dimentional work progress, in order to implement underground structure, the survey line major part is to lay along the subsurface geological structure direction, is called " main profile ", lay some surveys line in addition main profile is connected, be called " cross-track ".For three-dimensional, main profile all is the same uniform grid survey line with cross-track in principle, and just name is different.
Fig. 2 is according to layer of the present invention position connection layout.Fig. 2 connects with complete layer position, the three-dimensional work area of DKH, tomography, has 4 layer positions in Fig. 2.
At step S3, the levels position is carried out at position of fault in the new layer position that forms connect, form the fault block of sealing.When tomography penetrates adjacently when two-layer up and down, be that the basis connects with these two adjacent layer position with the line of rent, it is the fault block on border with the tomography that the layer position that connects will form at least one.
Fig. 3 is according to tomography block diagram of the present invention.Fig. 4 is the example illustration of the tomography piecemeal of making according to the present invention, wherein, different tomography piecemeals is indicated different colors.
At step S4, carry out section according to the fault block that forms with the layer position that is connected.Fig. 5 is the stratum section split position line with the stratum trisection between two layer positions in the three-dimensional work area of DKH according to exemplary embodiment of the present invention, as the slice position among Fig. 5 with the stratum trisection between two layer positions.
Calculating Deng time section needs the control of layer position, and the microsection manufacture method is fixed when waiting, and is not discussed at this.The making of horizon slice is that the layer position according to sealing is time (or degree of depth) reference mark along the main profile direction with newly-generated layer position (or this layer position adds, subtracts a certain fixed value), extract corresponding geological data, form the data volume of a new layer position type, this layer bit data body shows that with planar fashion the section that is this geological data shows.Each trap-up zone it should be noted that the section that this method is made is repetition containing the trap-up zone, so can be divided into the upper lower burrs mode and show.Method for making and the horizon slice of the section of interlayer stratum are similar, different is, the interlayer section is between certain two layer position, according to the number of slices that will make, calculate the Along ent of two interlayers, containing the trap-up zone, need calculate Along ent respectively according to upper lower burrs, after calculating all Along ents, extract corresponding geological data according to Along ent, form the data volume of at least one new layer position type, these layers bit data body shows that with planar fashion the section that is this geological data shows, equally, the section that this method is made is repetition containing the trap-up zone, and each trap-up zone can be divided into the upper lower burrs mode and show.Here need to prove, layer bit data body is shown as an irregular geometric curved surfaces at the spread in stereoscopic three-dimensional space in the space, 3-D display can see that the tectonic relief of layer position curved surface changes and the locus, can also see layer expression of bit attribute value; And " layer bit data a body show with planar fashion " refers to the projection in one plane of layer bit data, and what see is a panel data, can only see that layer bit attribute value and plane relative position change, and can't see the locus and structure changes.
Fig. 6 and Fig. 7 illustrate the stratum slice plane figure of the upper lower burrs that obtains according to the present invention respectively, specifically, Fig. 6 and Fig. 7 are dish and the lower wall demonstration results of the slice of data body that extracts at whole 3-D data volume along cut-off rule wherein, wherein, Fig. 6 is that last dish shows the result, and Fig. 7 shows the result for lower wall.As shown in Figure 6 and Figure 7, the scope of irregular tomography polygon delineation wherein is for containing the trap-up zone, and these two figure are different in the trap-up zone, and other zones are identical.
By using according to the method that contains the earthquake microsection manufacture of overthrust fault of the present invention, the stratum slice efficiency that contains overthrust fault significantly improves, and is reliable in the fault region result calculated, and the data of footwall are credible, can be used for carrying out reservoir prediction, therefore measurable area increases.Reservoir border performance with faults control is obvious.The present invention can solve the stratum section that contains overthrust fault and analyze reservoir integrated forecasting problem, and the reservoir prediction in the corresponding oil-gas exploration and development, oil and gas detection etc. have direct effect.
Although specifically shown with reference to exemplary embodiment of the present invention and described the present invention, but it should be appreciated by those skilled in the art, under the situation that does not break away from the spirit and scope of the present invention that are defined by the claims, can carry out various changes on form and the details to it.
Claims (5)
1. method that contains the earthquake microsection manufacture of overthrust fault comprises:
Geologic horizon and the tomography that contains overthrust fault made an explanation;
The geologic horizon of the explanation results of geologic horizon and tomography is connected with tomography forms new layer position;
The fault block that the levels position is connected to form sealing is carried out at position of fault in the new layer position that forms;
Carry out section according to the fault block that forms with the layer position that is connected,
Wherein, be connected with tomography at the geologic horizon with the explanation results of geologic horizon and tomography and form in the step of new layer position, with the explanation results of geologic horizon and tomography according to a certain end of main profile direction from survey line, with each section of the layer position of explaining with couple together with its crossing tomography, form a new layer position
Wherein, carrying out the levels position in the new layer position that will form at position of fault is connected to form in the step of fault block of sealing, when tomography penetrates adjacent when two-layer up and down, be that the basis connects with these two adjacent layer position with the line of rent, it is the fault block on border with the tomography that the layer position that connects forms at least one.
2. the method for claim 1, wherein adopt the seismic interpretation method that contains against covering fracture that geologic horizon and tomography are made an explanation, containing the overthrust fault zone, same geologic horizon was repeated to pick up and show in the different time.
3. the method for the method of claim 1, wherein carrying out section comprises: section when waiting, horizon slice and interlayer section.
4. method as claimed in claim 3, wherein, in horizon slice, adding, subtract a certain fixed value along the main profile direction with newly-generated layer position or this layer position according to the layer position of sealing is time or degree of depth reference mark, extract corresponding geological data, form the data volume of a new layer position type, this layer bit data body shows that with planar fashion the section that is this geological data shows.
5. method as claimed in claim 3, wherein, in the interlayer section, between certain two layer position, calculate the Along ent of two interlayers according to the number of slices that will make, in containing the trap-up zone, after calculating Along ent respectively and calculate all Along ents according to upper lower burrs, extract the data volume that corresponding geological data forms at least one new layer position type according to Along ent, these layers bit data body shows that with planar fashion the section that is this geological data shows.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110145744 CN102353990B (en) | 2011-06-01 | 2011-06-01 | Method for manufacturing seismic slice with overthrust fault |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110145744 CN102353990B (en) | 2011-06-01 | 2011-06-01 | Method for manufacturing seismic slice with overthrust fault |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102353990A CN102353990A (en) | 2012-02-15 |
| CN102353990B true CN102353990B (en) | 2013-08-14 |
Family
ID=45577582
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201110145744 Active CN102353990B (en) | 2011-06-01 | 2011-06-01 | Method for manufacturing seismic slice with overthrust fault |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102353990B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103941288B (en) * | 2014-04-16 | 2016-08-17 | 中国海洋石油总公司 | A kind of objective geologic body means of interpretation based on horizontal navigation |
| US9720131B2 (en) * | 2014-06-05 | 2017-08-01 | Chevron U.S.A. Inc. | System and method of building complex earth models |
| CN105301637B (en) * | 2015-10-10 | 2017-09-15 | 中国石油天然气集团公司 | The construction method and device of two-dimentional geological reservoir section |
| CN107703542B (en) * | 2017-08-31 | 2020-04-10 | 中国石油天然气集团公司 | Method and device for determining seismic horizon of reverse fault area |
| CN110308444B (en) * | 2019-08-08 | 2021-03-09 | 中国矿业大学(北京) | Road horizon intelligent identification and interference source elimination method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6018499A (en) * | 1997-11-04 | 2000-01-25 | 3Dgeo Development, Inc. | Three-dimensional seismic imaging of complex velocity structures |
| CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2010257595B2 (en) * | 2009-06-09 | 2013-12-05 | Shell Internationale Research Maatschappij B.V. | Method for analysing seismic data with attribute functions |
-
2011
- 2011-06-01 CN CN 201110145744 patent/CN102353990B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6018499A (en) * | 1997-11-04 | 2000-01-25 | 3Dgeo Development, Inc. | Three-dimensional seismic imaging of complex velocity structures |
| CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
Non-Patent Citations (4)
| Title |
|---|
| 三维地震勘探在复杂逆掩断层地区的应用及效果;吴博等;《西部探矿工程》;20050731;157-159页 * |
| 付晓飞等.逆掩断层垂向封闭性定量模拟实验及评价方法.《地质科学》.2004,(第02期),223-233页. |
| 吴博等.三维地震勘探在复杂逆掩断层地区的应用及效果.《西部探矿工程》.2005,157-159页. |
| 逆掩断层垂向封闭性定量模拟实验及评价方法;付晓飞等;《地质科学》;20040430(第02期);223-233页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102353990A (en) | 2012-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106934858B (en) | Three-dimensional geological modeling method and system for scale region of mining area | |
| Deng et al. | Two distinct strike-slip fault networks in the Shunbei area and its surroundings, Tarim Basin: Hydrocarbon accumulation, distribution, and controlling factors | |
| US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
| CN107808413A (en) | A kind of three-dimensional geological modeling method based on GOCAD | |
| CN102353990B (en) | Method for manufacturing seismic slice with overthrust fault | |
| CN103700141B (en) | Multistage dimension geologic modeling analysis method of coal seam floor water gushing passageway | |
| CN104375179A (en) | Method for looking for oil and gas reservoir based on TRAP-3D software | |
| US11953637B2 (en) | Methods and systems of generating simulation grids for different structural interpretations | |
| US11409024B2 (en) | Methods and systems for generating simulation grids via zone by zone mapping from design space | |
| CN102221707B (en) | Seismic attribute extracting method for geologic structure containing over-thrust fault | |
| CN105652341A (en) | Quick quantitative ancient landform recovering method with old structural effect considered | |
| KR101131517B1 (en) | Processing system and method for subsurface information in oil sands reservoir using statistical analysis of well logging data | |
| Al Rassas et al. | Application of 3d reservoir geological model on es1 formation, block nv32, shenvsi oilfield, China | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| Lescoffit et al. | Quantifying the impact of fault modeling parameters on production forecasting for clastic reservoirs | |
| Shao et al. | 3D Geological Modeling under Extremely Complex Geological Conditions. | |
| Ali et al. | Integrating 4D seismic data with production related effects at Enfield, North West Shelf, Australia | |
| Martin et al. | Three-dimensional modeling of Pennsylvanian sandstone units in the mature Dudley oil field, Illinois, USA | |
| CN107991702B (en) | The simple accurate determination method of normal fault timing | |
| Ma et al. | Development of a Full-Field Parallel Model to Design Pressure Maintenance Project in the Wara Reservoir, Greater Burgan Field, Kuwait | |
| Bezkhodarnov et al. | Prediction of Reservoir Properties from Seismic Data by Multivariate Geostatistics Analysis | |
| Shepherd | 3-D geocellular modeling | |
| Wang et al. | Super-large scale horizontal well position optimization method and infilling practice in high water-cut large complex fluvial reservoir based on multidisciplinary innovative techniques | |
| Podnebesnykh et al. | New approach to the evaluation of the structure of initial reserves in Ozhginskoe gas-oil field | |
| Gazar et al. | Adapting Field Development Plan of Compartmentalized Oil Rim Reservoir with Significant Gas Cap and Water Aquifer-Case Study |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20120215 Assignee: SICHUAN JISAITE TECHNOLOGY CO.,LTD. Assignor: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. Contract record no.: 2014510000163 Denomination of invention: Method for manufacturing seismic slice with overthrust fault Granted publication date: 20130814 License type: Exclusive License Record date: 20141016 |
|
| LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180129 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan 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. |