CN106842318A - Method and device for determining geophysical two-dimensional representation of microscopic erosion pore - Google Patents
Method and device for determining geophysical two-dimensional representation of microscopic erosion pore Download PDFInfo
- Publication number
- CN106842318A CN106842318A CN201611261876.7A CN201611261876A CN106842318A CN 106842318 A CN106842318 A CN 106842318A CN 201611261876 A CN201611261876 A CN 201611261876A CN 106842318 A CN106842318 A CN 106842318A
- Authority
- CN
- China
- Prior art keywords
- target area
- microcosmic erosion
- microcosmic
- erosion hole
- substratum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000003628 erosive effect Effects 0.000 title claims abstract description 241
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000011148 porous material Substances 0.000 title claims abstract description 45
- 239000011435 rock Substances 0.000 claims abstract description 51
- 238000009826 distribution Methods 0.000 claims abstract description 47
- 230000004044 response Effects 0.000 claims abstract description 42
- 238000005266 casting Methods 0.000 claims abstract description 16
- 238000011161 development Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims description 33
- 238000005481 NMR spectroscopy Methods 0.000 claims description 25
- 230000005311 nuclear magnetism Effects 0.000 claims description 22
- 238000012512 characterization method Methods 0.000 claims description 20
- 238000002389 environmental scanning electron microscopy Methods 0.000 claims description 13
- 238000004090 dissolution Methods 0.000 claims description 12
- 239000004744 fabric Substances 0.000 claims description 2
- 238000007796 conventional method Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 6
- 239000000463 material Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000000047 product Substances 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 238000011017 operating method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000009738 saturating Methods 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
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 embodiment of the application provides a method and a device for determining geophysical two-dimensional representation of microscopic erosion pores. According to the method, a microscopic erosion pore development area is determined through a casting body slice and a scanning electron microscope, then a logging response characteristic of the microscopic erosion pore is researched by combining a logging curve and a nuclear magnetic logging porosity curve, and then a rock physical quantity version is determined; and secondly, obtaining a micro erosion pore logging characteristic curve according to the rock physical gauge plate, and further determining the distribution characteristics of the micro erosion pores through waveform indication inversion. Therefore, the technical problem that the distribution characteristics of the microscopic erosion pores cannot be accurately determined in the conventional method is solved.
Description
Technical field
The application is related to petroleum exploration domain, more particularly to a kind of microcosmic erosion hole geophysics two-dimensional characterization determination side
Method and device.
Background technology
With the continuous improvement of oil-gas exploration technology, find to be visited in deep layer and hyposmosis-ultra-low penetration Sandstone Gas Reservoir
Bright oil gas rate more and more higher, therefore be increasingly taken seriously.But, the difficulties in exploration of this kind of reservoir is also very big.Wherein, it is above-mentioned
Hyposmosis-extra-low permeability Sandstone Gas Reservoir typically refers to meet porosity less than 10%, and sandstone of the permeability less than 1md is stored up
Layer.
Specifically, deep layer and hyposmosis-extra-low permeability Sandstone Gas Reservoir typically have larger resource potential, while
With poor storage and collection performance.Research finds that the development of secondary microcosmic erosion hole, can improve deep layer and hyposmosis-spy is hypotonic
The storage and collection performance of saturating Sandstone Gas Reservoir.But due to being influenceed by seismic data quality, it is currently based on poststack two-dimension earthquake money
Material goes the research of the development characteristics for characterizing microcosmic erosion hole little.Using existing method, found on well-log information secondary
The logging response character of microcosmic erosion hole is just highly difficult, and the earthquake of microcosmic erosion hole is found on poststack 2-d seismic data
Response characteristic is just more difficult.Therefore, when existing method is embodied, often exist can not accurately determine it is micro- in target area
See the technical problem of the distribution characteristics of dissolution porosity.
Regarding to the issue above, effective solution is not yet proposed at present.
The content of the invention
The embodiment of the present application determines method and apparatus there is provided a kind of microcosmic erosion hole geophysics two-dimensional characterization, to solve
Certainly existing microcosmic erosion hole determine method exist can not accurately determine that microcosmic erosion hole is distributed spy in two dimensional cross-section
The technical problem levied.
The embodiment of the present application determines method there is provided a kind of microcosmic erosion hole geophysics two-dimensional characterization, including：
Obtain sample, the log data of the target area, the target area of section sandstone of coring of target area
The nuclear-magnetism porosity logs in domain and the two-dimentional poststack seismic data of the target area；
Casting body flake and ESEM experimental analysis are carried out to the sample, analysis result is obtained；
According to the analysis result and the log data, by core Location, the microcosmic molten of the sample is determined
Pit gap development area；
Nuclear-magnetism porosity logs of log data, the target area according to the target area and described
The microcosmic erosion pores'growth area of sample, determines the logging response character of the microcosmic erosion hole of the target area；
The target area is divided into multiple continuous substratums according to depth, and according to the well logging of the microcosmic erosion hole
Response characteristic, is that each substratum in the multiple continuous substratum sets up rock physicses amount version respectively；
According to the rock physicses amount version of each substratum in multiple continuous substratums, determine that microcosmic erosion hole logging character is bent
Line；
The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, according to the microcosmic erosion hole
Logging character curve, inverting is indicated by seismic waveform, obtains the distribution characteristics of the microcosmic erosion hole of the target area.
In one embodiment, the logging response character of the microcosmic erosion hole includes：Velocity of longitudinal wave is more than default
Speed, and/or, nuclear magnetic resonance log porosity is less than default porosity.
In one embodiment, according to the target area log data, the nuclear-magnetism hole of the target area
The microcosmic erosion pores'growth area of porosity well-log information and the sample, determines the survey of the microcosmic erosion hole of the target area
Well response characteristic, including：
The interval transit time curve in log data according to the target area, determines the compressional wave of the target area
Speed；
Nuclear-magnetism porosity logs according to the target area, determine the nuclear magnetic resonance log hole of the target area
Degree；
Microcosmic erosion pores'growth area according to the sample, determines the pre-set velocity and the default porosity.
In one embodiment, it is the multiple continuous according to the logging response character of the microcosmic erosion hole
Each substratum in substratum sets up rock physicses amount version respectively, including：
According to the logging response character of the microcosmic erosion hole, to each substratum in the multiple continuous substratum,
The nuclear magnetic resonance log porosity with place substratum, with velocity of longitudinal wave as ordinate, sets up the rock of each substratum as abscissa respectively
Physical quantity version.
In one embodiment, according to the rock physicses amount version of each substratum in multiple continuous substratums, determine microcosmic
Dissolution porosity logging character curve, including：
From the multiple continuous substratum microcosmic erosion hole is obtained in the rock physicses amount version of each substratum by substratum
Information；
The information of the microcosmic erosion hole of each substratum that will be obtained, it is bent as the microcosmic erosion hole logging character
Line.
In one embodiment, the data volume of inverting, root are indicated as seismic waveform with the two-dimentional poststack seismic data
According to the microcosmic erosion hole logging character curve, inverting is indicated by seismic waveform, obtain the microcosmic molten of the target area
The distribution characteristics of pit gap, including：
The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, and with the microcosmic erosion hole
Logging character curve as constraints, preferably radio-frequency component, to set up inverse model；
According to the inverse model, carry out seismic waveform and indicate inverting, obtain the microcosmic erosion hole of the target area
Distribution characteristics.
In one embodiment, the data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, and
Using the microcosmic erosion hole logging character curve as constraints, preferably radio-frequency component, to set up inverse model, including：
The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, and with the microcosmic erosion hole
Logging character curve as constraints, preferably radio-frequency component, to set up initial model；
According to the initial model, inverting obtains the target area and crosses microcosmic erosion in the two-dimentional stacked seismic data of well
Interstitial space distribution characteristics；
Microcosmic erosion interstitial space distribution characteristics and the target area in well two dimension stacked seismic data are crossed by described
Microcosmic erosion hole logging character curve on well is compared, and obtains relative error；
If the relative error is less than predetermined threshold value, using the initial model as the inverse model.
In one embodiment, microcosmic erosion interstitial space distribution spy in well two dimension stacked seismic data is being crossed by described
The microcosmic erosion hole logging character curve levied on the well with the target area is compared, after obtaining relative error, institute
Stating method also includes：
If the relative error is more than or equal to the predetermined threshold value, radio-frequency component is reselected, and according to selecting again
The radio-frequency component selected sets up initial model, until according to the initial model set up, the well excessively of the target area that inverting is obtained
Microcosmic erosion hole on microcosmic erosion interstitial space distribution characteristics in two-dimentional stacked seismic data and the well of the target area
Relative error between gap logging character curve is less than the predetermined threshold value, then using current initial model as the inverting mould
Type.
Based on identical inventive concept, the embodiment of the present application additionally provides a kind of microcosmic erosion hole geophysics bivariate table
Determining device is levied, including：
Acquisition module, the sample of section sandstone of coring for obtaining target area, the log money of the target area
The two-dimentional poststack seismic data of material, the nuclear-magnetism porosity logs of the target area and the target area；
Analysis module, for carrying out casting body flake and ESEM experimental analysis to the sample, obtains analysis result；
First determining module, for according to the analysis result and the log data, by core Location, it is determined that
The microcosmic erosion pores'growth area of the sample；
Second determining module, for according to the analysis result and the log data, by core Location, it is determined that
The microcosmic erosion pores'growth area of the sample；
3rd determining module, for the target area to be divided into multiple continuous substratums according to depth, and according to described
The logging response character of microcosmic erosion hole, is that each substratum in the multiple continuous substratum sets up rock physicses amount respectively
Version；
4th determining module, for the rock physicses amount version according to each substratum in multiple continuous substratums, determines microcosmic
Dissolution porosity logging character curve；
5th determining module, the data volume for indicating inverting as seismic waveform with the two-dimentional poststack seismic data, root
According to the microcosmic erosion hole logging character curve, inverting is indicated by seismic waveform, obtain the microcosmic molten of the target area
The distribution characteristics of pit gap.
In one embodiment, second determining module includes：
First determining unit, for the interval transit time curve in the log data according to the target area, it is determined that
The velocity of longitudinal wave of the target area；
Second determining unit, for the nuclear-magnetism porosity logs according to the target area, determines the target area
The nuclear magnetic resonance log porosity in domain；
3rd determining unit, for the microcosmic erosion pores'growth area according to the sample, determine the pre-set velocity and
The default porosity.
In one embodiment, the 4th determining module includes：
Information acquisition unit, for each substratum from the multiple continuous substratum rock physicses amount version in by substratum
Obtain the information of microcosmic erosion hole；
4th determining unit, for will obtain each substratum microcosmic erosion hole information, as described microcosmic molten
Pit gap logging character curve.
In one embodiment, the 5th determining module includes：
Model determining unit, the data volume for indicating inverting as seismic waveform with the two-dimentional poststack seismic data, and
Using the microcosmic erosion hole logging character curve as constraints, preferably radio-frequency component, to set up inverse model；
5th determining unit, inverting is indicated for according to the inverse model, carrying out seismic waveform, obtains the target area
The distribution characteristics of the microcosmic erosion hole in domain.
In the embodiment of the present application, determine microcosmic erosion pores'growth area by using casting body flake and ESEM, tie
Close log and nuclear magnetic resonance log porosity curve obtains the logging response character of microcosmic erosion hole, with reference to the rock thing set up
Reason template and then the determination continuous logging character curve of microcosmic erosion hole；Further according to two-dimentional poststack seismic data and microcosmic erosion
Hole logging character curve negotiating seismic waveform indicates inverting to determine the distribution characteristics of microcosmic erosion hole, expands the model of prediction
Enclose, solve the problems, such as that the use of existing method presence is limited in scope and determine precision low technical with microcosmic erosion hole, reached
The accurate technique effect for determining microcosmic erosion distribution of pores in two-dimentional stacked seismic data.
Brief description of the drawings
In order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, are not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is the treatment that method is determined according to a kind of microcosmic erosion hole geophysics two-dimensional characterization of the embodiment of the present application
Flow chart；
Fig. 2 is a kind of composition of the microcosmic erosion hole geophysics two-dimensional characterization determining device according to the embodiment of the present application
Structure chart；
Fig. 3 is to provide a kind of microcosmic erosion hole geophysics two-dimensional characterization using the embodiment of the present application to determine method/dress
Put W1 borehole logging tool microcosmic erosion pores'growth characteristic synthetic block diagrams in the target area for obtaining；
Fig. 4 is to provide a kind of microcosmic erosion hole geophysics two-dimensional characterization using the embodiment of the present application to determine method/dress
Put the schematic diagram of the Sandstone Gas Reservoir microcosmic erosion porous rocks physics differentiation amount version for obtaining；
Fig. 5 is to provide a kind of microcosmic erosion hole geophysics two-dimensional characterization using the embodiment of the present application to determine method/dress
Put the schematic diagram for crossing W1 wells Crossline (road) direction poststack 2-d seismic data for obtaining；
Fig. 6 is to provide a kind of microcosmic erosion hole geophysics two-dimensional characterization using the embodiment of the present application to determine method/dress
Put the schematic diagram that W1 well Crossline directions microcosmic erosion hole seismic waveform indicates inverting section of crossing for obtaining.
Specific embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, below in conjunction with the application reality
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example is only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
In view of at present due to being influenceed by seismic data quality, sign microcosmic erosion is gone based on poststack 2-d seismic data
The research of the development characteristics of hole is little.The characteristics of due to microcosmic erosion hole itself, i.e.,：Microcosmic erosion pore scale first is non-
It is often small, typically it is difficult to find the logging response character of microcosmic erosion hole on the well curve routinely surveyed；Secondly as microcosmic molten
Erosion Pore Genesis are more, are difficult to find its seismic response in conventional two dimension stacked seismic data.Therefore, using existing side
Method, finds that the logging response character of secondary microcosmic erosion hole is highly difficult on well-log information, and in poststack 2-d seismic data
The upper seismic response features for finding microcosmic erosion hole are just more difficult.Therefore, existing method in the specific implementation, is often deposited
The technical problem of the distribution characteristics of the microcosmic erosion hole in it accurately can not determine target area.Asked for above-mentioned technology is produced
The basic reason of topic, the application considers that casting body flake and ESEM can be first passed through, and determines on the basis of core Location micro-
See dissolution porosity development area；Secondly log and nuclear magnetic resonance log porosity curve are combined, the well logging of microcosmic erosion hole is analyzed
Response characteristic, with reference to rock physicses template, and then obtains microcosmic erosion hole logging character curve；Finally according to microcosmic erosion hole
Gap logging character curve and two-dimentional poststack seismic data, inverting is indicated by seismic waveform, determines the distribution of microcosmic erosion hole
Feature.So as to solve existing microcosmic erosion hole determine method exist can not accurately determine it is microcosmic molten in target area
The technical problem of distribution of pores is lost, the accurate technique effect for determining microcosmic erosion distribution of pores in two dimensional cross-section has been reached.
Based on above-mentioned thinking thinking, this application provides a kind of microcosmic erosion hole geophysics two-dimensional characterization determination side
Method.Refer to Fig. 1.A kind of microcosmic erosion hole geophysics two-dimensional characterization that the application is provided determines method, can specifically wrap
Include following steps (step 101 to step 107).
Step 101：Obtain sample, the log data of the target area, the institute of section sandstone of coring of target area
State the nuclear-magnetism porosity logs of target area and the two-dimentional poststack seismic data of the target area.
In one embodiment, general prospect pit all has log data, and nuclear-magnetism porosity logs can be with
There is what is be directed to measure, above-mentioned log data and nuclear-magnetism porosity logs can often be directly obtained.Target area takes
The sample of heart section sandstone, typically will be sampled acquisition to section sample of sandstone of coring.
Step 102：Casting body flake and ESEM experimental analysis are carried out to the sample, analysis result is obtained.
Have in an implementation method, in order to subsequently accurately determine the microcosmic erosion pores'growth area of the sample, examine again
The logging response character for considering microcosmic erosion hole is not generally relatively obvious, can first determine the analysis result of sample.I.e., it is possible to
First pass through at the same carry out casting body flake to sample and scanning electron microscope analysis obtain sample casting body flake and ESEM analysis
As a result, subsequently to use.
Step 103：According to the analysis result and the log data, by core Location, the sample is determined
Microcosmic erosion pores'growth area.
In one embodiment, in order to accurately determine the microcosmic erosion pores'growth area of sample, specifically can be according to upper
Analysis result and well-log information are stated, core Location is carried out；On the basis of core Location, the microcosmic erosion hole hair of sample is determined
Area is educated, so as to the follow-up logging response character for determining microcosmic erosion hole.During specific implementation, can according to log data and
The rock core depth of well logging and section of coring carries out core Location jointly.
Step 104：The nuclear-magnetism porosity logging money of log data, the target area according to the target area
The microcosmic erosion pores'growth area of material and the sample, determines that the log response of the microcosmic erosion hole of the target area is special
Levy.
In one embodiment, by analyzing the microcosmic erosion hole in microcosmic erosion pores'growth area, it may be determined that
The logging response character of the microcosmic erosion hole can specifically include：Velocity of longitudinal wave is relatively large, nuclear magnetic resonance log porosity phase
To smaller.I.e. velocity of longitudinal wave is more than pre-set velocity, and/or, nuclear magnetic resonance log porosity is less than default porosity.
In one embodiment, in order to determine the logging response character of the microcosmic erosion hole of target area, specifically may be used
Performed (S1 to S3) with according to following steps.
S1：The interval transit time curve in log data according to the target area, determines the target area
Velocity of longitudinal wave.
S2：Nuclear-magnetism porosity logs according to the target area, determine the nuclear magnetic resonance log hole of the target area
Porosity.
S3：Microcosmic erosion pores'growth area according to the sample, determines the pre-set velocity and the default porosity.
And then can be according to the velocity of longitudinal wave of above-mentioned acquisition, nuclear magnetic resonance log porosity, pre-set velocity and the default hole
Whether degree, determine whether each velocity of longitudinal wave and nuclear magnetic resonance log porosity meet the logging response character of microcosmic erosion hole, i.e.,
Velocity of longitudinal wave is met more than pre-set velocity, and/or, nuclear magnetic resonance log porosity is less than default porosity.
Step 105：The target area is divided into multiple continuous substratums according to depth, and according to the microcosmic erosion hole
The logging response character of gap, is that each substratum in the multiple continuous substratum sets up rock physicses amount version respectively.
In one embodiment, in order to ensure the precision of result and subsequently determine microcosmic erosion hole logging character curve
Accuracy and continuity, target area can be divided into multiple substratums by depth according to default depth, carry out respectively
Treatment.For example, the formation thickness of target area is 100m, can be according to the preset thickness of every substratum 10m, by the stratum average mark
Into 10 substratums.It is of course also possible to as the case may be, the stratum after 100m is divided into 5 substratums of not equal thickness, wherein,
First substratum 10m, the second substratum 20m, the 3rd substratum 20m, the 4th substratum 40, the 5th substratum 10m.As long as ensureing follow-up acquisition
The rock physicses amount version of each substratum is continuous accurate.In this regard, the application is not construed as limiting.
It should be noted that subsequently obtain rock physicses amount version as far as possible accurately, continuously in order that obtaining, it is embodied
When, can also as the case may be or implementation condition, by the several substratums in many points of target area, to cause the thickness of each substratum
It is small as far as possible, you can to reduce with by the follow-up depth bounds for obtaining rock physicses amount version, so can further reduce follow-up
The error of the microcosmic erosion hole logging character curve of determination.
In one embodiment, in order to determine the rock physicses amount version of each substratum, during specific implementation, can be with：According to
The logging response character of the microcosmic erosion hole is small with place respectively to each substratum in the multiple continuous substratum
The nuclear magnetic resonance log porosity of layer is abscissa, with velocity of longitudinal wave as ordinate, sets up the rock physicses amount version of each substratum.
Step 106：According to the rock physicses amount version of each substratum in multiple continuous substratums, determine that microcosmic erosion hole is surveyed
Well indicatrix.
In one embodiment, in order to determine microcosmic erosion hole logging character curve, specifically can be according to following step
It is rapid to perform (S1 to S2).
S1：From the multiple continuous substratum microcosmic erosion hole is obtained in the rock physicses amount version of each substratum by substratum
The information of gap.
S2：The information of the microcosmic erosion hole of each substratum that will be obtained, as the microcosmic erosion hole logging character
Curve.
In one embodiment, in order to obtain each rock physicses amount version in microcosmic erosion hole information, specifically
Can be microcosmic erosion pores'growth region by the analysis upper left corner area in each rock physicses amount version that falls during implementation,
Obtain the information of the microcosmic erosion hole in rock physicses amount version.Because, it is generally, quick due to microcosmic erosion hole
Sense curvilinear characteristic, so in the ideal situation, the microcosmic erosion hole information of rock physicses amount version can mainly be distributed in the rock
The upper left corner of physical quantity version.Certainly, different research blocks are due to the difference of depositional environment, diagenesis etc., microcosmic erosion hole
Data differ to establish a capital and fully fall in the upper left corner, during specific implementation, can be adjusted accordingly according to actual conditions.
Step 107：The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, according to described microcosmic
Dissolution porosity logging character curve, by seismic waveform indicate inverting, obtain the target area microcosmic erosion hole point
Cloth feature.
In one embodiment, in order to accurately determine the distribution characteristics of the microcosmic erosion hole of target area, specifically
During implementation, can be performed according to following steps (S1 to S2).
S1：The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, and with the microcosmic erosion
Hole logging character curve as constraints, preferably radio-frequency component, to set up inverse model.
S2：According to the inverse model, carry out seismic waveform and indicate inverting, obtain the microcosmic erosion hole of the target area
The distribution characteristics of gap.
In one embodiment, in order to set up accurate inverse model, can be according to following steps (S1 during specific implementation
To S4) implement.
S1：The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, and with the microcosmic erosion
Hole logging character curve as constraints, preferably radio-frequency component, to set up initial model, that is, find microcosmic erosion hole with
Seismic waveform relation.
S2：According to the initial model, inverting obtains the target area and crosses microcosmic in the two-dimentional stacked seismic data of well
Dissolution porosity spatial distribution characteristic.
S3：Microcosmic erosion interstitial space distribution characteristics and the target area in well two dimension stacked seismic data are crossed by described
Well on microcosmic erosion hole logging character curve be compared, obtain relative error.
S4：If the relative error is less than predetermined threshold value, using the initial model as the inverse model.
In one embodiment, in order to be corrected more than or equal to the initial model of predetermined threshold value to relative error.
Cross micro- on well of the microcosmic erosion interstitial space distribution characteristics with the target area in well two dimension stacked seismic data by described
See dissolution porosity logging character curve to be compared, obtain after relative error, methods described can also include：If the ratio
It is more than or equal to the predetermined threshold value compared with error, then reselects radio-frequency component, and set up just according to the radio-frequency component for reselecting
Beginning model, until on the initial model according to foundation, the two dimension stacked seismic data of well excessively of the target area that inverting is obtained
Microcosmic erosion interstitial space distribution characteristics and the target area well on microcosmic erosion hole logging character curve between
Relative error be less than the predetermined threshold value, then using current initial model as the inverse model.
In one embodiment, in order to correct initial model, the radio-frequency component that reselects can specifically include：Adjust
The parameters such as the parameter values of whole initial model high frequency components, such as peak frequency, dominant frequency.
In the embodiment of the present application, determine microcosmic erosion pores'growth area by using casting body flake and ESEM, join
Examine log and nuclear magnetic resonance log porosity curve determines the logging response character of microcosmic erosion hole, with reference to rock physicses template
And then determine microcosmic erosion hole logging character curve；Finally according to microcosmic erosion hole logging character curve and two-dimentional poststack ground
Shake data, inverting is indicated by seismic waveform, determines microcosmic erosion porosity distribution, cannot be accurate so as to solve existing method
The technical problem of the microcosmic erosion porosity distribution in target area really is determined in two dimensional cross-section.
A kind of microcosmic erosion hole geophysics bivariate table is additionally provided based on same inventive concept, in the embodiment of the present invention
Determining device is levied, as described in the following examples.Due to a kind of microcosmic erosion hole geophysics two-dimensional characterization determining device solution
Certainly the principle of problem determines that method is similar to a kind of microcosmic erosion hole geophysics two-dimensional characterization, therefore a kind of microcosmic erosion hole
The implementation of gap geophysics two-dimensional characterization determining device may refer to a kind of microcosmic erosion hole geophysics two-dimensional characterization and determine
The implementation of method, repeats part and repeats no more.Used below, term " unit " or " module " can realize predetermined function
Software and/or hardware combination.Although the device described by following examples is preferably realized with software, hardware,
Or the realization of the combination of software and hardware is also that may and be contemplated.Fig. 2 is referred to, is a kind of micro- of the embodiment of the present invention
A kind of composition structure chart of dissolution porosity geophysics two-dimensional characterization determining device is seen, the device can include：Acquisition module
201st, analysis module 202, the first determining module 203, the second determining module 204, the 3rd determining module 205, the 4th determining module
206 and the 5th determining module 207.The structure is specifically described below.
Acquisition module 201, specifically can be used for obtaining the sample of section sandstone of coring of target area, the target area
The two-dimentional poststack earthquake money of log data, the nuclear-magnetism porosity logs of the target area and the target area
Material.
Analysis module 202, for carrying out casting body flake and ESEM experimental analysis to the sample, obtains analysis knot
Really.
First determining module 203, for according to the analysis result and the log data, by core Location,
Determine the microcosmic erosion pores'growth area of the sample.
Second determining module 204, for according to the analysis result and the log data, by core Location,
Determine the microcosmic erosion pores'growth area of the sample.
3rd determining module 205, specifically can be used for for the target area being divided into multiple continuous substratums according to depth,
And be that each substratum in the multiple continuous substratum is set up respectively according to the logging response character of the microcosmic erosion hole
Rock physicses amount version.
4th determining module 206, specifically can be used for the rock physicses amount according to each substratum in multiple continuous substratums
Version, determines microcosmic erosion hole logging character curve.
5th determining module 207, specifically can be used for indicating inverting by seismic waveform of the two-dimentional poststack seismic data
Data volume, according to the microcosmic erosion hole logging character curve, inverting is indicated by seismic waveform, obtain the target area
The distribution characteristics of the microcosmic erosion hole in domain.
In one embodiment, the logging response character of the microcosmic erosion hole can specifically include：Velocity of longitudinal wave
More than pre-set velocity, and/or, nuclear magnetic resonance log porosity is less than default porosity.
In one embodiment, in order to determine the logging response character of microcosmic erosion hole, second determining module
204 can specifically include：
First determining unit, for the interval transit time curve in the log data according to the target area, it is determined that
The velocity of longitudinal wave of the target area；
Second determining unit, for the nuclear-magnetism porosity logs according to the target area, determines the target area
The nuclear magnetic resonance log porosity in domain；
3rd determining unit, for the microcosmic erosion pores'growth area according to the sample, determine the pre-set velocity and
The default porosity.
It should be noted that the core Location, can refer to specifically according to log data and well logging and section of coring
Rock core depth carry out core Location jointly.
In one embodiment, in order to determine rock physicses amount version, the 3rd determining module 205 specifically can be used for
According to the logging response character of the microcosmic erosion hole, to each substratum in the multiple continuous substratum, respectively with institute
It is abscissa in the nuclear magnetic resonance log porosity of substratum, with velocity of longitudinal wave as ordinate, sets up the rock physicses amount version of each substratum.
In one embodiment, in order to determine microcosmic erosion hole logging character curve, the 4th determining module 206
Can specifically include：
Information acquisition unit, for each substratum from the multiple continuous substratum rock physicses amount version in by substratum
Obtain the information of microcosmic erosion hole；
4th determining unit, for will obtain each substratum microcosmic erosion hole information, as described microcosmic molten
Pit gap logging character curve.
In one embodiment, in order to determine the distribution characteristics of the microcosmic erosion hole of target area, the described 5th is true
Cover half block 207 can specifically include：
Model determining unit, the data volume for indicating inverting as seismic waveform with the two-dimentional poststack seismic data, and
Using the microcosmic erosion hole logging character curve as constraints, preferably radio-frequency component, to set up inverse model；
5th determining unit, inverting is indicated for according to the inverse model, carrying out seismic waveform, obtains the target area
The distribution characteristics of the microcosmic erosion hole in domain.
It should be noted that the above-mentioned model determining unit inverse model good in order to set up precision, during specific implementation, can be with
The data volume of inverting is first indicated as seismic waveform with the two-dimentional poststack seismic data, and it is special to be logged well with the microcosmic erosion hole
Curve is levied as constraints, preferably radio-frequency component, to set up initial model；Further according to the initial model, inverting obtains institute
State target area and cross microcosmic erosion interstitial space distribution characteristics in the two-dimentional stacked seismic data of well；And then cross well two dimension by described
Microcosmic erosion interstitial space distribution characteristics is logged well with the microcosmic erosion hole on the well of the target area in stacked seismic data
Indicatrix is compared, and obtains relative error；At this moment, if the relative error is less than predetermined threshold value, will be described initial
Model is used as the inverse model；If the relative error is more than or equal to the predetermined threshold value, radio-frequency component is reselected,
And initial model is set up according to the radio-frequency component for reselecting, until according to the initial model set up, the mesh that inverting is obtained
On microcosmic erosion interstitial space distribution characteristics and the well of the target area crossed in well two dimension stacked seismic data in mark region
Microcosmic erosion hole logging character curve between relative error be less than the predetermined threshold value, then current initial model is made
It is the inverse model.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.Especially for system reality
Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
It should be noted that system, device, module or unit that above-mentioned implementation method is illustrated, specifically can be by computer
Chip or entity are realized, or are realized by the product with certain function.For convenience of description, in this manual, retouch
It is divided into various units with function when stating apparatus above to describe respectively.Certainly, can be the function of each unit when the application is implemented
Realized in same or multiple softwares and/or hardware.
Additionally, in this manual, adjective as such as first and second can be only used for an element or dynamic
Make to be made a distinction with another element or action, without requiring or implying any actual this relation or order.Permit in environment
Perhaps in the case of, in only element, part or step is should not be interpreted as limited to reference to element or part or step (s)
It is individual, and can be in element, part or step one or more etc..
As can be seen from the above description, a kind of microcosmic erosion hole geophysics two that the embodiment of the present application is provided
Dimension table levies determination method and apparatus.Microcosmic erosion pores'growth area is determined by using casting body flake and ESEM, with reference to survey
Well curve and nuclear magnetic resonance log porosity curve determine the logging response character of microcosmic erosion hole, with reference to rock physicses template and then
Determine microcosmic erosion hole logging character curve；Finally according to microcosmic erosion hole logging character curve and two-dimentional poststack earthquake money
Material, inverting is indicated by seismic waveform, microcosmic erosion porosity distribution is determined, so as to solve existing microcosmic erosion hole
Determine determine the low technical problem of precision with microcosmic erosion hole using to be limited in scope present in method, reached and cutd open in two dimension
The accurate technique effect for determining microcosmic erosion distribution of pores on face；Sample is carried out further through using casting body flake and ESEM
Analysis, analysis result, log data and nuclear-magnetism porosity logs further according to sample determine microcosmic erosion hole
Logging response character, further increases the precision of microcosmic erosion hole determination；Also inverting is indicated to determine by seismic waveform micro-
The distribution characteristics of dissolution porosity is seen, the distribution and precision for determining microcosmic erosion hole is further improved.
During scene is embodied at one, a kind of microcosmic erosion hole geophysics two dimension is provided using the embodiment of the present application
Characterize and determine that method, device is predicted to the microcosmic erosion hole in somewhere.During specific implementation, can be according to following steps
(step 1 to step 3) is performed.
Step 1, determine microcosmic erosion hole logging response character.
Require prospect pit it is necessary to have interval transit time curve, (i.e. in the application included by log data one first
Plant profile) and nuclear magnetic resonance log porosity data (the nuclear-magnetism porosity logs i.e. in the application).Secondly, to section of coring
Sample of sandstone will carry out casting body flake and scanning electron microscope analysis.And then core Location is carried out, on the basis of core Location, find out micro-
See dissolution porosity Growth Depth scope (the microcosmic erosion pores'growth area i.e. in the application).Finally, with reference to velocity of longitudinal wave and core
It can be found that the place of microcosmic erosion pores'growth, velocity of longitudinal wave is larger, and porosity is smaller (to be obtained final product for magnetic logging porosity curve
The logging response character of microcosmic erosion hole in the application).The W1 borehole logging tool microcosmic erosions hole that Fig. 3 specifically is can be referred to sent out
Educate characteristic synthetic block diagram.
Step 2, by rock physicses amount version, obtain microcosmic erosion hole logging character curve.
For velocity of longitudinal wave and nuclear-magnetism porosity in the range of above-mentioned microcosmic erosion pores'growth logging response character (i.e. this
The logging response character of the microcosmic erosion hole in application), set up with nuclear magnetic resonance log porosity as abscissa, be with velocity of longitudinal wave
The rock physicses amount version of ordinate, it can be found that microcosmic erosion hole is mainly distributed on a left side for rock physicses amount version from amount version
Upper angle, you can to determine that microcosmic erosion hole has larger velocity of longitudinal wave and less nuclear magnetic resonance log porosity.Divide substratum again
This rock physicses amount version is successively set up, microcosmic erosion hole information is read from each amount version, and then can obtain continuous micro-
See dissolution porosity logging character curve.The Sandstone Gas Reservoir microcosmic erosion porous rocks physical quantity version of Fig. 4 can specifically be referred to.
Step 3, by logging well, joint seismic waveform indicates inverting, characterizes and obtains microcosmic erosion hole in two dimensional cross-section
Distribution characteristics.
The microcosmic erosion pore character curve obtained with step 2 is constraints, preferably goes out suitable radio-frequency component, is built
Vertical rational initial model (finding microcosmic erosion hole and seismic waveform relation), indicates inversion prediction to cross well and cuts open with waveform
Microcosmic erosion porosity distribution on face.Contrasted with microcosmic erosion pore character curve on well, if on waveform inversion and well
The rate of coincideing is relatively low, then initial model is re-established, untill efficiency of inverse process is preferable.Fig. 5 can specifically be referred to crosses W1 wells
Crossline (road) direction poststack 2-d seismic data schematic diagram and Fig. 6's is to cross W1 well Crossline directions microcosmic erosion hole
Gap seismic waveform indicates inverting generalized section.Wherein, well side curve is microcosmic erosion pore character curve in Fig. 6.From this
Cross well inverting section and can see the microcosmic erosion hole obtained by inverting and the statistics (signature logging beside well pillar on well
Curve) it is consistent substantially, illustrate the microcosmic erosion hole set up with reference to casting body flake and ESEM and rock physicses amount version
Indicatrix, inverting is indicated by seismic waveform, can relatively accurately characterize development of the microcosmic erosion hole in two-dimensional space
Degree.
A kind of microcosmic erosion hole geophysics two-dimensional characterization provided by above-mentioned concrete application the embodiment of the present application is true
Determine the example of method and apparatus, a kind of microcosmic erosion hole geophysics bivariate table provided by the embodiment of the present application is provided
Levy and determine that method and apparatus can solve the microcosmic erosion that can not accurately determine target area that existing method is present really
Technical problem of the hole in the distribution characteristics of two-dimensional space.
Although different microcosmic erosion holes are mentioned in teachings herein determines method or apparatus, the application is not
It is confined to be professional standard or the situation described by embodiment etc., some professional standards or uses self-defined mode or reality
Apply example description practice processes on embodiment amended slightly can also realize above-described embodiment it is identical, equivalent or close,
Or the implementation result being anticipated that after deformation.Using data acquisition, treatment, output, judgment mode after these modifications or deformation etc.
Embodiment, within the scope of still may belong to the optional embodiment of the application.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on conventional or noninvasive
The means of the property made can include more or less operating procedures.The step of being enumerated in embodiment order is only numerous steps
A kind of mode in execution sequence, unique execution sequence is not represented.When device or client production in practice is performed, can
Performed or executed in parallel (such as at parallel processor or multithreading with according to embodiment or method shown in the drawings order
The environment of reason, even distributed data processing environment).Term " including ", "comprising" or its any other variant be intended to contain
Lid nonexcludability is included, so that process, method, product or equipment including a series of key elements not only will including those
Element, but also other key elements including being not expressly set out, or also include being this process, method, product or equipment
Intrinsic key element.In the absence of more restrictions, be not precluded from the process including the key element, method, product or
Also there are other identical or equivalent elements in person's equipment.
Device that above-described embodiment is illustrated or module etc., can specifically be realized by computer chip or entity, or by having
There is the product of certain function to realize.For convenience of description, it is divided into various modules with function during description apparatus above to retouch respectively
State.Certainly, the function of each module can be realized in same or multiple softwares and/or hardware when the application is implemented,
Can will realize that the module of same function is realized by the combination of multiple submodule.Device embodiment described above is only
Schematically, for example, the division of the module, only a kind of division of logic function, can there is other drawing when actually realizing
The mode of dividing, such as multiple module or components can be combined or be desirably integrated into another system, or some features can be ignored,
Or do not perform.
It is also known in the art that in addition to realizing controller in pure computer readable program code mode, it is complete
Entirely can by by method and step carry out programming in logic come cause controller with gate, switch, application specific integrated circuit, may be programmed
Logic controller realizes identical function with the form of embedded microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the device for realizing various functions included to its inside can also be considered as the structure in hardware component.Or
Person even, can be used to realizing that the device of various functions is considered as not only being the software module of implementation method but also can be hardware
Structure in part.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure, class etc..The application can also be in a distributed computing environment put into practice, in these DCEs,
Task is performed by the remote processing devices connected by communication network.In a distributed computing environment, program module can
With in the local and remote computer-readable storage medium including including storage device.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On the part that is contributed to prior art in other words can be embodied in the form of software product, the computer software product
Can store in storage medium, such as ROM/RAM, magnetic disc, CD, including some instructions are used to so that a computer equipment
(can be personal computer, mobile terminal, server, or network equipment etc.) performs each embodiment of the application or implementation
Method described in some parts of example.
Each embodiment in this specification is described by the way of progressive, same or analogous portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.The application can be used for crowd
In more general or special purpose computing system environments or configuration.For example：Personal computer, server computer, handheld device or
Portable set, laptop device, multicomputer system, the system based on microprocessor, set top box, programmable electronics set
Standby, network PC, minicom, mainframe computer, the DCE including any of the above system or equipment etc..
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without deviating from spirit herein, it is desirable to which appended claim includes these deformations and changes without deviating from the application.
Claims (12)
1. a kind of microcosmic erosion hole geophysics two-dimensional characterization determines method, it is characterised in that including：
Obtain the sample of section sandstone of coring of target area, the log data of the target area, the target area
Nuclear-magnetism porosity logs and the two-dimentional poststack seismic data of the target area；
Casting body flake and ESEM experimental analysis are carried out to the sample, analysis result is obtained；
According to the analysis result and the log data, by core Location, the microcosmic erosion hole of the sample is determined
Gap development area；
The nuclear-magnetism porosity logs and the sample of log data, the target area according to the target area
Microcosmic erosion pores'growth area, determine the logging response character of the microcosmic erosion hole of the target area；
The target area is divided into multiple continuous substratums according to depth, and according to the log response of the microcosmic erosion hole
Feature, is that each substratum in the multiple continuous substratum sets up rock physicses amount version respectively；
According to the rock physicses amount version of each substratum in multiple substratums, microcosmic erosion hole logging character curve is determined；
The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, is logged well according to the microcosmic erosion hole
Indicatrix, inverting is indicated by seismic waveform, obtains the distribution characteristics of the microcosmic erosion hole of the target area.
2. method according to claim 1, it is characterised in that
The logging response character of the microcosmic erosion hole includes：Velocity of longitudinal wave is more than pre-set velocity, and/or, nuclear magnetic resonance log hole
Porosity is less than default porosity.
3. method according to claim 2, it is characterised in that log data according to the target area, described
The nuclear-magnetism porosity logs of target area and the microcosmic erosion pores'growth area of the sample, determine the target area
The logging response character of microcosmic erosion hole, including：
The interval transit time curve in log data according to the target area, determines the compressional wave speed of the target area
Degree；
Nuclear-magnetism porosity logs according to the target area, determine the nuclear magnetic resonance log porosity of the target area；
Microcosmic erosion pores'growth area according to the sample, determines the pre-set velocity and the default porosity.
4. method according to claim 1, it is characterised in that according to the logging response character of the microcosmic erosion hole,
For each substratum in the multiple continuous substratum sets up rock physicses amount version respectively, including：
According to the logging response character of the microcosmic erosion hole, to each substratum in the multiple continuous substratum, difference
Nuclear magnetic resonance log porosity with place substratum, with velocity of longitudinal wave as ordinate, sets up the rock physicses of each substratum as abscissa
Amount version.
5. method according to claim 1, it is characterised in that according to the rock thing of each substratum in multiple continuous substratums
Reason amount version, determines the continuous logging character curve of microcosmic erosion hole, including：
Obtain the letter of microcosmic erosion hole in the rock physicses amount version of each substratum by substratum from the multiple continuous substratum
Breath；
The information of the microcosmic erosion hole of each substratum that will be obtained, as the microcosmic erosion hole logging character curve.
6. method according to claim 1, it is characterised in that with the two-dimentional poststack seismic data be that seismic waveform is indicated
The data volume of inverting, according to the microcosmic erosion hole logging character curve, inverting is indicated by seismic waveform, obtains the mesh
The distribution characteristics of the microcosmic erosion hole in region is marked, including：
The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, and is logged well with the microcosmic erosion hole
Indicatrix as constraints, preferably radio-frequency component, to set up inverse model；
According to the inverse model, carry out seismic waveform indicate inverting, obtain the target area microcosmic erosion hole point
Cloth feature.
7. method according to claim 6, it is characterised in that with the two-dimentional poststack seismic data be that seismic waveform is indicated
The data volume of inverting, and using the microcosmic erosion hole logging character curve as constraints, preferably radio-frequency component, to set up
Inverse model, including：
The data volume of inverting is indicated as seismic waveform with the two-dimentional poststack seismic data, and is logged well with the microcosmic erosion hole
Indicatrix as constraints, preferably radio-frequency component, to set up initial model；
According to the initial model, inverting obtains the target area and crosses microcosmic erosion hole in the two-dimentional stacked seismic data of well
Spatial distribution characteristic；
By on well of the microcosmic erosion interstitial space distribution characteristics with the target area in the two dimension stacked seismic data of well excessively
Microcosmic erosion hole logging character curve be compared, obtain relative error；
If the relative error is less than predetermined threshold value, using the initial model as the inverse model.
8. method according to claim 7, it is characterised in that microcosmic molten on it described will cross well two dimension stacked seismic data
Pit gap spatial distribution characteristic is compared with the microcosmic erosion hole logging character curve on the well of the target area, obtains
After relative error, methods described also includes：
If the relative error is more than or equal to the predetermined threshold value, reselect radio-frequency component, and according to reselecting
Radio-frequency component sets up initial model, until according to the initial model set up, the two dimension of well excessively of the target area that inverting is obtained
Microcosmic erosion interstitial space distribution characteristics in stacked seismic data is surveyed with the microcosmic erosion hole on the well of the target area
Relative error between well indicatrix is less than the predetermined threshold value, then using current initial model as the inverse model.
9. a kind of microcosmic erosion hole geophysics two-dimensional characterization determining device, it is characterised in that including：
Acquisition module, sample, the log data of the target area, the institute of section sandstone of coring for obtaining target area
State the nuclear-magnetism porosity logs of target area and the two-dimentional poststack seismic data of the target area；
Analysis module, for carrying out casting body flake and ESEM experimental analysis to the sample, obtains analysis result；
First determining module, for according to the analysis result and the log data, by core Location, it is determined that described
The microcosmic erosion pores'growth area of sample；
Second determining module, for according to the analysis result and the log data, by core Location, it is determined that described
The microcosmic erosion pores'growth area of sample；
3rd determining module, for the target area to be divided into multiple continuous substratums according to depth, and according to described microcosmic
The logging response character of dissolution porosity, is that each substratum in the multiple continuous substratum sets up rock physicses amount version respectively；
4th determining module, for the rock physicses amount version according to each substratum in multiple continuous substratums, determines microcosmic erosion
Hole logging character curve；
5th determining module, the data volume for indicating inverting as seismic waveform with the two-dimentional poststack seismic data, according to institute
Microcosmic erosion hole logging character curve is stated, inverting is indicated by seismic waveform, obtain the microcosmic erosion hole of the target area
The distribution characteristics of gap.
10. device according to claim 9, it is characterised in that second determining module includes：
First determining unit, for the interval transit time curve in the log data according to the target area, it is determined that described
The velocity of longitudinal wave of target area；
Second determining unit, for the nuclear-magnetism porosity logs according to the target area, determines the target area
Nuclear magnetic resonance log porosity；
3rd determining unit, for the microcosmic erosion pores'growth area according to the sample, determines the pre-set velocity and described
Default porosity.
11. devices according to claim 9, it is characterised in that the 4th determining module includes：
Information acquisition unit, for being obtained by substratum in the rock physicses amount version of each substratum from the multiple continuous substratum
The information of microcosmic erosion hole；
4th determining unit, for will obtain each substratum microcosmic erosion hole information, as the microcosmic erosion hole
Gap logging character curve.
12. devices according to claim 9, it is characterised in that the 5th determining module includes：
Model determining unit, the data volume for indicating inverting as seismic waveform with the two-dimentional poststack seismic data, and with institute
Microcosmic erosion hole logging character curve is stated as constraints, preferably radio-frequency component, to set up inverse model；
5th determining unit, inverting is indicated for according to the inverse model, carrying out seismic waveform, obtains the target area
The distribution characteristics of microcosmic erosion hole.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611261876.7A CN106842318B (en) | 2016-12-30 | 2016-12-30 | Method and device for determining geophysical two-dimensional representation of microscopic erosion pore |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611261876.7A CN106842318B (en) | 2016-12-30 | 2016-12-30 | Method and device for determining geophysical two-dimensional representation of microscopic erosion pore |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106842318A true CN106842318A (en) | 2017-06-13 |
| CN106842318B CN106842318B (en) | 2019-01-18 |
Family
ID=59114990
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611261876.7A Active CN106842318B (en) | 2016-12-30 | 2016-12-30 | Method and device for determining geophysical two-dimensional representation of microscopic erosion pore |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106842318B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108897051A (en) * | 2018-05-15 | 2018-11-27 | 吉林大学 | A kind of high-resolution ground nuclear magnetic resonance imaging method |
| CN111694054A (en) * | 2019-03-15 | 2020-09-22 | 中国石油化工股份有限公司 | Sandstone reservoir mineral corrosion strength calculation method and system based on logging curve |
| CN118052993A (en) * | 2024-04-12 | 2024-05-17 | 山东省地震工程研究院 | Method, system and equipment for determining characteristics of sliding surface of earthquake fault based on scanning electron microscope |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102650701A (en) * | 2012-04-27 | 2012-08-29 | 中国石油天然气股份有限公司 | Method for predicting reservoir wave impedance based on dipole wave |
| CN103077548A (en) * | 2012-05-14 | 2013-05-01 | 中国石油化工股份有限公司 | Method for establishing solution vug reservoir body distribution model of fractured-vuggy carbonate rock reservoir |
| CN103645519A (en) * | 2013-12-17 | 2014-03-19 | 中国石油大学(华东) | Volcanic rock weathering crust identification and classification standard based oil and gas exploration method |
| CN104853822A (en) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | Method for evaluating shale gas reservoir and searching sweet spot region |
| CN104898161A (en) * | 2014-03-05 | 2015-09-09 | 中国石油集团东方地球物理勘探有限责任公司 | Effective sandstone predicting method based on logging response simulator |
| CN105093308A (en) * | 2014-05-19 | 2015-11-25 | 中国石油化工股份有限公司 | Well trajectory design method and system |
| KR101620506B1 (en) * | 2015-10-15 | 2016-05-12 | 한국지질자원연구원 | Stochastic optimization technique for in situ stresses from the probability distributions of rock strength and borehole breakout width |
| CN105626058A (en) * | 2015-12-30 | 2016-06-01 | 中国石油天然气股份有限公司 | Method and device for determining development degree of reservoir karst |
-
2016
- 2016-12-30 CN CN201611261876.7A patent/CN106842318B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102650701A (en) * | 2012-04-27 | 2012-08-29 | 中国石油天然气股份有限公司 | Method for predicting reservoir wave impedance based on dipole wave |
| CN103077548A (en) * | 2012-05-14 | 2013-05-01 | 中国石油化工股份有限公司 | Method for establishing solution vug reservoir body distribution model of fractured-vuggy carbonate rock reservoir |
| CN103645519A (en) * | 2013-12-17 | 2014-03-19 | 中国石油大学(华东) | Volcanic rock weathering crust identification and classification standard based oil and gas exploration method |
| CN104898161A (en) * | 2014-03-05 | 2015-09-09 | 中国石油集团东方地球物理勘探有限责任公司 | Effective sandstone predicting method based on logging response simulator |
| CN105093308A (en) * | 2014-05-19 | 2015-11-25 | 中国石油化工股份有限公司 | Well trajectory design method and system |
| CN104853822A (en) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | Method for evaluating shale gas reservoir and searching sweet spot region |
| KR101620506B1 (en) * | 2015-10-15 | 2016-05-12 | 한국지질자원연구원 | Stochastic optimization technique for in situ stresses from the probability distributions of rock strength and borehole breakout width |
| CN105626058A (en) * | 2015-12-30 | 2016-06-01 | 中国石油天然气股份有限公司 | Method and device for determining development degree of reservoir karst |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108897051A (en) * | 2018-05-15 | 2018-11-27 | 吉林大学 | A kind of high-resolution ground nuclear magnetic resonance imaging method |
| CN108897051B (en) * | 2018-05-15 | 2019-10-01 | 吉林大学 | A kind of high-resolution ground nuclear magnetic resonance imaging method |
| CN111694054A (en) * | 2019-03-15 | 2020-09-22 | 中国石油化工股份有限公司 | Sandstone reservoir mineral corrosion strength calculation method and system based on logging curve |
| CN111694054B (en) * | 2019-03-15 | 2021-12-24 | 中国石油化工股份有限公司 | Sandstone reservoir mineral corrosion strength calculation method and system based on logging curve |
| CN118052993A (en) * | 2024-04-12 | 2024-05-17 | 山东省地震工程研究院 | Method, system and equipment for determining characteristics of sliding surface of earthquake fault based on scanning electron microscope |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106842318B (en) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105334536B (en) | Tight sandstone reservoir chicken-wire cracking system effectiveness evaluation method | |
| CN104992468B (en) | Three-dimensional geological modeling method for fracture-cave carbonate reservoir | |
| CN103775057B (en) | Method and device for identifying effective reservoir of tight oil and gas reservoir | |
| CN106772680B (en) | Two-dimensional post-stack seismic data-based micro-crack characterization determination method and device | |
| CN107402176B (en) | method and device for determining porosity of crack | |
| CN109388817A (en) | A kind of Reservoir Fracture three-dimensional modeling method | |
| CN105701319B (en) | A kind of microfacies modeling method under horizontal well rule well pattern | |
| CN108397130B (en) | Boring method and device | |
| CN106355571B (en) | Method and device for determining quality of dolomite reservoir | |
| CN110414723A (en) | The method, apparatus and system of fractured hydrocarbon reservoir history matching based on microseismic event | |
| Oeth et al. | Characterization of small scale heterogeneity to predict acid fracture performance | |
| CN106842318B (en) | Method and device for determining geophysical two-dimensional representation of microscopic erosion pore | |
| CN106249300B (en) | Method and device for determining TOC content and heterogeneity of source rock based on well-seismic combination | |
| CN103412335B (en) | A kind of method utilizing earthquake thing phase body predicting reservoir | |
| Reeves et al. | Development and Calibration of Dual-Permeability Flow Models with DiscontinuousFault Networks | |
| CN108319743A (en) | Method and device for predicting abundance of oil-gas resources of ancient marine carbonate rock | |
| CN112085305B (en) | Method and device for evaluating reservoir seam network performance | |
| CN106842316B (en) | Crack determination method and apparatus | |
| CN110910267B (en) | Method and device for determining sea-phase carbonate rock interlayer | |
| CN111239809B (en) | Dominant reservoir determination method and system based on attribute fusion principle | |
| Zhang et al. | Two-parameter prestack seismic inversion of porosity and pore-structure parameter of fractured carbonate reservoirs: Part 2—Applications | |
| CN104914463B (en) | The imaging method of small scale FRACTURE PREDICTION | |
| Al-Ameri et al. | Investigating geomechanical considerations on suitable layer selection for hydraulically fractured horizontal wells placement in tight reservoirs | |
| Rustad et al. | Pore scale estimation, up scaling and uncertainty modeling for multiphase properties | |
| CN109581498A (en) | Reservoir thickness distribution determination method, system, equipment, readable medium |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |