CN104977618B - A kind of method evaluated shale gas reservoir and find dessert area - Google Patents
A kind of method evaluated shale gas reservoir and find dessert area Download PDFInfo
- Publication number
- CN104977618B CN104977618B CN201410140366.9A CN201410140366A CN104977618B CN 104977618 B CN104977618 B CN 104977618B CN 201410140366 A CN201410140366 A CN 201410140366A CN 104977618 B CN104977618 B CN 104977618B
- Authority
- CN
- China
- Prior art keywords
- data
- seismic
- parameter
- reservoir
- log
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The present invention is a kind of integrated application rock physicses, log data, comprehensive or wide-azimuth 3D seismic data synthesis progress shale gas evaluating reservoir and the method for finding shale gas exploration and development " dessert area ".The rock geophysical survey of application experiment room, log data comprehensive analysis, rock physicses are modeled, vertical seismic profiling (VSP) data, the high-precision top layer comprehensive modeling static corrections processing of comprehensive or wide-azimuth 3D seismic data, prestack denoising, amplitude is with respect to fidelity processing, geological data drives the raising resolution processes of surface seismic data in borehole restraint and well, fine excision, iteration speed is analyzed, interval velocity field analysis, pre-stack depth migration processing, High-resolution Processing, divide orientation processing, high-precision anisotropy processing, the inversion procedure of prestack poststack data, the comprehensive geophysical survey technology such as analysis of neural network, analyze the feature in " dessert area " in shale reservoir, the gas-bearing property prospect of accurate evaluation shale gas reservoir and " the dessert area " for drawing a circle to approve shale gas exploration and development.
Description
Technical field
It is a kind of integrated application rock physicses, log data, comprehensive the invention belongs to applied geophysics exploitation method
Or the comprehensive geophysical survey such as wide-azimuth 3D seismic data and vertical seismic profiling (VSP) data technology carries out shale gas reservoir and commented
Valency and the method for finding shale gas exploration and development dessert area.
Background technology
Shale gas aboundresources, shale gas exploration and development is expected to alleviate the energy crisis faced, but shale gas is used as one kind
Petroleum resources, although its Reservoir model is different from conventional oil gas reservoir, exploration and development is in the exploratory stage, is concentrated mainly on shale gas
Reservoir model, study in terms of geologic feature, effect of the geophysical techniques in shale gas exploration and development is also to be developed.
Although rock geophysics, geophysical log, seismic exploration technique play vital effect in conventional gas and oil exploration,
But the Reservoir model in shale is focusing only on to shale gas reservoir, in terms of unique geology characteristic, it is rare for the comprehensive of shale gas
Geophysical exploration method is closed, geophysical techniques are in marginal state in shale gas exploration and development research.
The research to shale gas is focused mostly in basic theory at present, and applied geophysics data is studied shale gas
Also in the exploratory stage.Li Zhirong etc. exists《South of Sichuan Basin, China shale gas seismic prospecting new development》(Gas industry, 2011,
31(4):40-43)In one text, on the basis of to south of Sichuan Basin, China shale interval geology, geophysical response characteristic analysis,
Tackled key problems by seismic data acquisition, processing and interpretation technique, form a set of more complete shale gas geophysical exploration and think
Road and techniqueflow, achieve the new development of shale gas seismic prospecting;Qi Baoquan etc. exists《Using the Logging Data To Evaluate Sichuan Basin
Southern shale gas reservoir》(Gas industry, 2011,31 (4):44-47)In one text, by Δ logR approach applications to the Sichuan Basin
In southern shale gas evaluating reservoir, in view of overlapping baseline when recognizing shale gas with porosity and resistivity curve overlay method
The influence of selection, the change of lithology etc., explores the well logging interpretation model of shale gas;Luo Rong etc. exists《Shale gas logging evaluation and ground
Shake prediction, monitoring technology are inquired into》(Gas industry, 2011,31 (4):34-39)In one text, for shale gas and conventional reservoir
Difference, inquired into application of the geophysical exploration technology in shale gas exploration and development, and propose development specifically designed for shale
Three-dimensional geophysical exploration, monitoring and the development technique of gas;Liu Shuanlian and Lu Huangsheng exist《Shale gas logging evaluation technical characterstic and
Evaluation method is inquired into》(Logging technique, 35 (2):113-116)In one text, from investigation North America shale gas successful exploration developing example
Start with, on the basis of reservoir geology background research, analyze the main difference of shale gas and conventional oil gas-bearing formation Logging Evaluation Method
It is different.According to shale gas exploration and development demand, the selection gist and logging evaluation technology of Chinese shale gas logging program have been inquired into.Carry
Go out shale minerals composition and reservoir structure evaluation, shale reservoir establishment of standard, types of fractures identification are commented with rock mechanics parameters
Research in terms of valency, the emphasis that can be evaluated as shale gas logging technique;Pay forever strong wait《Shale gas reservoir reservoir fracturing
Experimental evaluation key technology》(Gas industry, 2011,31 (4):51-54)In one text, from Rock Elastic Parameters angle,
Analysis compared for compact sandstone gas and shale gas reservoir mechanical property feature, for shale rock brittleness feature and reservoir core
The experimental evaluation key technology such as sensitiveness, has carried out substantial amounts of experimental evaluation research, and with the high tracer monitoring of live pressure-break,
Ground micro-seismic Fracturing Monitoring result has carried out comparative analysis, and the exploitation to shale gas is significant;Liu Zhenwu etc. exists《Page
Rock gas exploration develops the demand to geophysical techniques》(Geophysical prospecting for oil, 2011,46 (5):810-818)In one text,
By the demand analysis to shale gas geophysical techniques and the prospect to future development, geophysical techniques conduct is explicitly pointed out
The key technology of shale gas evaluating reservoir and storey increase design, will play an important role in shale gas exploration and development;Nie Xin etc.
People exists《Application of the logging technique in the evaluation of shale gas reservoir mechanical property》(Journal of Engineering Geophysics, 2012,9 (4):
433-439)In one text, several logging methods such as acoustic imaging, resistivity imaging, array acoustic are summarized in shale gas reservoir power
Application and meaning in terms of properties evaluations, and the limitation and applicable elements of various logging methods are analyzed, illustrate to combine
This several logging method can effectively evaluate the mechanical property of shale gas reservoir;Hao build it is winged et al.《Shale gas geophysics
Logging evaluation is summarized》(Advances in Geophysics, 2012,27 (4):1624-1632)In one text, text is for the external especially U.S.
The present situation of recent shale gas exploration and development has carried out extensive literature survey, the current external shale gas geophysical log skill of summary
The current situation of art, for the logging program of the conventional gas bearing shale of the different phase introduction of exploration and development, then summarizes shale
Mud logging response characteristic, and the important parameter of shale gas reservoir evaluation methods and evaluating reservoir has been discussed in detail, include exchange
Carbon content, rock forming mineral component and content, porosity, air content and rock mechanics parameters, finally propose shale gas geophysics
The problem of well logging study is present and development trend.
In summary, in the shale gas exploration carried out at present, only sound out logged well or seismic exploration technique application
Test, not yet application comprehensive geophysical survey technology evaluate the gas-bearing property prospect of shale gas reservoir and find exploration and development shale
" the dessert area " of gas, also not publicly in shale gas evaluating reservoir how the detailed description of integrated application geophysical exploration technology and
Detail.
The content of the invention
For problems of the prior art, the present invention provides a kind of integrated application rock physicses, log data, Quan Fang
Position or wide-azimuth 3D seismic data and vertical seismic profiling (VSP) data evaluation shale gas reservoir and the method for finding dessert area.
The present invention is realized by following steps:
1)Different directions core column is drilled through in the core column of all drilling well different buried depths in exploratory area, core column is vacuumized simultaneously
Pressurization saturation is carried out to it with rock stratum mineralized water resistivity identical mineralized water;
Described different directions are, levels vertical with attitude of stratum and in angle of 45 degrees.
Described core column is the cm of diameter 2.5, the cm of length 5.
2)Under the conditions of laboratory simulation underground confined pressure and pore pressure, the dynamic and static state of the core column after measurement saturation
Elastic parameter, attenuation of elastic wave coefficient, frequency dispersion effect and p-and s-wave velocity anisotropy coefficient, obtain rock core dynamic and static
The conversion relational expression of modulus of elasticity, carries out anisotropic rock physical analogy and elastic parameter is calculated with intersecting;
According to intersection result, combination and shale gas " dessert area " parameter of sensibility elasticity parameter or sensibility elasticity parameter are obtained
Corresponding dependency relation, ask for and predict the parameter or parameter combination in shale gas " dessert area ";
3)The log data of all drillings in exploratory area is obtained, processing is corrected to the log data for surveying all in area,
Eliminate shadow of the factors such as wellhole environment, hole deviation change, well liquid change, the change of well temperature and logger error to log
Ring, the optimal log of stratum physical property change can truly be reflected by obtaining;
Using many ensaying methods and core test analysis method, calculate subterranean minerals composition and content, density of earth formations,
Compressional wave and shear wave velocity and porosity, and the rock physicses from earth's surface to shaft bottom are set up according to full well section geophysical log curve
Model;
Described optimal log is to eliminate the change of drilling internal diameter, hole deviation change, well liquid change, the change of well temperature, well logging
Speed is uneven, downhole instrument is stuck, after the non-at the uniform velocity rotation of instrument and the factor such as logger error, can truly reflect stratum
The optimal log of change in physical.
4)The replacement perturbation analysis of the attributes such as fluid, porosity, lithology data is carried out to the log after correction process;
Described perturbation analysis is the corresponding log obtained by changing after formation fluid, porosity or lithology, is looked for
Go out changing rule of the correspondence log with different fluid, porosity or lithology.
5)Mineral constituent analysis is done using logging principle associate(d) matrix method for solving is optimized to optimal log, obtained
The content and its regularity of distribution of mineral in full well section, and calculate the total saturation degree of mineralogical composition and stratum;
Described mineral are clay, calcite, quartz, pyrite, total content of organic carbon(TOC)With the mineral such as dolomite.
Described optimal log is clay mineral curve, bulk density curve, stratum uranium content in log data
Curve, neutron porosity curve, resistivity curve, compressional wave time difference curve and shear wave slowness curve.
6)Full well section petrophysical model is set up, by the velocity of longitudinal wave predicted according to petrophysical model, shear wave velocity, close
Degree, vertical S-wave impedance and Poisson's ratio curve are contrasted with the log surveyed, to predict the identical journey with measured curve
Spend to verify the reliability and reasonability of petrophysical model;
7)With step 2)Core column measurement dynamic and static elastic parameter, attenuation of elastic wave coefficient, frequency dispersion effect and
P-and s-wave velocity anisotropy coefficient demarcation by log calculate or predict come result;
8)The rock constituents perturbation analysis of total content of organic carbon, quartz, clay mineral etc. is carried out to log data;
Described rock constituents perturbation analysis is the percentage composition by changing different minerals in petrophysical model, meter
Corresponding log is calculated, according to the size of the log variable quantity calculated, corresponding Mineral change is found out most sensitive
Property parameters or Sensitive Attributes parameter combination.
9)Many attribute intersections are carried out to various reservoir attribute parameters, favourable shale interval are obtained according to cross plot result each
Attributive character, it is determined that for predicting shale gas " dessert area " associated parameter or parameter combination;
Described parameter or parameter combination is modulus of elasticity, Young's modulus of elasticity, density, the coefficient of rigidity, springform
The sum of products Young's modulus of elasticity and the product of density of amount and the product, the coefficient of rigidity and density of density.
10)Utilize step 6)The full well section petrophysical model set up, obtains the artificial earthquake synthesis of petrophysical model
Record or trace gather, well shake demarcation is carried out with log data and artificial earthquake composite traces or trace gather and is handled, in shale reservoir depth
Nearby carry out AVO(Amplitude changes with geophone offset)And AVA(Amplitude is with azimuthal variation)Analysis;
11)Comprehensive or wide-azimuth 3D seismic data is gathered in exploratory area;
12)Two dimension Walkaway VSP are gathered in the well in exploratory area(Mobile geophone offset vertical seismic profiling (VSP))Or three-dimensional VSP
(Vertical seismic profiling (VSP))Data;Or with Three Dimensional Ground geological data synchronous acquisition two dimension Walkaway VSP(Mobile geophone offset
Vertical seismic profiling (VSP))Or three-dimensional VSP(Vertical seismic profiling (VSP))Data;
13)To the two dimension or three-dimensional VSP in exploratory area(Vertical seismic profiling (VSP))Data are according to the depth and ground of well geophone
Seismic wave carries out velocity analysis, migration imaging and inverting from what ground reached well geophone when walking, obtain accurate formation velocity,
Earth-attenuation coefficient(Q values)With the anisotropic parameters of each formation velocity;
14)Or wide-azimuth 3D seismic data comprehensive to ground carries out high-precision top layer comprehensive modeling, calculates static correction
Amount, carries out static corrections processing;Surface seismic data is handled with vertical seismic profiling (VSP) data-driven in borehole restraint and well, ground is improved
The resolution ratio and precision of geological data, then carry out fine excision and iteration speed is calculated, then complete velocity modeling and three-dimensional
Pre-stack time migration and three-dimensional pre-stack depth migration imaging processing;
Described top layer comprehensive modeling static correction is:Static corrections processing, prestack denoising, amplitude compensation, Q values(Earth-attenuation)
Compensation, surface consistent deconvolution and predictive deconvolution amplitude are with respect to fidelity processing.
15)Raising resolution processes are carried out to the data after the processing of three-dimensional pre-stack depth migration imaging;
16)With the seismic channel high resolution processing method and tool of the imparametrization analysis of spectrum based on statistical adaptive signal theory
There is the high-resolution subsurface reflective information estimating method of fidelity, high-resolution is carried out to the data after the processing of three-dimensional pre-stack depth migration
Rate processing.
Described reflective information method of estimation is to be based on statistical signal self-adaptive processing, using nonparametric spectral analysis method and
High-resolution subsurface reflective information estimating method with fidelity, is keeping original seismic data information and is not damaging to greatest extent
Lose in original data on the premise of small geological information, obtain high-resolution complex seismic trace collection.
Described reflective information method of estimation is based on statistical signal self-adaptive processing nonparametric spectrum analysis theory, passes through simulation
With respect to the statistical nature of interference, the adaptively reflection amplitudes progress to different time position is stably estimated exactly, so as to carry
High section resolution ratio, widens frequency band, and original seismic data information can be kept to greatest extent and does not lose micro- in original data
Small geological information, obtains the high-resolution complex seismic trace collection of fidelity.
17)Accurate buried depth, thickness, occurrence and the planar distribution of shale reservoir are extracted from three-dimension high-resolution seismic data;
18)Inversion of three dimensional high-resolution post-stack seismic data is to obtain the seismic attributes data body of post-stack inversion, for solving
Release tomography and crack;
19)Utilize relevant and association attributes(Similitude, characteristic value similitude)Inclination angle and inclination angle orientation characteristics, minimax
Curvature, positive camber and negative cruvature attribute describe and characterized the Distribution Characteristics of subsurface fault, crack crack and tectonic boundary;
20)Utilize KSOM(Unsupervised adaptive statistical model)Neural computing method, it is automatic by nonlinear way
To coherence, minimum and maximum curvature, curvature morphological indices, 6 attribute such as instantaneous inclination angle and inclination angle orientation are classified, according to
The distribution characteristics of fracture spacing determines seismic facies body, sets up Earthquake Faulting phase, draws tomography and fracture belt distributed data body, uses
To characterize seismic facies anomalous body and slit band;
21)Automatic tomography pickup is carried out using poststack attribute data;
Described tomography pickup is to calculate section automatically based on coherent body, characteristic value similitude or curvature body, it is determined that macroscopical
Crack and craven fault.
22)Processing is corrected and evened up in the optimization, denoising, stretching for carrying out prestack seismic gather;
23)The elliptical velocity inverting of earthquake data before superposition is carried out, while according to the change of interval velocity in shale reservoir and difference
It is different, determine that the higher-pressure region in shale reservoir is determined in strata pressure doubling-up;
Described elliptical velocity inverting is to RMS(Root-mean-square value)The bearing data body of speed carries out elliptical velocity point
Analysis, obtains fracture strike orientation and compressional wave anisotropic parameters.
24)Carry out the AVO of three-dimensional earthquake data before superposition(Amplitude changes with geophone offset)Synchronous wave impedance is anti-with compressional wave shear wave
Drill;
The synchronous wave impedance inversion of described compressional wave shear wave is to calculate AVO(Amplitude changes with geophone offset)Gradient attribute, and
Inverting angular stack seismic data, synchronously obtains p-wave impedance, S-wave impedance and other derivation resilient properties, particularly λ ρ
(The product of modulus of elasticity and density, μ ρ(The product of the coefficient of rigidity and density)、Eρ(Young's modulus of elasticity and density multiply
Product).
25)Carry out the elliptic inverse of the anisotropic parameters of three-dimensional earthquake data before superposition;
Described elliptic inverse is that azimuthal gradient and speed do elliptic inverse, to obtain Thomson(Thomsen)Ginseng
Number.Converted by rock physicses, by the geomechanics anisotropic parameters of layer for the purpose of Thomson Parameter Switch, such as Young elasticity
Modulus, Poisson's ratio etc.;
26)Carry out the modulus of elasticity λ ρ of earthquake data before superposition(The product of modulus of elasticity and density)、μρ(The coefficient of rigidity
With the product of density)、Eρ(The product of Young's modulus of elasticity and density)Elliptic inverse, obtain anisotropic elasticity modulus, lead to
Rock physical analysis is crossed, anisotropic elasticity modulus is converted to the reservoir parameter of target zone;
Described reservoir parameter is rock brittleness, lithology, porosity, fluid, total organic carbon(TOC)Content etc..
27)To various sign tomographies and the joint geologic interpretation of the seismic properties in crack and demarcation;
Described joint geologic interpretation and demarcation reservoir petrologic characteristic parameter body log calibration, crack pit shaft into
As data and/or core analyzing data are demarcated, large scale tomography and microcosmic tomography with pressure break micro-seismic monitoring achievement and pit shaft into
As data demarcation, stress anisotropy carries out local demarcation with pressure break micro-seismic monitoring achievement.Calibration process i.e. with calculated value with
Measured result is contrasted, and finds out difference value or coefficient correlation between the two, then to calculated value carry out system correction or
Correction, to ensure that the calculated value in the local eyeball in underground is consistent with measurement result.
28)According to rammell fracture development status, determine that possible completion formation damage area and fracturing fluid disturb offset well
Possibility;
29)According to step 2)Rock core dynamic and static modulus of elasticity conversion relational expression, by three-dimensional earthquake data before superposition
Anisotropic elastic wave Simultaneous Retrieving obtain the dynamic modulus of elasticity be converted to static modulus of elasticity;
30)Using the correlation of static modulus of elasticity and rock brittleness, the fragility of shale reservoir is determined(Rupturable property)Point
Cloth rule and feature, the completion and design of hydraulic fracturing scheme of optimum level well;
The completion and design of hydraulic fracturing scheme of described optimum level well be by horizontal well be laid in fragility it is higher and be easy to pressure
In the shale containing high total organic carbon split, and the spacing of optimization design each fracturing section.
31)Using the regularity of distribution of static modulus of elasticity or derivation static modulus of elasticity in shale reservoir, shale is determined
The Brittleness of reservoir, obtains the orientation and intensity of local crustal stress, determines shale reservoir interrupting layer, crack and the orientation in crack
Trend and dense degree, draw a circle to approve and predict the high total organic carbon in shale reservoir(TOC)High stratum in content and shale reservoir
Pressure area;
32)The various favourable parameters of the comprehensive shale gas reservoir obtained, with reference to the accurate buried depth, thickness, production of shale reservoir
Shape and planar distribution, obtain the gas-bearing property prospect of shale gas reservoir and draw a circle to approve " the dessert area " of shale gas exploration and development.
The high totals content of organic carbon of described favourable parameters, including but not limited to shale, the fragility of shale reservoir, tomography,
The orientation and density in crack and crack, the orientation of local crustal stress and intensity, partial high pressure area and porosity distribution.
The present invention can analyze the relation between reservoir parameter and rock geophysical property, accurate determination shale reservoir
Accurate buried depth, thickness, occurrence and planar distribution, evaluate total content of organic carbon in shale gas reservoir or abundance of organic matter exactly
The development degree in Paraclase crack, the both macro and micro intensity azimuth distribution rule of crustal stress, calculating in distribution, prediction exploratory area
Local pressure exceptions area and porosity distribution, overall merit shale gas in the fragility and toughness characteristics on stratum, prediction shale reservoir
The gas-bearing property prospect of reservoir and " the dessert area " for drawing a circle to approve shale gas exploration and development, horizontal well is carried out using Comprehensive Geophysics achievement
Track design and Fracturing Project optimization, for shale gas extensive exploration and successfully exploitation provide important geophysics into
Really.
Accurate buried depth, thickness, occurrence, planar distribution, TOC of the invention according to shale reservoir(Total content of organic carbon)Or have
The features such as the distribution of machine matter abundance, the development degree equal strength azimuth distribution rule in Paraclase crack, can be evaluated shale gas storage
The gas-bearing property prospect of floor and prediction " dessert area " distribution, instruct the design and Fracturing Project optimization of shale gas horizontal well path, are
The extensive exploration and exploitation of shale gas provide important geophysical techniques guarantee.
Brief description of the drawings
Fig. 1 is schematic flow sheet of the present invention.
Embodiment
The present invention is described in detail below in conjunction with accompanying drawing.
The present invention is mainly made up of step in detail below:
It is of the invention to be realized by following steps, as shown in Figure 1:
1)Different directions core column is drilled through in the core column of all drilling well different buried depths in exploratory area, core column is vacuumized simultaneously
Pressurization saturation is carried out to it with rock stratum mineralized water resistivity identical mineralized water.Different directions are, water vertical with attitude of stratum
It is gentle in angle of 45 degrees, core column is the cm of diameter 2.5, the cm of length 5.
2)Under the conditions of laboratory simulation underground confined pressure and pore pressure, the dynamic and static state of the core column after measurement saturation
Elastic parameter, attenuation of elastic wave coefficient, frequency dispersion effect and compressional wave shear wave velocity anisotropy coefficient, obtain rock core dynamic and static
The conversion relational expression of modulus of elasticity, carries out anisotropic rock physical analogy and elastic parameter is calculated with intersecting.According to intersection
As a result, the combination of sensibility elasticity parameter or sensibility elasticity parameter and the corresponding dependency relation of shale gas " dessert area " parameter are obtained,
Ask for and predict the parameter or parameter combination in shale gas " dessert area ".
Step 1)With 2)It is that the rock core dynamic and the measure of static elastic parameter in left side in Fig. 1 and analysis are calculated.
3)The log data of all drillings in exploratory area is obtained, processing is corrected to surveying all log datas in area, disappears
Except influence of the factor to log such as wellhole environment, hole deviation change, well liquid change, the change of well temperature and logger errors,
The optimal log of stratum physical property change can truly be reflected by obtaining.Using many ensaying methods and core test point
Analysis method, calculates subterranean minerals composition and content, density of earth formations, compressional wave shear wave velocity and porosity, and according to the full well section earth
Physical logging curve sets up the petrophysical model from earth's surface to shaft bottom.Optimal log is to eliminate the change of drilling internal diameter, well
Oblique change, well liquid, which change, well temperature changes, logging speed is uneven, downhole instrument is stuck, instrument is non-at the uniform velocity rotates and well logging
After instrument error factor, it can truly reflect the optimal log of stratum physical property change.
4)The replacement disturbance of fluid, porosity, the attribute of lithology data is carried out to the log after correction processAnalysis.
Perturbation analysis is the corresponding log obtained by changing after formation fluid, porosity or lithology, finds out correspondence log
With the changing rule of different fluid, porosity or lithology.
5)Mineral constituent analysis is done using logging principle associate(d) matrix method for solving is optimized to optimal log, obtained
The content and its regularity of distribution of mineral in full well section, and calculate the total saturation degree of mineralogical composition and stratum.Mineral are clay, Fang Xie
Stone, quartz, pyrite, total content of organic carbon(TOC)With the mineral such as dolomite.Optimal log is the clay in log data
Mineral curve, bulk density curve, stratum uranium content curve, neutron porosity curve, resistivity curve, compressional wave time difference curve and
Shear wave slowness curve.
6)Full well section petrophysical model is set up, by the velocity of longitudinal wave predicted according to petrophysical model, shear wave velocity, close
Degree, vertical S-wave impedance and Poisson's ratio curve are contrasted with the log surveyed, to predict the identical journey with measured curve
Spend to verify the reliability and reasonability of petrophysical model.
7)With step 2)Core column measurement dynamic and static elastic parameter, attenuation of elastic wave coefficient, frequency dispersion effect and
P-and s-wave velocity anisotropy coefficient demarcation by log calculate or predict come result.
8)The rock constituents perturbation analysis of total content of organic carbon, quartz, clay mineral etc. is carried out to log data.Rock group
It is the percentage composition by changing different minerals in petrophysical model to divide perturbation analysis, calculates corresponding log, root
According to the size of the log variable quantity calculated, corresponding Mineral change most sensitive property parameters or Sensitive Attributes are found out
The combination of parameter.
9)Many attribute intersections are carried out to various reservoir attribute parameters, favourable shale interval are obtained according to cross plot result each
Attributive character, it is determined that for predicting shale gas " dessert area " associated parameter or parameter combination.Parameter or parameter combination are bullets
Property modulus, Young's modulus of elasticity, density, the coefficient of rigidity, product, the coefficient of rigidity and the density of modulus of elasticity and density
Sum of products Young's modulus of elasticity and density product.
10)Utilize step 6)The full well section petrophysical model set up, obtains the artificial earthquake synthesis of petrophysical model
Record or trace gather, well shake demarcation is carried out with log data and artificial earthquake composite traces or trace gather and is handled, in shale reservoir depth
Nearby carry out AVO(Amplitude changes with geophone offset)And AVA(Amplitude is with azimuthal variation)Analysis.
Step 3)To step 10)Be log data is corrected in Fig. 1, mineral constituent calculating, geophysical log number
Perturbation analysis, artificial earthquake composite traces and AVO/AVA trace gathers are replaced according to analysis and rock physicses modeling, rock constituents and attribute
The work such as analysis.
11)Comprehensive or wide-azimuth 3D seismic data is gathered in exploratory area.
12)Two dimension Walkaway VSP are gathered in the well in exploratory area(Mobile geophone offset vertical seismic profiling (VSP))Or three-dimensional VSP
(Vertical seismic profiling (VSP))Data;Or with Three Dimensional Ground geological data synchronous acquisition two dimension Walkaway VSP(Mobile geophone offset
Vertical seismic profiling (VSP))Or three-dimensional VSP(Vertical seismic profiling (VSP))Data.
13)To the two dimension or three-dimensional VSP in exploratory area(Vertical seismic profiling (VSP))Data are according to the depth and ground of well geophone
Seismic wave carries out velocity analysis, migration imaging and inverting from what ground reached well geophone when walking, obtain accurate formation velocity,
Earth-attenuation coefficient(Q values)With the anisotropic parameters of each formation velocity.
14)Or wide-azimuth 3D seismic data comprehensive to ground carries out high-precision top layer comprehensive modeling, calculates static correction
Amount, carries out static corrections processing;Surface seismic data is handled with vertical seismic profiling (VSP) data-driven in borehole restraint and well, ground is improved
The resolution ratio and precision of geological data, then carry out fine excision and iteration speed is calculated, then complete velocity modeling and three-dimensional
Pre-stack time migration and three-dimensional pre-stack depth migration imaging processing.Top layer comprehensive modeling static correction is:Static corrections processing, prestack are gone
Make an uproar, amplitude compensation, Q values(Earth-attenuation)Compensation, surface consistent deconvolution and predictive deconvolution amplitude are with respect to fidelity processing.
Step 11)To step 14)It is that the comprehensive or wide-azimuth 3D seismic data of collection and two-dimensional movement geophone offset are vertical
Seismic profile or three-dimensional perpendicular seismic profile data, and carry out the processing of vertical seismic profiling (VSP) data and hung down with borehole restraint and well
Straight seismic profile data-driven processing surface seismic data processing.
15)Raising resolution processes are carried out to the data after the processing of three-dimensional pre-stack depth migration imaging.
16)With the seismic channel high resolution processing method and tool of the imparametrization analysis of spectrum based on statistical adaptive signal theory
There is the high-resolution subsurface reflective information estimating method of fidelity, high-resolution is carried out to the data after the processing of three-dimensional pre-stack depth migration
Rate processing.Reflective information method of estimation is to be based on statistical signal self-adaptive processing, using nonparametric spectral analysis method and with guarantor
The high-resolution subsurface reflective information estimating method really spent, is keeping original seismic data information and is not losing original to greatest extent
In data on the premise of small geological information, high-resolution complex seismic trace collection is obtained.Reflective information method of estimation is based on statistics and believed
Number self-adaptive processing nonparametric spectrum analysis theory, by simulating the statistical nature of relative interference, adaptively to different time position
The reflection amplitudes put carry out stabilization and estimated exactly, so as to improve section resolution ratio, widen frequency band, can keep to greatest extent
Original seismic data information and small geological information in original data is not lost, obtain the high-resolution complex seismic trace collection of fidelity.
17)Accurate buried depth, thickness, occurrence and the planar distribution of shale reservoir are extracted from three-dimension high-resolution seismic data.
Step 15)To step 17)It is that the data after the processing of three-dimensional pre-stack depth migration imaging improve at resolution ratio
The structure interpretation of shale reservoir is managed and carried out, the information such as accurate buried depth, thickness, occurrence and the planar distribution of shale reservoir are extracted.
18)Inversion of three dimensional high-resolution post-stack seismic data is to obtain the seismic attributes data body of post-stack inversion, for solving
Release tomography and crack.
19)Utilize relevant and association attributes(Similitude, characteristic value similitude)Inclination angle and inclination angle orientation characteristics, minimax
Curvature, positive camber and negative cruvature attribute describe and characterized the Distribution Characteristics of subsurface fault, crack crack and tectonic boundary.
20)Utilize KSOM(Unsupervised adaptive statistical model)Neural computing method, it is automatic by nonlinear way
To coherence, minimum and maximum curvature, curvature morphological indices, 6 attribute such as instantaneous inclination angle and inclination angle orientation are classified, according to
The distribution characteristics of fracture spacing determines seismic facies body, sets up Earthquake Faulting phase, draws tomography and fracture belt distributed data body, uses
To characterize seismic facies anomalous body and slit band.
21)Automatic tomography pickup is carried out using poststack attribute data.Tomography pickup is to be based on coherent body, characteristic value similitude
Or curvature body calculates section automatically, determines macrocrack and craven fault.
Step 18)To step 21)It is that inversion procedure, neutral net meter are carried out to three-dimension high-resolution post-stack seismic data
Calculate, then obtain the Distribution Characteristics of subsurface fault, crack crack and tectonic boundary.
22)Processing is corrected and evened up in the optimization, denoising, stretching for carrying out prestack seismic gather.Including a point azimuthal velocity analysis,
Divide the process steps such as orientation, subangle and full angle superposition.
23)The elliptical velocity inverting of earthquake data before superposition is carried out, while according to the change of interval velocity in shale reservoir and difference
It is different, determine that the higher-pressure region in shale reservoir is determined in strata pressure doubling-up.Described elliptical velocity inverting is to RMS(Root-mean-square value)Speed
The bearing data body of degree carries out elliptical velocity analysis, obtains fracture strike orientation and compressional wave anisotropic parameters.
24)Carry out the AVO of three-dimensional earthquake data before superposition(Amplitude changes with geophone offset)The synchronous wave impedance inversion of ripple in length and breadth.
Compressional wave shear wave synchronization wave impedance inversion is to calculate AVO(Amplitude changes with geophone offset)Gradient attribute, and inverting angular stack
Shake data, synchronously obtains p-wave impedance, S-wave impedance and other derivation resilient properties, particularly λ ρ(Modulus of elasticity and density
Product, μ ρ(The product of the coefficient of rigidity and density)、Eρ(The product of Young's modulus of elasticity and density).
25)Carry out the elliptic inverse of the anisotropic parameters of three-dimensional earthquake data before superposition.Elliptic inverse is azimuthal ladder
Degree and speed do elliptic inverse, to obtain Thomson(Thomsen)Parameter, is converted by rock physicses, and Thomson parameter is turned
It is changed to the geomechanics anisotropy parameter of target zone, such as Young's modulus of elasticity, Poisson's ratio.
26)Carry out modulus of elasticity, the λ ρ of earthquake data before superposition(The product of modulus of elasticity and density)、μρ(Shearing elasticity mould
The product of amount and density)、Eρ(The product of Young's modulus of elasticity and density)Elliptic inverse, obtain anisotropic elasticity modulus,
By Rock physical analysis, anisotropic elasticity modulus is converted to the reservoir parameter of target zone.Reservoir parameter be rock brittleness,
Lithology, porosity, fluid, total organic carbon(TOC)Content etc..
Step 22)To step 26)Be prestack seismic gather is optimized, denoising, stretching correction and even up processing, then
Inversion procedure is carried out, the modulus of elasticity that inverting is obtained is converted to the reservoir parameter of target zone, such as rock brittleness, lithology, hole
Degree, fluid, total organic carbon(TOC)Content etc..
27)To various sign tomographies and the joint geologic interpretation of the seismic properties in crack and demarcation.Joint geologic interpretation with
Reservoir petrologic characteristic parameter body log calibration is demarcated, crack is demarcated with pit shaft imaging data and/or core analyzing data,
Large scale tomography and microcosmic tomography are demarcated with pressure break micro-seismic monitoring achievement and pit shaft imaging data, stress anisotropy pressure break
Micro-seismic monitoring achievement carries out local demarcation.Calibration process is contrasted with calculated value with measured result, is found out between the two
Difference value or coefficient correlation, the correction or correction of system are then carried out to calculated value, to ensure in the local eyeball in underground
Calculated value is consistent with measurement result.
28)According to rammell fracture development status, determine that possible completion formation damage area and fracturing fluid disturb offset well
Possibility.
29)According to step 2)Rock core dynamic and static modulus of elasticity conversion relational expression, by three-dimensional earthquake data before superposition
Anisotropic elastic wave Simultaneous Retrieving obtain the dynamic modulus of elasticity be converted to static modulus of elasticity.
30)Using the correlation of static modulus of elasticity and rock brittleness, the fragility of shale reservoir is determined(Rupturable property)Point
Cloth rule and feature, the completion and design of hydraulic fracturing scheme of optimum level well.The completion and design of hydraulic fracturing scheme of optimum level well be
Horizontal well is laid in fragility is higher and be easy in the shale containing high total organic carbon of pressure break, and optimization design each fracturing section
Spacing.
31)Using the regularity of distribution of static modulus of elasticity or derivation static modulus of elasticity in shale reservoir, determine shale
The Brittleness of reservoir, obtains the orientation and intensity of local crustal stress, determines shale reservoir interrupting layer, crack and the orientation in crack
Trend and dense degree, draw a circle to approve and predict the high total organic carbon in shale reservoir(TOC)High stratum in content and shale reservoir
Pressure area.
32)The various favourable parameters of the comprehensive shale gas reservoir obtained, with reference to the accurate buried depth, thickness, production of shale reservoir
Shape and planar distribution, obtain the gas-bearing property prospect of shale gas reservoir and draw a circle to approve " the dessert area " of shale gas exploration and development.
Step 27)To step 32)It is the joint geologic interpretation to the various seismic properties for characterizing tomographies and crack and demarcation.
And the various favourable parameters of the shale gas reservoir of shale gas reservoir are obtained by integrated interpretation, finally determine gas-bearing property prospect doubling-up
Determine " the dessert area " of shale gas exploration and development(The quantitative analysis flow seen below Fig. 1).
Claims (19)
1. a kind of method evaluated shale gas reservoir and find dessert area, feature are through the following steps that realize:
1) drill through different directions core column in the core column of all drilling well different buried depths in exploratory area, by core column vacuumize and with
Rock stratum mineralized water resistivity identical mineralized water carries out pressurization saturation to it;
2) under the conditions of laboratory simulation underground confined pressure and pore pressure, the dynamic and static elastic of the core column after measurement saturation
Parameter, attenuation of elastic wave coefficient, frequency dispersion effect and p-and s-wave velocity anisotropy coefficient, obtain rock core dynamic and static elastic mould
The conversion relational expression of amount, carries out anisotropic rock physical analogy and elastic parameter is calculated with intersecting;
According to intersection result, the combination for obtaining sensibility elasticity parameter or sensibility elasticity parameter is corresponding with shale gas dessert area parameter
Dependency relation, asks for and predicts the parameter or parameter combination in shale gas dessert area;
3) log data of all drillings in exploratory area is obtained, processing is corrected to surveying the log data in area, eliminates wellhole
The influence of environment, hole deviation change, well liquid change, the change of well temperature and logger error to log, acquisition can be true
Reflect the optimal log of stratum physical property change;
Using many ensaying methods and core test analysis method, subterranean minerals composition and content, density of earth formations, in length and breadth are calculated
Wave velocity and porosity, and lithology/rock physicses mould from earth's surface to shaft bottom is set up according to full well section geophysical log curve
Type;
4) fluid, porosity, lithology data are carried out to the log after correction process and replaces perturbation analysis;
5) mineral constituent analysis is done using logging principle associate(d) matrix method for solving is optimized to optimal log, obtains full well
The content and its regularity of distribution of mineral in section, and calculate the total saturation degree of mineralogical composition and stratum;
6) set up full well section lithology/petrophysical model, by the velocity of longitudinal wave predicted according to petrophysical model, shear wave velocity,
Density, compressional wave and S-wave impedance and Poisson's ratio curve and the log of actual measurement are contrasted, with prediction and measured curve
Degree of agreement verify the reliability and reasonability of lithology/petrophysical model;
7) use step 2) the core column dynamic and static elastic parameter, attenuation of elastic wave coefficient, frequency dispersion effect and in length and breadth that measure
Wave velocity anisotropy coefficient demarcation by log calculate or predict come result;
8) the rock constituents perturbation analysis of total content of organic carbon, quartz, clay mineral is carried out to log data;
9) many attribute intersections are carried out to reservoir attribute parameter, each attributive character of shale interval is obtained according to cross plot result, really
The fixed parameter or parameter combination for being used to predict that shale gas dessert area is associated;
10) utilize step 6) set up full well section petrophysical model, obtain petrophysical model artificial earthquake composite traces
Or trace gather, well shake demarcation is carried out with log data and artificial earthquake composite traces or trace gather and is handled, near shale reservoir depth
Enter row amplitude to analyze with azimuthal variation with geophone offset change and amplitude;
11) comprehensive or wide-azimuth 3D seismic data is gathered in exploratory area;
12) two-dimensional movement geophone offset vertical seismic profiling (VSP) or three-dimensional perpendicular seismic profile data are gathered in the well in exploratory area;Or
With Three Dimensional Ground geological data synchronous acquisition two-dimensional movement geophone offset vertical seismic profiling (VSP) or three-dimensional perpendicular seismic profile data;
13) to the two dimension or three-dimensional perpendicular seismic profile data in exploratory area according to the depth and seismic wave of well geophone from ground
Reach well geophone carries out velocity analysis, migration imaging and inverting when walking, obtain accurate formation velocity, earth-attenuation system
The anisotropic parameters of number and each formation velocity;
14) comprehensive to ground or wide-azimuth 3D seismic data carries out high-precision top layer comprehensive modeling, and deriving static correction values are entered
Row static corrections processing;Surface seismic data is handled with vertical seismic profiling (VSP) data-driven in borehole restraint and well, surface seismic is improved
The resolution ratio and precision of data, then carry out fine excision and iteration speed is calculated, then complete velocity modeling and three-dimensional prestack
Time migration and three-dimensional pre-stack depth migration imaging processing;
15) raising resolution processes are carried out to the data after the processing of three-dimensional pre-stack depth migration imaging;
16) with the seismic channel high resolution processing method of the imparametrization analysis of spectrum based on statistical adaptive signal theory and with guarantor
The high-resolution subsurface reflective information estimating method really spent, is carried out at high-resolution to the data after the processing of three-dimensional pre-stack depth migration
Reason;
17) accurate buried depth, thickness, occurrence and the planar distribution of shale reservoir are extracted from three-dimension high-resolution seismic data;
18) inversion of three dimensional high-resolution post-stack seismic data is disconnected for explaining to obtain the seismic attributes data body of post-stack inversion
Layer and crack;
19) using relevant and association attributes inclination angle and inclination angle orientation characteristics, minimax curvature, positive camber and negative cruvature attribute come
Describe and characterize subsurface fault, the Distribution Characteristics of crack crack and tectonic boundary;
20) unsupervised adaptive statistical model neural computing method is utilized, by nonlinear way automatic phasing dryness, most
Small and maximum curvature, curvature morphological indices, instantaneous inclination angle and inclination angle orientation is classified, according to the distribution characteristics of fracture spacing come
Determine seismic facies body, set up Earthquake Faulting phase, draw tomography and fracture belt distributed data body, for characterize seismic facies anomalous body and
Slit band;
21) automatic tomography pickup is carried out using poststack attribute data;
22) processing is corrected and evened up in the optimization, denoising, stretching for carrying out prestack seismic gather;
23) the elliptical velocity inverting of earthquake data before superposition is carried out, while according to the change of interval velocity in shale reservoir and difference, really
Determine strata pressure and draw a circle to approve the higher-pressure region in shale reservoir;
24) amplitude for carrying out three-dimensional earthquake data before superposition changes the synchronous wave impedance inversion with ripple in length and breadth with geophone offset;
25) elliptic inverse of the anisotropic parameters of three-dimensional earthquake data before superposition is carried out;
26) elliptic inverse of the modulus of elasticity of earthquake data before superposition is carried out, anisotropic elasticity modulus is obtained, passes through rock physicses
Analysis, anisotropic elasticity modulus is converted to the reservoir parameter of target zone;
27) to sign tomography and the joint geologic interpretation of the seismic properties in crack and demarcation;
28) according to rammell fracture development status, determine that possible completion formation damage area and fracturing fluid disturb the possibility of offset well
Property;
29) according to step 2) rock core dynamic and static modulus of elasticity conversion relational expression, by each of three-dimensional earthquake data before superposition
The dynamic modulus of elasticity that anisotropy elastic wave Simultaneous Retrieving is obtained is converted to static modulus of elasticity;
30) using the correlation of static modulus of elasticity and rock brittleness, the fragility regularity of distribution and feature of shale reservoir are determined, it is excellent
Change the completion and design of hydraulic fracturing scheme of horizontal well;
31) regularity of distribution using static modulus of elasticity or derivation static modulus of elasticity in shale reservoir, determines shale reservoir
Brittleness, obtain the orientation and intensity of local crustal stress, determine shale reservoir interrupting layer, crack and the azimuth tendency in crack
And dense degree, predict and draw a circle to approve the high formation pressure area in high total content of organic carbon and shale reservoir in shale reservoir;
32) parameter of the comprehensive shale gas reservoir obtained, with reference to accurate buried depth, thickness, occurrence and the planar distribution of shale reservoir,
Obtain the gas-bearing property prospect of shale gas reservoir and draw a circle to approve the dessert area of shale gas exploration and development.
2. method according to claim 1, feature is step 1) described in different directions be, level vertical with attitude of stratum or into
45 degree of angles.
3. method according to claim 1, feature is step 1) described in core column be 2.5 centimetres of diameter, 5 centimetres of length.
4. method according to claim 1, feature is step 3) described in optimal log be to eliminate the change of drilling internal diameter, well
Tiltedly change, well liquid changes, well temperature changes, logging speed is uneven, downhole instrument is stuck, instrument non-at the uniform velocity rotation and logging instrument
After device error component, it can truly reflect the optimal log of stratum physical property change.
5. method according to claim 1, feature is step 4) described in perturbation analysis be by changing formation fluid, porosity
Or the corresponding log obtained after lithology, find out changing rule of the correspondence log with different fluid, porosity or lithology.
6. method according to claim 1, feature is step 5) described in mineral be clay, it is calcite, quartz, pyrite, total
Organic carbon content and dolomite.
7. method according to claim 1, feature is step 5) described in optimal log be clay pit in log data
Thing curve, bulk density curve, stratum uranium content curve, neutron porosity curve, resistivity curve, compressional wave time difference curve and horizontal stroke
Ripple deviation curve.
8. method according to claim 1, feature is step 8) described in rock constituents perturbation analysis be by changing rock thing
The percentage composition of different minerals in model is managed, corresponding log is calculated, according to the log variable quantity calculated
Size, finds out the combination of corresponding Mineral change most sensitive property parameters or Sensitive Attributes parameter.
9. method according to claim 1, feature is step 9) described in parameter or parameter combination be modulus of elasticity, Young elasticity
Modulus, density, the coefficient of rigidity, the product of modulus of elasticity and density, the sum of products Young bullet of the coefficient of rigidity and density
The product of property modulus and density.
10. method according to claim 1, feature is step 14) described in top layer comprehensive modeling static correction be:At static correction
Reason, prestack denoising, amplitude compensation, Q value complements are repaid, surface consistent deconvolution and predictive deconvolution amplitude are with respect to fidelity processing.
11. method according to claim 1, feature is step 16) described in reflective information method of estimation be to be based on statistical signal
Self-adaptive processing, the high-resolution subsurface reflective information estimating method using nonparametric spectral analysis method and with fidelity, most
Keep to limits original seismic data information and do not lose in original data on the premise of small geological information, obtain high-resolution
Complex seismic trace collection.
12. method according to claim 1, feature is step 21) described in tomography pickup be similar based on coherent body, characteristic value
Property or curvature body calculate section automatically, determine macrocrack and craven fault.
13. method according to claim 1, feature is step 23) described in elliptical velocity inverting be to root-mean-square value speed
Bearing data body carries out elliptical velocity analysis, obtains fracture strike orientation and compressional wave anisotropic parameters.
14. method according to claim 1, feature is step 24) described in the synchronous wave impedance inversion of ripple in length and breadth be calculated amplitude
The gradient attribute changed with geophone offset, and inverting angular stack seismic data, synchronously obtain p-wave impedance, S-wave impedance and bullet
The product of property modulus and density, the product of the coefficient of rigidity and density, the product of Young's modulus of elasticity and density.
15. method according to claim 1, feature is step 25) described in elliptic inverse be that azimuthal gradient and speed are done
Elliptic inverse, to obtain Thomson parameter, is converted by rock physicses, by the geomechanics of layer for the purpose of Thomson Parameter Switch
Anisotropy parameter.
16. method according to claim 1, feature is step 26) described in reservoir parameter be rock brittleness, lithology, porosity,
Fluid, total content of organic carbon.
17. method according to claim 1, feature is step 27) described in joint geologic interpretation and demarcation reservoir petrologic characteristic
Parameter body log calibration, crack is demarcated with pit shaft imaging data and/or core analyzing data, large scale tomography and microcosmic
Tomography is demarcated with pressure break micro-seismic monitoring achievement and pit shaft imaging data, and stress anisotropy pressure break micro-seismic monitoring achievement is entered
The local demarcation of row, calibration process is contrasted with calculated value with measured result, finds out difference value between the two or phase relation
Number, then carries out the correction or correction of system to calculated value, to ensure calculated value and measurement result in the local eyeball in underground
Unanimously.
18. method according to claim 1, feature is step 30) described in optimum level well completion and design of hydraulic fracturing scheme
It is that horizontal well is laid in into fragility is higher and be easy in the shale containing high total organic carbon of pressure break, and optimization design each fracturing section
Spacing.
19. method according to claim 1, feature is step 32) described in parameter, include total content of organic carbon, the page of shale
Fragility, tomography, crack and the orientation in crack and density of rock reservoir, the orientation of local crustal stress and intensity, partial high pressure area and
Porosity is distributed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410140366.9A CN104977618B (en) | 2014-04-09 | 2014-04-09 | A kind of method evaluated shale gas reservoir and find dessert area |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410140366.9A CN104977618B (en) | 2014-04-09 | 2014-04-09 | A kind of method evaluated shale gas reservoir and find dessert area |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104977618A CN104977618A (en) | 2015-10-14 |
CN104977618B true CN104977618B (en) | 2017-09-22 |
Family
ID=54274293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410140366.9A Active CN104977618B (en) | 2014-04-09 | 2014-04-09 | A kind of method evaluated shale gas reservoir and find dessert area |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104977618B (en) |
Families Citing this family (76)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105651966A (en) * | 2016-01-18 | 2016-06-08 | 山东科技大学 | Shale oil and gas high-quality reservoir stratum evaluation method and parameter determination method |
CN105822295A (en) * | 2016-03-17 | 2016-08-03 | 成都创源油气技术开发有限公司 | Stratum fracture and cave development three-dimensional modeling method |
CN105822294A (en) * | 2016-03-17 | 2016-08-03 | 成都创源油气技术开发有限公司 | Shale gas reservoir geological evaluation single well modeling method |
CN105840175A (en) * | 2016-03-17 | 2016-08-10 | 成都创源油气技术开发有限公司 | Three-dimensional fracture modeling method |
CN105842751B (en) * | 2016-03-17 | 2019-02-22 | 成都创源油气技术开发有限公司 | Shale reservoir fracture evaluation method |
CN105700013B (en) * | 2016-04-26 | 2018-06-22 | 中国石油天然气集团有限公司 | The method for building shale gas dessert productivity model |
CN107436452A (en) * | 2016-05-27 | 2017-12-05 | 中国石油化工股份有限公司 | Hydrocarbon source rock Forecasting Methodology and device based on probabilistic neural network algorithm |
CN107451310B (en) * | 2016-05-31 | 2020-09-04 | 中国石油化工股份有限公司 | Classification evaluation method and device based on shale source-storage correlation |
CN107655789B (en) * | 2016-07-25 | 2020-03-20 | 中国石油化工股份有限公司 | Method for determining density of organic matter in shale |
CN106526669B (en) * | 2016-09-19 | 2018-10-23 | 中国石油化工股份有限公司 | A kind of Seismic Reservoir Prediction method of shale oil-gas reservoir |
CN106503834A (en) * | 2016-09-30 | 2017-03-15 | 中国石油天然气股份有限公司 | A kind of Forecasting Methodology in the fine and close oil dessert area of the ultralow porosity permeability reservoir of lacustrine facies |
CN106501858A (en) * | 2016-12-20 | 2017-03-15 | 南华大学 | Geophysics evaluation methodology is managed layer by layer in a kind of shale gas ground |
CN106855485B (en) | 2016-12-20 | 2019-08-06 | 中国石油天然气股份有限公司 | A kind of conversion method of sound state elastic parameter |
CN106772607A (en) * | 2017-01-19 | 2017-05-31 | 山西省煤炭地质物探测绘院 | A kind of method for predicting coal bed gas dessert |
CN106908854B (en) * | 2017-02-08 | 2019-01-18 | 中国石油天然气股份有限公司 | The determination method and apparatus in the dessert area of target area |
CN106761728B (en) * | 2017-02-14 | 2019-10-01 | 中国石油大学(北京) | A kind of recognition methods of the advantageous interval of marine facies shale formation |
CN107218032A (en) * | 2017-05-02 | 2017-09-29 | 中国石油大学(北京) | Compact reservoir sugariness evaluation method and device |
CN107143330B (en) * | 2017-05-25 | 2020-06-23 | 中石化石油工程技术服务有限公司 | Shale gas reservoir quality logging evaluation method |
CN109143368B (en) * | 2017-06-28 | 2020-04-07 | 中国石油化工股份有限公司 | Pre-mountain area earthquake denoising method |
CN109425896B (en) * | 2017-08-25 | 2021-04-30 | 中国石油天然气股份有限公司 | Dolomite oil and gas reservoir distribution prediction method and device |
US10928536B2 (en) * | 2017-12-07 | 2021-02-23 | Saudi Arabian Oil Company | Mapping chemostratigraphic signatures of a reservoir with rock physics and seismic inversion |
CN108171376B (en) * | 2017-12-27 | 2020-09-18 | 中国石油化工股份有限公司 | Total organic carbon prediction method and device, electronic equipment and storage medium |
CN108363100B (en) * | 2018-01-16 | 2020-05-08 | 太原理工大学 | Seismic geological identification method for coalbed methane dessert area based on sequestration unit and rough set |
CN108397130B (en) * | 2018-02-02 | 2019-11-08 | 中国石油天然气集团有限公司 | Boring method and device |
CN110486005B (en) * | 2018-05-14 | 2021-10-15 | 中国石油化工股份有限公司 | Method for identifying dessert layer of shale gas well |
CN108957534B (en) * | 2018-06-12 | 2020-02-14 | 中国石油天然气股份有限公司 | Method and device for predicting gas saturation |
CN109063232A (en) * | 2018-06-15 | 2018-12-21 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | One kind being used for shale gas well reservoir "sweet spot" evaluation method |
CN109113730B (en) * | 2018-07-12 | 2021-11-30 | 中国石油天然气股份有限公司 | Method, device and system for determining dessert region for shale oil in-situ conversion development |
CN108956952A (en) * | 2018-08-01 | 2018-12-07 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | The prediction technique and exploitation method of mud shale series of strata reservoir geology dessert between salt |
CN108594305A (en) * | 2018-08-13 | 2018-09-28 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Shale reservoir gas-bearing property analysis method and device |
CN109298448B (en) * | 2018-09-06 | 2020-08-21 | 中国海洋石油集团有限公司 | Prediction method and device for compact gas fracturing engineering dessert |
CN111123357B (en) * | 2018-10-31 | 2022-03-29 | 中国石油天然气股份有限公司 | Method and device for determining evaluation index of unconventional oil and gas reservoir |
CN109581531A (en) * | 2018-11-02 | 2019-04-05 | 中国石油天然气股份有限公司大港油田分公司 | A kind of unconventional oil and gas dessert quantitative evaluation method |
CN109492938B (en) * | 2018-12-04 | 2022-06-14 | 同济大学 | Dessert indicator factor-based deep carbonate reservoir quality evaluation method |
CN109887614B (en) * | 2019-01-22 | 2021-03-09 | 中国石油天然气股份有限公司 | Hydraulic fracture analysis method and device |
CN111679318B (en) * | 2019-03-11 | 2022-12-02 | 中国石油天然气股份有限公司 | Method and system for identifying fracture system in reservoir |
CN111898785A (en) * | 2019-05-05 | 2020-11-06 | 中国石油天然气集团有限公司 | Fracture toughness spatial distribution characteristic prediction method and system based on shale |
CN110068864B (en) * | 2019-06-03 | 2024-02-06 | 中铁西南科学研究院有限公司 | Method for detecting stratum cavity and non-compaction of tunnel |
CN110295882A (en) * | 2019-06-24 | 2019-10-01 | 中国石油化工股份有限公司 | Basin edge complicated structure band normal pressure shale gas horizontal well integrated design method |
CN112180443B (en) * | 2019-07-04 | 2024-03-01 | 中国石油天然气集团有限公司 | Shale gas two-dimensional seismic dessert area optimization method and device |
CN110618460A (en) * | 2019-07-22 | 2019-12-27 | 中国石油化工股份有限公司 | Micro-logging azimuth weighted interpolation modeling method combined with horizon information |
CN112305598A (en) * | 2019-07-29 | 2021-02-02 | 中国石油化工股份有限公司 | Complex geological special-shaped body reservoir prediction method, storage medium and computing equipment |
CN110488386B (en) * | 2019-09-20 | 2022-03-25 | 西南石油大学 | Anisotropic rock physical calibration method based on shale crystal geometric factor orientation function |
CN112578433A (en) * | 2019-09-27 | 2021-03-30 | 中国石油天然气集团有限公司 | Seismic data processing method and device |
CN112711065A (en) * | 2019-10-25 | 2021-04-27 | 中国石油天然气集团有限公司 | Pre-stack seismic inversion method and device |
CN113138412A (en) * | 2020-01-20 | 2021-07-20 | 中国石油天然气集团有限公司 | Deep shale porosity earthquake prediction method and device |
CN113138407A (en) * | 2020-01-20 | 2021-07-20 | 中国石油天然气集团有限公司 | Deep shale gas multi-scale fracture earthquake prediction method and system |
CN113267809A (en) * | 2020-02-17 | 2021-08-17 | 中国石油天然气集团有限公司 | Method and device for predicting I-type shale reservoir |
CN111624652B (en) * | 2020-04-17 | 2023-04-25 | 中国石油天然气股份有限公司 | Method for determining geological parameters in a formation |
CN113552621A (en) * | 2020-04-23 | 2021-10-26 | 中国石油天然气集团有限公司 | Shale gas ground stress determination method and device |
CN111506861B (en) * | 2020-04-24 | 2020-11-13 | 成都捷科思石油天然气技术发展有限公司 | Method for calculating crack strength of favorable region of target layer |
CN111411947B (en) * | 2020-05-14 | 2023-09-15 | 中国石油天然气集团有限公司 | Method for predicting formation pressure based on HTL parameters |
CN111608649B (en) * | 2020-05-21 | 2023-04-11 | 中国地质调查局成都地质调查中心 | Method for predicting beneficial area of exogenous supply type shale gas exploration |
CN113700466B (en) * | 2020-05-22 | 2023-10-31 | 中国石油天然气股份有限公司 | Method, device, equipment and medium for detecting oil gas in deep carbonate reservoir |
CN112051624B (en) * | 2020-06-23 | 2023-05-23 | 江苏华勘环保科技有限公司 | Three-dimensional visualization method for determining geological survey of solid waste site |
CN113866824A (en) * | 2020-06-30 | 2021-12-31 | 中国石油化工股份有限公司 | Method for determining elastic parameters of stratum between salts, storage medium and computer equipment |
CN111948718B (en) * | 2020-07-21 | 2023-11-28 | 中国石油天然气集团有限公司 | Method and device for predicting total organic carbon content of shale gas reservoir |
CN111852467B (en) * | 2020-07-28 | 2023-04-07 | 核工业北京地质研究院 | Method and system for delineating extension range of sandstone uranium ore body |
CN111897012A (en) * | 2020-08-12 | 2020-11-06 | 重庆地质矿产研究院 | Shale gas single-well geological comprehensive evaluation method |
CN112100906B (en) * | 2020-08-21 | 2021-06-22 | 中海石油(中国)有限公司 | Data-driven large-scale density modeling method, computing device and storage medium |
CN112505769B (en) * | 2020-11-25 | 2024-03-26 | 重庆地质矿产研究院 | Shale gas earthquake monitoring intelligent evaluation method based on dynamic geological engineering big data |
CN112505778B (en) * | 2020-12-01 | 2022-03-11 | 西南石油大学 | Three-dimensional in-situ characterization method for heterogeneity of shale storage and generation performance |
CN112859197B (en) * | 2021-03-31 | 2024-03-22 | 中国石油天然气集团有限公司 | Digital wellbore resistivity simulation method based on homogenized stratum electromagnetic field theory |
CN113253342A (en) * | 2021-05-11 | 2021-08-13 | 电子科技大学 | Method for constructing complex three-dimensional seismic model label by combining logging data |
CN113419284B (en) * | 2021-06-30 | 2022-02-08 | 东北石油大学 | Method for identifying physical facies double desserts of well logging rock based on cluster analysis |
CN113775330B (en) * | 2021-09-23 | 2023-12-19 | 中国石油集团东方地球物理勘探有限责任公司 | Comprehensive geophysical exploration system for high-temperature geothermal field and geothermal dessert area evaluation method |
CN114910499B (en) * | 2022-02-23 | 2023-09-26 | 河海大学 | Method for quantitatively predicting shale gas reservoir based on stratum difference diagenetic effect |
CN115079261B (en) * | 2022-06-06 | 2023-07-11 | 吉林大学 | Compact sandstone gas reservoir evaluation method based on multiparameter quantitative interpretation template |
CN115308796B (en) * | 2022-06-23 | 2023-07-07 | 吉林大学 | Shale gas-bearing and bedding structure prediction method based on anisotropic dispersion properties |
CN114994759B (en) * | 2022-08-02 | 2022-12-02 | 中国科学院地质与地球物理研究所 | Intelligent carbon seal storage box identification method and system based on GAN network |
CN115577616A (en) * | 2022-09-15 | 2023-01-06 | 福瑞升(成都)科技有限公司 | Carbonatite fracture-cave earthquake depicting method and device based on deep learning |
CN115266782B (en) * | 2022-09-27 | 2023-03-24 | 中国科学院地质与地球物理研究所 | Method for evaluating unconventional reservoir dual desserts based on dual-energy CT technology |
CN115291300B (en) * | 2022-09-28 | 2022-12-30 | 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) | Three-dimensional exploration analysis system and evaluation method based on complex geological conditions |
CN115877463B (en) * | 2022-10-27 | 2024-02-23 | 西南石油大学 | Shale gas reservoir TOC prediction method |
CN116068616B (en) * | 2023-03-06 | 2023-06-20 | 西安石油大学 | Reservoir abnormal formation pressure research method and system |
CN116663430B (en) * | 2023-07-27 | 2023-09-22 | 中国科学院地质与地球物理研究所 | Method and system for quantitatively representing porosity of high-temperature geothermal field reservoir structure |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1595201A (en) * | 2004-06-25 | 2005-03-16 | 大庆油田有限责任公司 | A method for improving seismic resolution |
WO2012036784A1 (en) * | 2010-09-13 | 2012-03-22 | Chevron U.S.A. Inc. | System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir |
CN102606151A (en) * | 2012-04-01 | 2012-07-25 | 中国石油大学(北京) | Method and device for predicting rock drillability of wildcat well before drilling |
CN102759748A (en) * | 2012-07-04 | 2012-10-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic reservoir prediction method by cross analysis based on target analysis |
-
2014
- 2014-04-09 CN CN201410140366.9A patent/CN104977618B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1595201A (en) * | 2004-06-25 | 2005-03-16 | 大庆油田有限责任公司 | A method for improving seismic resolution |
WO2012036784A1 (en) * | 2010-09-13 | 2012-03-22 | Chevron U.S.A. Inc. | System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir |
CN102606151A (en) * | 2012-04-01 | 2012-07-25 | 中国石油大学(北京) | Method and device for predicting rock drillability of wildcat well before drilling |
CN102759748A (en) * | 2012-07-04 | 2012-10-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic reservoir prediction method by cross analysis based on target analysis |
Non-Patent Citations (3)
Title |
---|
基于地球化学-地球物理的松辽盆地上白垩统油页岩识别与资源评价;贾建亮;《中国博士学位论文全文数据库 基础科学辑》;20120915(第9期);第24-110页 * |
非常规油气勘探、评价和开发新方法;王拓等;《岩性油气藏》;20131231;第25卷(第6期);第35-39、61页 * |
页岩气"甜点"评价与预测;周德华等;《石油实验地质》;20120331;第34卷(第2期);第109-114页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104977618A (en) | 2015-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104977618B (en) | A kind of method evaluated shale gas reservoir and find dessert area | |
Liu et al. | 3D geomechanical modeling and numerical simulation of in-situ stress fields in shale reservoirs: a case study of the lower Cambrian Niutitang formation in the Cen'gong block, South China | |
CN104853822A (en) | Method for evaluating shale gas reservoir and searching sweet spot region | |
Liu et al. | Simulation of paleotectonic stress fields and quantitative prediction of multi-period fractures in shale reservoirs: a case study of the Niutitang Formation in the Lower Cambrian in the Cen'gong block, South China | |
CN103454685B (en) | Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness | |
CN103256046B (en) | Unconventionaloil pool hides method and the device that horizontal well stitches the simulation of long fracturing parameter entirely | |
Castagna et al. | Relationships between compressional-wave and shear-wave velocities in clastic silicate rocks | |
Yaghoubi | Hydraulic fracturing modeling using a discrete fracture network in the Barnett Shale | |
Angus et al. | Integrated hydro-mechanical and seismic modelling of the Valhall reservoir: a case study of predicting subsidence, AVOA and microseismicity | |
AU2012260680A1 (en) | A method to aid in the exploration, mine design, evaluation and/or extraction of metalliferous mineral and/or diamond deposits | |
CN101231346A (en) | Method for estimating coal, rock mass physical mechanics parameter through seismic wave velocity | |
Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
Liu et al. | The in situ stress field and microscale controlling factors in the Ordos Basin, central China | |
CN106249300B (en) | It shakes to combine based on well and determines source rock TOC contents and anisotropism method and apparatus | |
Xia et al. | Application of 3D fine seismic interpretation technique in Dawangzhuang area, Bohai Bay Basin, Northeast China | |
Feng et al. | 3D numerical simulation of heterogeneous in situ stress field in low-permeability reservoirs | |
Li et al. | An Integrated quantitative modeling approach for fault-related fractures in tight sandstone reservoirs | |
Li et al. | Quantitative prediction of multi-period tectonic fractures based on integrated geological-geophysical and geomechanics data in deep carbonate reservoirs of Halahatang oilfield in northern Tarim Basin | |
Feng et al. | Development characteristics and quantitative prediction of multiperiod fractures in superdeep thrust-fold belt | |
CN104820237B (en) | The method for calculating the differential body on stratum | |
Chan | Subsurface geophysical characterization of the crystalline Canadian Shield in northeastern Alberta: implications for geothermal development | |
Chen et al. | Description technology of fractured-vuggy carbonate reservoir in halahatang oil field, tarim basin—take the ha 7 test area as an example | |
Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
Gu et al. | Investigation of fractures using seismic computerized crosshole tomography | |
Velez et al. | Horizontal Well Evaluation to Determine Geological Facies, Mechanical Properties, and Natural Fracture Changes Using Slim Through-the-Bit Dipole Sonic and Oil-Based Microimaging Tools. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |