CN106960076A - A kind of sedimentary facies prediction method based on hydrodynamics - Google Patents
A kind of sedimentary facies prediction method based on hydrodynamics Download PDFInfo
- Publication number
- CN106960076A CN106960076A CN201710111834.3A CN201710111834A CN106960076A CN 106960076 A CN106960076 A CN 106960076A CN 201710111834 A CN201710111834 A CN 201710111834A CN 106960076 A CN106960076 A CN 106960076A
- Authority
- CN
- China
- Prior art keywords
- section
- deposition
- sedimentary facies
- data
- deposit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of sedimentary facies prediction method based on hydrodynamics.It includes determining main control sandy ground quality factor;With stratum buried history, basin is recovered;According to the difference of particle mode of transport, deposit is divided into two classes;The deposition for calculating deposit is characterized according to the simulation of deposition component, simulation is repeated until simulation terminates;The panel data in 1 4400 moment is chosen, is analyzed while choosing a section, determines section sedimentary facies;Suitable source direction optionally takes a section perpendicular to source direction, and the data screening of section is extracted, the sedimentary facies of section is determined with reference to the feature of map;Determine the corresponding relation of each component ratio relation of section and sedimentary facies:According to the corresponding relation in the section of selection and the section SEDIMENTARY FACIES ANALYSIS section of determination between sand mud ratio and sedimentary facies;Determine planar depositions phase.Have the advantages that crystallizing field sedimentary facies can be quantitatively described.
Description
Technical field
The present invention relates to oil-gas exploration and development technical field, more specifically it is a kind of sedimentary facies based on hydrodynamics
Forecasting Methodology.
Background technology
In traditional geological research, often it is be combined with each other, is mutually authenticated using a variety of methods for sedimentary facies prediction, than
Such as geology, geophysics, geochemistry, ichnology method and facies marker quantitative analysis, so as to judge area's sedimentary facies classification;
But these methods can not isolate progress, the precision and accuracy of SEDIMENTARY FACIES ANALYSIS are otherwise reduced.Its cost of prior art is big,
Experimentation is big by human factor, it is impossible to press close to truth.
The content of the invention
It is an object of the invention to provide a kind of sedimentary facies prediction method based on hydrodynamics, crystallizing field sedimentary facies is carried out
Quantitative description, makes crystallizing field sedimentary facies more conform to actual geological process, meets the demand of geologist.
To achieve these goals, the technical scheme is that:A kind of sedimentary facies prediction method based on hydrodynamics,
It is characterized in that:Comprise the following steps:
Step 1:It is determined that mainly controlling sandy ground quality factor:From the sediment dynamics feature of sedimentary sand bodies, carry out sedimentology
Investigation, it is determined that the research fluctuations including landforms of area's sedimentary sand bodies, the tectonic geology feature of slope change and including
Flow, flow velocity, the Sediment dynamics governing factor of horizontal plane fluctuation;
Step 2:With stratum buried history, basin is recovered;On the basis of basin reconstruction, basin fill sand is set up
Body system dynamics simulation system, the shadow that source analyte, hydrodynamic condition parameter and construction characteristic parameter develop to sand-body distribution
Ring;
The Flow Field Distribution in the parallel finite element method analog study area of Navier-Stocks equations is first depending on, by flow field
Distributed data is analyzed with geologic information, adjusts hydrodynamic condition parameter, meets geological conditions;
Governing equation includes continuity equation, the equation of momentum;
Navier-Stocks equations are to accelerating including viscosity, density, gravity on the basis of continuity and the equation of momentum
The sign of the hydrodynamics influence factor of degree;
Continuity equation formula is as follows:
Wherein, ζ:Free Surface is to the distance of reference plane, m;
D:The depth of water below reference plane, m;
U、V:Average speed, m/s;
Q:Flow, m3/s;
Equation of momentum formula is as follows:
Wherein:ζ:Free Surface is to the distance of reference plane, m;
D:The depth of water below reference plane, m;
U、V:Average speed, m/s;
ρ:Fluid density, kg/m3;
P:Pressure, Pa;
H:The depth of water H=d+ ζ, m altogether;
E:The evaporation capacity of system, mm;
vV:Vertical vortex viscosity, m2/s;
fv、fu:Vertically, the coefficient of friction of horizontal direction;
Step 3:According to the difference of particle mode of transport, deposit is divided into Cohesive Sediment and the class of non-sticky silt two,
Deposition simulation is carried out using diffusion equation for Cohesive Sediment, diffusion equation formula is as follows:
Wherein:c:The concentration of deposition component, kg/m3;
u,v,ω:The velocity component of current, m/s;
εs,x, εs,y, εs,z:The diffusion coefficient of deposition component, m2/s;
ωs:The sinking speed of deposition component, m/s;
Deposition simulation is carried out using the settlement function of Van Rijn deposit for non-sticky silt, Van Rijn's is heavy
The settlement function formula of product thing is as follows:
Wherein:s:The relative density ρ of deposition components/ρw;
Ds:Deposition component representative diameter, μm;
υ:The kinematic viscosity of water body, m2/s;
Cohesive Sediment is set as a fixed value, only the deposition and the condition of degrading to Cohesive Sediment use restraint, and (this glues
The deposition of the silt of property degrades the value that academia does not determine, experience when using according to geologist is set);
Step 4:Characterized according to the simulation of deposition component in above-mentioned steps 3, deposition in analysis deposit, degrade, pass by
Do not stay the judgement of three kinds of sedimentary conditions, calculate the deposition of deposit, (ρ is density, m with density formula ρ=m/v of rock
It is volume for quality, v) volume of deposit is calculated, according to calculating the size of crystallizing field so as to calculate the deposit thickness of deposit,
Deposit thickness is added in initial bottom shape, the deposition bottom shape of subsequent time is used as;2 are repeated the above steps to step 4, until mould
Plan terminates, and realizes the quantitatively characterizing to sedimentary basin;
Step 5:The panel data in the 1-4400 moment is chosen, is analyzed while choosing a section, determines section
Sedimentary facies;Suitable source direction optionally takes a section perpendicular to source direction, and the data screening of section is extracted, made
The clay content broken line graph of section, sand content broken line graph, geomorphic feature distribution map, the thickness difference broken line graph at two moment, with reference to map
Curve fluctuations, the clay content corresponding to same distance, sand content, landforms elevation, thickness difference numerical value change feature determine
The sedimentary facies of section;
Step 6:Determine the corresponding relation of each (sand mud) component ratio relation of section and sedimentary facies:According to the section of selection
And the corresponding relation in the section SEDIMENTARY FACIES ANALYSIS section of determination between sand mud ratio and sedimentary facies, show that sand mud is with sinking in section
Corresponding relation between product phase;
Step 7:Determine planar depositions phase:Clay content is made in panel data at the time of by the above-mentioned steps 5 of selection
Distribution map, further according to the section sand mud proportionate relationship and the corresponding relation of sedimentary facies that are drawn in above-mentioned steps 6 and each map institute
The feature prediction planar depositions phase reflected.
In the above-mentioned technical solutions, in the step 1, from the sediment dynamics feature of sedimentary sand bodies, deposition is carried out
Investigation is learned, it is determined that the fluctuations including landforms, the tectonic geology feature of slope change and the bag of research area's sedimentary sand bodies
Include flow, flow velocity, the Sediment dynamics governing factor of horizontal plane fluctuation, obtain studying the geomorphologic map in area, the depositional model of deposit,
Hydrodynamic flow, flow velocity, research area's sedimentation time.
In the above-mentioned technical solutions, in step 3, on the basis of gridding is carried out to research area, according to point of individual well
Layer data, individual well SEDIMENTARY FACIES ANALYSIS performance data, individual well reservoir structure parameter, the reservoir cyclicity analyze data of individual well, individual well
Sediment composition analyze data, the sedimentary evolution history data of individual well, the paleostructure of individual well recover data, with the side of difference
Method, obtains the change number of terrain data, the component data in riverbed, the rugosity data in riverbed and the reference plane in analog study area
According to.
In the above-mentioned technical solutions, in step 3, the boundary condition of setting models selects suitable boundary condition, simulation water
Dynamic condition and petrofacies data with individual well, component analysis data carry out checking adjustment, finally obtains suitable hydrodynamic condition.
In the above-mentioned technical solutions, in step 3, according to the analyze data of the rock stratum of individual well, each coordinate points position is calculated
Silt content, the supply ratio to deposit in border enters row constraint, with deposition transmission equation, sedimentation equation simulation again
The deposition of existing deposit.
A kind of sedimentary facies prediction method based on hydrodynamics of the present invention includes hydrodynamic model and deposition mould
Three models of type and landform model,
Hydrodynamic model is passed according to the diffusion of convective-diffusion equation (N-S equations) simulated flow and the diffusion of deposit
It is defeated, carrying is transmitted to deposit on the basis of hydrodynamics,
Deposit transmission is mainly according to the different situations of silt, for Cohesive Sediment and the Bu Tong difference of non-sticky silt
The sedimentation of deposit is simulated with diffusion equation and sedimentation equation;
The type for determining sedimentary facies is the relative position of judgment standard face and deposition surface on the basis of consideration is hydrodynamic
Put, come judge deposit deposition whether, reach geologic simulation effect.
The invention has the advantages that:
(1) quantitative description is carried out to crystallizing field sedimentary facies, crystallizing field sedimentary facies is more conformed to actual geological process,
Meet the demand of geologist;
(2) different depositional environments can be simulated and analyzed;
(3) do not limited in simulation process by space, it is few by human factor in simulation process;
(4) result of simulation is the parameter of quantization, can be very good to carry out the Treatment Analysis in later stage;
(5) low cost of simulation, result can be debugged repeatedly.
Brief description of the drawings
Fig. 1 is embodiment of the present invention T=3500 moment simulation regions clay content broken line graph.
Fig. 2 is embodiment of the present invention T=3500 moment simulation regions amount broken line graph containing sand 1.
Fig. 3 is embodiment of the present invention T=3500 moment simulation regions amount broken line graph containing sand 2.
Fig. 4 is embodiment of the present invention T=3500 moment simulation regions section geomorphic feature distribution map.
Fig. 5 is the profile that the embodiment of the present invention is chosen along source direction.
Fig. 6 is embodiment of the present invention T=3500 moment simulation regions plane clay content distribution map.
Fig. 7 is the inventive method process chart.
Fig. 8 is the difference figure of embodiment of the present invention T=3500 moment and T=4000 moment sand thickness.
In Fig. 1, abscissa represents distance, and unit is km;Ordinate represents hundred shared by the mud of T=3500 moment simulation regions
Divide ratio;
In Fig. 2, abscissa represents distance, and unit is km;Ordinate is represented shared by the sand 1 of T=3500 moment simulation regions
Percentage;
In Fig. 3, abscissa represents distance, and unit is km;Ordinate is represented shared by the sand 2 of T=3500 moment simulation regions
Percentage;
In Fig. 4, abscissa represents distance, and unit is km;Ordinate represents elevation, and unit is m;
In Fig. 5, abscissa represents distance, and unit is km;Ordinate represents distance, and unit is km;
In Fig. 6, abscissa represents distance, and unit is m;Ordinate represents distance, and unit is m;
In Fig. 8, abscissa represents distance, and unit is km;Ordinate represents distance, and unit is m.
Embodiment
The performance that the invention will now be described in detail with reference to the accompanying drawings, but they do not constitute limitation of the invention, only
It is for example.Make advantages of the present invention more clear by explanation simultaneously and be readily appreciated that.
Understood refering to accompanying drawing:A kind of sedimentary facies prediction method based on hydrodynamics, it is characterised in that:Including following step
Suddenly:
Step 1:It is determined that mainly controlling sandy ground quality factor:From the sediment dynamics feature of sedimentary sand bodies, carry out sedimentology
Investigation, it is determined that the research fluctuations including landforms of area's sedimentary sand bodies, the tectonic geology feature of slope change and including
Flow, flow velocity, the Sediment dynamics governing factor of horizontal plane fluctuation;
Step 2:With stratum buried history, basin is recovered;On the basis of basin reconstruction, basin fill sand is set up
Body system dynamics simulation system, the shadow that source analyte, hydrodynamic condition parameter and construction characteristic parameter develop to sand-body distribution
Ring;
The Flow Field Distribution in the parallel finite element method analog study area of Navier-Stocks equations is first depending on, by flow field
Distributed data is analyzed with geologic information, adjusts hydrodynamic condition parameter, meets geological conditions;
Governing equation includes continuity equation, the equation of momentum;
Navier-Stocks equations are to accelerating including viscosity, density, gravity on the basis of continuity and the equation of momentum
The sign of the hydrodynamics influence factor of degree;
Continuity equation formula is as follows:
Wherein, ζ:Free Surface is to the distance of reference plane, m;
D:The depth of water below reference plane, m;
U、V:Average speed, m/s;
Q:Flow, m3/s;
Equation of momentum formula is as follows:
Wherein:ζ:Free Surface is to the distance of reference plane, m;
D:The depth of water below reference plane, m;
U、V:Average speed, m/s;
ρ:Fluid density, kg/m3;
P:Pressure, Pa;
H:The depth of water H=d+ ζ, m altogether;
E:The evaporation capacity of system, mm;
vV:Vertical vortex viscosity, m2/s;
fv、fu:Vertically, the coefficient of friction of horizontal direction;
Step 3:According to the difference of particle mode of transport, deposit is divided into Cohesive Sediment and the class of non-sticky silt two,
Deposition simulation is carried out using diffusion equation for Cohesive Sediment, diffusion equation formula is as follows:
Wherein:c:The concentration of deposition component, kg/m3;
u,v,ω:The velocity component of current, m/s;
εs,x, εs,y, εs,z:The diffusion coefficient of deposition component, m2/s;
ωs:The sinking speed of deposition component, m/s;
Deposition simulation is carried out using the settlement function of Van Rijn deposit for non-sticky silt, Van Rijn's is heavy
The settlement function formula of product thing is as follows:
Wherein:s:The relative density ρ of deposition components/ρw;
Ds:Deposition component representative diameter, μm;
υ:The kinematic viscosity of water body, m2/s;
Cohesive Sediment is set as a fixed value, only the deposition and the condition of degrading to Cohesive Sediment use restraint;
Step 4:Characterized according to the simulation of deposition component in above-mentioned steps 3, deposition in analysis deposit, degrade, pass by
The judgement of three kinds of sedimentary conditions is not stayed, the deposition of deposit is calculated, and deposit is calculated with density formula ρ=m/v of rock
Volume, according to calculating the size of crystallizing field so as to calculate the deposit thickness of deposit, deposit thickness is added to initial bottom shape
On, it is used as the deposition bottom shape of subsequent time;2 are repeated the above steps to step 4, until simulation terminates, is realized to sedimentary basin
Quantitatively characterizing;
Step 5:The panel data in the 1-4400 moment is chosen, is analyzed while choosing a section, determines section
Sedimentary facies;Suitable source direction optionally takes a section perpendicular to source direction, and the data screening of section is extracted, made
The clay content broken line graph of section, sand content broken line graph, geomorphic feature distribution map, the thickness difference broken line graph at two moment, with reference to map
Curve fluctuations, the clay content corresponding to same distance, sand content, landforms elevation, thickness difference numerical value change feature determine
The sedimentary facies of section;
Step 6:Determine the corresponding relation of each component ratio relation of section and sedimentary facies:According to the section of selection and determination
Section SEDIMENTARY FACIES ANALYSIS section in corresponding relation between sand mud ratio and sedimentary facies, draw in section sand mud and sedimentary facies it
Between corresponding relation;
Step 7:Determine planar depositions phase:Clay content is made in panel data at the time of by the above-mentioned steps 5 of selection
Distribution map, further according to the section sand mud proportionate relationship and the corresponding relation of sedimentary facies that are drawn in above-mentioned steps 6 and each map institute
The feature prediction planar depositions phase reflected.
In the step 1, from the sediment dynamics feature of sedimentary sand bodies, carry out sedimentology investigation, it is determined that grinding
Study carefully the fluctuations including landforms of area's sedimentary sand bodies, the tectonic geology feature of slope change and including flow, flow velocity, horizontal plane
The Sediment dynamics governing factor of fluctuation, obtain studying the geomorphologic map in area, the depositional model of deposit, hydrodynamic flow, flow velocity,
Study area's sedimentation time.
In step 3, on the basis of gridding is carried out to research area, according to the individual-layer data of individual well, individual well sedimentary facies
Analysis results data, individual well reservoir structure parameter, the reservoir cyclicity analyze data of individual well, the sediment composition of individual well analysis number
According to, the sedimentary evolution history data of individual well, the paleostructure of individual well recover data, with the method for difference, obtain analog study area
Terrain data, the component data in riverbed, the rugosity data in riverbed and the delta data of reference plane.
In step 3, the boundary condition of setting models selects suitable boundary condition, simulation hydrodynamic condition and and individual well
Petrofacies data, component analysis data carry out checking adjustment, finally obtain suitable hydrodynamic condition.
In step 3, according to the analyze data of the rock stratum of individual well, the content of the silt of each coordinate points position, opposite side are calculated
The supply ratio of deposit enters row constraint in boundary, and the deposition of deposit is reproduced with deposition transmission equation, sedimentation equation simulation.
Embodiment
Setting simulation region is 100m × 100m (long × wide), and the major parameter of simulation region is as shown in table 1 below:
The parameter setting of the simulated experiment of table 1
It can be seen from table 1:Simulate sedimentation time 1 year, calculate 6-7 days time, made using Delft3D QUICKPLOT
Go out T=3500 moment sedimentary section geomorphic feature distribution maps (as shown in Fig. 4) and profile (such as Fig. 5 chosen along source direction
It is shown);The related cross-sectional data at T=3500 moment is exported, the clay content broken line graph at T=3500 moment is made using EXCEL
(as shown in Figure 1), T=3500 moment amount broken line graph (as shown in Figure 2) containing sand 1, the T=3500 moment contains the amount broken line graph of sand 2 (such as
Shown in Fig. 3);The related cross-sectional data at T=4000 moment is exported, makes that T=3500 and T=4000 is same to be cutd open using Excel
The thick difference figure (as shown in Figure 8) of facing sand;The correlation plane data at T=3500 moment are exported, T=3500 is made using surfer
Moment plane clay content distribution map (as shown in Figure 6);The expression landforms that numerical value is more than 0 are drawn according to section geomorphic feature distribution map
On horizontal plane, i.e. delta leading edge parfacies or delta plain facies, numerical value are less than 0 expression under horizontal plane, i.e.,
Prodelta parfacies;Found out according to the thick poor figure of T=3500 and T=4000 same profiles sand when being represented when numerical value is more than 0 two
Sand body deposition is carved, sand body is degrading during less than 0 two moment of expression, then it is three that numerical value, which is more than the 0 sand body area of the main force for being delta,
Angle continent leading edge;This sedimentary subfacies is determined in conjunction with the distribution of T=3500 moment clay contents plans, it is known that more than clay content is
Prodelta parfacies, clay content is few for delta front sub-facies;Section and plan analyze the result come and combined i.e. respectively
It can determine that final sedimentary subfacies.
A kind of sedimentary facies prediction method based on hydrodynamics of the present invention compared with physical deposition simulated experiment, into
This (work station) it is smaller, experiment parameter can be adjusted repeatedly, experimentation is small by human factor, experimental situation closer to
Truth, the result of simulation is the parameter of quantization, can be very good to carry out the Treatment Analysis in later stage, the cost of simulation compared with
Low, result can be debugged repeatedly.
Other unaccounted parts belong to prior art.
Claims (5)
1. a kind of sedimentary facies prediction method based on hydrodynamics, it is characterised in that:Comprise the following steps:
Step 1:It is determined that mainly controlling sandy ground quality factor:From the sediment dynamics feature of sedimentary sand bodies, carry out sedimentology investigation
Research, it is determined that the research fluctuations including landforms of area's sedimentary sand bodies, the tectonic geology feature of slope change and including flow,
Flow velocity, the Sediment dynamics governing factor of horizontal plane fluctuation;
Step 2:With stratum buried history, basin is recovered;On the basis of basin reconstruction, basin fill sand body body is set up
It is dynamics simulation system, the influence that source analyte, hydrodynamic condition parameter and construction characteristic parameter develop to sand-body distribution;
The Flow Field Distribution in the parallel finite element method analog study area of Navier-Stocks equations is first depending on, by Flow Field Distribution
Data are analyzed with geologic information, adjust hydrodynamic condition parameter, meet geological conditions;
Governing equation includes continuity equation, the equation of momentum;
Navier-Stocks equations are to including viscosity, density, acceleration of gravity on the basis of continuity and the equation of momentum
The sign of hydrodynamics influence factor;
Continuity equation formula is as follows:
Wherein, ζ:Free Surface is to the distance of reference plane, m;
D:The depth of water below reference plane, m;
U、V:Average speed, m/s;
Q:Flow, m3/s;
Equation of momentum formula is as follows:
Wherein:ζ:Free Surface is to the distance of reference plane, m;
D:The depth of water below reference plane, m;
U、V:Average speed, m/s;
ρ:Fluid density, kg/m3;
P:Pressure, Pa;
H:The depth of water H=d+ ζ, m altogether;
E:The evaporation capacity of system, mm;
vV:Vertical vortex viscosity, m2/s;
fv、fu:Vertically, the coefficient of friction of horizontal direction;
Step 3:According to the difference of particle mode of transport, deposit is divided into Cohesive Sediment and the class of non-sticky silt two,
Deposition simulation is carried out using diffusion equation for Cohesive Sediment, diffusion equation formula is as follows:
Wherein:c:The concentration of deposition component, kg/m3;
u,v,ω:The velocity component of current, m/s;
εs,x, εs,y, εs,z:The diffusion coefficient of deposition component, m2/s;
ωs:The sinking speed of deposition component, m/s;
Deposition simulation, Van Rijn deposit are carried out using the settlement function of Van Rijn deposit for non-sticky silt
Settlement function formula it is as follows:
Wherein:s:The relative density ρ of deposition components/ρw;
Ds:Deposition component representative diameter, μm;
υ:The kinematic viscosity of water body, m2/s;
Cohesive Sediment is set as a fixed value, only the deposition and the condition of degrading to Cohesive Sediment use restraint;
Step 4:Characterized according to the simulation of deposition component in above-mentioned steps 3, deposition in analysis deposit, degrade, pass by and do not stay
The judgement of three kinds of sedimentary conditions, calculates the deposition of deposit, and the body of deposit is calculated with density formula ρ=m/v of rock
Product, according to the size of crystallizing field is calculated so as to calculate the deposit thickness of deposit, deposit thickness is added in initial bottom shape, made
For the deposition bottom shape of subsequent time;2 are repeated the above steps to step 4, until simulation terminates, the quantitative table to sedimentary basin is realized
Levy;
Step 5:The panel data in the 1-4400 moment is chosen, is analyzed while choosing a section, determines that section is deposited
Phase;Suitable source direction optionally takes a section perpendicular to source direction, and the data screening of section is extracted, section is made
Clay content broken line graph, sand content broken line graph, geomorphic feature distribution map, the thickness difference broken line graph at two moment, with reference to the song of map
Line fluctuations, the clay content corresponding to same distance, sand content, landforms elevation, thickness difference numerical value change feature determines section
Sedimentary facies;
Step 6:Determine the corresponding relation of each component ratio relation of section and sedimentary facies:According to the section of selection and cuing open for determining
Corresponding relation in the SEDIMENTARY FACIES ANALYSIS section of face between sand mud ratio and sedimentary facies, draws in section between sand mud and sedimentary facies
Corresponding relation;
Step 7:Determine planar depositions phase:Point of clay content is made in panel data at the time of by the above-mentioned steps 5 of selection
Butut, is reflected further according to the section sand mud proportionate relationship and the corresponding relation and each map of sedimentary facies drawn in above-mentioned steps 6
The feature prediction planar depositions phase gone out.
2. a kind of sedimentary facies prediction method based on hydrodynamics according to claim 1, it is characterised in that:The step
In 1, from the sediment dynamics feature of sedimentary sand bodies, carry out sedimentology investigation, it is determined that the bag of research area's sedimentary sand bodies
Include the fluctuations of landforms, the tectonic geology feature of slope change and the Sediment dynamics control fluctuated including flow, flow velocity, horizontal plane
Factor processed, obtains studying the geomorphologic map in area, the depositional model of deposit, hydrodynamic flow, flow velocity, research area's sedimentation time.
3. a kind of sedimentary facies prediction method based on hydrodynamics according to claim 1 or 2, it is characterised in that:In step
In rapid 3, to research area carry out gridding on the basis of, according to the individual-layer data of individual well, individual well SEDIMENTARY FACIES ANALYSIS performance data,
Individual well reservoir structure parameter, the reservoir cyclicity analyze data of individual well, the sediment composition analyze data of individual well, the deposition of individual well
Evolutionary History data, the paleostructure of individual well recover data, with the method for difference, obtain terrain data, the riverbed in analog study area
Component data, the rugosity data in riverbed and the delta data of reference plane.
4. a kind of sedimentary facies prediction method based on hydrodynamics according to claim 3, it is characterised in that:In step 3,
The boundary condition of setting models, selects suitable boundary condition, simulates hydrodynamic condition and petrofacies data with individual well, component point
Analyse data and carry out checking adjustment, finally obtain suitable hydrodynamic condition.
5. a kind of sedimentary facies prediction method based on hydrodynamics according to claim 4, it is characterised in that:In step 3,
According to the analyze data of the rock stratum of individual well, the content of the silt of each coordinate points position, the supply to deposit in border are calculated
Ratio enters row constraint, and the deposition of deposit is reproduced with deposition transmission equation, sedimentation equation simulation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111834.3A CN106960076A (en) | 2017-02-28 | 2017-02-28 | A kind of sedimentary facies prediction method based on hydrodynamics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111834.3A CN106960076A (en) | 2017-02-28 | 2017-02-28 | A kind of sedimentary facies prediction method based on hydrodynamics |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106960076A true CN106960076A (en) | 2017-07-18 |
Family
ID=59469960
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710111834.3A Pending CN106960076A (en) | 2017-02-28 | 2017-02-28 | A kind of sedimentary facies prediction method based on hydrodynamics |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106960076A (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107577877A (en) * | 2017-09-07 | 2018-01-12 | 长江大学 | Deposition method for numerical simulation based on natural material resource |
CN107609253A (en) * | 2017-09-07 | 2018-01-19 | 长江大学 | The deposition method for numerical simulation of carbonate |
CN107607434A (en) * | 2017-09-15 | 2018-01-19 | 湖南军信环保股份有限公司 | A kind of sludge sand content detection method |
CN108229043A (en) * | 2018-01-12 | 2018-06-29 | 中国海洋大学 | Consider the deep-sea SPAR type floating wind turbine Analysis of Fatigue methods of vortex-induced effect |
CN109242952A (en) * | 2018-08-14 | 2019-01-18 | 长江大学 | A kind of husky body three-dimensional modeling method and system based on double material resource sink sedimentation experiments |
CN110609319A (en) * | 2019-08-20 | 2019-12-24 | 中国石油大学(华东) | Method for establishing time domain sequence stratigraphic section by using forward modeling of sedimentation |
CN111985124A (en) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | Single sand body model construction method, system, equipment and readable storage medium |
CN111985125A (en) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | Single sand body modeling method, system, equipment and readable storage medium |
CN112149312A (en) * | 2020-10-13 | 2020-12-29 | 中国石油天然气股份有限公司 | Method and device for determining deposition microphase, electronic equipment and storage medium |
CN112782759A (en) * | 2020-01-20 | 2021-05-11 | 中国石油天然气集团有限公司 | Prediction method of complex reservoir |
CN113295849A (en) * | 2021-05-21 | 2021-08-24 | 新疆大学 | Method for determining ancient wind direction of continental lake basin through distribution pattern of deposition system |
CN113309511A (en) * | 2021-07-13 | 2021-08-27 | 西安石油大学 | Method for acquiring compact oil accumulation power |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5844799A (en) * | 1996-01-26 | 1998-12-01 | Institut Francais Du Petrole | Method for simulating the filling of a sedimentary basin |
CN105701319A (en) * | 2016-03-08 | 2016-06-22 | 中国石油大学(华东) | Sedimentary microfacies modeling method under horizontal well rule development well net |
CN106202749A (en) * | 2016-07-15 | 2016-12-07 | 长江大学 | Deposition analogy method deposition process controlled based on datum level |
-
2017
- 2017-02-28 CN CN201710111834.3A patent/CN106960076A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5844799A (en) * | 1996-01-26 | 1998-12-01 | Institut Francais Du Petrole | Method for simulating the filling of a sedimentary basin |
CN105701319A (en) * | 2016-03-08 | 2016-06-22 | 中国石油大学(华东) | Sedimentary microfacies modeling method under horizontal well rule development well net |
CN106202749A (en) * | 2016-07-15 | 2016-12-07 | 长江大学 | Deposition analogy method deposition process controlled based on datum level |
Non-Patent Citations (1)
Title |
---|
王杨君 等: "水动力数值模拟的河控三角洲分支河道演化研究", 《地质科技情报》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107577877A (en) * | 2017-09-07 | 2018-01-12 | 长江大学 | Deposition method for numerical simulation based on natural material resource |
CN107609253A (en) * | 2017-09-07 | 2018-01-19 | 长江大学 | The deposition method for numerical simulation of carbonate |
CN107609253B (en) * | 2017-09-07 | 2021-03-09 | 长江大学 | Carbonate rock deposition numerical simulation method |
CN107607434A (en) * | 2017-09-15 | 2018-01-19 | 湖南军信环保股份有限公司 | A kind of sludge sand content detection method |
CN108229043A (en) * | 2018-01-12 | 2018-06-29 | 中国海洋大学 | Consider the deep-sea SPAR type floating wind turbine Analysis of Fatigue methods of vortex-induced effect |
CN108229043B (en) * | 2018-01-12 | 2021-05-11 | 中国海洋大学 | Deep sea SPAR type floating fan fatigue damage analysis method considering vortex-induced effect |
CN109242952A (en) * | 2018-08-14 | 2019-01-18 | 长江大学 | A kind of husky body three-dimensional modeling method and system based on double material resource sink sedimentation experiments |
CN109242952B (en) * | 2018-08-14 | 2024-04-16 | 长江大学 | Sand body three-dimensional modeling method and system based on double-object-source water tank deposition experiment |
CN110609319A (en) * | 2019-08-20 | 2019-12-24 | 中国石油大学(华东) | Method for establishing time domain sequence stratigraphic section by using forward modeling of sedimentation |
CN110609319B (en) * | 2019-08-20 | 2021-12-10 | 中国石油大学(华东) | Method for establishing time domain sequence stratigraphic section by using forward modeling of sedimentation |
CN112782759A (en) * | 2020-01-20 | 2021-05-11 | 中国石油天然气集团有限公司 | Prediction method of complex reservoir |
CN111985125A (en) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | Single sand body modeling method, system, equipment and readable storage medium |
CN111985124B (en) * | 2020-07-15 | 2023-07-07 | 北京金阳普泰石油技术股份有限公司 | Single sand body model construction method, system, equipment and readable storage medium |
CN111985125B (en) * | 2020-07-15 | 2023-08-11 | 北京金阳普泰石油技术股份有限公司 | Single sand modeling method, system, equipment and readable storage medium |
CN111985124A (en) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | Single sand body model construction method, system, equipment and readable storage medium |
CN112149312A (en) * | 2020-10-13 | 2020-12-29 | 中国石油天然气股份有限公司 | Method and device for determining deposition microphase, electronic equipment and storage medium |
CN112149312B (en) * | 2020-10-13 | 2024-01-30 | 中国石油天然气股份有限公司 | Method and device for determining deposition microphase, electronic equipment and storage medium |
CN113295849A (en) * | 2021-05-21 | 2021-08-24 | 新疆大学 | Method for determining ancient wind direction of continental lake basin through distribution pattern of deposition system |
CN113295849B (en) * | 2021-05-21 | 2022-12-13 | 新疆大学 | Method for determining ancient wind direction of continental lake basin through distribution pattern of deposition system |
CN113309511A (en) * | 2021-07-13 | 2021-08-27 | 西安石油大学 | Method for acquiring compact oil accumulation power |
CN113309511B (en) * | 2021-07-13 | 2023-08-25 | 西安石油大学 | Method for obtaining tight oil reservoir power |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106960076A (en) | A kind of sedimentary facies prediction method based on hydrodynamics | |
CN106202749B (en) | The deposition analogy method that deposition process is controlled based on datum level | |
US9514096B2 (en) | Method of developing a sedimentary basin from a stratigraphic simulation of multilithologic filling taking account of fine sediment transport | |
Lunt et al. | Deposits of the sandy braided South Saskatchewan River: Implications for the use of modern analogs in reconstructing channel dimensions in reservoir characterization | |
Issautier et al. | Impacts of fluvial reservoir heterogeneity on connectivity: implications in estimating geological storage capacity for CO2 | |
CN106959261A (en) | It is a kind of to predict sedimentary particle distribution and the method for proportioning | |
Balázs et al. | Tectonic and climatic controls on asymmetric half‐graben sedimentation: Inferences from 3‐D numerical modeling | |
CN104698492A (en) | Method for calculating abnormal stratum pressure | |
Trabelsi et al. | GIS-based subsurface databases and 3-D geological modeling as a tool for the set up of hydrogeological framework: Nabeul–Hammamet coastal aquifer case study (Northeast Tunisia) | |
CN105137482A (en) | Method for calculating ancient gradient of sedimentary body | |
Ige et al. | Evaluation of aquifer hydraulic characteristics using geoelectrical sounding, pumping and laboratory tests: A case study of Lokoja and Patti Formations, Southern Bida Basin, Nigeria | |
Koneshloo et al. | The impact of geological uncertainty on primary production from a fluvial reservoir | |
Lastras et al. | Conceptual and numerical modeling of the BIG'95 debris flow, western Mediterranean Sea | |
Kurah et al. | Reservoir characterization and volumetric estimation of reservoir fluids using simulation and analytical methods: a case study of the coastal swamp depobelt, Niger Delta Basin, Nigeria | |
Changmin et al. | Hierarchical architectural element analysis for sandy conglomerate deposits of Baikouquan Formation, Mahu Sag | |
CN104698494A (en) | Method for calculating abnormal stratum pressure | |
Yin et al. | Characterization of marine shale in Western Hubei Province based on unmanned aerial vehicle oblique photographic data. | |
Yuan et al. | Utilizing Integrated Prediction Error Filter Analysis (INPEFA) to divide base-level cycle of fan-deltas: a case study of the Triassic Baikouquan Formation in Mabei slope area, Mahu Depression, Junggar Basin, China | |
Qi et al. | Hierarchy and subsurface correlation of muddy baffles in lacustrine delta fronts: a case study in the X Oilfield, Subei Basin, China | |
Yong et al. | Application of stratigraphic-sedimentological forward modeling of sedimentary processes to predict high-quality reservoirs within tight sandstone | |
George et al. | Estimation of aquifer hydraulic parameters via complementing surfacial geophysical measurement by laboratory measurements on the aquifer core samples | |
Farhat et al. | Structural control of hydrogeological setting revealed by gravity and Time Domain Electromagnetic data analysis: the Moulares basin (Southwestern Tunisia) | |
CN106919752A (en) | Based on the lake basin landforms quantitative forecasting technique that flow field characterizes | |
Zhang et al. | Sedimentology and sequence stratigraphy of the second member of Shuangyang Formation, Y45 Block, Moliqing oilfield, Yitong Basin, China | |
Abu-Hashish et al. | 3D geological modeling of the Upper Cretaceous reservoirs in GPT oil field, Abu Sennan area, Western Desert, Egypt |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170718 |