CN107817519A - Basement lithology recognition methods and system based on weight magnetic shake Comprehensive interpretation - Google Patents
Basement lithology recognition methods and system based on weight magnetic shake Comprehensive interpretation Download PDFInfo
- Publication number
- CN107817519A CN107817519A CN201610824523.7A CN201610824523A CN107817519A CN 107817519 A CN107817519 A CN 107817519A CN 201610824523 A CN201610824523 A CN 201610824523A CN 107817519 A CN107817519 A CN 107817519A
- Authority
- CN
- China
- Prior art keywords
- anomaly
- gravity
- magnetic
- seismic
- lithology
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000005484 gravity Effects 0.000 claims abstract description 59
- 239000000758 substrate Substances 0.000 claims abstract description 57
- 230000005389 magnetism Effects 0.000 claims abstract description 23
- 238000009826 distribution Methods 0.000 claims abstract description 22
- 238000011160 research Methods 0.000 claims abstract description 21
- 238000005553 drilling Methods 0.000 claims abstract description 16
- 239000011435 rock Substances 0.000 claims description 39
- 230000002159 abnormal effect Effects 0.000 claims description 22
- 238000001228 spectrum Methods 0.000 claims description 20
- 238000004458 analytical method Methods 0.000 claims description 13
- 238000007689 inspection Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 230000002547 anomalous effect Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 238000000205 computational method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 241000923606 Schistes Species 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000010438 granite Substances 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 2
- 241001026509 Kata Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 240000000249 Morus alba Species 0.000 description 1
- 235000008708 Morus alba Nutrition 0.000 description 1
- 240000001987 Pyrus communis Species 0.000 description 1
- 235000014443 Pyrus communis Nutrition 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010224 classification analysis Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 239000011044 quartzite Substances 0.000 description 1
- 239000010454 slate Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of basement lithology recognition methods based on weight magnetic shake Comprehensive interpretation and system, this method includes the Seismic reflection character of Study of recognition area basement lithology, lithology and the pattern of seismic reflection are established according to recognition result, and attribute is carried out to seismic data and calculated, the result calculated according to attribute generates seismic properties distribution map；Gravity earthquake peeling is carried out to Bouguer anomaly data to calculate to obtain the residual gravity anomaly of survey region substrate；Magnetic anomalies continuation calculating is carried out to magnetic anomalies data to obtain the residual magnetism of survey region substrate exception；Gravity-magnetic anomaly block plan is established according to residual gravity anomaly and residual magnetism extremely；Overlapping research is carried out based on seismic properties distribution map and gravity-magnetic anomaly block plan, and combines drilling well explanation and the basement lithology of survey region is identified.This method realizes the identification of basement lithology in the plane, improves the accuracy rate of Lithology Discrimination, and to Eastern China, other basin basement lithology identifications have important reference.
Description
Technical field
The invention belongs to field of geophysical exploration, more particularly to a kind of fundamental rock based on weight magnetic shake Comprehensive interpretation
Property recognition methods and system.
Background technology
The basement hydrocarbon reservoir Exploration Domain new as one, contains huge Exploration Potential.Basement hydrocarbon reservoir reservoir is commented
Valency is the emphasis of basement hydrocarbon reservoir description, and the identification of the basement lithology of basement rock is the basis of basement rock evaluating reservoir.Lithology can not be known
Not or identification mistake will mislead a series of result of later stage geological prospectings.
By determining the lithology distributions feature of basement rock, storage and collection performance and the formation of different lithology basement rock could be implemented
Mechanism, the pool forming rules of basement hydrocarbon reservoir are further disclosed, theory support are provided for basement hydrocarbon reservoir Exploration Potential, to continue to seek
Large-medium oil-and-gas fields are looked for provide foundation.The physical property and production capacity difference of different lithology reservoir are larger, and lithology is to shadows such as well logging, earthquakes
The influence rung often than fluid, crack is bigger, therefore identifies that basement lithology is the basic and crucial of development basement rock evaluating reservoir.
The Main Means of prior art identification basement lithology include being identified based on well-log information, based on aero magnetic data entering
Row identification is identified using the method for seismic prospecting.There are the following problems for the above method, and method of seismic prospecting has high score
Resolution, but because deep seismic reflection data quality is poor, effective wave energy is weaker, can not be obtained completely according to seismic data
Infrastructure form.And although aero magnetic data can effectively reflect fracture, igneous rock, metamorphic rock isodensity, magnetic anomaly body,
But due to target seeker gyro it is intrinsic the characteristics of, it serious multi-solution is present.
To sum up, the present invention propose it is a kind of based on aero magnetic data and seismic prospecting data carry out the method for comprehensive analysis explanation come
Improve the reliability of basement lithology identification.
The content of the invention
One of technical problems to be solved by the invention are to need to provide a kind of reliability for improving basement lithology identification
Method.
In order to solve the above-mentioned technical problem, embodiments herein provide firstly a kind of based on weight magnetic shake integrated interpretation
The basement lithology recognition methods of technology, including：The Seismic reflection character of Study of recognition area basement lithology, is established according to recognition result
Lithology and the pattern of seismic reflection, and attribute is carried out to seismic data and calculated, the result calculated according to attribute generates seismic properties
Distribution map；Gravity-earthquake peeling is carried out to Bouguer anomaly data to calculate to obtain the residual gravity anomaly of survey region substrate；It is right
Magnetic anomalies data carries out magnetic anomalies continuation calculating to obtain the residual magnetism of survey region substrate exception；According to the residue
Gravity anomaly and residual magnetism establish gravity-magnetic anomaly block plan extremely；Based on the seismic properties distribution map and gravity-magnetic anomaly subregion
Figure carries out overlapping research, and combines drilling well and explain, the basement lithology of survey region is identified.
Preferably, include before attribute calculating is carried out to seismic data：According to collection research area substrate data and
Collection appear and drill cores sample determine need carry out Lithology Discrimination lithology classification, determine the earthquake of each reflecting interface
Reflectance signature, horizon calibration is carried out to the important seismic interface of survey region；Seismic interpretation, institute are carried out according to the result of horizon calibration
Stating seismic interpretation includes the explanation of layer position and fault interpretation；Gridding processing is carried out to the seismic interpretation result in basement top boundary line,
Result based on gridding processing carries out attribute calculating.
Preferably, the attribute, which calculates, includes the calculating of RMS amplitude attribute, waveform clustering attribute calculates and relative wave resistance
Family dependents of military personel in the liberated areas's property calculates.
Preferably, it is described that lithology and the pattern of seismic reflection are established according to recognition result, including：Based on the seismic reflection
Feature, lithology and the pattern of seismic reflection are established by single well analysis and well shake contrast, wherein, the Seismic reflection character refers to
Seismic reflection border, inside and out feature.
Preferably, it is described that Bouguer anomaly data progress gravity-earthquake peeling is calculated and remained with obtaining survey region substrate
Remaining gravity anomaly, including：Bouguer anomaly data is carried out superficial part sedimentary calculate extremely with obtain survey region substrate rock mass with
The gravity anomaly of upper sedimentary；Deep regional field is carried out to Bouguer anomaly data to calculate extremely to obtain survey region substrate rock mass
The regional field in following deep is abnormal；The superficial part sedimentary exception and deep regional field exception are subtracted with Bouguer anomaly data to obtain
To the residual gravity anomaly of survey region substrate.
Preferably, the superficial part sedimentary exception that carried out to Bouguer anomaly data calculates, including：Seismic data is solved
Release, establish Stratigraphic framework, to extract main density interface；Sedimentary density mould is established based on the main density interface
Type；Forward modelling is carried out based on the sedimentary density model, to obtain more than the substrate rock mass caused by shallow-layer density interface
The gravity anomaly of sedimentary.
Preferably, the deep regional field exception that carried out to Bouguer anomaly data calculates, including：To from Bouguer anomaly data
In the gravity anomaly upward continuation that observes, obtain deep regional field；The gravity anomaly is calculated by wavelet analysis
Frequency spectrum；Deep diopter filtering is carried out to the frequency spectrum, it is abnormal that deep regional field is chosen from filter result.
Preferably, it is described that magnetic anomalies continuation calculating is carried out to magnetic anomalies data to obtain the residue of survey region substrate
Magnetic anomalies, including：Magnetic anomalies upward continuation to being observed from magnetic anomalies data, the magnetic anomalies are converted to
Space more than inspection surface；The frequency spectrum of the magnetic anomalies is calculated by wavelet analysis；Remove the superficial part in the frequency spectrum
Magnetic anomalies；Moving average method is used to be filtered the frequency spectrum to remove the deep magnetic anomalies in the frequency spectrum.
Preferably, it is described that gravity-magnetic anomaly block plan is established according to the residual gravity anomaly and residual magnetism extremely, including：
The residual gravity anomaly and residual magnetism are filtered extremely using moving average method；To the residual gravity anomaly and remain
Remaining magnetic anomalies carry out vertical second dervative calculating to obtain gravity-magnetic anomaly block plan.
Embodiments herein additionally provides a kind of basement lithology identifying system based on weight magnetic shake Comprehensive interpretation, bag
Include：Seismic properties distribution map establishes module, the Seismic reflection character of its Study of recognition area basement lithology, is established according to recognition result
Lithology and the pattern of seismic reflection, and attribute is carried out to seismic data and calculated, the result calculated according to attribute generates seismic properties
Distribution map；Residual gravity anomaly acquisition module, it carries out gravity-earthquake peeling to Bouguer anomaly data and calculated to obtain research area
The residual gravity anomaly of domain substrate；Residual magnetism exception acquisition module, it carries out magnetic anomalies continuation meter to magnetic anomalies data
Calculate to obtain the residual magnetism of survey region substrate exception；Gravity-magnetic anomaly block plan establishes module, and it is according to the residual gravity
Abnormal and residual magnetism establishes gravity-magnetic anomaly block plan extremely；Overlap identification module, its be based on the seismic properties distribution map and
Gravity-magnetic anomaly block plan carries out overlapping research, and combines drilling well and explain, the basement lithology of survey region is identified.
Compared with prior art, one or more of such scheme embodiment can have the following advantages that or beneficial to effect
Fruit：
By combining drilling data and well log interpretation, there is provided a kind of technical method for identifying sedimentary basin basement lithology,
The identification of basement lithology is realized in plane, improves the accuracy rate of Lithology Discrimination, to other basin basement lithologies of Eastern China
Identification has important reference.
Other advantages, the target of the present invention, and feature will be illustrated in the following description to a certain extent, and
And to a certain extent, based on will be apparent to those skilled in the art to investigating hereafter, Huo Zheke
To be instructed from the practice of the present invention.The target and other advantages of the present invention can be wanted by following specification, right
Specifically noted structure is sought in book, and accompanying drawing to realize and obtain.
Brief description of the drawings
Accompanying drawing is used for providing to the technical scheme of the application or further understanding for prior art, and constitution instruction
A part.Wherein, the accompanying drawing for expressing the embodiment of the present application is used for the technical side for explaining the application together with embodiments herein
Case, but do not form the limitation to technical scheme.
Fig. 1 is the flow according to the basement lithology recognition methods that Comprehensive interpretation is shaken based on weight magnetic of the embodiment of the present invention
Schematic diagram；
Fig. 2 is the hexahedron coordinate schematic diagram of hexahedron forward modelling；
Fig. 3 is the explanation according to the basement lithology recognition methods that Comprehensive interpretation is shaken based on weight magnetic of the embodiment of the present invention
Route schematic diagram；
Fig. 4 is the LS areas basement lithology distribution schematic diagram according to an example of the invention；
Fig. 5 is the structure according to the basement lithology identifying system that Comprehensive interpretation is shaken based on weight magnetic of the embodiment of the present invention
Schematic diagram.
Embodiment
Embodiments of the present invention are described in detail below with reference to drawings and Examples, and how the present invention is applied whereby
Technological means solves technical problem, and the implementation process for reaching relevant art effect can fully understand and implement according to this.This Shen
Each feature that please be in embodiment and embodiment, can be combined with each other under the premise of not colliding, the technical scheme formed
Within protection scope of the present invention.
Fig. 1 is the flow according to the basement lithology recognition methods that Comprehensive interpretation is shaken based on weight magnetic of the embodiment of the present invention
Schematic diagram, as illustrated, this method comprises the following steps：
Step S110, the Seismic reflection character of Study of recognition area basement lithology, lithology and earthquake are established according to recognition result
The pattern of reflection, and attribute is carried out to seismic data and calculated, the result calculated according to attribute generates seismic properties distribution map.
Step S120, gravity-earthquake peeling is carried out to Bouguer anomaly data to calculate to obtain the residue of survey region substrate
Gravity anomaly.
Step S130, magnetic anomalies continuation calculating is carried out to magnetic anomalies data to obtain the remanence of survey region substrate
Power is abnormal.
Step S140, gravity-magnetic anomaly block plan is established according to residual gravity anomaly and residual magnetism extremely.
Step S150, overlapping research is carried out based on seismic properties distribution map and gravity-magnetic anomaly block plan, and combines drilling well solution
Release and the basement lithology of survey region is identified.
Specifically, before being identified, the substrate data of survey region is first collected, for example, being appeared by gathering and analyzing
Determine to need the lithology classification for carrying out Lithology Discrimination with drill cores sample.The data of collection includes having bored the drilling well for meeting substrate
Data, drilling well substrate well section curvilinear characteristic to be analyzed, target area and its adjacent area petrophysics property are summarized in research, including
Density, magnetic etc..Thin section identification under core observation and mirror is carried out to the sample of collection, understands the basement lithology type in research area.
In step s 110, before attribute calculating is carried out to seismic data, horizon calibration is carried out to survey region first,
The concrete meaning of horizon calibration includes boring research area the well development composite traces work for meeting substrate, while determines each reflection circle
The Seismic reflection character in face.
Then seismic interpretation is carried out according to the result of horizon calibration, and basement top circle is determined by the result of seismic interpretation
Line.Seismic interpretation work is carried out to research area's sedimentary, seismic interpretation work includes the explanation of layer position and fault interpretation, emphasized
Basement top boundary line, while by explaining key section, build the Stratigraphic framework in research area.
Gridding processing is carried out to basement top boundary line again, attribute calculating is carried out with the result handled based on gridding.It is right
The grid that at least 10*10 is carried out in basement top boundary line is explained, then its gridding is handled, and is calculated for ensuing seismic properties
Carry out data preparation.
Attribute calculating is carried out to seismic data includes the calculating of RMS amplitude attribute, the calculating of waveform clustering attribute and relative ripple
Impedance attribute calculates.
Specifically, amplitude attribute namely RMS amplitude are that the average value of amplitude attribute square extracts square root again, due to shaking
Amplitude is at averagely preceding square, and therefore, it is very sensitive to king-sized amplitude, therefore can be used for identifying lithology.
Waveform clustering attribute is spread situation of the geological information in transverse direction of same nature.Based on waveform similarity
Seismic facies neural net method carries out classification analysis in window when waiting to the variation characteristic of seismic waveform, and analysis waveform is in the plane
Changing rule.
Relative Wave Impedance attribute is a kind of wave impedance inversion attribute of no well modeling, its inversion formula such as expression formula (1) institute
Show：
In formula, R is reflectance factor, ρ1、ρ2And v1、v2Be respectively reflecting interface up and down rock density and seismic wave at it
Between spread speed.
If ρ v are ρ2v2With ρ1v1Average value, then have：
That is the integration of reflectance factor is proportional to wave impedance ρ v natural logrithm.Its advantages of is that no well control can also be anti-
Drill, it is not necessary to the demarcation of reflectance factor, and also recursion cumulative errors are small.It is less that the well of substrate is met due to boring, some well substrates are bent
Line is not also complete, effectively can not model constraint inverting using drilling well.And the relative impedances inverting of no restriction from borehole data is then especially suitable for base
The inverting at bottom.
Further in step s 110, the Seismic reflection character of Study of recognition area basement lithology, it is anti-with earthquake to establish lithology
Emission mode.Specifically, identifying Seismic reflection character first, the Seismic reflection character that identification obtains is then based on, passes through individual well point
The relation that analysis and well shake contrast are established between the lithology of survey region substrate and seismic reflection pattern.Wherein, Seismic reflection character
Including the composite traces for combining drilling well, profile obtains.
In the step s 120, gravity earthquake peeling calculating is carried out to Bouguer anomaly data to remain to obtain survey region substrate
Remaining gravity anomaly, specifically includes following steps：
Step S1201, superficial part sedimentary is carried out to Bouguer anomaly data to calculate extremely to obtain survey region substrate rock mass
The gravity anomaly of above sedimentary
Step S1202, deep regional field is carried out to Bouguer anomaly data to calculate extremely to obtain survey region substrate rock mass
The regional field in following deep is abnormal.
Step S1203, superficial part sedimentary exception and deep regional field exception are subtracted with Bouguer anomaly data to be studied
The residual gravity anomaly of substrate areas.
Wherein, superficial part sedimentary calculates extremely seeks to shell the gravity anomaly for studying area's substrate rock mass above sedimentary
From.When carrying out the calculating extremely of superficial part sedimentary to Bouguer anomaly data, seismic data is explained first main to extract
Density interface, be then based on main density interface and establish sedimentary density model, further according to the physical properties of rock data of statistics
Residual density is calculated, stratum Midst density interface and residual density are subjected to forward modelling, gone out according to these information forward modellings shallow
The gravity anomaly of substrate rock mass above sedimentary caused by layer density interface.
In one embodiment of the invention, forward modelling uses hexahedron forward modelling method, and detailed process is as follows：
It is a series of upright hexahedron that length, width and height are respectively △ x, △ y and △ z first by STRATIGRAPHIC DIVISION,
As shown in Figure 2.
To calculating point i, i=1 ..., M (M is a calculating point number), hexahedron j, j=1 ..., N (N is hexahedron number) have
Following relation,
In formula：
x1=xj-xi- △ x/2,
x2=xj-xi+ △ x/2,
y1=yj-yi- △ y/2,
y2=yj-yi+ △ y/2,
Z1=zj-zi。
In formula, △ g (i) are the GRAVITY ANOMALIES for calculating point, and σ (j) is that hexahedral density (is commonly used for ambient field
2.67kg/m3) density contrast, G=6.67 × 10-8CGS is gravitational constant.X, y, z is hexahedral length, width and height.
Deep regional field calculates extremely to be sought to peel off the regional field for studying deep below area's target substrate rock mass.Right
When Bouguer anomaly data carries out deep regional field and calculated extremely, first to the gravity anomaly that is observed from Bouguer anomaly data to
Upper continuation, gravity anomaly is converted to space more than inspection surface, the frequency of gravity anomaly is then calculated by wavelet analysis
Spectrum, then deep diopter filtering is carried out to the frequency spectrum of gravity anomaly, it is abnormal that deep regional field is chosen from filter result.
Specifically, gravity anomaly is stacking reaction of the earth's surface to all anomalous body information in underground, bulk effect is gravity
Principal character.It is abnormal cube to be inversely proportional to distance between anomalous body with observation station, i.e., exceptional value to inspection surface to anomalous body it
Between distance reaction it is sensitive.The speed that gravitational field is decayed with the increase of distance and tool residual density and the geologic body body of gravity
Product, buried depth are relevant.The gravitational field that volume is big, buries corresponding to deep geologic body is decayed slow with the increase for above prolonging distance；Small volume,
The gravitational field corresponding to shallow geologic body is buried with the increase for above prolonging distance to decay soon.Therefore, prolonging upwards using gravity anomaly
Deep variable density feature can be protruded by opening up, extremely can be close with qualitative analysis underground different depth by the upward continuation of different height
Spend the Regional Characteristics of body.By analyzing the change for the inspection surface exception that observed anomaly is converted to different height, can analyze
The distribution situation of different depth and scale anomalous body.
Upward continuation is exactly the gravity anomaly on inspection surface to be converted to space more than inspection surface.Frequency domain continuation calculates
It is conventional computational methods.Observed anomaly △ g are carried out first to pay formula conversion, abnormal frequency spectrum G (u, v)=FFT is calculated
{△g(x,y)}；Thereafter frequency spectrum is multiplied by continuation conversion factor Wherein z is the height of continuation
Degree, then frequency domain transformation result is carried out instead to pay formula conversion △ gu(x, y)=FFT-1{ F (u, v) }, just obtain different after continuation
Constant value.
Obtaining the abnormal (G of superficial part sedimentaryq) and the abnormal (G in deep regional fields) after, subtract superficial part with Bouguer anomaly and deposit
Layer is abnormal and deep regional field is abnormal, to obtain the residual gravity anomaly (G of survey region substratej)。
In step s 130, magnetic anomalies are carried out to magnetic anomalies data to calculate to obtain the remanence of survey region substrate
Power is abnormal, specifically includes following steps：
Step S1301, the magnetic anomalies upward continuation to being observed from magnetic anomalies data, magnetic anomalies are converted
Space more than to inspection surface.
Step S1302, the frequency spectrum of magnetic anomalies is calculated by wavelet analysis.
Step S1303, the superficial part magnetic anomalies in frequency spectrum are removed.
Step S1304, moving average method is used to be filtered above-mentioned frequency spectrum to remove the deep magnetic anomalies in frequency spectrum.
Specifically, magnetic force continuation computational methods are identical with gravity continuation computational methods, i.e., being prolonged using the upward of magnetic anomalies
Deep magnetic force change feature can be protruded by opening up, extremely can be with qualitative analysis underground different depth magnetic by the upward continuation of different height
The Regional Characteristics of power body.
Moving average method refers on geophysical anomaly figure, selectes the window of a certain size, will be all different in window
Constant value is done sums averagely, the exceptional value using required average value as window center point.By point away from or line-spacing moving window, repeat
This averaging method, until completing said process to view picture figure, this process is referred to as moving average.Moving average is equivalent to low pass filtered
Ripple.
Next in step S140, magnetic anomaly again is established according to the residual gravity anomaly and residual magnetism being calculated extremely
Normal block plan, is specifically included, and residual gravity anomaly and residual magnetism are filtered extremely using moving average method, to residue weight
Power is abnormal and residual magnetism carries out vertical second dervative calculating to obtain border gravity-magnetic anomaly block plan extremely, wherein, abnormal two
The Zero value line of rank vertical derivatives is substantially corresponding with the border of anomalous body.
Finally in step S150, obtain the substrate gravity-magnetic anomaly block plan obtained in step S140 and in step S110
Based seismic property distribution figure overlapping research, with reference to drilling well explain, to research area's substrate expansion Lithology Discrimination, finally give this
Area's basement lithology distribution map.
The entire flow of the embodiment of the present invention is as shown in Figure 3.
It is generally only individually to carry out lithology solution using such as well logging, weight magnetoelectricity or earthquake geophysics means in the prior art
Release, although some have used heavy magnetic shake technology, the utilization of earthquake still based on the profile analysis of subjectivity, does not use more
Add objective seismic properties and cluster analysis seismic facies technology.
And combination of embodiment of the present invention drilling data and well log interpretation, there is provided a kind of skill for identifying sedimentary basin basement lithology
Art method, the identification of basement lithology is realized in the plane, improve the accuracy rate of Lithology Discrimination, to other basins of Eastern China
Basement lithology identification has important reference.
The method of the embodiment of the present invention can give full play to the advantage of various data, make up respective deficiency；In combination with brill
Well data and well log interpretation, make interpretation results have more objectivity, improve the reliability of explanation.
The method and step and beneficial effect of the present invention are described further below by an example.
The example is located at northeast LS areas, collects LS areas substrate data first, gathers and analyzes and appears and drill cores
Sample, statistical analysis LS ground region petrophysical parameter.
By analyzing drill log feature, rock core and thin section identification are observed, it is volcanic rock to determine local area substrate principal lithologic
(including rhyolite, andesite, basalt and tufa stone etc.), intrusive rock (predominantly granite), carbonate rock, epimetamorphic rock
(phyllite, slate, relict sand shale), mesometamorphism rock (schist class) and kata metamorphism rock (gneiss).
To LS Earthquakes reflecting interfaces T3, T4 1And T5(these three numerals are survey region stratigraphic boundaries, are also apparent
Geology reflecting interface, that is, the object of horizon calibration, wherein T3, T4 1It is the boundary line of sedimentary, and T5It is the top of substrate rock mass
Face boundary line) carry out horizon calibration, rift layer carries out seismic interpretation work, emphasized basement top boundary line T5, forms 10*10 nets
Lattice, and to its gridding.13 key sections of emphasized, the Stratigraphic framework in research area is built, it is remaining to establish area below
Density model is prepared.
Shaken and combined by well, studied the characteristics of seismic of each lithology, establish lithology and the pattern of seismic reflection.Development is ground
Study carefully area's based seismic attribute to calculate, mainly using RMS amplitude, waveform clustering and relative impedances attribute.Wherein, the choosing of window when
Select it is very crucial, when the too small validity that target interval can be made to exceed extraction scope, can also influence some parameters of window；When window it is excessive
The off-note of target interval can then be made can not be made a distinction from Anomalies of Backgrounds.The present invention is by comparing, final choice T5-
5ms to T5+50ms this when window.Obtain LS areas substrate RMS amplitude attributed graph, waveform clustering attributed graph and relative impedances
Figure, while by reading the property value of well point position, make seismic properties scatter diagram, establish lithology-seismic properties relation.
Carry out gravity-earthquake peeling to LS areas Bouguer anomaly data to calculate.Superficial part sedimentary is carried out first to calculate extremely
With deep regional field computation.
By seismic data, the main density interface T in the area is extracted3、T4 1And T5, establish Depth Domain residual density model；
According to the physical properties of rock material computation residual density of statistics, result of calculation：T3It is 2.3kg/m to earth's surface residual density3, T3To T4 1
Density is 2.47kg/m3,T4 1To T5Density is 2.57kg/m3,T5Following residual density is 2.64kg/m3。
It is a series of upright hexahedron that length, width and height are respectively △ x, △ y, △ z by STRATIGRAPHIC DIVISION, using six
Face body forward modelling method, calculates T3To earth's surface, T3To T4 1And T4 1To T5The gravity anomaly sum of these three shallow layer stratum sections
(Gq)。
Then carry out deep regional field to calculate extremely, frequency domain continuation computational methods, upward continuation are used to Bouguer anomaly
It is abnormal (Gs) to obtain deep regional field.
Bouguer anomaly is subtracted into superficial part sedimentary exception and deep regional field is abnormal, obtains remaining substrate gravity anomaly
(Gj), i.e. Gj=G0-Gq-Gs, wherein G0For original Bouguer anomaly.
Ask for remaining substrate magnetic anomalies.LS areas T5Above sedimentary magnetic force is weaker, and only flint ridge group development has
The volcanic rock of magnetic, therefore the main exception used magnetic anomalies continuation computational methods, remove superficial part localized magnetization body.Again to prolonging
Open up result and carry out moving average removal deep regional exception, obtain remaining substrate magnetic anomalies.
To substrate gravity anomaly GjWith substrate magnetic anomalies MjAsk for moving average and vertical second dervative calculates, obtain anti-
Reflect gravity anomaly and the magnetic anomalies distributed area of substrate destination layer.
The based seismic attributed graph and substrate gravity-magnetic anomaly block plan of the above-mentioned acquisition of comprehensive analysis, according to basement lithology weight magnetic
Mark sheet is shaken, is explained with reference to drilling well, Lithology Discrimination is deployed to LS areas substrate, finally gives area's basement lithology distribution map, such as
Shown in table 1 and Fig. 4.
Table 1LS areas basement lithology weight magnetic shake mark sheet
According to basement lithology seismic properties feature, gravity anomaly and magnetic anomalies feature, with reference to data such as drilling wells, by pear tree
Rift basement lithology dividing elements are granite phase region, relict sand shale phase region, epimetamorphic rock phase region, mesometamorphism petrofacies area, depth
Metamorphic rock phase region, carbonate rock phase region and volcanic rock phase region.
Granite petrofacies basement lithology shown as in RMS amplitude attribute-and it is weak, value of relative impedance is higher, waveform change
Greatly；And poor continuity.Low-density, low-middle magnetic characteristic, negative G exception, substrate magnetic are shown as in basement lithology gravity anomaly
Shown as in power exception it is low-in positive magnetic anomalies.The surrounding of rift is mainly distributed on, ten rooms south and elm are distributed in inside rift
The Taibei.
Limestone petrofacies show as high RMS amplitude value in basement lithology, and medium value of relative impedance, waveform is complete, continuity
It is good；High density, without magnetic-weak magnetic feature.It is distributed near the room of the north eight and Yang great Cheng of rift.
Volcanic Lithofacies basement lithology shows as RMS amplitude and relative impedances change greatly, it is impossible to determines scope, waveform
Imperfect, poor continuity, local continuous；Midst density, high magnetic characteristic.It is distributed near elm platform, is largely distributed in rift
To the north of northern eight rooms；
Kata metamorphism rock petrofacies basement lithology shows as medium RMS amplitude value and value of relative impedance, and waveform changes greatly, even
In continuous property-poor；High density, high magnetic characteristic.Predominantly gneiss, scope very little, it is distributed in the room of rift ten area.
Weak RMS amplitude value and value of relative impedance during mesometamorphism rock petrofacies basement lithology is shown as, waveform change greatly, even
Continuous property is poor；High density, in-high density magnetic.Be distributed in mulberry tree platform nearby and small city be southwestern, including quartzite, phyllite and
Proterozoic erathem schist etc..
Epimetamorphic rock petrofacies basement lithology weak RMS amplitude value and value of relative impedance in showing as, waveform change greatly, even
Continuous property is poor；High density, middle magnetic characteristic.For the main body of basement lithology in rift, large area is distributed in rift, pottery predominantly difficult to understand
The schist of system.
Rotten glutenite falls within epimetamorphic rock, but it individually is attributed into one kind here.Its petrofacies basement lithology is shown as
Low-middle RMS amplitude value, relative impedances value changes are big；Low-density, nonmagnetic feature.Big elm is mainly distributed on to Qin Jiatun
Between.
Tested using known drilling well, accuracy is more than 90%.By the research of the present invention, substrate is substantially increased
The efficiency of Lithology Discrimination, the distribution of lithology can be accurately determined, further clear and definite buried hills Exploration targets.Pair plus
The paces of fast oil-production capacity construction have important practical significance.The achievement in research of the technology will to other basins of Eastern China and
The identification of block basement lithology has important reference function.
Fig. 5 is the structure according to the basement lithology identifying system that Comprehensive interpretation is shaken based on weight magnetic of the embodiment of the present invention
Schematic diagram, as illustrated, the identifying system includes：
Seismic properties distribution map establishes module 51, its execution step S110 operation, residual gravity anomaly acquisition module 52,
It performs step S120 operation, and residual magnetism exception acquisition module 53, it performs step S130 operation, gravity-magnetic anomaly point
Area's figure establishes module 54, and it performs step S140 operation, overlaps identification module 55, and it performs step S150 operation.Herein
Repeat no more.
Those skilled in the art should be understood that above-mentioned each module of the invention or each step can use general calculating
Device realizes that they can be concentrated on single computing device, or be distributed in the network that multiple computing devices are formed
On, alternatively, they can be realized with the program code that computing device can perform, it is thus possible to be stored in storage
Performed in device by computing device, they are either fabricated to each integrated circuit modules respectively or will be more in them
Individual module or step are fabricated to single integrated circuit module to realize.So, the present invention be not restricted to any specific hardware and
Software combines.
Although disclosed herein embodiment as above, described content only to facilitate understand the present invention and adopt
Embodiment, it is not limited to the present invention.Any those skilled in the art to which this invention pertains, this is not being departed from
On the premise of the disclosed spirit and scope of invention, any modification and change can be made in the implementing form and in details,
But the scope of patent protection of the present invention, still should be subject to the scope of the claims as defined in the appended claims.
Claims (10)
1. a kind of basement lithology recognition methods based on weight magnetic shake Comprehensive interpretation, including：
The Seismic reflection character of Study of recognition area basement lithology, lithology and the pattern of seismic reflection are established according to recognition result, and
Carry out attribute to seismic data to calculate, the result calculated according to attribute generates seismic properties distribution map；
Gravity-earthquake peeling is carried out to Bouguer anomaly data to calculate to obtain the residual gravity anomaly of survey region substrate；
Magnetic anomalies continuation calculating is carried out to magnetic anomalies data to obtain the residual magnetism of survey region substrate exception；
Gravity-magnetic anomaly block plan is established according to the residual gravity anomaly and residual magnetism extremely；
Overlapping research is carried out based on the seismic properties distribution map and gravity-magnetic anomaly block plan, and combines drilling well and explains, to research
The basement lithology in region is identified.
2. according to the method for claim 1, it is characterised in that include carrying out to seismic data before attribute calculates：
According to the substrate data in research area of collection and appearing and drill cores sample for collection, it is determined that needing to carry out Lithology Discrimination
Lithology classification；
The Seismic reflection character of each reflecting interface is determined, horizon calibration is carried out to the important seismic interface in research area；
Seismic interpretation is carried out according to the result of horizon calibration, the seismic interpretation includes the explanation of layer position and fault interpretation；
Gridding processing is carried out to the seismic interpretation result in basement top boundary line, the result based on gridding processing is belonged to
Property calculate.
3. according to the method for claim 2, it is characterised in that the attribute calculate include the calculating of RMS amplitude attribute,
Waveform clustering attribute calculates and Relative Wave Impedance attribute calculates.
4. according to the method in any one of claims 1 to 3, it is characterised in that described that lithology is established according to recognition result
With the pattern of seismic reflection, including：
Based on the Seismic reflection character, lithology and the pattern of seismic reflection are established by single well analysis and well shake contrast, wherein,
The Seismic reflection character refer to seismic reflection border, inside and out feature.
5. according to the method for claim 1, it is characterised in that described that gravity-earthquake peeling is carried out to Bouguer anomaly data
Calculate to obtain the residual gravity anomaly of survey region substrate, including：
Superficial part sedimentary is carried out to Bouguer anomaly data to calculate extremely to obtain the weight of survey region substrate rock mass above sedimentary
Power is abnormal；
Deep regional field is carried out to Bouguer anomaly data to calculate extremely to obtain the region in deep below survey region substrate rock mass
Field is abnormal；
The superficial part sedimentary exception and deep regional field exception are subtracted with Bouguer anomaly data to obtain survey region substrate
Residual gravity anomaly.
6. according to the method for claim 5, it is characterised in that described that superficial part sedimentary exception is carried out to Bouguer anomaly data
Calculate, including：
Seismic data is explained, establishes Stratigraphic framework, to extract main density interface；
Sedimentary density model is established based on the main density interface；
Forward modelling is carried out based on the sedimentary density model, to obtain more than the substrate rock mass caused by shallow-layer density interface
The gravity anomaly of sedimentary.
7. the method according to claim 5 or 6, it is characterised in that described that deep regional field is carried out to Bouguer anomaly data
It is abnormal to calculate, including：
Gravity anomaly upward continuation to being observed from Bouguer anomaly data, obtains deep regional field；
The frequency spectrum of the gravity anomaly is calculated by wavelet analysis；
Deep diopter filtering is carried out to the frequency spectrum, it is abnormal that deep regional field is chosen from filter result.
8. according to the method for claim 1, it is characterised in that described that magnetic anomalies continuation meter is carried out to magnetic anomalies data
Calculate to obtain the residual magnetism of survey region substrate exception, including：
Magnetic anomalies upward continuation to being observed from magnetic anomalies data, the magnetic anomalies are converted to more than inspection surface
Space；
The frequency spectrum of the magnetic anomalies is calculated by wavelet analysis；
Remove the superficial part magnetic anomalies in the frequency spectrum；
Moving average method is used to be filtered the frequency spectrum to remove the deep magnetic anomalies in the frequency spectrum.
9. according to the method for claim 1, it is characterised in that described different according to the residual gravity anomaly and residual magnetism
Gravity-magnetic anomaly block plan often is established, including：
The residual gravity anomaly and residual magnetism are filtered extremely using moving average method；
Vertical second dervative calculating is carried out extremely to the residual gravity anomaly and residual magnetism to obtain gravity-magnetic anomaly block plan.
10. a kind of basement lithology identifying system based on weight magnetic shake Comprehensive interpretation, including：
Seismic properties distribution map establishes module, the Seismic reflection character of its Study of recognition area basement lithology, is built according to recognition result
Vertical lithology and the pattern of seismic reflection, and attribute is carried out to seismic data and calculated, the result calculated according to attribute generates earthquake category
Property distribution map；
Residual gravity anomaly acquisition module, it carries out gravity-earthquake peeling to Bouguer anomaly data and calculated to obtain survey region
The residual gravity anomaly of substrate；
Residual magnetism exception acquisition module, it carries out magnetic anomalies continuation calculating to obtain survey region base to magnetic anomalies data
The residual magnetism at bottom is abnormal；
Gravity-magnetic anomaly block plan establishes module, and it establishes gravity-magnetic anomaly point extremely according to the residual gravity anomaly and residual magnetism
Qu Tu；
Identification module is overlapped, it is based on the seismic properties distribution map and gravity-magnetic anomaly block plan carries out overlapping research, and combines
Drilling well is explained, the basement lithology of survey region is identified.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610824523.7A CN107817519A (en) | 2016-09-14 | 2016-09-14 | Basement lithology recognition methods and system based on weight magnetic shake Comprehensive interpretation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610824523.7A CN107817519A (en) | 2016-09-14 | 2016-09-14 | Basement lithology recognition methods and system based on weight magnetic shake Comprehensive interpretation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107817519A true CN107817519A (en) | 2018-03-20 |
Family
ID=61601342
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610824523.7A Pending CN107817519A (en) | 2016-09-14 | 2016-09-14 | Basement lithology recognition methods and system based on weight magnetic shake Comprehensive interpretation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107817519A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111323811A (en) * | 2018-12-13 | 2020-06-23 | 中国石油天然气集团有限公司 | Multi-well seismic data imaging method and system |
| CN112363246A (en) * | 2020-10-29 | 2021-02-12 | 中国石油天然气集团有限公司 | Igneous rock lithology identification method and device |
| CN112485841A (en) * | 2019-09-11 | 2021-03-12 | 中国石油天然气集团有限公司 | Deep stratum lithology identification method and device |
| CN113406713A (en) * | 2021-03-31 | 2021-09-17 | 北京中科地物能源技术有限公司 | Method for evaluating residual gas development potential of developed well of microgravity monitoring gas reservoir |
| CN114048636A (en) * | 2022-01-11 | 2022-02-15 | 中国测绘科学研究院 | Gravity anomaly calculation method and device based on wavelet transformation |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6502037B1 (en) * | 1999-04-02 | 2002-12-31 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
| CN101520518A (en) * | 2008-02-25 | 2009-09-02 | 中国石油集团东方地球物理勘探有限责任公司 | Method for recognizing the lithology of petrosilex by using the combined characteristics of gravity-magnetic-electronic anomaly |
| CN103675944A (en) * | 2013-11-27 | 2014-03-26 | 核工业北京地质研究院 | Comprehensive geophysical method for sedimentary basin uranium resource exploration |
| CN104635261A (en) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | Gravity, magnetism, electrical method and earthquake combined split modeling method for mountain fronts |
| CN105005097A (en) * | 2015-07-17 | 2015-10-28 | 中国石油化工股份有限公司 | Method for comprehensive recognition of igneous rocks by employing gravity, magnetism, electromagnetism, and earthquake data |
-
2016
- 2016-09-14 CN CN201610824523.7A patent/CN107817519A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6502037B1 (en) * | 1999-04-02 | 2002-12-31 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
| CN101520518A (en) * | 2008-02-25 | 2009-09-02 | 中国石油集团东方地球物理勘探有限责任公司 | Method for recognizing the lithology of petrosilex by using the combined characteristics of gravity-magnetic-electronic anomaly |
| CN104635261A (en) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | Gravity, magnetism, electrical method and earthquake combined split modeling method for mountain fronts |
| CN103675944A (en) * | 2013-11-27 | 2014-03-26 | 核工业北京地质研究院 | Comprehensive geophysical method for sedimentary basin uranium resource exploration |
| CN105005097A (en) * | 2015-07-17 | 2015-10-28 | 中国石油化工股份有限公司 | Method for comprehensive recognition of igneous rocks by employing gravity, magnetism, electromagnetism, and earthquake data |
Non-Patent Citations (4)
| Title |
|---|
| 孟令顺 等: "《勘探重力学与地磁学》", 31 December 2008, 地质出版社 * |
| 张元高: "松辽盆地北部古中央隆起带深层构造研究及基底岩性预测", 《万方》 * |
| 王谦身 等: "《重力学》", 30 April 2003, 地震出版社 * |
| 鲁宝亮 等: "南海北部盆地基底岩性地震-重磁响应特征与识别", 《地球物理学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111323811A (en) * | 2018-12-13 | 2020-06-23 | 中国石油天然气集团有限公司 | Multi-well seismic data imaging method and system |
| CN112485841A (en) * | 2019-09-11 | 2021-03-12 | 中国石油天然气集团有限公司 | Deep stratum lithology identification method and device |
| CN112363246A (en) * | 2020-10-29 | 2021-02-12 | 中国石油天然气集团有限公司 | Igneous rock lithology identification method and device |
| CN113406713A (en) * | 2021-03-31 | 2021-09-17 | 北京中科地物能源技术有限公司 | Method for evaluating residual gas development potential of developed well of microgravity monitoring gas reservoir |
| CN113406713B (en) * | 2021-03-31 | 2024-03-22 | 北京中科地物能源技术有限公司 | Method for evaluating residual gas development potential of microgravity monitoring gas reservoir developed well |
| CN114048636A (en) * | 2022-01-11 | 2022-02-15 | 中国测绘科学研究院 | Gravity anomaly calculation method and device based on wavelet transformation |
| CN114048636B (en) * | 2022-01-11 | 2022-05-03 | 中国测绘科学研究院 | Gravity anomaly calculation method and device based on wavelet transformation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105510993B (en) | The buried extrusion pressing type complexity gypseous-salt rock identification of foreland basin and distribution forecasting method | |
| CN107817519A (en) | Basement lithology recognition methods and system based on weight magnetic shake Comprehensive interpretation | |
| CN105675635B (en) | Tight rock component relative amount and brittleness index determine method and apparatus | |
| CN107894615B (en) | A kind of method of quantitative evaluation 3-D seismics attribute forecast reservoir parameter validity | |
| CN108931814A (en) | A method of the basement rock FRACTURE PREDICTION based on the fusion of more attributes | |
| CN106772579B (en) | Earthquake prestack inversion method and apparatus in a kind of girdle | |
| Caers et al. | Stochastic integration of seismic data and geologic scenarios: A West Africa submarine channel saga | |
| CN105044770A (en) | Compact glutenite gas reservoir quantificational prediction method | |
| CN107748399B (en) | Method for identifying deep tectonic layer of mountain front zone by utilizing gravity interface inversion | |
| CN110231652A (en) | A kind of seismic facies extracting method of the Noise application space cluster based on density | |
| CN109655894A (en) | Carbonate rock ancient stream channel seismic inversion low frequency model construction method and system | |
| CN105842733A (en) | Shale reservoir earthquake identification method | |
| CN106033127B (en) | Crustal stress azimuthal seismic Forecasting Methodology based on shear wave velocity rate of change | |
| CN106371139B (en) | A kind of method and device using interval amplitude attribute identification lithological pool | |
| Thachaparambil | Discrete 3D fracture network extraction and characterization from 3D seismic data—A case study at Teapot Dome | |
| Alizadeh et al. | Using stacked generalization ensemble method to estimate shear wave velocity based on downhole seismic data: A case study of Sarab-e-Zahab, Iran | |
| Ran et al. | Volcanic gas reservoir characterization | |
| Méndez et al. | Rock type prediction and 3D modeling of clastic paleokarst fillings in deeply-buried carbonates using the Democratic Neural Networks Association technique | |
| Hickman | Parasequence-scale stratigraphic variability of lithology and porosity of Mississippian Meramec reservoirs and the relationships to production characteristics, STACK trend, Oklahoma | |
| CN106990449A (en) | A kind of terrestrial face interlayer screen work based on multidimensional data builds characterizing method | |
| Zhang et al. | Hierarchical modeling of carbonate fault-controlled Paleokarst systems: A case study of the Ordovician reservoir in the Tahe Oilfield, Tarim Basin | |
| CN114185083A (en) | Quantitative evaluation method for fault sealing in clastic rock stratum | |
| Méndez et al. | Predicting and 3D modeling of karst zones using seismic facies analysis in Ordovician carbonates of the Tahe oilfield, China | |
| Erzeybek Balan | Characterization and modeling of paleokarst reservoirs using multiple-point statistics on a non-gridded basis | |
| CN113514884A (en) | Compact sandstone reservoir prediction 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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180320 |
|
| RJ01 | Rejection of invention patent application after publication |