CN109858701A - The quantitative identification method and system of Fractured edge water gas reservoir water influx - Google Patents
The quantitative identification method and system of Fractured edge water gas reservoir water influx Download PDFInfo
- Publication number
- CN109858701A CN109858701A CN201910111218.7A CN201910111218A CN109858701A CN 109858701 A CN109858701 A CN 109858701A CN 201910111218 A CN201910111218 A CN 201910111218A CN 109858701 A CN109858701 A CN 109858701A
- Authority
- CN
- China
- Prior art keywords
- water
- gas reservoir
- influx
- linear regression
- production
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The present invention provides the quantitative identification methods and system of a kind of Fractured edge water gas reservoir water influx, and the method includes: obtaining geology and creation data, obtain geologic parameter and the means of production according to the geology and creation data and water influx quantitative judge algorithm；According to the geologic parameter and the means of production and default gas reservoir and water body radius value, is calculated by Fractured edge water gas reservoir water influx calculation formula and obtain water enchroachment (invasion) magnitude；According to water drive gas reservoir matter balance equation and the water enchroachment (invasion) magnitude, performance graph is obtained by linear regression method, is calculated according to the performance graph and obtains the corresponding linear regression coeffficient of the performance graph；When the linear regression coeffficient is greater than predetermined threshold, water influx recognition result is obtained according to the performance graph and the geologic parameter and the means of production.
Description
Technical field
The present invention relates to oil-gas field development technical fields, more particularly to a kind of quantifying for Fractured edge water gas reservoir water influx
Recognition methods and system.
Background technique
Fractured side water water drive gas reservoir occupies critical positions in China's gas field development, and for such gas reservoir, side water is altered
Into may make many reservoirs by water seal every reducing the control area of gas well, seriously affect gas well deliverability and production stability.When
When yield is lower than critical liquid carrying flow rate, it is likely that lead to the serious hydrops in gas well shaft bottom and stop production.Water influx can be calculated, according to water
It invades situation and takes certain drainage measure, discharge formation gas to resume production.The calculating of water influx is also gas reservoir production management work
The basic content made.The dynamic water enchroachment (invasion) numerical quantity for rapidly calculating gas reservoir in time, facilitates the pipe of gas reservoir development process
Reason.It according to water enchroachment (invasion) magnitude, can not only determine that water body parameter can also carry out gas reservoir protection prediction, be mentioned for the production plan in later period
For foundation.
It however for a long time, is unfolded for Single Medium gas reservoir mostly to the recognition methods research work of water influx,
The research of fracture-hole dual media gas reservoir water influx quantitative judge is less.It is difficult fracture side according to existing method
The water influx of Gas Reservoirs, which is made, accurately to be identified, and then influences the adjustment of the formulation and production measure of subsequent development scheme.For more
Fractured edge water gas reservoir is cost-effectively developed, quantitative judge water influx is just particularly important.
Summary of the invention
It is an object of that present invention to provide the quantitative identification methods and system of a kind of Fractured edge water gas reservoir water influx, to water
The amount of invading (also including water body radius, gas reservoir radius, dynamic geological reserves) quantitative assessment, is supervised conducive to the dynamic of Fractured edge water gas reservoir
Survey and production management.
In order to achieve the above object, the quantitative identification method of Fractured edge water gas reservoir water influx provided by the present invention, specific to wrap
Contain: obtaining geology and creation data, geologic parameter is obtained according to the geology and creation data and water influx quantitative judge algorithm
And the means of production；According to the geologic parameter and the means of production and default gas reservoir and water body radius value, pass through Fractured side aqueous vapor
It hides water influx calculation formula and calculates acquisition water enchroachment (invasion) magnitude；According to water drive gas reservoir matter balance equation and the water enchroachment (invasion) magnitude, pass through
Linear regression method obtains performance graph, is calculated according to the performance graph and obtains the corresponding linear regression system of the performance graph
Number；When the linear regression coeffficient is greater than predetermined threshold, according to the performance graph and the geologic parameter and the means of production
Obtain water influx recognition result.
In above-mentioned Fractured edge water gas reservoir water influx quantitative identification method, it is preferred that the geology and creation data packet
Containing gas reservoir, gas and parameters of formation water and Production development data.
In above-mentioned Fractured edge water gas reservoir water influx quantitative identification method, it is preferred that the geologic parameter and production money
Material includes water influx quantitative judge algorithm input parameter and water influx quantitative judge algorithm output parameter；The water influx is quantitatively known
Other algorithm input parameter includes: reservoir pressure, deviation factor for gas, fracture permeabgility, gas reservoir porosity, system compressibility,
When gas viscosity, reservoir effective thickness, elastic storativity ratio, interporosity flow coefficient, irreducible water saturation, formation water volume factor, production
Between, surface gas rate, gas volume factor, tired gas production, tired water yield；The water influx quantitative judge algorithm output parameter packet
Contain: gas reservoir radius, water body radius, water enchroachment (invasion) magnitude.
In above-mentioned Fractured edge water gas reservoir water influx quantitative identification method, it is preferred that according to the geologic parameter and in advance
If gas reservoir and water body radius value, calculating acquisition water enchroachment (invasion) magnitude by Fractured edge water gas reservoir water influx calculation formula includes: by institute
Water body radius, time and pressure parameter zero dimension in default gas reservoir and water body radius value are stated, and brings Fractured edge water gas reservoir into
Zero dimension water influx is obtained in water influx calculation formula；It is public by stehfest numerical inversion according to the zero dimension water influx
Formula obtains the nondimensional water influx under real space；Water influx is obtained according to level pressure formula and the nondimensional water influx
Value.
In above-mentioned Fractured edge water gas reservoir water influx quantitative identification method, it is preferred that according to water drive gas reservoir material balance
Equation and the water enchroachment (invasion) magnitude, obtain performance graph by linear regression method, obtain fitting a straight line according to the performance graph
Linear regression coeffficient include: the water enchroachment (invasion) magnitude is brought into water drive gas reservoir matter balance equation, and combine geologic parameter and
The means of production establish the relationship scatter plot between gas reservoir pressure and oil in place, pass through linear regression according to the relationship scatter plot
Algorithm obtains performance graph, is calculated according to the performance graph and obtains the corresponding linear regression coeffficient of the performance graph.
In above-mentioned Fractured edge water gas reservoir water influx quantitative identification method, it is preferred that the method further includes: when described
When linear regression coeffficient is less than or equal to predetermined threshold, the default gas reservoir and water body radius are adjusted according to the relationship scatter plot
Value；Acquisition linear regression coeffficient is recalculated according to the adjustment adjusted default gas reservoir and water body radius value；Until described
When linear regression coeffficient is greater than predetermined threshold, according to the corresponding performance graph of the linear regression coeffficient and the geologic parameter and
The means of production obtain water influx recognition result.
The present invention also provides a kind of quantitative judge system of Fractured edge water gas reservoir water influx, the system is obtained comprising data
Modulus block, computing module, building module and identification module；The data acquisition module is for obtaining geology and creation data, root
Geologic parameter and the means of production are obtained according to the geology and creation data and water influx quantitative judge algorithm；The computing module is used
According to the geologic parameter and the means of production and default gas reservoir and water body radius value, pass through Fractured edge water gas reservoir water enchroachment (invasion) meter
It calculates formula and calculates acquisition water enchroachment (invasion) magnitude；The building module is used for according to water drive gas reservoir matter balance equation and the water influx
Value obtains performance graph by linear regression method, is calculated according to the performance graph and obtains the corresponding line of the performance graph
Property regression coefficient；The identification module is used for when the linear regression coeffficient is greater than predetermined threshold, according to the performance graph
Water influx recognition result is obtained with the geologic parameter and the means of production.
In the quantitative judge system of above-mentioned Fractured edge water gas reservoir water influx, it is preferred that the computing module includes place
Manage unit, inverting unit and level pressure unit；The processing unit is used for water body in the default gas reservoir and water body radius value half
Diameter, time and pressure parameter zero dimension, and bring acquisition zero dimension water enchroachment (invasion) in Fractured edge water gas reservoir water influx calculation formula into
Amount；The inverting unit is used to obtain real space by stehfest numerical inversion formula according to the zero dimension water influx
Under nondimensional water influx；The level pressure unit is used to obtain water influx according to level pressure formula and the nondimensional water influx
Value.
In the quantitative judge system of above-mentioned Fractured edge water gas reservoir water influx, it is preferred that the building module also includes:
The water enchroachment (invasion) magnitude is brought into water drive gas reservoir matter balance equation, and geologic parameter and the means of production is combined to establish gas reservoir pressure
Relationship scatter plot between oil in place obtains performance graph, root by linear regression algorithm according to the relationship scatter plot
It is calculated according to the performance graph and obtains the corresponding linear regression coeffficient of the performance graph.
In the quantitative judge system of above-mentioned Fractured edge water gas reservoir water influx, it is preferred that the identification module also includes:
When the linear regression coeffficient is less than or equal to predetermined threshold, the default gas reservoir and water are adjusted according to the relationship scatter plot
Body radius value；The quantitative judge system recalculates acquisition according to the adjustment adjusted default gas reservoir and water body radius value
Linear regression coeffficient；It is corresponding dynamic according to the linear regression coeffficient when linear regression coeffficient is greater than predetermined threshold
State curve and the geologic parameter and the means of production obtain water influx recognition result.
The present invention also provides a kind of computer equipment, including memory, processor and storage on a memory and can located
The computer program run on reason device, the processor realize the above method when executing the computer program.
The present invention also provides a kind of computer readable storage medium, the computer-readable recording medium storage has in execution
State the computer program of method.
The quantitative identification method of Fractured edge water gas reservoir water influx provided by the present invention and system identification go out Fractured side
Gas Reservoirs water influx, and gas reservoir radius, water body radius and Production development curve are obtained, it is that subsequent production dynamic monitoring and formulation are closed
The development plan of reason lays the foundation.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present invention, constitutes part of this application, not
Constitute limitation of the invention.In the accompanying drawings:
Fig. 1 is the method flow diagram according to one embodiment of the present of invention；
Fig. 2 is that water influx calculates step schematic diagram；
Fig. 3 is the schematic diagram calculation that water influx is calculated according to level pressure formula；
Fig. 4 is production performance graph schematic diagram；
Fig. 5 is practical Fractured edge water gas reservoir water influx calculated result proof diagram；
Fig. 6 is the quantitative identification method explanation results schematic diagram of Fractured edge water gas reservoir water influx；
Fig. 7 is the quantitative judge system explanation results schematic diagram of Fractured edge water gas reservoir water influx.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to apply to the present invention whereby
Technological means solves technical problem, and the realization process for reaching technical effect can fully understand and implement.It needs to illustrate
As long as not constituting conflict, each feature in each embodiment and each embodiment in the present invention can be combined with each other, institute
The technical solution of formation is within the scope of the present invention.
In addition, step shown in the flowchart of the accompanying drawings can be in the department of computer science of such as a group of computer-executable instructions
It is executed in system, although also, logical order is shown in flow charts, and it in some cases, can be to be different from herein
Sequence execute shown or described step.
It please refers to shown in Fig. 1, the quantitative identification method of Fractured edge water gas reservoir water influx provided by the present invention, it is specific to wrap
Contain: S101 obtains geology and creation data, obtains geology according to the geology and creation data and water influx quantitative judge algorithm
Parameter and the means of production；S102 passes through crack according to the geologic parameter and the means of production and default gas reservoir and water body radius value
Property edge water gas reservoir water influx calculation formula calculate obtain water enchroachment (invasion) magnitude；S103 is according to water drive gas reservoir matter balance equation and the water
Magnitude is invaded, performance graph is obtained by linear regression method, it is corresponding to calculate the acquisition performance graph according to the performance graph
Linear regression coeffficient；S104 is when the linear regression coeffficient is greater than predetermined threshold, with described according to the performance graph
Matter parameter and the means of production obtain water influx recognition result.Wherein, the geology and creation data may include gas reservoir, gas and ground
Layer water parameter and Production development data；The geologic parameter and the means of production may include water influx quantitative judge algorithm input parameter
With water influx quantitative judge algorithm output parameter；The water influx quantitative judge algorithm input parameter includes: reservoir pressure, gas
Deviation factors, fracture permeabgility, gas reservoir porosity, system compressibility, gas viscosity, reservoir effective thickness, elastic storativity ratio,
Interporosity flow coefficient, irreducible water saturation, formation water volume factor, production time, surface gas rate, gas volume factor, tired production gas
Amount, tired water yield；The water influx quantitative judge algorithm output parameter includes: gas reservoir radius, water body radius, water enchroachment (invasion) magnitude.
Due to being difficult to obtain gas reservoir radius and water body radius parameter value, in an embodiment of the present invention, according to Fractured side
When Gas Reservoirs water influx calculation formula calculates water influx, the mode that default gas reservoir and water body radius value are verified can be used；Tool
Body, it please refers to shown in Fig. 2, in this embodiment, according to the geologic parameter and default gas reservoir and water body radius value, by splitting
Seam property edge water gas reservoir water influx calculation formula, which calculates, to be obtained water enchroachment (invasion) magnitude and include: S201 is by the default gas reservoir and water body radius value
Middle water body radius, time and pressure parameter zero dimension, and bring into Fractured edge water gas reservoir water influx calculation formula and obtain nothing
Dimension water influx；S202 is obtained under real space according to the zero dimension water influx by stehfest numerical inversion formula
Nondimensional water influx；S203 obtains water enchroachment (invasion) magnitude according to level pressure formula and the nondimensional water influx.Wherein, the level pressure
The level pressure formula of Everdingen and Hurst proposition can be used in formula.
By the above-mentioned means, being carried out after calculating acquisition water influx using default gas reservoir and water body radius value, it need to be carried out
Effectively verifying, for this purpose, in an embodiment of the present invention, mainly making Production development song according to water drive gas reservoir matter balance equation
Line is prepared for verifying water influx calculated value, specifically, being led to according to water drive gas reservoir matter balance equation and the water enchroachment (invasion) magnitude
It crosses linear regression method and obtains performance graph, include according to the linear regression coeffficient that the performance graph obtains fitting a straight line: will
The water enchroachment (invasion) magnitude is brought into water drive gas reservoir matter balance equation, and combine geologic parameter and the means of production establish gas reservoir pressure with
Relationship scatter plot between oil in place obtains performance graph by linear regression algorithm according to the relationship scatter plot, according to
The performance graph, which calculates, obtains the corresponding linear regression coeffficient of the performance graph.In this embodiment, meter is mainly utilized
For obtained water influx as a result, according to water drive gas reservoir matter balance equation and in conjunction with creation data, linear regression goes out Production development
Curve, and obtain the related coefficient of fitting a straight line；Thereafter, calculated water enchroachment (invasion) magnitude is substituted into water drive gas reservoir matter balance equation
In, in conjunction with creation data and actual parameter, make about PH~GPThe scatter plot of relationship, linear regression, which goes out, produces performance graph, and
The correlation coefficient value of regression straight line is obtained, correlation coefficient value can reflect the accuracy that water influx calculates；With this, Production development is utilized
The correlation coefficient value of curve matching is constrained, final output Fractured edge water gas reservoir water influx quantitative judge result.
In the above-described embodiments, parameter zero dimension formula:
In formula, rDFor dimensionless radius, zero dimension；tDFor non dimensional time, zero dimension；PDFor zero dimension pressure, it is no because
It is secondary；R is any side water body radius, reFor side water radius, rgFor oil reservoir radius, m；KfFor fracture permeabgility, μm2；T is time, h；
φmFor matrix porosity, dimensionless；φfFor fracture porosity, dimensionless；ctfFor crack system compressibility, MPa-1；ctmFor
Matrix system compressibility, MPa-1；P is strata pressure, MPa；H is gas reservoir effective thickness, m；qgFor surface gas rate, m3/d；Bg
For gas volume factor, m3/m3；μgFor gas viscosity, mPas.
Zero dimension water influx calculation formula under the space Laplace are as follows:
In formula,For in the transformed zero dimension water influx of Laplace, dimensionless；U is Laplace Space Operators, nothing
Dimension；WithIt is the single order of Second Type and the Bessel function of zeroth order, dimensionless respectively.
Stehfest numerical inversion obtains the zero dimension water influx under real space:
In formula, WeDFor zero dimension water influx, dimensionless.
The level pressure solution using diffusion equation for recycling Everdingen and Hurst to propose, it is assumed that one at gas-water interface
It is a small that pressure drop, the water that this corresponding pressure drop of calculating enters from water layer are lauched each pressure drop fallen under gross pressure
The amount of invading is cumulative to obtain total water influx.And zero dimension water influx under real space is had into dimension, obtain the following institute of water influx formula
Show:
......
In formula, t0Indicate initial time, h；t1Indicate first moment, h；t2Indicate second moment, h；tj-1Indicate the
J-1 moment, h；PiIndicate initial time t0Corresponding pressure, MPa；P1Indicate first moment t1Corresponding pressure, MPa；P2
Indicate second moment t2Corresponding pressure, MPa；Pj-1Indicate -1 moment t of jthj-1Corresponding pressure, MPa；ΔP0It indicates just
Beginning section pressure difference, MPa；ΔP1Indicate first segment pressure difference, MPa；ΔP2Indicate second segment pressure difference, MPa；ΔPj-1Indicate -1 section of jth pressure
Difference, MPa；Indicate first segment average pressure, MPa；Indicate second segment average pressure, MPa；Indicate -1 section of jth pressure
Difference, MPa.
According to one embodiment of present invention, in step 4, according to water drive gas reservoir matter balance equation, in conjunction with step 3
The water enchroachment (invasion) magnitude of middle output and practical geology manufacturing parameter, make PH~GPRelationship scatter plot, linear regression go out Production development song
Line, and obtain linear regression related coefficient.Wherein, water drive gas reservoir matter balance equation is as follows:
In formula, cpFor rock compressibility, MPa-1；cwFor the water flooding compressed coefficient, MPa-1；swcFor irreducible water saturation,
Dimensionless.WeFor water influx, m3；WpFor surface water rate, m3；BwFor formation water volume factor, m3/m3；AgFor gas reservoir area, m2。
It enables:
In formula, PHFor gas reservoir H pressure, MPa；PHiFor the gas reservoir H pressure under initial condition, MPa.Then water drive gas reservoir substance is flat
Weigh equation are as follows:
Then water drive gas reservoir Production development curve is PH~GPRelation curve.
In the above-described embodiments, it when the linear regression coeffficient is less than or equal to predetermined threshold, is dissipated according to the relationship
Point diagram adjusts the default gas reservoir and water body radius value；Again according to the adjustment adjusted default gas reservoir and water body radius value
It calculates and obtains linear regression coeffficient；When the linear regression coeffficient is greater than predetermined threshold, according to the linear regression coeffficient
Corresponding performance graph and the geologic parameter and the means of production obtain water influx recognition result.In real work, according to substance
Equilibrium equation is it is found that PH~GPRelation curve is straight line.If the Production development curve correlation coefficient value of linear regression is greater than
90%, then it is assumed that the gas reservoir radius value and water body radius value of setting are accurate, and the water enchroachment (invasion) magnitude of calculating is reliable, can directly export knot
Fruit；If not satisfied, then returning to step 3, water body radius and gas reservoir radius value are adjusted, water enchroachment (invasion) magnitude is recalculated.Water body radius
Value can be adjusted according to the scatter plot in step 4, if scatterplot is in nose-up tendency, illustrate to assume that radius is greater than real radius value,
Assumed value should be reduced；If scatterplot falls trend under being in, illustrate that the radius assumed is less than real radius, Ying Zeng great assumed value.Constantly adjust
Whole gas reservoir radius and water body radius value, until the Production development curve correlation coefficient value of linear regression is greater than 90%.
Below by way of come for certain Fractured edge water gas reservoir to method of the present invention carry out overall description, the gas reservoir
Underlying parameter is as shown in Table 1 and Table 2.
Table 1
Table 2
Since a possibility that gas reservoir is Fractured edge water gas reservoir, and water body alters air inlet hiding from outside, is very big.From the means of production
From the point of view of, which there is production water phenomenon, and water yield is in increase tendency, shows to be likely to occur water enchroachment (invasion) phenomenon.Therefore, may be used
Using the quantitative identification method of Fractured edge water gas reservoir water influx of the present invention.For Fractured edge water gas reservoir water influx
Quantitative identification method for, top priority is that geologic parameter and creation data are collected and are arranged, including gas reservoir, gas
And parameters of formation water and Production development data, and analyzed it according to the reliability of data.It is arranged to data
After analysis, determine water influx quantitative identification method outputs and inputs parameter.Inputting parameter includes: the gas reservoir radius of hypothesis, vacation
Fixed water body radius, reservoir pressure, deviation factor for gas, fracture permeabgility, gas reservoir porosity, system compressibility, gas are viscous
Degree, reservoir effective thickness, elastic storativity ratio, interporosity flow coefficient, irreducible water saturation, formation water volume factor, production time, ground
Gas production, gas volume factor, tired gas production, tired water yield.Quantitative identification method output parameter includes: gas reservoir radius, water body
Radius, water enchroachment (invasion) magnitude.It please refers to shown in Fig. 2 and Fig. 3, t in Fig. 30Indicate initial time, h；t1Indicate first moment, h；t2Table
Show second moment, h；tj-1Indicate -1 moment of jth, h；PiIndicate initial time t0Corresponding pressure, MPa；P1Indicate first
A moment t1Corresponding pressure, MPa；P2Indicate second moment t2Corresponding pressure, MPa；Pj-1Indicate -1 moment t of jthj-1
Corresponding pressure, MPa；ΔP0Indicate initial segment pressure difference, MPa；ΔP1Indicate first segment pressure difference, MPa；ΔP2Indicate second segment pressure
Difference, MPa；ΔPj-1Indicate -1 section of pressure difference of jth, MPa.Due to being difficult to obtain gas reservoir radius and water body radius parameter value, according to crack
Property edge water gas reservoir water influx calculation formula when calculating water influx, it is necessary first to it is assumed that gas reservoir radius and water body radius, secondly by water
Then body radius, time and pressure parameter zero dimension substitute into formula and calculate the zero dimension water influx under the space Laplace
Solution, obtains the nondimensional water influx solution under real space using stehfest numerical inversion formula, finally combines
The level pressure formula that Everdingen and Hurst is proposed, has dimension for water influx, finally obtains water enchroachment (invasion) magnitude.
Thereafter it is please referred to shown in Fig. 4 again, according to output parameter as a result, calculated water enchroachment (invasion) magnitude is substituted into water drive gas reservoir object
In matter equilibrium equation, in conjunction with creation data and actual parameter, make about PH~GPThe scatter plot of relationship, linear regression go out to produce
Performance graph, and obtain the correlation coefficient value of regression straight line.If scatterplot is in nose-up tendency, illustrate that the radius assumed is greater than practical half
Diameter should reduce assumed value；If scatterplot falls trend under being in, illustrate to assume that radius is less than real radius, Ying Zeng great assumed value.It please refers to
Shown in Fig. 5, according to matter balance equation it is found that PH~GPRelation curve is straight line.In the process, if linear regression go out
Production development curve correlation coefficient is greater than 90%., then it is assumed that it is assumed that water body radius and gas reservoir radius it is reliable, the water influx of calculating
Value is accurate；If not satisfied, then returning in step 3, water body radius and gas reservoir radius value are assumed again, and carry out to water influx
It recalculates, until meeting fitting Production development linear correlation and being greater than 90%.The related coefficient of linear regression herein is
0.9856, i.e. correlation is 98.56%, meets the requirements, according to gained Production development curve, can calculate dynamic geological reserves, G
=a/b=1487 (108m3).Finally output is as a result, obtain gas reservoir radius rgFor 4400m, water body radius reFor 5280m, dynamically
Matter reserves G is 1487 × 108m3And water enchroachment (invasion) magnitude W under present conditioneIt is 2.27 × 104m3.Specific Fractured edge water gas reservoir signal
Figure is as shown in fig. 6, corresponding water enchroachment (invasion) magnitude W of each momente, gas reservoir radius rg, water body radius reAnd dynamic geological reserves G such as the following table 3
It is shown.
Table 3
It please refers to shown in Fig. 7, it is described the present invention also provides a kind of quantitative judge system of Fractured edge water gas reservoir water influx
System includes data acquisition module, computing module, building module and identification module；The data acquisition module is for obtaining geology
And creation data, geologic parameter and the means of production are obtained according to the geology and creation data and water influx quantitative judge algorithm；
The computing module is used to pass through Fractured side according to the geologic parameter and the means of production and default gas reservoir and water body radius value
Gas Reservoirs water influx calculation formula, which calculates, obtains water enchroachment (invasion) magnitude；The building module is used for according to water drive gas reservoir matter balance equation
And the water enchroachment (invasion) magnitude, performance graph is obtained by linear regression method, is calculated according to the performance graph and obtains the dynamic
The corresponding linear regression coeffficient of curve；The identification module is used for when the linear regression coeffficient is greater than predetermined threshold, according to
The performance graph and the geologic parameter and the means of production obtain water influx recognition result.
In the above-described embodiments, the computing module includes processing unit, inverting unit and level pressure unit；The processing is single
Member is used for water body radius, time and pressure parameter zero dimension in the default gas reservoir and water body radius value, and brings crack into
Property edge water gas reservoir water influx calculation formula in obtain zero dimension water influx；The inverting unit is used for according to the zero dimension water enchroachment (invasion)
Amount, the nondimensional water influx under real space is obtained by stehfest numerical inversion formula；The level pressure unit is used for root
Water enchroachment (invasion) magnitude is obtained according to level pressure formula and the nondimensional water influx.
In the above-described embodiments, the building module also includes: bringing the water enchroachment (invasion) magnitude into water drive gas reservoir material balance
In equation, and geologic parameter and the means of production is combined to establish the relationship scatter plot between gas reservoir pressure and oil in place, according to institute
It states relationship scatter plot and performance graph is obtained by linear regression algorithm, calculated according to the performance graph and obtain the performance graph
Corresponding linear regression coeffficient.
In the above-described embodiments, the identification module also includes: when the linear regression coeffficient is less than or equal to predetermined threshold
When value, the default gas reservoir and water body radius value are adjusted according to the relationship scatter plot；The quantitative judge system is according to adjustment
The default gas reservoir of adjustment afterwards and water body radius value recalculate acquisition linear regression coeffficient；Until the linear regression coeffficient
When greater than predetermined threshold, obtained according to the corresponding performance graph of the linear regression coeffficient and the geologic parameter and the means of production
Water influx recognition result.
The present invention also provides a kind of computer equipment, including memory, processor and storage on a memory and can located
The computer program run on reason device, the processor realize the above method when executing the computer program.
The present invention also provides a kind of computer readable storage medium, the computer-readable recording medium storage has in execution
State the computer program of method.
The quantitative identification method of Fractured edge water gas reservoir water influx provided by the present invention and system identification go out Fractured side
Gas Reservoirs water influx, and gas reservoir radius, water body radius and Production development curve are obtained, it is that subsequent production dynamic monitoring and formulation are closed
The development plan of reason lays the foundation.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the present invention, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or 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 counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection scope of invention.
Claims (12)
1. a kind of quantitative identification method of Fractured edge water gas reservoir water influx, which is characterized in that the method includes:
Geology and creation data are obtained, geologic parameter is obtained according to the geology and creation data and water influx quantitative judge algorithm
And the means of production；
According to the geologic parameter and the means of production and default gas reservoir and water body radius value, pass through Fractured edge water gas reservoir water influx
Calculation formula, which calculates, obtains water enchroachment (invasion) magnitude；
According to water drive gas reservoir matter balance equation and the water enchroachment (invasion) magnitude, performance graph is obtained by linear regression method, according to
The performance graph, which calculates, obtains the corresponding linear regression coeffficient of the performance graph；
When the linear regression coeffficient is greater than predetermined threshold, according to the performance graph and the geologic parameter and the means of production
Obtain water influx recognition result.
2. the quantitative identification method of Fractured edge water gas reservoir water influx according to claim 1, which is characterized in that describedly
Matter and creation data include gas reservoir, gas and parameters of formation water and Production development data.
3. the quantitative identification method of Fractured edge water gas reservoir water influx according to claim 1, which is characterized in that describedly
Matter parameter and the means of production include water influx quantitative judge algorithm input parameter and water influx quantitative judge algorithm output parameter；
The water influx quantitative judge algorithm input parameter includes: reservoir pressure, deviation factor for gas, fracture permeabgility, gas reservoir
Porosity, system compressibility, gas viscosity, reservoir effective thickness, elastic storativity ratio, interporosity flow coefficient, irreducible water saturation,
Layer water volume coefficient, production time, surface gas rate, gas volume factor, tired gas production, tired water yield；
The water influx quantitative judge algorithm output parameter includes: gas reservoir radius, water body radius, water enchroachment (invasion) magnitude.
4. the quantitative identification method of Fractured edge water gas reservoir water influx according to claim 1, which is characterized in that according to institute
Geologic parameter and default gas reservoir and water body radius value are stated, is calculated by Fractured edge water gas reservoir water influx calculation formula and obtains water enchroachment (invasion)
Magnitude includes:
By water body radius, time and pressure parameter zero dimension in the default gas reservoir and water body radius value, and bring Fractured into
Zero dimension water influx is obtained in edge water gas reservoir water influx calculation formula；
According to the zero dimension water influx, the nondimensional water enchroachment (invasion) under real space is obtained by stehfest numerical inversion formula
Amount；
Water enchroachment (invasion) magnitude is obtained according to level pressure formula and the nondimensional water influx.
5. the quantitative identification method of Fractured edge water gas reservoir water influx according to claim 1, which is characterized in that according to water
Material balance equation for gas reservoir and the water enchroachment (invasion) magnitude are driven, performance graph is obtained by linear regression method, it is bent according to the dynamic
The linear regression coeffficient that line obtains fitting a straight line includes: the water enchroachment (invasion) magnitude brought into water drive gas reservoir matter balance equation, and
The relationship scatter plot between gas reservoir pressure and oil in place is established in conjunction with geologic parameter and the means of production, according to the relationship scatterplot
Figure obtains performance graph by linear regression algorithm, and it is corresponding linear to calculate the acquisition performance graph according to the performance graph
Regression coefficient.
6. the quantitative identification method of Fractured edge water gas reservoir water influx according to claim 5, which is characterized in that the side
Method also includes:
When the linear regression coeffficient is less than or equal to predetermined threshold, the default gas reservoir is adjusted according to the relationship scatter plot
And water body radius value；
Acquisition linear regression coeffficient is recalculated according to the adjustment adjusted default gas reservoir and water body radius value；
When the linear regression coeffficient is greater than predetermined threshold, according to the corresponding performance graph of the linear regression coeffficient and institute
It states geologic parameter and the means of production obtains water influx recognition result.
7. a kind of quantitative judge system of Fractured edge water gas reservoir water influx, which is characterized in that the system includes data acquisition
Module, computing module, building module and identification module；
The data acquisition module is quantitative according to the geology and creation data and water influx for obtaining geology and creation data
Recognizer obtains geologic parameter and the means of production；
The computing module is used to pass through crack according to the geologic parameter and the means of production and default gas reservoir and water body radius value
Property edge water gas reservoir water influx calculation formula calculate obtain water enchroachment (invasion) magnitude；
The building module is used to be obtained according to water drive gas reservoir matter balance equation and the water enchroachment (invasion) magnitude by linear regression method
Performance graph is obtained, is calculated according to the performance graph and obtains the corresponding linear regression coeffficient of the performance graph；
The identification module is used for when the linear regression coeffficient is greater than predetermined threshold, with described according to the performance graph
Matter parameter and the means of production obtain water influx recognition result.
8. the quantitative judge system of Fractured edge water gas reservoir water influx according to claim 7, which is characterized in that the meter
Calculating module includes processing unit, inverting unit and level pressure unit；
The processing unit is used for water body radius, time and pressure parameter zero dimension in the default gas reservoir and water body radius value
Change, and brings acquisition zero dimension water influx in Fractured edge water gas reservoir water influx calculation formula into；
The inverting unit is used to obtain real space by stehfest numerical inversion formula according to the zero dimension water influx
Under nondimensional water influx；
The level pressure unit is used to obtain water enchroachment (invasion) magnitude according to level pressure formula and the nondimensional water influx.
9. the quantitative judge system of Fractured edge water gas reservoir water influx according to claim 7, which is characterized in that the structure
Modeling block also includes: the water enchroachment (invasion) magnitude being brought into water drive gas reservoir matter balance equation, and combines geologic parameter and production money
Material establishes the relationship scatter plot between gas reservoir pressure and oil in place, is obtained according to the relationship scatter plot by linear regression algorithm
Performance graph is obtained, is calculated according to the performance graph and obtains the corresponding linear regression coeffficient of the performance graph.
10. the quantitative judge system of Fractured edge water gas reservoir water influx according to claim 9, which is characterized in that described
Identification module also includes: when the linear regression coeffficient is less than or equal to predetermined threshold, being adjusted according to the relationship scatter plot
The default gas reservoir and water body radius value；
The quantitative judge system recalculates acquisition linearly according to the adjustment adjusted default gas reservoir and water body radius value
Regression coefficient；It is bent according to the corresponding dynamic of the linear regression coeffficient when linear regression coeffficient is greater than predetermined threshold
Line and the geologic parameter and the means of production obtain water influx recognition result.
11. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor
Calculation machine program, which is characterized in that the processor realizes any side of claim 1 to 6 when executing the computer program
Method.
12. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage has perform claim
It is required that the computer program of 1 to 6 any the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910111218.7A CN109858701B (en) | 2019-02-12 | 2019-02-12 | Quantitative identification method and system for water invasion of fractured edge water-gas reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910111218.7A CN109858701B (en) | 2019-02-12 | 2019-02-12 | Quantitative identification method and system for water invasion of fractured edge water-gas reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109858701A true CN109858701A (en) | 2019-06-07 |
| CN109858701B CN109858701B (en) | 2021-04-16 |
Family
ID=66897676
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910111218.7A Active CN109858701B (en) | 2019-02-12 | 2019-02-12 | Quantitative identification method and system for water invasion of fractured edge water-gas reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109858701B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110688764A (en) * | 2019-10-09 | 2020-01-14 | 西南石油大学 | Calculation method for water body multiple of water drive gas reservoir |
| CN111506865A (en) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | Shale gas well production fitting and predicting method |
| CN112329243A (en) * | 2020-11-09 | 2021-02-05 | 中国海洋石油集团有限公司 | Inversion method and device for water invasion process of heterogeneous reservoir high-water-cut horizontal well |
| CN112328953A (en) * | 2019-09-16 | 2021-02-05 | 中国石油天然气股份有限公司 | Water invasion identification method, device and equipment for gas well and readable storage medium |
| CN112487594A (en) * | 2019-08-23 | 2021-03-12 | 中国石油天然气股份有限公司 | Oil reservoir water body multiple calculation method and device |
| CN113268846A (en) * | 2020-02-17 | 2021-08-17 | 中国石油天然气股份有限公司 | Method, device, terminal and storage medium for determining water control strategy with water-gas reservoir |
| CN113818873A (en) * | 2021-10-20 | 2021-12-21 | 中国石油大学(北京) | Method and device for calculating dynamic reserves of gas reservoir under action of hidden gas-containing area |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760081A (en) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | Gas reservoir prediction method and system for carbonate reservoir based on pore structure characteristics |
| CN104484550A (en) * | 2014-11-06 | 2015-04-01 | 中国石油天然气股份有限公司 | Method and device for acquiring reservoir capacity of gas storage reservoir reconstructed from flooded gas reservoir |
| CN109184676A (en) * | 2018-09-21 | 2019-01-11 | 中国地质大学（武汉） | Volume evaluation method is effectively transformed in a kind of shale gas reservoir |
-
2019
- 2019-02-12 CN CN201910111218.7A patent/CN109858701B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760081A (en) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | Gas reservoir prediction method and system for carbonate reservoir based on pore structure characteristics |
| CN104484550A (en) * | 2014-11-06 | 2015-04-01 | 中国石油天然气股份有限公司 | Method and device for acquiring reservoir capacity of gas storage reservoir reconstructed from flooded gas reservoir |
| CN109184676A (en) * | 2018-09-21 | 2019-01-11 | 中国地质大学（武汉） | Volume evaluation method is effectively transformed in a kind of shale gas reservoir |
Non-Patent Citations (2)
| Title |
|---|
| 祝浪涛等: "致密油藏直井体积压裂压力分析模型", 《大庆石油地质与开发》 * |
| 窦祥骥等: "应力敏感微纳米孔隙储层压裂裂缝反演方法", 《油气地质与采收率》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112487594A (en) * | 2019-08-23 | 2021-03-12 | 中国石油天然气股份有限公司 | Oil reservoir water body multiple calculation method and device |
| CN112487594B (en) * | 2019-08-23 | 2022-11-01 | 中国石油天然气股份有限公司 | Oil reservoir water body multiple calculation method and device |
| CN112328953A (en) * | 2019-09-16 | 2021-02-05 | 中国石油天然气股份有限公司 | Water invasion identification method, device and equipment for gas well and readable storage medium |
| CN112328953B (en) * | 2019-09-16 | 2023-12-22 | 中国石油天然气股份有限公司 | Water invasion identification method, device and equipment for gas well and readable storage medium |
| CN110688764A (en) * | 2019-10-09 | 2020-01-14 | 西南石油大学 | Calculation method for water body multiple of water drive gas reservoir |
| CN113268846A (en) * | 2020-02-17 | 2021-08-17 | 中国石油天然气股份有限公司 | Method, device, terminal and storage medium for determining water control strategy with water-gas reservoir |
| CN113268846B (en) * | 2020-02-17 | 2022-11-04 | 中国石油天然气股份有限公司 | Method, device, terminal and storage medium for determining water control strategy with water-gas reservoir |
| CN111506865A (en) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | Shale gas well production fitting and predicting method |
| CN112329243A (en) * | 2020-11-09 | 2021-02-05 | 中国海洋石油集团有限公司 | Inversion method and device for water invasion process of heterogeneous reservoir high-water-cut horizontal well |
| CN112329243B (en) * | 2020-11-09 | 2023-08-25 | 中国海洋石油集团有限公司 | Heterogeneous reservoir high-water-content horizontal well water invasion process inversion method and device |
| CN113818873A (en) * | 2021-10-20 | 2021-12-21 | 中国石油大学(北京) | Method and device for calculating dynamic reserves of gas reservoir under action of hidden gas-containing area |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109858701B (en) | 2021-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109858701A (en) | The quantitative identification method and system of Fractured edge water gas reservoir water influx | |
| Eshraghi et al. | Optimization of miscible CO2 EOR and storage using heuristic methods combined with capacitance/resistance and Gentil fractional flow models | |
| Chen et al. | Robust constrained optimization of short-and long-term net present value for closed-loop reservoir management | |
| CA3079352A1 (en) | Enhancing reservoir production optimization through integrating inter-well tracers | |
| CA2997608C (en) | History matching of hydrocarbon production from heterogenous reservoirs | |
| RU2590265C2 (en) | Systems and methods for assessment of moments of penetration of fluid in locations of production wells | |
| US20080015832A1 (en) | Method for Predicting Acid Placement in Carbonate Reservoirs | |
| CN104879103A (en) | Layered water injection effect analysis method | |
| CN106460493A (en) | Method for improved design of hydraulic fracture height in a subterranean laminated rock formation | |
| GB2563241A (en) | Hydraulic fracturing simulation | |
| CN109655903B (en) | Shale layer transverse wave velocity prediction method and system | |
| Kapangaziwiri et al. | Towards revised physically based parameter estimation methods for the Pitman monthly rainfall-runoff model | |
| CN102536200A (en) | Method for predicting primary capacity of compact carbonate rock gas bearing formations | |
| CN104504472B (en) | Method and device for predicting well productivity of gas storage gas-water interaction area | |
| CN116306385B (en) | Oil reservoir fracturing imbibition energy increasing numerical simulation method, system, equipment and medium | |
| CN110397425A (en) | Coal seam gas production well flowing bottomhole pressure (FBHP) control system and control method | |
| Bui et al. | Formation damage simulation of a multi-fractured horizontal well in a tight gas/shale oil formation | |
| CN104636818B (en) | Method and device for predicting gas storage gas-oil interaction area well productivity | |
| Abdideh et al. | Stress field analysis and its effect on selection of optimal well trajectory in directional drilling (case study: southwest of Iran) | |
| CN105844011B (en) | A kind of calculation of permeability based on capillary model | |
| Qin et al. | Analysis of spatiotemporal variations of excess pore water pressure during mechanized tunneling using genetic programming | |
| Feng et al. | Effects of porous properties of rock on near-wellbore hydraulic fracture complexity | |
| US11585192B2 (en) | Method and system for reactively defining valve settings | |
| Jin | An integrated geomechanics and petrophysics study of hydraulic fracturing in naturally fractured reservoirs | |
| Liu et al. | Research on stress-sensitivity of fractured porous media |
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 |