CN108595875A - A kind of method, apparatus and system of determining side water gas condensate reservoir water breakthrough time - Google Patents
A kind of method, apparatus and system of determining side water gas condensate reservoir water breakthrough time Download PDFInfo
- Publication number
- CN108595875A CN108595875A CN201810430878.7A CN201810430878A CN108595875A CN 108595875 A CN108595875 A CN 108595875A CN 201810430878 A CN201810430878 A CN 201810430878A CN 108595875 A CN108595875 A CN 108595875A
- Authority
- CN
- China
- Prior art keywords
- water
- gas
- phase
- condensate
- flow velocity
- 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
Abstract
The embodiment of the present application discloses a kind of method, apparatus and system of determining side water gas condensate reservoir water breakthrough time, and the method includes the percolation flow velocities that the Darcy effect of non-Darcy effects and water phase based on gas phase determines water phase；The Darcy effect of non-Darcy effects and water phase based on gas phase determines the percolation flow velocity of water phase；Water breakthrough time and gas well are established away from the functional relation between at a distance from the water of side according to the percolation flow velocity of the water phase；Processing is modified to the functional relation using reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree；The water breakthrough time in target work area is determined according to the functional relation after correcting process.Using each embodiment of the application, the water breakthrough time of side water gas condensate reservoir can be more accurately predicted.
Description
Technical field
The present invention relates to petroleum gas seismic exploration technique fields, particularly, are related to a kind of determining side water gas condensate reservoir
The method, apparatus and system of water breakthrough time.
Background technology
Water breakthrough time prediction is edge water gas reservoir and a major issue in bottomwater gas field development process, studies people both at home and abroad
Member has made extensive work to it, derives different water breakthrough time computation models, and has edge closing water-setting gassing to hide the water breakthrough time
Research is still in the starting stage at initial stage.
During gas reservoir development, since the intrusion of side water can be such that gas reservoir section constantly changes.With the exploitation of gas reservoir, Bian Shui
Constantly moved from home position to gas well with the interface of gas-bearing formation, gas well starts water breakthrough after gas and water interface reaches shaft bottom, this when
Between be known as side water water breakthrough time.Currently, water gas condensate reservoir water breakthrough time computation model in side usually considers retrograde condensation effect
It influences, but the prediction result of water breakthrough time and practical compare still have certain deviation.Therefore, side water-setting gassing how is improved
The precision of prediction for hiding the water breakthrough time is those skilled in the art's technical problem urgently to be resolved hurrily.
Invention content
A kind of method, apparatus for being designed to provide determining side water gas condensate reservoir water breakthrough time of the embodiment of the present application and it is
System, can be more accurately predicted the water breakthrough time of side water gas condensate reservoir.
The method, apparatus and system of a kind of determining side water gas condensate reservoir water breakthrough time provided by the present application be by include with
What under type was realized：
A kind of method of determining side water gas condensate reservoir water breakthrough time, including：
The Darcy effect of non-Darcy effects and water phase based on gas phase determines the percolation flow velocity of water phase；
Water breakthrough time and gas well are established away from the functional relation between at a distance from the water of side according to the percolation flow velocity of the water phase；
Using reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree to the function
Relationship is modified processing；
The water breakthrough time in target work area is determined according to the functional relation after correcting process.
The method of the embodiment of the present application deckle water gas condensate reservoir water breakthrough time really, the non-Darcy effects based on gas phase
And the Darcy effect of water phase determines the percolation flow velocity of water phase, including：
The Darcy effect of non-Darcy effects and water phase based on gas phase determines the mapping between aqueous vapor two phase fluid flow speed
Relationship；
The percolation flow velocity of gas phase is determined based on the volume factor of condensate gas；
The percolation flow velocity of water phase is determined using the percolation flow velocity and the mapping relations of the gas phase.
The method of the embodiment of the present application deckle water gas condensate reservoir water breakthrough time really, the volume factor based on condensate gas
Determine the percolation flow velocity of gas phase, including：
The percolation flow velocity of gas phase is determined according to following formula：
Wherein, vgIndicate the percolation flow velocity of gas phase, qgIndicate the yield of condensate gas, BgIndicate the volume factor of condensate gas, h
Indicate that reservoir thickness, r indicate tonguing point to the distance in shaft bottom.
The method of the embodiment of the present application deckle water gas condensate reservoir water breakthrough time really, the percolation flow velocity of the determining water phase,
Including：
The percolation flow velocity of water phase is determined according to following formula：
Wherein, vg、vwGas phase, the percolation flow velocity of water phase, M are indicated respectivelywgIndicate that aqueous vapor mobility ratio, β indicate Non-Darcy's flow
Dynamic coefficient, μwIndicate the viscosity of gas phase, kwIndicate the permeability of water phase, vgIndicate the percolation flow velocity of gas phase.
The method of the embodiment of the present application deckle water gas condensate reservoir water breakthrough time really, it is described that the functional relation is repaiied
Positive processing, including：
Remaining condensate saturation degree is determined based on the volume factor of condensate gas；
Correspondingly, described be saturated using reservoir pore space initial water saturation, residual gas saturation and remaining condensate
Degree is modified processing to the functional relation and includes：Using reservoir pore space initial water saturation, residual gas saturation and
Processing is modified to the functional relation based on the remaining condensate saturation degree after the volume factor determination of condensate gas.
The method of the embodiment of the present application deckle water gas condensate reservoir water breakthrough time really, the volume factor based on condensate gas
Determine remnants condensate saturation degrees, including：
Remaining condensate saturation degree is determined according to following formula：
Wherein, SocIndicate that remaining condensate saturation degree, t indicate time, qgIndicate the yield of condensate gas, BgIndicate condensate gas
Volume factor, h indicate reservoir thickness, dC indicate unit volume condensate gas in condensate increment (reduction amount), rbIt indicates
The obstruction radius of critical condensate saturation degree in reservoir gaps.
On the other hand, the embodiment of the present application also provides a kind of device of determining side water gas condensate reservoir water breakthrough time, including：
Percolation flow velocity determining module, the Darcy effect for non-Darcy effects and water phase based on gas phase determine water phase
Percolation flow velocity；
Module is built, for establishing water breakthrough time and gas well away between at a distance from the water of side according to the percolation flow velocity of the water phase
Functional relation；
Correcting module, for full using reservoir pore space initial water saturation, residual gas saturation and remaining condensate
Processing is modified to the functional relation with degree；
Water breakthrough time determining module, the water breakthrough time for determining target work area according to the functional relation after correcting process.
The device of the embodiment of the present application deckle water gas condensate reservoir water breakthrough time really, the percolation flow velocity determining module packet
It includes：
Percolation flow velocity relation determination unit, the Darcy effect for non-Darcy effects and water phase based on gas phase determine water
Mapping relations between gas two phase fluid flow speed；
Gas phase percolation flow velocity determination unit, the percolation flow velocity for determining gas phase based on the volume factor of condensate gas；
Water phase percolation flow velocity determination unit, for using the gas phase percolation flow velocity and the mapping relations determine water
The percolation flow velocity of phase.
The device of the embodiment of the present application deckle water gas condensate reservoir water breakthrough time really, including processor and handled for storage
The memory of device executable instruction is realized when described instruction is executed by the processor and is included the following steps：
The Darcy effect of non-Darcy effects and water phase based on gas phase determines the percolation flow velocity of water phase；
Water breakthrough time and gas well are established away from the functional relation between at a distance from the water of side according to the percolation flow velocity of the water phase；
Using reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree to the function
Relationship is modified processing；
The water breakthrough time in target work area is determined according to the functional relation after correcting process.
On the other hand, the embodiment of the present application also provides a kind of system of determining side water gas condensate reservoir water breakthrough time, including extremely
A few processor and the memory for storing computer executable instructions, the processor is realized above-mentioned when executing described instruction
Described in any one embodiment the step of method.
A kind of method, apparatus for determining side water gas condensate reservoir water breakthrough time that this specification one or more embodiment provides
And system, it can further consider the shadow of nearly well area gas phase non-Darcy flow effect on the basis of considering retrograde condensation effect
It rings, the opposite side water gas condensate reservoir water breakthrough time is predicted, to further improve the prediction of side water gas condensate reservoir water breakthrough time
Precision.
Description of the drawings
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments described in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, other drawings may also be obtained based on these drawings.In the accompanying drawings：
Fig. 1 is a kind of flow signal of the embodiment of the method for determining side water gas condensate reservoir water breakthrough time that this specification provides
Figure；
Fig. 2 is side water gas condensate reservoir gas liquid two-phase flow model schematic in one embodiment that this specification provides；
Fig. 3 is water breakthrough time prediction contrast schematic diagram in the example that this specification provides；
Fig. 4 is water breakthrough time prediction contrast schematic diagram in another example that this specification provides；
Fig. 5 is a kind of modular structure of the device embodiment for determining side water gas condensate reservoir water breakthrough time that this specification provides
Schematic diagram.
Specific implementation mode
In order to make those skilled in the art more fully understand the technical solution in this specification, below in conjunction with this explanation
Attached drawing in book one or more embodiment carries out the technical solution in this specification one or more embodiment clear, complete
Site preparation describes, it is clear that described embodiment is only specification a part of the embodiment, instead of all the embodiments.Based on saying
Bright book one or more embodiment, it is obtained by those of ordinary skill in the art without making creative efforts all
Other embodiment should all belong to the range of this specification example scheme protection.
Water breakthrough time prediction is edge water gas reservoir and a major issue in bottomwater gas field development process, studies people both at home and abroad
Member has made extensive work to it, derives different water breakthrough time computation models, and has edge closing water-setting gassing to hide the water breakthrough time
Research is still in the starting stage at initial stage.Currently, water gas condensate reservoir water breakthrough time computation model in side usually considers retrograde condensation effect
Influence, but the prediction result of water breakthrough time and practical compared to still having certain deviation.More especially for high-yield well
Obviously.For above-mentioned technical problem, the application is on the basis of considering retrograde condensation effect, by further considering nearly well area gas
The influence of phase non-Darcy flow effect, opposite side water gas condensate reservoir water breakthrough time are predicted so that the prediction result of acquisition with it is existing
The field operating mode goodness of fit is preferable, and prediction result is more accurate.Meanwhile also to can be suitably used for different working conditions following for the scheme of the application
The water breakthrough time of water gas condensate reservoir predicts there is versatility.
Fig. 1 is that a kind of embodiment of the method flow of determining side water gas condensate reservoir water breakthrough time that this specification provides is shown
It is intended to.Although present description provides such as following embodiments or method operating procedure shown in the drawings or apparatus structure, it is based on
Less behaviour after routine may include either more in the method or device without performing creative labour or part merging
Make step or modular unit.In the step of there is no necessary causalities in logicality or structure, the execution of these steps is suitable
Sequence or the modular structure of device are not limited to this specification embodiment or execution shown in the drawings sequence or modular structure.The side
Device in practice, server or the end product of method or modular structure are in application, can be according to embodiment or attached drawing institute
Method or modular structure carry out the sequence execution shown is either parallel to be executed (such as the ring of parallel processor or multiple threads
Border, the even implementation environment including distributed treatment, server cluster).
Specific one embodiment is as shown in Figure 1, this specification provides the side of deckle water gas condensate reservoir water breakthrough time really
In one embodiment of method, the method may include：
S2, the percolation flow velocity that water phase is determined based on the non-Darcy effects of gas phase and the Darcy effect of water phase.
In the present embodiment, gas phase flow model in porous media, and reaching based on water phase can be built based on the non-Darcy effects of gas phase
Western effect builds water phase flow model in porous media；Then, the seepage flow of water phase is determined further according to gas phase flow model in porous media and water phase flow model in porous media
Speed.The high speed seepage flow of near wellbore zone gas generates more apparent non-Darcy effects, considers the non-Darcy effects of gas phase
And Darcy effect structure gas phase flow model in porous media, the water phase flow model in porous media of water phase, it can further increase and obtain gas and water two-phase
The accuracy of percolation flow velocity.
As illustrated in fig. 2, it is assumed that the gas and water boundary perimeter of gas condensate reservoir has a bite gas-producing well B, the well and initial reservoir gas,
The distance between water boundaries are approximately equal to a.Gas and water boundary has a water particle A to do planar radial flowing to gas-producing well B.The t times
Afterwards, the distance between the well and water particle A are ra, by tbtAfter time, water particle A reaches B wells.In order to facilitate retouching for scheme
It states, under conditions of meeting the side water gas condensate reservoir gas and water characteristics of motion, can make the following assumptions：1. tilt reservoir uniform thickness and
Matter；2. during gas and water two phase fluid flow, with movement displacement, and residual gas is formed in water enchroachment (invasion) region；3. ignore capillary force with
The influence of gravity, viscous force control the flowing of gas and water two-phase in the earth formation；4. fluid is micro- compressible, gas phase, the density of liquid phase and
Viscosity is definite value；5. condensate critically connected graph is definite value.
Based under above-mentioned assumed condition, in one embodiment of this specification, the gas phase flow model in porous media and water phase are oozed
Flow model can be expressed as respectively：
Wherein, Pg、PwGas phase, the pressure of water phase are indicated respectively；μg、μwGas phase, the viscosity of water phase are indicated respectively；vg、vwPoint
It Biao Shi not gas phase, the percolation flow velocity of water phase；kg、kwGas phase, the permeability of water phase are indicated respectively；R indicates high gas rate well B and water quality
The distance between point A；β indicates non-Darcy flow coefficient；ρgIndicate the density of gas phase.
In one embodiment of this specification, in gas and water boundary, the pressure difference dp of gas phaseg/ dr also has an impact water phase,
I.e. gas and water two-phase is equal in the pressure difference of boundary.Based on the gas and water two-phase boundary pressure differential, by above-mentioned water, gas two
The flow model in porous media of phase is combined, it may be determined that between water phase percolation flow velocity and gas phase percolation flow velocity (aqueous vapor two phase fluid flow speed)
Mapping relations.The mapping relations may be used different mathematical way conversions and be established.As this specification provides
In one embodiment, gas and water two-phase meets the following conditions in the pressure difference of boundary：
Formula (1), formula (2) and formula (3) progress simultaneous can be obtained：
Wherein, aqueous vapor mobility ratio MwgIt can be expressed as：
Wherein, λw、λgThe mobility of water phase, gas phase is indicated respectively.
Thus can be according to the percolation flow velocity of the quantitative determination water phase of the percolation flow velocity of gas phase using formula (4).Certainly,
This specification other embodiment is not excluded for the mapping relations of other forms, or deformation, the transformation etc. of the formula (4) based on offer
Mapping relations form.
In one embodiment of this specification, the percolation flow velocity of gas phase can be determined based on the volume factor of condensate gas, so
Afterwards, the percolation flow velocity for further determining water phase, the accuracy that the percolation flow velocity to improve water phase determines.Wherein, condensate gas
Volume factor may include ratio in unit volume shared by condensate gas.In one or more embodiment of this specification,
The percolation flow velocity of gas phase can be determined according to following formula：
Wherein, qgIndicate the yield of condensate gas, BgIndicate that the volume factor of condensate gas, h indicate reservoir thickness.Analytical formula
(6) it is found that qgBgIndicate that the seepage discharge of condensate gas, 2 π hr indicate corresponding cross-sectional area, to be quantified based on formula (6)
Determine the percolation flow velocity of gas phase.Meanwhile also fully considering that condensate gas proportion is different in actual formation, further improves
The accuracy that the percolation flow velocity of gas phase is determined using above-mentioned computation model, so as to further increase the standard of water breakthrough time prediction
Exactness.
S4, water breakthrough time and gas well are established according to the percolation flow velocity of the water phase closed away from the function between at a distance from the water of side
System, using reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree to the functional relation
It is modified processing；
In the present embodiment, the water breakthrough time can be established based on the percolation flow velocity of the water phase of above-mentioned acquisition with gas well away from side water
Functional relation between distance (i.e. the distance between gas well and gas and water interface), then, then considers reservoir pore space and initially contains
The influence of water saturation, residual gas saturation and remaining condensate saturation degree to the water breakthrough time, carries out above-mentioned functional relation
It corrects, and determines the water breakthrough time of target work area side water gas condensate reservoir according to revised functional relation.Such as, for above-mentioned step
Water particle A in rapid S2 can be expressed as within the dt times to the distance of high gas rate well B tonguings：
Wherein, φ indicates reservoir porosity.
Further consider the influence that the irreducible water, residual gas and the condensate that are detained in reservoir pore space predict the water breakthrough time,
Above-mentioned formula (7) can be then modified for：
Wherein, SwiIndicate initial water saturation, SgrIndicate residual gas saturation, SocIndicate remaining condensate saturation degree.
In some embodiments of the application, initial water saturation, residual gas saturation and remaining condensate can be directly acquired
The data of saturation degree are brought into above-mentioned formula (8), are completed to water breakthrough time and gas well away from the functional relation between at a distance from the water of side
Amendment.In one embodiment that this specification provides, the volume factor for being also based on condensate gas determines that remaining condensate is full
And degree.For example, the speedup of remaining condensate saturation degree can be determined first according to following formula：
Wherein, qgIndicate the yield of condensate gas, BgIndicate that the volume factor of condensate gas, h indicate that reservoir thickness, dC indicate single
The increment (reduction amount) of condensate, r in the volume condensate gas of positionbIndicate the obstruction half of critical condensate saturation degree in reservoir gaps
Then diameter determines the saturation degree of condensate further according to the speedup of condensate saturation degree.In the present embodiment, remaining condensation is being determined
Condensate gas proportion in actual formation has been fully considered during oily saturation degree, it, can to more meet actual formation situation
To further increase the accuracy of follow-up water breakthrough time prediction.It is of course also possible to directly substitute into formula (9) in formula (8), it can
：
Wherein, obstruction radius rbIt can be expressed as：
Wherein, Z indicates that deviation factor, T indicate that formation temperature, Y indicate that retrograde condensation coefficient, K indicate gas reservoir effective permeability,
PRIndicate original formation pressure.
S6, the water breakthrough time that target work area is determined according to the functional relation after correcting process.
In the present embodiment, the number of initial water saturation, residual gas saturation, remaining condensate saturation degree can be obtained
According to, and the percolation flow velocity data that determining water phase is calculated in above-mentioned steps are obtained, it is true according to the functional relation after correcting process
The water breakthrough time of deckle water gas condensate reservoir.Such as the functional relation after the correcting process of above-mentioned acquisition, it is assumed that before gas well liquid loading
Original gas and water interface is r at a distance from high-yield wella, as t=0, then ra=a integrates formula (10), can obtain：
Formula (4), formula (6) and formula (11) are substituted into formula (12), integrating and arranging can obtain：
Wherein：
Further formula (13) can simplify arranging and be：
Then water breakthrough time tbtIt can be expressed as：
To which seeing for side water gas condensate reservoir can be quantitatively determined according to formula (14), (15), (16), (17) and (19)
The water time.The scheme provided using this specification said one or multiple embodiments, on the basis of considering retrograde condensation effect,
By further considering that the influence of gas phase non-Darcy flow effect, opposite side water gas condensate reservoir water breakthrough time are predicted, so as to
To improve the precision of side water gas condensate reservoir water breakthrough time prediction.Meanwhile also fully combining condensate gas proportion in actual formation
Difference, the water breakthrough time is quantitatively determined, to further improve the accuracy of water breakthrough time prediction, particularly with containing in stratum
The larger work area of water, the raising of precision of prediction become apparent.
In order to further verify the practicability and feasibility of this method, this specification additionally provides one using said program
A specific example.Northwestward gas field directional well A1 wells, A2 wells, A3 wells and A4 wells belong to two blocks in the gas field, accordingly
Basic parameter is shown in Table 1.
Fig. 3 and Fig. 4 are the gas well water breakthrough times of the related basic gain of parameter based on A1 wells with distance of the gas well away from side water
Prediction result.Fig. 3 is indicated under different gas and water boundary conditions, is considered non-Darcy effects and is not considered that non-Darcy effects are corresponding
Water breakthrough time compares.As can be seen from Figure 3, gas phase non-Darcy flow effect is whether considered, with gas and water frontier distance
Constantly increase, the water breakthrough time gradually increases, and amplitude gradually increases with the increase of frontier distance, therefore gas well water breakthrough pushes away
Late.When frontier distance is identical, consider that the water breakthrough time of gas phase non-Darcy flow effect is shorter, gas well water breakthrough becomes early；On boundary
When smaller with a distance from gas well, the water breakthrough time influenced by gas phase non-Darcy effects it is smaller, with the continuous increase of distance, water breakthrough
Time is influenced gradually to increase by gas phase non-Darcy effects.
Under the conditions of Fig. 4 indicates different gas and water frontier distances, while ignoring retrograde condensation effect and non-Darcy effects, does not consider
Retrograde condensation acts on, and ignores the side water gas condensate reservoir water breakthrough time of retrograde condensation effect and the effect prediction of non-darcy simultaneously.By Fig. 4
It is found that ignoring the water breakthrough time longest of retrograde condensation effect and non-Darcy effects simultaneously, the water breakthrough time is the latest.When frontier distance is identical
When, consider that the water breakthrough time ratio of retrograde condensation effect ignores shorter when the factor, gas well water breakthrough change morning；When frontier distance is smaller,
Water breakthrough time is influenced little by retrograde condensation, but with the gradual increase of frontier distance, the water breakthrough time is influenced not by retrograde condensation
It is disconnected to increase.
14 mouthfuls of well basic parameter tables of table
Meanwhile in another example of this specification, also to the program prediction using the offer of this specification above-described embodiment
The precision of prediction of model and several common prediction models compares and analyzes, and analysis result is shown in Table 2.Wherein, model 1 is to ignore
Gas phase non-Darcy effects model corresponding with retrograde condensation effect, model 2 are imitated to consider that the non-darcy of gas phase is ignored in retrograde condensation effect
It is the corresponding model of scheme that this specification embodiment provides to answer corresponding model, model 3.
2 each the model calculation of table compares
As seen from Table 2, the precision of prediction of model 1 is minimum, and relative error is up to 45.25%；And model 2 compares model 1
For, precision of prediction slightly improves, but relative error is still very big, between 18.01%~27.15%；And the knot that model 3 calculates
Fruit relative error is smaller, is less than 8%, and compared with above-mentioned common calculation model, precision improves 14.64%~38.92%, and existing
The field operating mode goodness of fit is preferable.By comparative analysis as can be seen that model 3 by considering retrograde condensation effect and the non-darcy of gas phase simultaneously
The influence of effect, and the influence of the volume factor of condensate gas has also been fully considered during calculating, so that prediction result
Relative error is smaller.Scheme to demonstrate the offer of this specification above-described embodiment can more Accurate Prediction side water-setting gassing
The water breakthrough time of Tibetan.Particularly with the high-yield well that nearly well area gas flow rate is larger, gas phase non-Darcy effects are more apparent, utilize
The program prediction water breakthrough time in this specification embodiment, more meet work area actual conditions.
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.Specifically it is referred to
The description of aforementioned relevant treatment related embodiment, does not do repeat one by one herein.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the action recorded in detail in the claims or step can be come according to different from the sequence in embodiment
It executes and desired result still may be implemented.In addition, the process described in the accompanying drawings not necessarily require show it is specific suitable
Sequence or consecutive order could realize desired result.In some embodiments, multitasking and parallel processing be also can
With or it may be advantageous.
A kind of method for determining side water gas condensate reservoir water breakthrough time that this specification one or more embodiment provides, can be with
On the basis of considering retrograde condensation effect, the influence of nearly well area gas phase non-Darcy flow effect, opposite side water-setting are further considered
Gassing is hidden the water breakthrough time and is predicted, to further improve the precision of side water gas condensate reservoir water breakthrough time prediction.
Based on the method for the water of deckle really gas condensate reservoir water breakthrough time described above, this specification one or more embodiment
A kind of device of determining side water gas condensate reservoir water breakthrough time is also provided.The device may include having used this specification implementation
System, software (application), module, component, server of example the method etc. simultaneously combine the necessary device for implementing hardware.It is based on
Same innovation thinking, this specification embodiment provide one or more embodiments in device as described in the following examples.
Since the implementation that device solves the problems, such as is similar to method, the implementation of the specific device of this specification embodiment can be joined
See the implementation of preceding method, overlaps will not be repeated.Used below, term " unit " or " module " may be implemented pre-
Determine the combination of the software and/or hardware of function.Although device described in following embodiment is preferably realized with software,
The realization of the combination of hardware or software and hardware is also that may and be contemplated.Specifically, one kind that Fig. 5 specifications provide
Determine that the modular structure schematic diagram of the device embodiment of side water gas condensate reservoir water breakthrough time, such as Fig. 5, described device may include：
Percolation flow velocity determining module 102, the Darcy effect that can be used for non-Darcy effects and water phase based on gas phase are true
Determine the percolation flow velocity of water phase；
Module 104 is built, can be used for establishing the water breakthrough time according to the percolation flow velocity of the water phase with gas well away from side water
Functional relation between distance；
Correcting module 106 can be used for solidifying using reservoir pore space initial water saturation, residual gas saturation and remnants
Condensate oil saturation degree is modified processing to the functional relation；
Water breakthrough time determining module 108 can be used for determining seeing for target work area according to the functional relation after correcting process
The water time.
In another embodiment of this specification, the percolation flow velocity determining module 102 may include percolation flow velocity relationship
Determination unit, gas phase percolation flow velocity determination unit and water phase percolation flow velocity determination unit, wherein
The percolation flow velocity relation determination unit can be used for the darcy effect of non-Darcy effects and water phase based on gas phase
It should determine the mapping relations between aqueous vapor two phase fluid flow speed；
The gas phase percolation flow velocity determination unit can be used for determining the seepage flow speed of gas phase based on the volume factor of condensate gas
Degree；
The water phase percolation flow velocity determination unit can be used for closing using the percolation flow velocity of the gas phase and the mapping
System determines the percolation flow velocity of water phase.
In another embodiment of this specification, the correcting module 106 may include condensate saturation degree determination unit
And amending unit, wherein
The condensate saturation degree determination unit can be used for determining that remaining condensate is full based on the volume factor of condensate gas
And degree；
The amending unit can be used for using reservoir pore space initial water saturation, residual gas saturation and be based on
Remaining condensate saturation degree after the volume factor of condensate gas determines is modified processing to the functional relation.
It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
A kind of device for determining side water gas condensate reservoir water breakthrough time that this specification one or more embodiment provides, can be with
On the basis of considering retrograde condensation effect, the influence of nearly well area gas phase non-Darcy flow effect, opposite side water-setting are further considered
Gassing is hidden the water breakthrough time and is predicted, to further improve the precision of side water gas condensate reservoir water breakthrough time prediction.
The method or apparatus described in above-described embodiment that this specification provides can realize that business is patrolled by computer program
It collects and records on a storage medium, the storage medium can be read and be executed with computer, realize this specification embodiment institute
The effect of description scheme.Therefore, this specification also provides a kind of device of determining side water gas condensate reservoir water breakthrough time, including processing
Device and the memory for storing processor-executable instruction are realized when described instruction is executed by the processor and are included the following steps：
The Darcy effect of non-Darcy effects and water phase based on gas phase determines the percolation flow velocity of water phase；
Water breakthrough time and gas well are established away from the functional relation between at a distance from the water of side according to the percolation flow velocity of the water phase；
Using reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree to the function
Relationship is modified processing；
The water breakthrough time in target work area is determined according to the functional relation after correcting process.
The storage medium may include the physical unit for storing information, typically by after information digitalization again with profit
The media of the modes such as electricity consumption, magnetic or optics are stored.The storage medium, which has, may include：Letter is stored in the way of electric energy
The device of breath such as, various memory, such as RAM, ROM；The device of information is stored in the way of magnetic energy such as, hard disk, floppy disk, magnetic
Band, core memory, magnetic bubble memory, USB flash disk；Using optical mode store information device such as, CD or DVD.Certainly, also its
Readable storage medium storing program for executing of his mode, such as quantum memory, graphene memory etc..
It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
The device of determining side water gas condensate reservoir water breakthrough time described in above-described embodiment a kind of can be made considering retrograde condensation
On the basis of, the influence of nearly well area gas phase non-Darcy flow effect, opposite side water gas condensate reservoir water breakthrough time are further considered
It is predicted, to further improve the precision of side water gas condensate reservoir water breakthrough time prediction.
This specification also provides a kind of system of determining side water gas condensate reservoir water breakthrough time, and the system can be individual
Water breakthrough time forecasting system can also be applied in a plurality of types of seismic prospectings processing, solution release system.The system can be with
For individual server, one or more the methods of this specification or one or more embodiments can also be include the use of
The server cluster of device, system (including distributed system), software (application), practical operation device, logic gates device,
Quantum computer etc. simultaneously combines the necessary terminal installation for implementing hardware.The two-dimensional intersection survey line static corrections processing system can be with
Memory including at least one processor and storage computer executable instructions, the processor execute real when described instruction
Described in existing above-mentioned any one or multiple embodiments the step of method.
It should be noted that system described above can also include others according to the description of method or device embodiment
Embodiment, concrete implementation mode are referred to the description of related method embodiment, do not repeat one by one herein.
The system of determining side water gas condensate reservoir water breakthrough time described in above-described embodiment a kind of can be made considering retrograde condensation
On the basis of, the influence of nearly well area gas phase non-Darcy flow effect, opposite side water gas condensate reservoir water breakthrough time are further considered
It is predicted, to further improve the precision of side water gas condensate reservoir water breakthrough time prediction.
It should be noted that this specification device or system described above according to the description of related method embodiment also
May include other embodiments, concrete implementation mode is referred to the description of embodiment of the method, does not go to live in the household of one's in-laws on getting married one by one herein
It states.Each embodiment in this specification is described in a progressive manner, and identical similar part is mutual between each embodiment
Mutually referring to each embodiment focuses on the differences from other embodiments.Especially for hardware+program
For class, storage medium+program embodiment, since it is substantially similar to the method embodiment, so description is fairly simple, it is related
Place illustrates referring to the part of embodiment of the method.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the action recorded in detail in the claims or step can be come according to different from the sequence in embodiment
It executes and desired result still may be implemented.In addition, the process described in the accompanying drawings not necessarily require show it is specific suitable
Sequence or consecutive order could realize desired result.In some embodiments, multitasking and parallel processing be also can
With or it may be advantageous.
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, vehicle-mounted human-computer interaction device, cellular phone, camera phone, smart phone, individual
Digital assistants, media player, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or
The combination of any equipment in these equipment of person.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module is realized can in the same or multiple software and or hardware when specification one or more, it can also
The module for realizing same function is realized by the combination of multiple submodule or subelement etc..Device embodiment described above is only
It is only illustrative, for example, in addition the division of the unit, only a kind of division of logic function can have in actual implementation
Dividing mode, such as multiple units or component can be combined or can be integrated into another system or some features can be with
Ignore, or does not execute.Another point, shown or discussed mutual coupling, direct-coupling or communication connection can be logical
Some interfaces are crossed, the INDIRECT COUPLING or communication connection of device or unit can be electrical, machinery or other forms.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, application-specific integrated circuit, programmable
The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions that its inside includes can also be considered as in hardware component.Or
Person even, can will be considered as either the software module of implementation method can be hardware again for realizing the device of various functions
Structure in component.
The present invention be with reference to according to the method for the embodiment of the present invention, the flow of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions every first-class in flowchart and/or the block diagram
The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided
Instruct the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine so that the instruction executed by computer or the processor of other programmable data processing devices is generated for real
The device for the function of being specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that instruction generation stored in the computer readable memory includes referring to
Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device so that count
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, in computer or
The instruction executed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, computing device includes one or more processors (CPU), input/output interface, net
Network interface and memory.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
Including so that process, method, commodity or equipment including a series of elements include not only those elements, but also wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that wanted including described
There is also other identical elements in the process of element, method or equipment.
It will be understood by those skilled in the art that this specification one or more embodiment can be provided as method, system or calculating
Machine program product.Therefore, this specification one or more embodiment can be used complete hardware embodiment, complete software embodiment or
The form of embodiment combining software and hardware aspects.Moreover, this specification one or more embodiment can be used at one or
It is multiple wherein include computer usable program code computer-usable storage medium (include but not limited to magnetic disk storage,
CD-ROM, optical memory etc.) on the form of computer program product implemented.
This specification one or more embodiment can computer executable instructions it is general on
Described in hereafter, such as program module.Usually, program module includes executing particular task or realization particular abstract data type
Routine, program, object, component, data structure etc..This this specification one can also be put into practice in a distributed computing environment
Or multiple embodiments, in these distributed computing environments, by being held by the connected remote processing devices of communication network
Row task.In a distributed computing environment, program module can be located at the local and remote computer including storage device
In storage medium.
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.In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ",
The description of " specific example " or " some examples " etc. means specific features described in conjunction with this embodiment or example, structure, material
Or feature is contained at least one embodiment or example of this specification.In the present specification, to the signal of above-mentioned term
Property statement must not necessarily be directed to identical embodiment or example.Moreover, specific features, structure, material or the spy of description
Point may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other,
Those skilled in the art can be by different embodiments or examples described in this specification and different embodiments or examples
Feature is combined.
Above is only an example of the present application, it is not intended to limit this application.For those skilled in the art
For, the application can have various modifications and variations.It is all within spirit herein and principle made by any modification, equivalent
Replace, improve etc., it should be included within the scope of claims hereof.
Claims (10)
1. a kind of method of determining side water gas condensate reservoir water breakthrough time, which is characterized in that including：
The Darcy effect of non-Darcy effects and water phase based on gas phase determines the percolation flow velocity of water phase；
Water breakthrough time and gas well are established away from the functional relation between at a distance from the water of side according to the percolation flow velocity of the water phase；
Using reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree to the functional relation
It is modified processing；
The water breakthrough time in target work area is determined according to the functional relation after correcting process.
2. the method for determining side water gas condensate reservoir water breakthrough time according to claim 1, which is characterized in that described to be based on gas
The non-Darcy effects of phase and the Darcy effect of water phase determine the percolation flow velocity of water phase, including：
The Darcy effect of non-Darcy effects and water phase based on gas phase determines the mapping relations between aqueous vapor two phase fluid flow speed；
The percolation flow velocity of gas phase is determined based on the volume factor of condensate gas；
The percolation flow velocity of water phase is determined using the percolation flow velocity and the mapping relations of the gas phase.
3. the method for determining side water gas condensate reservoir water breakthrough time according to claim 2, which is characterized in that described based on solidifying
The volume factor of gassing determines the percolation flow velocity of gas phase, including：
The percolation flow velocity of gas phase is determined according to following formula：
Wherein, vgIndicate the percolation flow velocity of gas phase, qgIndicate the yield of condensate gas, BgIndicate that the volume factor of condensate gas, h indicate
Reservoir thickness, r indicate tonguing point to the distance in shaft bottom.
4. the method for determining the side water gas condensate reservoir water breakthrough time according to claim 1-3 any one of them, which is characterized in that institute
The percolation flow velocity of determining water phase is stated, including：
The percolation flow velocity of water phase is determined according to following formula：
Wherein, vg、vwGas phase, the percolation flow velocity of water phase, M are indicated respectivelywgIndicate that aqueous vapor mobility ratio, β indicate non-Darcy flow system
Number, μwIndicate the viscosity of gas phase, kwIndicate the permeability of water phase, vgIndicate the percolation flow velocity of gas phase.
5. the method for determining side water gas condensate reservoir water breakthrough time according to claim 1, which is characterized in that described to described
Functional relation is modified processing, including：
Remaining condensate saturation degree is determined based on the volume factor of condensate gas；
Correspondingly, described utilize reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree pair
The functional relation is modified processing：Using reservoir pore space initial water saturation, residual gas saturation and it is based on
Remaining condensate saturation degree after the volume factor of condensate gas determines is modified processing to the functional relation.
6. the method for determining side water gas condensate reservoir water breakthrough time according to claim 5, which is characterized in that described based on solidifying
The volume factor of gassing determines remaining condensate saturation degree, including：
Remaining condensate saturation degree is determined according to following formula：
Wherein, SocIndicate that remaining condensate saturation degree, t indicate time, qgIndicate the yield of condensate gas, BgIndicate the body of condensate gas
Product coefficient, h indicate that reservoir thickness, dC indicate the increment (reduction amount) of condensate in unit volume condensate gas, rbIndicate reservoir
The obstruction radius of critical condensate saturation degree in gap.
7. a kind of device of determining side water gas condensate reservoir water breakthrough time, which is characterized in that including：
Percolation flow velocity determining module, the Darcy effect for non-Darcy effects and water phase based on gas phase determine the seepage flow of water phase
Speed；
Module is built, for establishing water breakthrough time and gas well away from the letter between at a distance from the water of side according to the percolation flow velocity of the water phase
Number relationship；
Correcting module, for utilizing reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree
Processing is modified to the functional relation；
Water breakthrough time determining module, the water breakthrough time for determining target work area according to the functional relation after correcting process.
8. the device of determining side water gas condensate reservoir water breakthrough time according to claim 7, which is characterized in that the seepage flow speed
Spending determining module includes：
Percolation flow velocity relation determination unit, the Darcy effect for non-Darcy effects and water phase based on gas phase determine aqueous vapor two
Mapping relations between phase percolation flow velocity；
Gas phase percolation flow velocity determination unit, the percolation flow velocity for determining gas phase based on the volume factor of condensate gas；
Water phase percolation flow velocity determination unit, for using the gas phase percolation flow velocity and the mapping relations determine water phase
Percolation flow velocity.
9. a kind of device of determining side water gas condensate reservoir water breakthrough time, which is characterized in that handled including processor and for storage
The memory of device executable instruction is realized when described instruction is executed by the processor and is included the following steps：
The Darcy effect of non-Darcy effects and water phase based on gas phase determines the percolation flow velocity of water phase；
Water breakthrough time and gas well are established away from the functional relation between at a distance from the water of side according to the percolation flow velocity of the water phase；
Using reservoir pore space initial water saturation, residual gas saturation and remaining condensate saturation degree to the functional relation
It is modified processing；
The water breakthrough time in target work area is determined according to the functional relation after correcting process.
10. a kind of system of determining side water gas condensate reservoir water breakthrough time, which is characterized in that including at least one processor and deposit
The memory of computer executable instructions is stored up, the processor realizes any one of claim 1-6 when executing described instruction
The step of the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810430878.7A CN108595875A (en) | 2018-05-08 | 2018-05-08 | A kind of method, apparatus and system of determining side water gas condensate reservoir water breakthrough time |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810430878.7A CN108595875A (en) | 2018-05-08 | 2018-05-08 | A kind of method, apparatus and system of determining side water gas condensate reservoir water breakthrough time |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108595875A true CN108595875A (en) | 2018-09-28 |
Family
ID=63635880
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810430878.7A Pending CN108595875A (en) | 2018-05-08 | 2018-05-08 | A kind of method, apparatus and system of determining side water gas condensate reservoir water breakthrough time |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108595875A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111255414A (en) * | 2020-01-15 | 2020-06-09 | 清华大学 | Method and device for inhibiting condensate gas reservoir reverse condensation |
| CN113356817A (en) * | 2021-07-19 | 2021-09-07 | 中海油田服务股份有限公司 | Perforation optimization design method for bottom water reservoir horizontal well based on prolonging anhydrous oil production period |
| CN114062647A (en) * | 2021-06-07 | 2022-02-18 | 中国石油大学(北京) | Method and device for determining pressure range of reverse condensation four-zone of sandstone condensate gas reservoir |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104376420A (en) * | 2014-11-20 | 2015-02-25 | 中国石油天然气股份有限公司 | Water breakthrough risk evaluation method and evaluation device for water-carrying gas reservoir gas well |
| CN107045671A (en) * | 2017-03-22 | 2017-08-15 | 重庆科技学院 | Water-producing gas well hydrops Risk Forecast Method |
-
2018
- 2018-05-08 CN CN201810430878.7A patent/CN108595875A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104376420A (en) * | 2014-11-20 | 2015-02-25 | 中国石油天然气股份有限公司 | Water breakthrough risk evaluation method and evaluation device for water-carrying gas reservoir gas well |
| CN107045671A (en) * | 2017-03-22 | 2017-08-15 | 重庆科技学院 | Water-producing gas well hydrops Risk Forecast Method |
Non-Patent Citations (1)
| Title |
|---|
| 明瑞卿 等: "边水凝析气藏高产井见水时间预测新模型", 《特种油气藏》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111255414A (en) * | 2020-01-15 | 2020-06-09 | 清华大学 | Method and device for inhibiting condensate gas reservoir reverse condensation |
| CN114062647A (en) * | 2021-06-07 | 2022-02-18 | 中国石油大学(北京) | Method and device for determining pressure range of reverse condensation four-zone of sandstone condensate gas reservoir |
| CN113356817A (en) * | 2021-07-19 | 2021-09-07 | 中海油田服务股份有限公司 | Perforation optimization design method for bottom water reservoir horizontal well based on prolonging anhydrous oil production period |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | Water coning calculations for vertical and horizontal wells | |
| CN108595875A (en) | A kind of method, apparatus and system of determining side water gas condensate reservoir water breakthrough time | |
| CN110998370B (en) | Parallel processing of invasion and seepage for large-scale high-resolution simulation of oil and gas secondary migration | |
| CN108131127B (en) | Method and device for obtaining gas-oil ratio of production of foam oil type extra heavy oil field | |
| Zalavadia et al. | A hybrid modeling approach to production control optimization using dynamic mode decomposition | |
| CN111323814B (en) | Method for quantitatively determining water saturation of sandstone reservoir based on rock physical template | |
| Li et al. | Local–global upscaling for compositional subsurface flow simulation | |
| Filippov et al. | Reservoir modeling of complex structure reservoirs on dynamic adaptive 3D Pebi-grid | |
| CN109978303B (en) | Method and device for determining exploitation measures of heavy oil reservoir and storage medium | |
| CN107153895A (en) | Superimposed Basins lithologic deposit Beneficial Zones of Exploring quantitative forecasting technique and device | |
| Rotzien et al. | Deepwater depositional environments | |
| Zeboudj et al. | Horizontal well performance flow simulation CFD-application | |
| CN113946931A (en) | Method and device for determining dominant factors of logging encounter card | |
| Chaki et al. | A Proxy Flow Modelling Workflow to Estimate Gridded Dynamic Properties and Well Production Rates by Deep Learning Algorithms | |
| Absalyamov et al. | Coupled static and dynamic high resolution modelling on high performance computer clusters | |
| Grassi et al. | Comments on Babu and Odeh's Model for Productivity Prediction of Partially Penetrated Horizontal Wells | |
| Qoqandi et al. | Production Forecasting in Tight Gas Reservoirs Using Long Short-Term Memory Methods (LSTM) | |
| US20220205359A1 (en) | Systems and methods for back-allocation of oil produced by waterflooding | |
| Falola et al. | Rapid High-Fidelity Forecasting for Geological Carbon Storage Using Neural Operator and Transfer Learning | |
| Bandara et al. | Effect of wettability alteration on long-term behavior of fluids in subsurface | |
| Nunna | Application of Pressure Transient Concepts for Improved Upscaling of Geologic Models | |
| CN108536892B (en) | Method, device and system for determining downward migration depth of petroleum | |
| Khait et al. | A GPU-based integrated simulation framework for modelling of complex subsurface applications | |
| US20230012733A1 (en) | Machine learning workflow for predicting hydraulic fracture initiation | |
| Haider | Advanced Production Diagnostics and Optimized Drilling Strategy: Unlocking the True Potential of Production Data Using Physics and Machine Learning |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180928 |
|
| RJ01 | Rejection of invention patent application after publication |