CN108629463A - Crustal stress prediction technique and device - Google Patents
Crustal stress prediction technique and device Download PDFInfo
- Publication number
- CN108629463A CN108629463A CN201810501796.7A CN201810501796A CN108629463A CN 108629463 A CN108629463 A CN 108629463A CN 201810501796 A CN201810501796 A CN 201810501796A CN 108629463 A CN108629463 A CN 108629463A
- Authority
- CN
- China
- Prior art keywords
- crustal stress
- well
- work area
- prediction time
- stress
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000012530 fluid Substances 0.000 claims abstract description 43
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 36
- 239000011148 porous material Substances 0.000 claims abstract description 36
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 34
- 238000002347 injection Methods 0.000 claims description 25
- 239000007924 injection Substances 0.000 claims description 25
- 238000009826 distribution Methods 0.000 claims description 12
- 239000007788 liquid Substances 0.000 claims description 9
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 3
- 239000003129 oil well Substances 0.000 claims 1
- 238000011161 development Methods 0.000 description 20
- 238000005516 engineering process Methods 0.000 description 7
- 230000006872 improvement Effects 0.000 description 7
- 230000009467 reduction Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 239000004744 fabric Substances 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/003—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by analysing drilling variables or conditions
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Business, Economics & Management (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Geochemistry & Mineralogy (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Marketing (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Analytical Chemistry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of crustal stress prediction technique of this specification embodiment offer and device.The method includes：The initial three-dimensional crustal stress distributed data in work area is provided；Receive prediction time and creation data；Based on the creation data, determine the well point position in the work area the prediction time formation pore fluid pressure；Based on the formation pore fluid pressure and the initial three-dimensional crustal stress distributed data, determine the well point position in the work area the prediction time crustal stress；Based on crustal stress of the well point position in the prediction time in the work area, determine crustal stress of the position in the prediction time between the well in the work area, so obtain the work area the prediction time three-dimensional ground stress distributed data.
Description
Technical field
This specification embodiment is related to technical field of geophysical exploration, more particularly to a kind of crustal stress prediction technique and dress
It sets.
Background technology
Crustal stress is that earth solid dielectric is additional by gravity, a variety of earth structure power, astronomical factor and spy pick engineering
The effect of power causes the mechanics parameter that response deforms in media interior unit.Accurate three-dimensional ground stress distributed data is oil
The important geologic basis of gas field development management and development plan adjustment.
Currently, being typically based on original crustal stress distributed data to instruct oil-gas field development and management, the original ground
Stress distribution data are usually crustal stress distributed data of the oil gas field before development of injection-production.But in the development of injection-production of oil gas field
In the process, water filling and oil recovery all can cause the crustal stress of oil gas field to change.Oil gas field the different development times dimensionally
Stress distribution is can be changed.Current technology does not account for reservoir engineering factor (such as water filling is recovered the oil) to crustal stress
The influence of changes in distribution, and then oil-gas field development management and the adjustment of development plan can not be instructed well.
Invention content
The purpose of this specification embodiment is to provide a kind of crustal stress prediction technique and device, with accurately predicting oil/gas field
It is distributed in the three-dimensional ground stress of different development times.
In order to achieve the above object, this specification embodiment provides a kind of crustal stress prediction technique, including：The original in work area is provided
Beginning three-dimensional ground stress distributed data；Receive prediction time and creation data；Based on the creation data, determine in the work area
Formation pore fluid pressure of the well point position in the prediction time；Based on the formation pore fluid pressure and described original three
Tie up crustal stress distributed data, determine the well point position in the work area the prediction time crustal stress；Based on the work area
Crustal stress of the interior well point position in the prediction time, position is on the ground of the prediction time between determining the well in the work area
Stress, so obtain the work area the prediction time three-dimensional ground stress distributed data.
In order to achieve the above object, this specification embodiment provides a kind of crustal stress prediction meanss, including：Unit is provided, is used
In the initial three-dimensional crustal stress distributed data for providing work area；Receiving unit, for receiving prediction time and creation data；First really
Order member, for be based on the creation data, determine the well point position in the work area the prediction time formation pore
Fluid pressure；Second determination unit is used for the formation pore fluid pressure and the initial three-dimensional crustal stress distributed data, really
Crustal stress of the well point position in the work area in the prediction time calmly；Third determination unit, for based in the work area
Crustal stress of the well point position in the prediction time, position is answered on the ground of the prediction time between determining the well in the work area
Power, so obtain the work area the prediction time three-dimensional ground stress distributed data.
The technical solution provided by above this specification embodiment can provide work area as it can be seen that in this specification embodiment
Initial three-dimensional crustal stress distributed data；Can determine the well point position in work area prediction time formation pore fluid pressure
Power；Can be based on the formation pore fluid pressure, determine the well point position in work area the prediction time crustal stress；It can
With based on crustal stress of the well point position in the prediction time in the work area, position is in the prediction time between determining well
Crustal stress, so obtain the work area the prediction time three-dimensional ground stress distributed data.The ground of this specification embodiment
Stress prediction method, it is contemplated that influence of the factor of production to crustal stress so that the crustal stress of prediction is distributed more accurate, Jin Erke
To instruct the adjustment of oil-gas field development management and development plan well.
Description of the drawings
In order to illustrate more clearly of this specification embodiment or technical solution in the prior art, below will to embodiment or
Attached drawing needed to be used in the description of the prior art is briefly described, it should be apparent that, the accompanying drawings in the following description is only
Some embodiments described in this specification, for those of ordinary skill in the art, in not making the creative labor property
Under the premise of, other drawings may also be obtained based on these drawings.
Fig. 1 is a kind of flow chart of crustal stress prediction technique of this specification embodiment；
Relational graphs of the Fig. 2 between a kind of crustal stress of this specification embodiment and water injection rate；
Relational graphs of the Fig. 3 between a kind of crustal stress of this specification embodiment and Liquid output；
Fig. 4 is a kind of primitively stress distribution schematic diagram of oil gas field of this specification embodiment；
Fig. 5 is the crustal stress distribution schematic diagram after 20 years of a kind of oil gas field of this specification embodiment；
Fig. 6 is a kind of illustrative view of functional configuration of crustal stress prediction meanss of this specification embodiment.
Specific implementation mode
Below in conjunction with the attached drawing in this specification embodiment, the technical solution in this specification embodiment is carried out clear
Chu is fully described by, it is clear that described embodiment is only this specification a part of the embodiment, rather than whole implementation
Example.The embodiment of base in this manual, those of ordinary skill in the art are obtained without creative efforts
Every other embodiment, should all belong to this specification protection range.
It please refers to Fig.1.This specification embodiment provides a kind of crustal stress prediction technique.The crustal stress prediction technique is to take
Business device is executive agent, be may comprise steps of.
Step S10：The initial three-dimensional crustal stress distributed data in work area is provided.
In the present embodiment, the region that the work area can be constructed with fingering row, specifically such as can be oil gas field.It is described
Crustal stress includes but not limited to vertical stress, horizontal maximum principal stress, horizontal minimum principal stress etc..The initial three-dimensional crustal stress
Distributed data can be the three-dimensional ground stress distributed data before note is adopted.
In the present embodiment, the initial three-dimensional crustal stress distributed data may include three-dimensional ground stress distributed model.Institute
It may include original crustal stress and the correspondence of spatial position to state three-dimensional ground stress distributed model.The original crustal stress can be with
Including crustal stress between original well and the crustal stress at original shaft location.The pit shaft may include water injection well and producing well etc..
In the present embodiment, the server can establish the three-dimensional on stratum using drilling well, well logging and three dimensional seismic data
Skeleton Model；Three dimensional seismic data inverting and well-log information can be utilized to establish three-dimensional petrophysical parameter model；Can carry out
The three-dimensional finite element simulation of crustal stress space distribution rule calculates；It can be based on the three-dimensional Skeleton Model, the three-dimensional
Petrophysical parameter model and three-dimensional finite element simulation result of calculation build three-dimensional ground stress distributed model.
Specifically, the server can utilize drilling well, well logging and three dimensional seismic data, compare and construct in detail stratigraphic division
On the basis of explanation, the three-dimensional Skeleton Model on stratum is established；It can be in the test of core sample static state petrophysical parameter and dynamic
On the basis of petrophysical parameter is tested and its compared and corrects, in conjunction with the three dimensional seismic data inverting under Log-constrained, utilize
Three dimensional seismic data calculates the distributed in three dimensions of rock mechanics parameters (such as elasticity modulus, Poisson's ratio etc.), obtains three-dimensional rock object
Manage parameter model；Can be carried out using well-log information individual well crustal stress (such as vertical stress, horizontal maximum principal stress, it is horizontal most
Minor principal stress etc.) it calculates and corrects, and then using the individual well crustal stress result of calculation after correction as constraint, utilize three-dimensional finite element
Numerical simulation technology calculates the crustal stress three-dimensional distribution rule on stratum；It can be on the basis of three-dimensional Skeleton Model, to log well ground
Stress calculating results are controlled as well point, and the three-dimensional ground stress calculated using three-dimensional finite element simulation technology is distributed as between well
Constraint, the method being combined using Decided modelling and stochastic modeling establish three-dimensional ground stress distributed model.
Step S12：Receive prediction time and creation data.
In the present embodiment, the prediction time can be any time of the oil gas field in the development phase.In view of oil gas field
Original crustal stress can just change during development of injection-production, the prediction time be typically larger than note adopts the moment.The work
May include at least one water injection well and at least one producing well in area.The creation data may include the note of each water injection well
The Liquid output of water and each producing well.
In the present embodiment, user directly can input the prediction time and the creation data in the server.
The server can receive the prediction time and the creation data.Certainly, user can also use other manner to input
The prediction time and the creation data, the present embodiment is to this and is not specifically limited.
Step S14：Based on the creation data, determine the well point position in the work area on the stratum of the prediction time
Pore fluid pressure.
In the present embodiment, the quantity of well point position can be one or more in the work area.Well in the work area
Point position can specifically include the well point position of water injection well and the well point position of producing well.
In the present embodiment, formation pore fluid pressure would generally be adopted with the note of oil gas field and be changed.In this way, institute
State server can use following formula calculate the well point position in the work area the prediction time formation pore fluid
Pressure.
In above formula (1), P can be the formation pore fluid pressure of prediction time；PwIt can be flowing bottomhole pressure (FBHP)；Q can be with
For the water injection rate of water injection well or the Liquid output of producing well；μ can be viscosity；B can be volume factor；K can be that stratum is effective
Permeability；H can be Effective thickness of formation；Ф can be stratum effecive porosity；C can be system compressibility；A can be
Supply area；Δ t can adopt the time interval between the moment for prediction time and note.
For example, for each water injection well in the work area, the server can calculate this using above-mentioned formula (1)
Formation pore fluid pressure of the well point position of water injection well in the prediction time.It should be noted that calculating water injection well
In the formation pore fluid pressure of the prediction time, the Q in above-mentioned formula (1) can be the water filling of water injection well for well point position
Amount.
As a further example, for each producing well in the work area, the server can be counted using above-mentioned formula (1)
Calculate formation pore fluid pressure of the well point position in the prediction time of the producing well.It should be noted that recovering the oil calculating
In the formation pore fluid pressure of the prediction time, the Q in above-mentioned formula (1) can be producing well for the well point position of well
Liquid output.
Step S16：Based on the formation pore fluid pressure and the initial three-dimensional crustal stress distributed data, determine described in
The crustal stress of well point position in work area in the prediction time.
In the present embodiment, water filling and oil recovery can influence point of crustal stress by influencing formation pore fluid pressure
Cloth.In this way, the server can calculate the well point position in the work area on the ground of the prediction time using following formula
Stress.It should be noted that the well point may include the well point of water injection well and the well point of producing well；The well point position can be with
The well point position of well point position and producing well including water injection well.
In above formula (2), σ1Can be crustal stress of the well point position in the prediction time；σ0It can be the original of well point position
Beginning crustal stress can specifically be obtained by the initial three-dimensional crustal stress distributed data in step S10；P can be prediction time
Formation pore fluid pressure；P0It can be prime stratum pore fluid pressure；α can be Biot coefficients；υ can be Poisson's ratio.
For example, for each water injection well in the work area, the server can calculate this using above-mentioned formula (2)
Crustal stress of the well point position of water injection well in the prediction time.
As a further example, for each producing well in the work area, the server can be counted using above-mentioned formula (2)
Calculate crustal stress of the well point position in the prediction time of the producing well.
Step S18：Based on crustal stress of the well point position in the prediction time in the work area, determine in the work area
Well between crustal stress of the position in the prediction time, and then obtain the work area the prediction time three-dimensional ground stress point
Cloth data.
In the present embodiment, in oilfield exploitation procedure, water filling and oil recovery can cause the fluid in stratum to flow, and make
It obtains formation pore fluid pressure to change, so that stratum deformation.The stratum deformation is mainly reflected in stratum not Tongfang
Upward principal stress changes.After the crustal stress on stratum changes, thus it is possible to vary stratum filtration parameter (such as permeability,
Porosity, compressed coefficient etc.), to further influence the fluid flowing in stratum.In this way, the server can provide fluid
Coupling model between crustal stress；Can based on crustal stress of the well point position in the prediction time in the work area, with
And the coupling model between the fluid and crustal stress, position is answered on the ground of the prediction time between determining the well in the work area
Power.Coupling model between the fluid and crustal stress can meet following formula.
In above formula (3), P is formation pore fluid pressure；E can be elasticity modulus；kijIt can be seepage coefficient；υ can be with
For Poisson's ratio；σ2Can between well crustal stress of the position in the prediction time.
In the present embodiment, position can be understood as removing the position other than the position of well point in the work area between the well.
For position between each well, the server can determine the water injection well or producing well that position is nearest between the well；It can be with base
In the well point position of water injection well in the crustal stress of the prediction time or the well point position of producing well in the prediction time
Coupling model between crustal stress and the fluid and crustal stress, position is answered on the ground of the prediction time between determining the well
Power.
In the present embodiment, by step S16, the server can obtain the well point position in the work area described
The crustal stress of prediction time.By step S18, position is in the prediction between the server can obtain the well in the work area
The crustal stress at moment.In this way, based on the well point position in the work area between the crustal stress and well of the prediction time position
In the crustal stress of the prediction time, the server can obtain the work area and divide in the three-dimensional ground stress of the prediction time
Cloth data.
In the present embodiment, the server can provide the initial three-dimensional crustal stress distributed data in work area；It can determine
The formation pore fluid pressure of well point position in work area in prediction time；It can be based on the formation pore fluid pressure, really
Determine the well point position in work area the prediction time crustal stress；It can be based on the well point position in the work area described pre-
The crustal stress for surveying the moment, determines crustal stress of the position in the prediction time between well, and then obtain the work area in the prediction
The three-dimensional ground stress distributed data at moment.The crustal stress prediction technique of this specification embodiment, it is contemplated that the factor of production is answered over the ground
The influence of power so that the crustal stress distribution of prediction is more accurate, and then can instruct oil-gas field development management and exploitation well
The adjustment of scheme.
The crustal stress prediction technique of the present embodiment, the crustal stress that can predict and evaluate work area in different development phases are three-dimensional
The regularity of distribution establishes the four-dimensional crustal stress distributed model changed over time, can be the Oil and gas field development scheme of different development phases
Adjustment and management provide important geologic basis, and then are the efficient and rational of China's densification low permeability oil and gas field and unconventionaloil pool
Exploitation provides technical support.Wherein, the four-dimensional crustal stress can refer to the three-dimensional ground stress in the oil gas field different development times.
In a Sample Scenario of the present embodiment, in oilfield exploitation procedure, due to the difference of injection-production relation so that
The variation of formation pore pressure, so that the original crustal stress of oil gas field changes.
Specifically, Fig. 2 is please referred to.Near water injection well, original crustal stress becomes well from the substantially uniform distribution of part
Eye nearby increases；Increasing degree far from crustal stress after wellbore is in the trend reduced.It is advised from the variation of individual well crustal stress and water filling
From the point of view of rule, as water filling crustal stress is in increased trend, the increase of horizontal minimum and maximum main crustal stress is shown.Crustal stress
Increasing degree and water injection intensity also have apparent relationship.When water injection intensity increases, the increasing degree of crustal stress also relative increase.
It recovers the oil and water filling is opposite process.Please refer to Fig. 3.Near producing well, original crustal stress is by local base
Originally being uniformly distributed becomes near wellbore reduction；It is in reduction trend to reduce amplitude far from wellbore crustal stress.From individual well crustal stress with adopt
From the point of view of the changing rule of liquid measure, as oil mining crustal stress is in the variation tendency reduced, horizontal minimum and maximum master is shown
The reduction of crustal stress.The reduction amplitude and Liquid output of crustal stress have apparent relationship.As accumulation Discharge rate increases, crustal stress subtracts
It is small；Liquid output is bigger, and the amplitude of reduction is also relatively bigger.
In this Sample Scenario, Fig. 4 and Fig. 5 are please referred to.Use the crustal stress prediction technique of this specification embodiment, prediction
20 years later crustal stress (Fig. 5) of Bohai gulf basin oil field development had occurred compared with original crustal stress (Fig. 4) it is prodigious
Variation.Specifically, the stress value at some positions becomes larger, and the stress value at some positions becomes smaller.By the actual measurement of 5 mouthfuls of wells of different phase
It examines, the relative error of maximum principal stress is 2.1%-6.8%, mean error 4.0%；The relative error of minimum principal stress is
2.3%-14.7%, mean error 7.5%.To illustrate the present embodiment crustal stress prediction technique reliability.
Please refer to Fig. 6.This specification embodiment also provides a kind of crustal stress prediction meanss, including with lower unit.
Unit 20, the initial three-dimensional crustal stress distributed data for providing work area are provided；
Receiving unit 22, for receiving prediction time and creation data；
First determination unit 24 determines the well point position in the work area described pre- for being based on the creation data
Survey the formation pore fluid pressure at moment；
Second determination unit 26, for being based on the formation pore fluid pressure and the initial three-dimensional crustal stress distribution number
According to, determine the well point position in the work area the prediction time crustal stress；
Third determination unit 28 is used for based on crustal stress of the well point position in the prediction time in the work area, really
Crustal stress of the position in the prediction time between well in the fixed work area, and then the work area is obtained in the prediction time
Three-dimensional ground stress distributed data.
In the 1990s, the improvement of a technology can be distinguished clearly be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So
And with the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow to be programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, programmable logic device
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, logic function determines device programming by user.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, designs and makes without asking chip maker
Dedicated IC chip 2.Moreover, nowadays, substitution manually makes IC chip, and this programming is also used instead mostly
" logic compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development
Seemingly, and the source code before compiling also handy specific programming language is write, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but there are many kind, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed are most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only needing method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, so that it may to be readily available the hardware circuit for realizing the logical method flow.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer.Specifically, computer for example may be used
Think personal computer, laptop computer, cellular phone, camera phone, smart phone, personal digital assistant, media play
It is any in device, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or these equipment
The combination of equipment.
As seen through the above description of the embodiments, those skilled in the art can be understood that this specification
The mode of required general hardware platform can be added to realize by software.Based on this understanding, the technical solution of this specification
Substantially the part that contributes to existing technology can be expressed in the form of software products in other words, the computer software
Product can be stored in a storage medium, such as ROM/RAM, magnetic disc, CD, including some instructions are used so that a computer
Equipment (can be personal computer, server either network equipment etc.) executes each embodiment of this specification or embodiment
Certain parts described in method.
Each embodiment in this specification is described in a progressive manner, identical similar portion between each embodiment
Point just to refer each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so description is fairly simple, related place is referring to embodiment of the method
Part explanation.
This specification can be used in numerous general or special purpose computing system environments or configuration.Such as：Personal computer,
Server computer, handheld device or portable device, laptop device, multicomputer system, microprocessor-based system,
Set top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer including any of the above system are set
Standby distributed computing environment etc..
This specification can describe in the general context of computer-executable instructions executed by a computer, such as journey
Sequence module.Usually, program module include routines performing specific tasks or implementing specific abstract data types, program, object,
Component, data structure etc..This specification can also be put into practice in a distributed computing environment, in these distributed computing environment
In, by executing task by the connected remote processing devices of communication network.In a distributed computing environment, program module
It can be located in the local and remote computer storage media including storage device.
Although depicting this specification by embodiment, it will be appreciated by the skilled addressee that there are many become for this specification
Shape and the spirit changed without departing from this specification, it is desirable to which the attached claims include these deformations and change without departing from this
The spirit of specification.
Claims (10)
1. a kind of crustal stress prediction technique, which is characterized in that including：
The initial three-dimensional crustal stress distributed data in work area is provided；
Receive prediction time and creation data；
Based on the creation data, determine the well point position in the work area the prediction time formation pore fluid pressure
Power；
Based on the formation pore fluid pressure and the initial three-dimensional crustal stress distributed data, the well point in the work area is determined
Crustal stress of the position in the prediction time；
Based on crustal stress of the well point position in the prediction time in the work area, position exists between determining the well in the work area
The crustal stress of the prediction time, so obtain the work area the prediction time three-dimensional ground stress distributed data.
2. the method as described in claim 1, which is characterized in that the initial three-dimensional crustal stress distributed data includes dimensionally answering
Power distributed model；Wherein, the three-dimensional ground stress distributed model includes original crustal stress and the correspondence of spatial position.
3. method as claimed in claim 2, the three-dimensional ground stress distributed model is built in the following way：
The three-dimensional Skeleton Model on stratum is established using drilling well, well logging and three dimensional seismic data；
Three-dimensional petrophysical parameter model is established using three dimensional seismic data inverting and well-log information；
The three-dimensional finite element simulation for carrying out crustal stress space distribution rule calculates；
Based on the three-dimensional Skeleton Model, the three-dimensional petrophysical parameter model and numerical simulation calculation as a result, structure three
Tie up crustal stress distributed model.
4. the method as described in claim 1, which is characterized in that the work area includes at least one water injection well and at least one adopts
Oil well；The creation data includes the water injection rate of the water injection well and the Liquid output of the producing well.
5. the method as described in claim 1, which is characterized in that the well point position in the determination work area is in the prediction
The formation pore fluid pressure at moment, including：
Using following formula calculate the well point position in the work area the prediction time formation pore fluid pressure；
Wherein,
P is the formation pore fluid pressure of prediction time；PwFor flowing bottomhole pressure (FBHP)；Q is the water injection rate of water injection well or the production liquid of producing well
Amount；μ is viscosity；B is volume factor；K is formation effective permeability；H is Effective thickness of formation；Ф is stratum effecive porosity；C
For system compressibility；A is supply area；Δ t adopts the time interval between the moment for prediction time and note.
6. the method as described in claim 1, which is characterized in that the well point position in the determination work area is in the prediction
The crustal stress at moment, including：
Using following formula calculate the well point position in the work area the prediction time crustal stress；
Wherein,
σ1For well point position the prediction time crustal stress；σ0For original crustal stress；P is the earth bore clearance flow of prediction time
Body pressure；α is Biot coefficients；υ is Poisson's ratio；P0For prime stratum pore fluid pressure.
7. the method as described in claim 1, which is characterized in that position is in the prediction between the well in the determination work area
The crustal stress at moment, including：
Coupling based on the well point position in the work area between the crustal stress and fluid and crustal stress of the prediction time
Model determines that position is in the crustal stress of the prediction time between the well in the work area.
8. the method for claim 7, the coupling model between the fluid and crustal stress meets following formula；
Wherein,
P is formation pore fluid pressure；E is elasticity modulus；kijFor seepage coefficient；υ is Poisson's ratio；σ2Position is described between well
The crustal stress of prediction time.
9. the method as described in claim 1, which is characterized in that the crustal stress include vertical stress, horizontal maximum principal stress,
At least one of horizontal minimum principal stress.
10. a kind of crustal stress prediction meanss, which is characterized in that including：
Unit, the initial three-dimensional crustal stress distributed data for providing work area are provided；
Receiving unit, for receiving prediction time and creation data；
First determination unit determines the well point position in the work area in the prediction time for being based on the creation data
Formation pore fluid pressure；
Second determination unit is used for the formation pore fluid pressure and the initial three-dimensional crustal stress distributed data, determines institute
State the well point position in work area the prediction time crustal stress；
Third determination unit, described in based on crustal stress of the well point position in the prediction time in the work area, determining
Crustal stress of the position in the prediction time between well in work area, so obtain the work area the prediction time dimensionally
Stress distribution data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810501796.7A CN108629463B (en) | 2018-05-23 | 2018-05-23 | Ground stress change prediction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810501796.7A CN108629463B (en) | 2018-05-23 | 2018-05-23 | Ground stress change prediction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108629463A true CN108629463A (en) | 2018-10-09 |
| CN108629463B CN108629463B (en) | 2022-08-09 |
Family
ID=63694049
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810501796.7A Active CN108629463B (en) | 2018-05-23 | 2018-05-23 | Ground stress change prediction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108629463B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427688A (en) * | 2019-07-29 | 2019-11-08 | 三峡大学 | A kind of crustal stress size prediction technique based on actual measurement vibration |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100191511A1 (en) * | 2007-08-24 | 2010-07-29 | Sheng-Yuan Hsu | Method For Multi-Scale Geomechanical Model Analysis By Computer Simulation |
| CN106014399A (en) * | 2016-05-31 | 2016-10-12 | 中国石油大学(北京) | Method for establishing high-precision three-dimensional ground stress model for heterogeneous formation |
| CN106991236A (en) * | 2017-04-05 | 2017-07-28 | 西南石油大学 | A kind of refracturing well and stratum selection method based on four-dimensional crustal stress dynamic change |
| CN107066718A (en) * | 2017-04-05 | 2017-08-18 | 西南石油大学 | A kind of four-dimensional dynamically stress simulation method |
-
2018
- 2018-05-23 CN CN201810501796.7A patent/CN108629463B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100191511A1 (en) * | 2007-08-24 | 2010-07-29 | Sheng-Yuan Hsu | Method For Multi-Scale Geomechanical Model Analysis By Computer Simulation |
| CN106014399A (en) * | 2016-05-31 | 2016-10-12 | 中国石油大学(北京) | Method for establishing high-precision three-dimensional ground stress model for heterogeneous formation |
| CN106991236A (en) * | 2017-04-05 | 2017-07-28 | 西南石油大学 | A kind of refracturing well and stratum selection method based on four-dimensional crustal stress dynamic change |
| CN107066718A (en) * | 2017-04-05 | 2017-08-18 | 西南石油大学 | A kind of four-dimensional dynamically stress simulation method |
Non-Patent Citations (3)
| Title |
|---|
| 梁何生等: "孔隙压力变化对地应力的影响研究", 《石油钻探技术》 * |
| 王连捷 等: "ANSYS 软件在求解地应力与流体耦合作用中的应用", 《地质力学学报》 * |
| 赵国: "加密调整井钻井地层压力预测方法研究与应用", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427688A (en) * | 2019-07-29 | 2019-11-08 | 三峡大学 | A kind of crustal stress size prediction technique based on actual measurement vibration |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108629463B (en) | 2022-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9068448B2 (en) | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs | |
| US9383465B2 (en) | Quantitative analysis of time-lapse seismic data | |
| US20100250215A1 (en) | Methods of modeling flow of gas within a reservoir | |
| US10359542B2 (en) | Generating dynamically calibrated geo-models in green fields | |
| CN109446735A (en) | A kind of generation method, equipment and the system of modeling logging data | |
| CN106153495A (en) | Rock gas absorption phase volume measuring method and device | |
| CN111999227A (en) | Stress sensitivity testing method and device | |
| CN108227032B (en) | Determination method and device for ground stress | |
| CN106970423B (en) | A kind of method and device of determining shale gas stratigraphic structure model | |
| CN107246900B (en) | A kind of method and device of determining oil-water interface position | |
| US20170184761A1 (en) | Method of determining a map of height of liquid hydrocarbon in a reservoir | |
| WO2015038162A1 (en) | Pseudo phase production simulation: a signal processing approach to assess quasi-multiphase flow production via successive analogous step-function relative permeability controlled models in reservoir flow simulation in order to rank multiple petro-physical realizations | |
| CN108386185A (en) | A kind of method and device of determining sand body connectivity | |
| Audigane et al. | Dual mesh method for upscaling in waterflood simulation | |
| CN108629463A (en) | Crustal stress prediction technique and device | |
| CN106845086B (en) | Formation pressure calculation method and device | |
| US11320555B2 (en) | Systems and methods for calibration of indeterministic subsurface discrete fracture network models | |
| CN115270533A (en) | Repeated fracturing design method and device, storage medium and electronic equipment | |
| CN109870743B (en) | Page-based seam space spread prediction method and device | |
| Casciano et al. | Well test simulation through Discrete Fracture Network modelling in a fractured carbonate reservoir | |
| Stephen et al. | Upscaling uncertainty analysis in a shallow-marine environment | |
| US11682167B2 (en) | Seamless scaling geomodeling | |
| CN109358364B (en) | Method, device and system for establishing underground river reservoir body geological model | |
| Vejbæk et al. | Downflank hydrocarbon potential identified using seismic inversion and geostatistics: Upper Maastrichtian reservoir unit, Dan Field, Danish Central Graben | |
| Cao et al. | An Estimation Model for Hydraulic Conductivity of Low‐Permeability and Unfilled Fractured Granite in Underground Water‐Sealed Storage Caverns |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |