CN107703544B - Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics - Google Patents

Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics Download PDF

Info

Publication number
CN107703544B
CN107703544B CN201710885605.7A CN201710885605A CN107703544B CN 107703544 B CN107703544 B CN 107703544B CN 201710885605 A CN201710885605 A CN 201710885605A CN 107703544 B CN107703544 B CN 107703544B
Authority
CN
China
Prior art keywords
oil
gas
indication
grid node
oil gas
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
CN201710885605.7A
Other languages
Chinese (zh)
Other versions
CN107703544A (en
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 CN201710885605.7A priority Critical patent/CN107703544B/en
Publication of CN107703544A publication Critical patent/CN107703544A/en
Application granted granted Critical
Publication of CN107703544B publication Critical patent/CN107703544B/en
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/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/51Migration
    • G01V2210/512Pre-stack
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/66Subsurface modeling
    • G01V2210/665Subsurface modeling using geostatistical modeling
    • G01V2210/6652Kriging

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

Oil gas forecasting method is changed with offset distance based on the indication using prestack seismic amplitude of geostatistical analysis the present invention provides a kind of.The oil gas forecasting method includes: to establish oil and gas indication curve according to log data, and on oil and gas indication curve, oil gas section has the first oil and gas indication value, and non-oil gas section has the second oil and gas indication value;AVO attribute volume is calculated;Analysis obtains the oil and gas indication curve optimal AVO attribute volume best with well bypass road AVO attribute data correlation;Oil gas section part non-in oil and gas indication curve and oil gas section part are subjected to histogram analysis with optimal AVO attribute volume respectively, the subordinating degree function relationship of non-oil gas section part and oil gas section part in oil and gas indication curve is calculated;Geostatistical analysis is carried out to oil and gas indication curve and obtains variogram;Oil gas sunykatuib analysis, which is carried out, using subordinating degree function relationship and variogram obtains petroleum-gas prediction result.It can be improved the precision and resolution ratio of fluid prediction according to the method for the present invention.

Description

Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics
Technical field
The invention belongs to the fluid prediction methods in geophysical exploration technology, and more particularly, being related to one kind can be high Predict to precision, high-resolution the geostatistics AVO analysis method of thin reservoir fluid.
Background technique
AVO (Amplitude variation offset, amplitude variation with Offset) technology is that a kind of research earthquake is anti- Penetrate the technology that amplitude changes with geophone offset, amplitude is the letter of subsurface rock and pore-fluid elastic parameter with the variation of offset distance Number, therefore the subsurface lithologic and its property of pore fluid reflected according to AVO changing rule can be used to directly prediction oil Gas and estimation reservoir lithology parameter.Traditional AVO analysis is exactly using Zoeppritz equation or its approximate equation by actual seismic road A kind of inversion method of the seismologic parameter (such as density, p-and s-wave velocity) of collection record inverting estimation rock.Conventional AVO analysis method Usually nonlinear prestack inversion problem linearization, the stability of solution is affected by initial model, and is easily trapped into Local minimum, and AVO analyzes the poor prediction for being unfavorable for reservoir fluid of the obtained result resolution ratio of inversion method.
Summary of the invention
Aiming at the problems existing in the prior art, it is an object of the invention to solve in deficiencies of the prior art At least one of.For example, it is an object of the invention to solve the problems, such as that traditional AVO analysis resolution ratio is lower, precision is not high.
To achieve the goals above, the present invention provides a kind of based on the indication using prestack seismic amplitude of geostatistical analysis with offset distance Change oil gas forecasting method.The oil gas forecasting method includes: A, establishes oil and gas indication curve according to log data, in the oil On gas indicative curve, oil gas section has the first oil and gas indication value, and non-oil gas section has the second oil and gas indication value;B, it is calculated AVO attribute volume;C, analysis obtains the oil and gas indication curve optimal AVO category best with well bypass road AVO attribute data correlation Property data volume;D, by oil gas section part non-in oil and gas indication curve and oil gas section part respectively with the optimal AVO attribute data Body carries out histogram analysis, and the subordinating degree function of non-oil gas section part and oil gas section part in oil and gas indication curve is calculated and closes System;E, geostatistical analysis is carried out to oil and gas indication curve and obtains variogram;F, using the subordinating degree function relationship and The variogram carries out oil gas sunykatuib analysis, obtains petroleum-gas prediction result.In an exemplary embodiment of the present invention, institute State amplitude variation with Offset attribute volume include L1 (norm), L2 (two norms), P+G, P-G, P*G, L1P+G, L1P-G, L1P*G, L2P+G, L2P-G, L2P*G, wherein P represents intercept, and G represents gradient.
In an exemplary embodiment of the present invention, the step F may include: F01, in amplitude variation with Offset Data grids are created on attribute volume;Grid node on F02, lookup data grids, and judge whether grid node needs to simulate, Wherein, needing grid node to be simulated is the unknown grid node of uncertain oil and gas indication value, does not need the grid simulated Node is the known grid node for having determined that oil and gas indication value;F03, lookup need in unknown grid node field to be simulated Grid node is known, until known grid node reaches predetermined quantity;F04, pass through Kriging regression algorithm, utilize the predetermined number The known grid node of amount calculate unknown grid node gram in gold number, then using gram in gold number and the variogram calculate Golden probability value in obtaining gram;F05, using unknown grid node gram in golden probability value and corresponding amplitude variation with Offset category The value of property data volume, Combined Ration determines the oil and gas indication value of unknown grid node to the subordinating degree function relationship, and will determine The unknown grid node of oil and gas indication value is as known grid node;F06, step F02 to F05 is repeated, until determining institute The oil and gas indication value for the unknown grid node having.
In an exemplary embodiment of the present invention, it in the step F02, can sequentially be swum according to random walk It walks data grids and searches grid node;
In an exemplary embodiment of the present invention, it in the step F03, can be carried out in a manner of spiral search It searches.
In an exemplary embodiment of the present invention, in the step F04, using gram in Kriging regression algorithm In golden weight coefficient in golden equation calculation gram, calculate unknown grid node gram in golden probability value.
Compared with prior art, the beneficial effect comprise that improving the precision and resolution ratio of fluid prediction.
Detailed description of the invention
Fig. 1 shows the petroleum-gas prediction side prestack AVO based on geostatistics of an exemplary embodiment of the present invention The flow chart of method oil gas simulation steps.
Fig. 2 shows the petroleum-gas prediction sides prestack AVO based on geostatistics of an exemplary embodiment of the present invention The flow chart of method point simulation.
Specific embodiment
Hereinafter, the folding based on geostatistics that the present invention will be described in detail will be carried out in conjunction with exemplary embodiment and attached drawing Preceding AVO oil gas forecasting method.
Prestack AVO oil gas forecasting method according to an exemplary embodiment of the present invention based on geostatistics may include as Lower step:
A, oil and gas indication curve is established according to log data, on the oil and gas indication curve, oil gas section has the first oil Gas indicated value (such as 1), non-oil gas section have the second oil and gas indication value (such as 0).
B, AVO attribute volume is calculated.AVO attribute volume include L1 (norm), L2 (two norms), P+G (cut Away from adding gradient), P-G (intercept subtracts gradient), P*G (intercept boarding degree), L1P+G, L1P-G, L1P*G, L2P+G, L2P-G, L2P*G Equal attribute volumes.Wherein, P represents intercept, and G represents gradient.
C, analysis obtains the oil and gas indication curve optimal AVO attribute data best with well bypass road AVO attribute data correlation Body.
D, by oil gas section part non-in oil and gas indication curve and oil gas section part respectively with the optimal AVO attribute volume Histogram analysis is carried out, the subordinating degree function of non-oil gas section part and oil gas section part in oil and gas indication curve is calculated and closes System.
E, geostatistical analysis is carried out to oil and gas indication curve and obtains variogram.
F, oil gas sunykatuib analysis is carried out using the subordinating degree function relationship and the variogram, obtains petroleum-gas prediction knot Fruit.The step F is specifically included:
F01, data grids are created on amplitude variation with Offset attribute volume.
F02, according to random walk, sequentially migration data grids search grid node, and judge whether grid node needs Simulation, wherein needing grid node to be simulated be the unknown grid node for not knowing oil and gas indication value, does not need to be simulated Grid node is the known grid node for having determined that oil and gas indication value.
F03, the known grid node needed in unknown grid node field to be simulated is searched in a manner of spiral search, directly Reach predetermined quantity to known grid node.
F04, pass through Kriging regression algorithm, the known grid node of the predetermined quantity is utilized to calculate unknown grid node Gram in gold number, then using gram in gold number and the variogram be calculated gram in golden probability value.
F05, using unknown grid node gram in golden probability value and corresponding amplitude variation with Offset attribute volume Value, Combined Ration determines the oil and gas indication value of unknown grid node to the subordinating degree function relationship, and will be determined that oil gas refers to The unknown grid node of indicating value is as known grid node.
F06, step F02 to F05 is repeated, until determining the oil and gas indication value of all unknown grid nodes.
Major technique thinking of the invention is the attribute volume for calculating AVO, by series of computation, is converted to 0 and 1 composition Data volume, 1 Indication of Oil-Gas, the 0 non-oil gas of instruction realizes oil and gas detection.In an exemplary embodiment of the present invention, it is based on The prestack AVO oil gas forecasting method of geostatistics may include steps of:
(1) a HCI indicative curve is established, HCI curve is an oil gas section HCI value by well logging data analysis first The indicated value curve for being 0 for 1, non-oil gas segment value.The indicated value of however, the present invention is not limited thereto, oil gas section and non-oil gas can be it It is any both can to distinguish different values.Here, HCI is the abbreviation of Hydrocarbon Indicator, and HCI curve is exactly Oil and gas indication curve, accuracy are controlled by log, and HCI is the curve of an Indication of Oil-Gas and water being made of 0 and 1.
(2) by calculating that tradition AVO analyzes as a result, obtain a series of relevant AVO attribute volumes, for example, L1, L2, P+G, The attribute volumes such as P-G, P*G, L1P+G, L1P-G, L1P*G, L2P+G, L2P-G and L2P*G, wherein P is intercept, and G is gradient, L1 For a norm, L2 is two norms.
(3) correlation between analysis HCI curve and AVO attribute track data by well, analyzes and obtains HCI-AVO attribute phase The highest HCI-AVO attribute curve of closing property, and it is saved.
(4) histogram analysis is carried out to HCI-0 and HCI-1 curve to the HCI-AVO attribute curve of preservation, and be calculated The subordinating degree function relationship of HCI-0 and HCI-1.Here, HCI-AVO attribute curve is by HCI curve and well bypass road AVO attribute The best AVO attribute of data dependence, and by its AVO attribute data corresponding to 0 and 1 save respectively to HCI-0 and HCI-1 In, so that followed by using, specifically, obtained subordinating degree function is the 4. step use in (6) step.HCI-0 is Refer to non-oil gas section part in oil and gas indication curve, oil gas section part in HCI-1 oil and gas indication curve.
(5) geostatistical analysis is carried out to the HCI indicative curve that front inputs and obtains variogram, for HC simulation point It is used in analysis method.Specifically, variogram the 3. using in step in (6) step)
(6) HC simulation is carried out, process is as shown in Figure 1 and Figure 2, it is realized by following steps:
1. according to random walk sequentially migration by not no data value grid node.Sequentially migration is random walk Refer to that herein is in sequentially search grid node in the random walk of foundation by certain sequential search grid node.
2. the known node (i.e. known grid node) in neighborhood is searched in a manner of spiral search at the grid node. Known node has value, and practical is exactly the point on log, has had built up grid node here, and belonging to well logging The indicated value of the HCI in well logging is imparted on the point of position.
3. by ordinary kriging interpolation algorithm gram in golden weight coefficient in golden equation calculation gram, calculate the point (i.e. Unknown mesh point) indicated value probability value P*;I.e. the value of unknown point is sought by ordinary kriging interpolation method.
4. using the value of AVO attribute volume corresponding to the indicator Kriging probability value P* and the grid node sought, Combined Ration to degree of membership relational matrix, two o'clock (gold probability value P* and AVO attribute value in gram) determine HC (oil gas) indicated value (0 or 1) indicated value that, then this is calculated is added to known point.As soon as a unknown point has been calculated, the value of unknown point is obtained, so After add it to known points according in, become the known data point that can be used in subsequent cycle.
5. repeat step 1. -4., until calculate the indicated values of all unknown points, EP (end of program).
The model of random distribution is randomly generated using geostatistical analysis process by the present invention, and using solution global optimum The method of change problem can effectively solve the problems, such as that the nonuniqueness in inverting and tradition AVO inverting are affected by initial model The problem of, it is more preferable compared to traditional AVO inversion method effect, and the resolution ratio of reservoir fluid prediction can be greatly improved, it improves To the accuracy rate of thin-layer fluid identification.
Although having been combined attached drawing above and exemplary embodiment describing the present invention, those of ordinary skill in the art It will be apparent to the skilled artisan that can be carry out various modifications to above-described embodiment in the case where not departing from spirit and scope of the claims.

Claims (5)

1. a kind of change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics, which is characterized in that the oil Gas prediction technique the following steps are included:
A, oil and gas indication curve is established according to log data, on the oil and gas indication curve, there is oil gas section the first oil gas to refer to Indicating value, non-oil gas section have the second oil and gas indication value;
B, amplitude variation with Offset attribute volume is calculated;
C, analysis obtains the optimal vibration best with well bypass road amplitude variation with Offset attribute data correlation of oil and gas indication curve Width is with offset distance change to attributes data volume;
D, oil gas section part non-in oil and gas indication curve and oil gas section part are changed with the optimal amplitude with offset distance respectively Attribute volume carries out histogram analysis, and being subordinate to for non-oil gas section part and oil gas section part in oil and gas indication curve is calculated Spend functional relation;
E, geostatistical analysis is carried out to oil and gas indication curve and obtains variogram;
F, carry out oil gas sunykatuib analysis using the subordinating degree function relationship and the variogram, obtain petroleum-gas prediction as a result, Wherein, the step F includes:
F01, data grids are created on amplitude variation with Offset attribute volume;
Grid node on F02, lookup data grids, and judge whether grid node needs to simulate, wherein need net to be simulated Lattice node is the unknown grid node of uncertain oil and gas indication value, and not needing the grid node simulated is to have determined that oil gas refers to The known grid node of indicating value;
F03, the known grid node needed in unknown grid node field to be simulated is searched, until known grid node reaches pre- Fixed number amount;
F04, pass through Kriging regression algorithm, gram of unknown grid node is calculated using the known grid node of the predetermined quantity In gold number, then using gram in gold number and the variogram be calculated gram in golden probability value;
F05, using unknown grid node gram in golden probability value and corresponding amplitude variation with Offset attribute volume value, Combined Ration determines the oil and gas indication value of unknown grid node to the subordinating degree function relationship, and oil and gas indication value will be determined Unknown grid node as known grid node;
F06, step F02 to F05 is repeated, until determining the oil and gas indication value of all unknown grid nodes.
2. according to claim 1 change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics, Be characterized in that, the amplitude variation with Offset attribute volume include L1, L2, P+G, P-G, P*G, L1P+G, L1P-G, L1P*G, L2P+G, L2P-G and L2P*G, wherein P is intercept, and G is gradient, and L1 is a norm, and L2 is two norms.
3. according to claim 1 change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics, It is characterized in that, in the step F02, according to random walk, sequentially migration data grids search grid node.
4. according to claim 1 change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics, It is characterized in that, in the step F03, is searched in a manner of spiral search.
5. according to claim 1 change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics, Be characterized in that, in the step F04, using in Kriging regression algorithm gram in golden weight coefficient in golden equation calculation gram, Calculate unknown grid node gram in golden probability value.
CN201710885605.7A 2017-09-27 2017-09-27 Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics Active CN107703544B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710885605.7A CN107703544B (en) 2017-09-27 2017-09-27 Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710885605.7A CN107703544B (en) 2017-09-27 2017-09-27 Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics

Publications (2)

Publication Number Publication Date
CN107703544A CN107703544A (en) 2018-02-16
CN107703544B true CN107703544B (en) 2019-06-04

Family

ID=61174526

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710885605.7A Active CN107703544B (en) 2017-09-27 2017-09-27 Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics

Country Status (1)

Country Link
CN (1) CN107703544B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109633781B (en) * 2018-08-16 2021-03-02 清能艾科(深圳)能源技术有限公司 Geological property acquisition method and device, electronic equipment and storage medium

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5440525A (en) * 1993-10-25 1995-08-08 Atlantic Richfield Company Seismic data hydrocarbon indicator
US5784334A (en) * 1996-03-13 1998-07-21 Atlantic Richfield Company Method and system for detecting hydrocarbon reservoirs using amplitude versus offset analysis of seismic signals
CN101281253A (en) * 2007-04-06 2008-10-08 中国石油集团东方地球物理勘探有限责任公司 Method for enhancing oil gas detecting accuracy using vibration amplitude with off-set distance variation characteristic
CN102147478A (en) * 2010-12-29 2011-08-10 中国海洋大学 Pre-stack low frequency signal recognition method of complex oil pool
CN102819038A (en) * 2012-07-27 2012-12-12 中国石油天然气股份有限公司 Method and system for identifying oil and water in carbonate rocks
CN102866426A (en) * 2012-09-20 2013-01-09 成都晶石石油科技有限公司 Method for analyzing oil gas information of rock body by applying amplitude versus offset (AVO) large-angle trace gathers
CN103076630A (en) * 2011-10-26 2013-05-01 中国石油化工股份有限公司 Hydrocarbon detection method based on elastic impedance gradient
CN103257360A (en) * 2013-04-25 2013-08-21 刘立峰 Method for identifying carbonate rock fluid based on fuzzy C mean cluster
CN105425292A (en) * 2015-11-20 2016-03-23 中国石油天然气集团公司 Oil and gas prediction method and oil and gas prediction device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7095677B2 (en) * 2003-05-27 2006-08-22 Paradigm Geophysical Crossplot analysis of A.V.O. anomolies in seismic surveying

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5440525A (en) * 1993-10-25 1995-08-08 Atlantic Richfield Company Seismic data hydrocarbon indicator
US5784334A (en) * 1996-03-13 1998-07-21 Atlantic Richfield Company Method and system for detecting hydrocarbon reservoirs using amplitude versus offset analysis of seismic signals
CN101281253A (en) * 2007-04-06 2008-10-08 中国石油集团东方地球物理勘探有限责任公司 Method for enhancing oil gas detecting accuracy using vibration amplitude with off-set distance variation characteristic
CN102147478A (en) * 2010-12-29 2011-08-10 中国海洋大学 Pre-stack low frequency signal recognition method of complex oil pool
CN103076630A (en) * 2011-10-26 2013-05-01 中国石油化工股份有限公司 Hydrocarbon detection method based on elastic impedance gradient
CN102819038A (en) * 2012-07-27 2012-12-12 中国石油天然气股份有限公司 Method and system for identifying oil and water in carbonate rocks
CN102866426A (en) * 2012-09-20 2013-01-09 成都晶石石油科技有限公司 Method for analyzing oil gas information of rock body by applying amplitude versus offset (AVO) large-angle trace gathers
CN103257360A (en) * 2013-04-25 2013-08-21 刘立峰 Method for identifying carbonate rock fluid based on fuzzy C mean cluster
CN105425292A (en) * 2015-11-20 2016-03-23 中国石油天然气集团公司 Oil and gas prediction method and oil and gas prediction device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"A Comparison of Hydrocarbon indicators derived from AVO analysis";Hong Feng 等;《SEG Technical Program Expanded Abstracts 2007》;20071231;第279-282页
"A new Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset";Changcheng Liu 等;《AAPG/SEG International Conference&Exhibition2016》;20160125;第1-12页
"致密砂岩储层的非线性AVO分析方法";杜增利 等;《成都理工大学学报(自然科学版)》;20120630;第39卷(第3期);第251-256页

Also Published As

Publication number Publication date
CN107703544A (en) 2018-02-16

Similar Documents

Publication Publication Date Title
CN110031896B (en) Seismic random inversion method and device based on multi-point geostatistics prior information
CN106054244B (en) The LPF of window multiple dimensioned full waveform inversion method when blocking
CN105388518B (en) A kind of centroid frequency and earthquake inversion of quality factor method in the united well of Frequency spectrum ratio
CN109407151B (en) Time-domain full waveform inversion method based on wave field local correlation time shift
CN106842295A (en) The waveform inversion method of logging information constrained
CN104237937B (en) Pre-stack seismic inversion method and system thereof
CN104502997A (en) Method for using fracture density curve to forecast fracture density body
CN106405647B (en) A kind of tuning inversion method on sedimentary formation thickness
CN108645994A (en) A kind of geology stochastic inversion methods and device based on Multiple-Point Geostatistics
CN104570101A (en) AVO (amplitude versus offset) three-parameter inversion method based on particle swarm optimization
CN109425900A (en) A kind of Seismic Reservoir Prediction method
CN103245972A (en) Method for determining complex geologic structure in two-dimensional space
CN106338764B (en) The super wire stripping recognition methods of bio-inspired computing stratigraphic trap oil-gas reservoir
CN108957534B (en) Method and device for predicting gas saturation
RU2490677C2 (en) Method for complex processing of geophysical data "litoscan" system for realising said method
CN107703544B (en) Change oil gas forecasting method with offset distance based on the indication using prestack seismic amplitude of geostatistics
CN112147677B (en) Method and device for generating parameter tag data of oil and gas reservoir
CN103543478A (en) Geologic morphological interpolation KM (Kriging and Multiple-point geostatistics) method
CN113219531B (en) Dense sandstone gas-water distribution identification method and device
CN104199106A (en) Seismic data residual static correction method based on ant colony algorithm
CN105589096B (en) A kind of sedimentary facies belt division methods based on D-S evidence theory
EP3186665B1 (en) Method of characterising a subsurface volume
CN115880455A (en) Three-dimensional intelligent interpolation method based on deep learning
CN110927786B (en) Seismic lithofacies prediction method and system based on virtual well random simulation
CN104297781A (en) Ray parameter domain unconstrained elastic parameter inversion method

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
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: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

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

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

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

Effective date of registration: 20201112

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.