CN108664677A - A kind of oil/gas well creation data analysis method - Google Patents
A kind of oil/gas well creation data analysis method Download PDFInfo
- Publication number
- CN108664677A CN108664677A CN201710212037.4A CN201710212037A CN108664677A CN 108664677 A CN108664677 A CN 108664677A CN 201710212037 A CN201710212037 A CN 201710212037A CN 108664677 A CN108664677 A CN 108664677A
- Authority
- CN
- China
- Prior art keywords
- yield
- pressure
- function
- double
- yield regularization
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 84
- 238000007405 data analysis Methods 0.000 title claims abstract description 18
- 238000004519 manufacturing process Methods 0.000 claims abstract description 36
- 230000000694 effects Effects 0.000 claims abstract description 21
- 239000000463 material Substances 0.000 claims abstract description 20
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 17
- 230000001360 synchronised effect Effects 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims description 12
- 239000003129 oil well Substances 0.000 claims description 9
- 206010016256 fatigue Diseases 0.000 claims description 5
- 235000013399 edible fruits Nutrition 0.000 claims description 4
- 239000000126 substance Substances 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 description 19
- 238000004458 analytical method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 230000035699 permeability Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000003825 pressing Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 239000012530 fluid Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000011067 equilibration Methods 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012517 data analytics Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Tourism & Hospitality (AREA)
- Primary Health Care (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Human Resources & Organizations (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Agronomy & Crop Science (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
A kind of oil/gas well creation data analysis method comprising：Step 1: building well yield Regularization pressure function to be analyzed according to original formation pressure, flowing bottomhole pressure (FBHP) and the daily output；Step 2: determining double-log plate drafting function according to yield Regularization pressure function, yield Regularization pressure double-log plate is drawn according to double-log plate drafting function；Step 3: being fitted to yield Regularization pressure double-log plate gross data according to measured data in the yield Regularization pressure double-log plate got, and the formation parameter of well to be analyzed is determined with pressure history fitting result according to double-log plate fitting result and yield.The problem of this method can effectively avoid the substance material balance time and not be monotonically changed with actual time, it eliminates because creation data interruption, big fluctuation explain the negative effect brought to creation data, the fitting synchronous with production history data of double-log plate data can be realized simultaneously, to the creation data Explanation Accuracy for improving creation data (the especially daily output) interruption, fluctuating larger oil/gas well.
Description
Technical field
The present invention relates to oil/gas well exploration and development technical fields, specifically, being related to a kind of oil/gas well creation data analysis
Method.
Background technology
Well testing is using permeation fluid mechanics as theoretical foundation, using test instrumentation as means, by measuring oily gas and water, the pressure of well
And the Production developments data such as yield, it studies and determines testing well and oil, the production capacity of gas and water layer, physical parameter and distinguish
Between other well or the technology of interlayer connected relation, all have to the exploitation of oil gas field new exploration areas, old area dynamic adjustment particularly important
Effect.
It is different according to the data source utilized to understand the parameter of certain certain well and corresponding reservoir, mainly there is pressure extensive
Two kinds of retrial well (referred to as pressing extensive well testing), creation data well test analysis common methods.However, existing creation data analysis method
Not applicable to part oil/gas well, this makes the explanation results that creation data is analyzed have very strong uncertain and multi-solution.
Invention content
To solve the above problems, the present invention provides a kind of oil/gas well creation data analysis method, the method includes：
Step 1: building well yield Regularization pressure to be analyzed according to original formation pressure, flowing bottomhole pressure (FBHP) and the daily output
Function or yield Regularization pseudopressure function；
Step 2: according to the yield Regularization pressure function or yield Regularization pseudopressure function, it is true based on tired yield
Determine double-log plate drafting function, according to the double-log plate drafting function draw yield Regularization pressure double-log plate or
Yield Regularization pseudopressure double-log plate；
Step 3: according to the yield Regularization pressure double-log plate or yield Regularization pseudopressure log-log graph that get
The measured data of match point is to the yield Regularization pressure double-log plate or yield Regularization pseudopressure log-log graph in version
Version is fitted, and determines well to be analyzed with pressure history fitting result according to double-log plate fitting result and yield
Relevant parameter.
According to one embodiment of present invention, in the step 1,
If the well to be analyzed is oil well, the yield Regularization pressure function is built according to following expression：
If the well to be analyzed is gas well, the yield Regularization pseudopressure function is calculated according to following expression：
Wherein, RNP indicates yield Regularization pressure, piIndicate original formation pressure, pwfBottom pressure is indicated, when t is indicated
Between, q indicates yield, ψiIndicate gas well prime stratum pseudopressure, ψwfIndicate gas well shaft bottom pseudopressure.
According to one embodiment of present invention, in the step 2,
Yield Regularization pressure derivative function is determined according to the yield Regularization pressure function, it is regular according to the yield
Change pressure function and yield Regularization pressure derivative function builds the yield Regularization pressure double-log plate；Or,
Yield Regularization pressure integral function is determined according to the yield Regularization pressure function, and is advised according to the yield
Integralization pressure integral function determines yield Regularization pressure integral derivative function, according to the yield Regularization pressure integral function
The yield Regularization pressure double-log plate is built with yield Regularization pressure integral derivative function；Or,
Yield Regularization pressure derivative function and yield Regularization pressure are determined according to the yield Regularization pressure function
Integral function, and yield Regularization pressure integral derivative function is determined according to the yield Regularization pressure integral function, according to
The yield Regularization pressure function, yield Regularization pressure derivative function, yield Regularization pressure integral function and yield rule
Integralization pressure integral derivative function builds the yield Regularization pressure double-log plate.
According to one embodiment of present invention, the yield Regularization pressure derivative function is determined according to following expression：
Wherein, RNPdIndicate that yield Regularization pressure derivative, RNP indicate that yield Regularization pressure, V indicate to tire out yield.
According to one embodiment of present invention, the yield Regularization pressure integral function is determined according to following expression：
Wherein, RNPiIndicate yield Regularization pressure integral, VeIndicate that the tired yield in the material balance time, RNP indicate production
Gauge integralization pressure, V indicate to tire out yield.
According to one embodiment of present invention, the yield Regularization pressure integral derivative letter is determined according to following expression
Number：
Wherein, RNPiIndicate yield Regularization pressure integral, RNPidIndicate that yield Regularization pressure integral derivative, V indicate
Tired yield.
According to one embodiment of present invention, in the step 3, by adjusting tired production match value and pressure match value
To adjust the fitting effect of yield Regularization pressure double-log plate.
According to one embodiment of present invention, the tired production match value is determined according to following expression：
The pressure match value is determined according to following expression：
Wherein, VMAnd RNPMTired production match value and pressure match value, V are indicated respectivelyTBAnd VSJDouble logarithmic curve is indicated respectively
The plate of middle match point tires out yield data and yield data, RNP are tired out in actual measurementTBAnd RNPSJIt indicates to be fitted in double logarithmic curve respectively
The plate yield Regularization pressure data and actual measurement yield Regularization pressure data of point.
According to one embodiment of present invention, it in the step 3, is also produced based on double-log plate fitting result
The history matching of amount and pressure data, optimizes the double-log plate fitting result by history matching, realizes figure
Version data are synchronous with production history data to be fitted, to the relevant parameter of the well to be analyzed after being optimized.
Oil/gas well creation data analysis method provided by the present invention is based on yield Regularization (quasi-) pressure and tired yield
Double-log plate carries out reservoir ginseng using oil/gas well creation data by the history matching of chart fitting, yield and pressure data
Number (such as permeability, skin factor) and pit shaft and transformation parameter (such as half length, fracture condudtiviy, SRV volumes of fracturing fracture) solution
It releases, the problem of capable of effectively avoiding the substance material balance time in this way and be not monotonically changed with actual time, eliminates because producing number
The negative effect brought is explained to creation data according to interruption, big fluctuation, to improve creation data (the especially daily output)
It is interrupted, fluctuates the creation data Explanation Accuracy of larger oil/gas well, improves yield and pressure history fitting effect, greatly improve
Explanation Accuracy, reduces multi-solution.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that understand through the implementation of the invention.The purpose of the present invention and other advantages can be by specification, rights
Specifically noted structure is realized and is obtained in claim and attached drawing.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is required attached drawing in technology description to do simple introduction：
Fig. 1 is the implementation process schematic diagram of oil/gas well creation data analysis method according to an embodiment of the invention；
Fig. 2 is double-log plate schematic diagram according to an embodiment of the invention；
Fig. 3 is creation data double-log chart fitting schematic diagram according to an embodiment of the invention；
Fig. 4 is the creation data double-log chart fitting schematic diagram of existing method；
Fig. 5 is history matching schematic diagram according to an embodiment of the invention；
Fig. 6 is the history matching schematic diagram of existing method.
Specific implementation mode
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to be applied to the present invention whereby
Technological means solves technical problem, and the realization process for reaching technique effect can fully understand and implement.It needs to illustrate
As long as not constituting conflict, each embodiment in the present invention and each feature in each embodiment can be combined with each other,
It is formed by technical solution within protection scope of the present invention.
Meanwhile in the following description, for illustrative purposes and numerous specific details are set forth, to provide to of the invention real
Apply the thorough understanding of example.It will be apparent, however, to one skilled in the art, that the present invention can not have to tool here
Body details or described ad hoc fashion are implemented.
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 is different according to the data source utilized to understand the parameter of certain certain well and corresponding reservoir, mainly there is pressure extensive
Two kinds of retrial well (referred to as pressing extensive well testing), creation data well test analysis common methods.
It is a kind of traditional well-logging method to press extensive well testing, is recorded after shaft bottom or well head placement force well, closing well when applying
Pressure change.The advantages of this method is that the close precision of pressure gauge sampling site is high, and explanation results can reflect that oil/gas well opening well and making production early stage believes
Breath.But it since shale gas reservoir reservoir permeability is extremely low, presses extensive well testing that can only record short-term partial data, can not reflect shale gas
The remote well informations such as the interference of well crack, crack boundary.Therefore the information for pressing extensive well testing that can be provided in shale gas well has very much
Limit.
Creation data well test analysis is to carry out well test analysis using the dynamic data of opening well and making production, and main method has the side PNR
Method (pressure normalizes method of flow, also referred to as Blasingame methods), Agarwal-Gardner methods, RNP method (flows
Normalize pressure method).Creation data well-logging method advantage is not need closing well, is not affected the normal production, and the length of record is passed through
Phase creation data explains well and reservoir characteristic parameter, can reflect the remote well information of oil/gas well, therefore this method is highly suitable for
The oil/gas well of low permeability reservoir.
But there are still two problems for above-mentioned existing creation data analysis method.First, existing creation data analysis side
Method mainly carries out well test analysis by the fitting of (quasi-) pressure and (quasi-) pressure derivative and material balance time double-log plate,
It is the relationship in analysis (quasi-) pressure and (quasi-) pressure derivative and material balance time.When the material balance time is with true production
Between monotonic increase when, double-log chart fitting is preferable.But since gas well is frequently accompanied by frequent production system in process of production
Changing (such as switch well, replaces choke size, inter-well interference at drain) or artificial origin leads to production status suddenly change, therefore
Creation data quality is usually poor, and pressure and data on flows are frequently present of interruption, fluctuation even missing, lead to the material balance time
It is very big so that it cannot monotonic increase with true production time fluctuating range.Pressure and pressure derivative and substance in this case
The fitting effect of equilibration time double-log plate can not then ensure that then explanation results have very strong uncertain and multi-solution.
Second, when existing creation data analysis method is mainly by (quasi-) pressure and (quasi-) pressure derivative with material balance
Between the fitting of double-log plate carry out well test analysis, the fitting and the history matching of yield and pressure data of double-log plate be
It independently carries out, can not ensure the synchronous chart fitting ensured and history matching effect.
Above-mentioned two problems so that existing creation data analysis method is not applicable to part oil/gas well, double-log plate
Fitting and yield can not ensure that synchronous fitting, explanation results have very strong uncertain and multi-solution with pressure history fitting.
For the above problem in the presence of the prior art, the present invention provides a kind of new oil/gas well creation data analyses
Method, this method are analyzed using yield Regularization (quasi-) pressure and its derivative and the double-log relationship plate of tired yield, are protected
The abscissa (tired yield) in double-log relationship plate has been demonstrate,proved with the monotonic increase of actual time, to solve because frequent raw
Production system changes (such as switch well, drain, replacement choke size, inter-well interference) or creation data caused by artificial origin is interrupted,
Big fluctuation and existing creation data Well Test Data Analysis Method is difficult to the problem accurately explained, and ensure that double-log plate and yield
Fitting effect synchronous with pressure history.
Fig. 1 shows the implementation process schematic diagram of oil/gas well creation data analysis method in the present embodiment.
The oil/gas well creation data analysis method that the present embodiment is provided needs to use bottom pressure during realizing
With the relational model of the tired yield in the production time, therefore the tired production in bottom pressure and production time is obtained firstly the need of structure
The relational model of amount.
For the variable-flow and transformation force data changed over time, when can calculate production by following expression to t
The tired yield carved：
Wherein, V indicates to tire out yield, and t indicates that production time, q indicate the daily output.
By taking a bite straight well oil well in infinitely great homogeneous formation as an example, strata pressure is distributed and is represented by when varying flow rate produces：
Wherein, p indicates pressure, piIndicate that original formation pressure, r indicate the point in stratum at a distance from pit shaft, B indicates body
Product coefficient, h indicate that formation thickness, k indicate that permeability, μ indicate gas viscosity.
χ indicates piezometric conductivity, and following expression may be used and be calculated：
Wherein, φ indicates porosity, CtIndicate the compressed coefficient.
Exist when the time is larger：
Expression formula (2) can be expressed as at this time：
Wherein, t1Indicate one level off to 0 time point.
Expression formula (1) is substituted into expression formula (5), bottom pressure can be expressed as：
Wherein, pwfIndicate bottom pressure,Indicate average product, a, b and c indicate equation coefficient.
From expression formula (6) as can be seen that for a bite straight well oil well in infinitely great homogeneous formation, when flowing reaches radial flow
When, it is linear relationship between bottom pressure and the logarithm of yield.
By taking a bite vertically fractured well in infinitely great homogeneous formation as an example, strata pressure is distributed and can indicate when varying flow rate produces
For：
Wherein, xfIndicate that fracture half-length, x and y indicate that any point is at a distance from well in stratum.Erf (x) indicates error letter
Following expression expression may be used in number：
Its bottom pressure can be expressed as：
It is smaller when the time according to the property of error function, i.e.,When, exist at this time：
At this point, expression formula (9) can be approximated to be：
Expression formula (1) is substituted into expression formula (11), bottom pressure can be expressed as：
From expression formula (12) as can be seen that for vertically fractured well in infinitely great homogeneous formation, when flowing reaches linear flow
When, it is linear relationship between bottom pressure and the square root of tired yield.
For gas well, the pressure in expression formula (6) and (12) is replaced using the pseudopressure of expression formula (13), can also obtain as
The equation of expression formula (6) and expression formula (12), to indicate the relationship between pseudopressure and tired yield.
Wherein, ψ indicates that gas well pseudopressure, z indicate Gas Compression Factor, p0Indicate that reference pressure (is preferably air
Pressure).
As shown in Figure 1, the oil/gas well that the present embodiment is provided generates data analysing method basis first in step S101
Original formation pressure, flowing bottomhole pressure (FBHP) and the daily output build the quasi- pressure of well yield Regularization pressure function or yield Regularization to be analyzed
Force function.
In the present embodiment, for the oil/gas well that the daily output and pressure change over time, this method according to acquisition primitively
Stressor layer, flowing bottomhole pressure (FBHP), the daily output calculate yield Regularization pressure or yield Regularization pseudopressure.If unmeasured flowing bottomhole pressure (FBHP)
Data, party's rule can be converted the casing pressure or oil pressure data that well head measures to flowing bottomhole pressure (FBHP) by pressure convert.
Specifically, if well to be analyzed is oil well, this method builds yield rule advantageously according to following expression
Integralization pressure function：
If well to be analyzed is gas well, similarly, it is regular that our rule builds yield advantageously according to following expression
Change pseudopressure function：
Wherein, piIndicate original formation pressure, pwfIndicate that bottom pressure, t indicate that time, q indicate yield, ψiIndicate gas well
Prime stratum pseudopressure, ψwfIndicate gas well shaft bottom pseudopressure.For oil well, RNP indicates yield Regularization pressure；For gas
For well, RNP then indicates yield Regularization pseudopressure.
Since for oil well and gas well, the creation data analytic process of the two is similar, therefore with oil well in the present embodiment
For carry out subsequent explanation.
As shown in Figure 1, in the present embodiment, after obtaining well yield Regularization pressure function to be analyzed, this method can be
Double-log plate drafting function, and basis in step s 103 are determined according to above-mentioned yield Regularization pressure function in step S102
Above-mentioned double-log plate drafting function builds yield Regularization pressure double-log plate.Specifically, in the present embodiment, this method
The double-log plate drafting function determined preferably includes：Yield Regularization pressure derivative function, yield Regularization pressure
Integral function, yield Regularization pressure integral derivative function and yield Regularization pressure function itself.
Wherein, in the present embodiment, following expression carry out table preferably may be used in yield Regularization pressure derivative function
Show：
Yield Regularization pressure integral function preferably may be used following expression and be indicated：
Yield Regularization pressure integral derivative function preferably may be used following expression and be indicated：
Wherein, RNPdIndicate yield Regularization pressure derivative, RNPiIndicate yield Regularization pressure integral, RNPidIndicate production
Gauge integralization pressure integral derivative, V indicate to tire out yield, VeIndicate the tired yield in the material balance time.
In the present embodiment, for oil well, material balance time teFollowing expression may be used to be calculated：
And for gas well, material balance time teFollowing expression then may be used to be calculated：
In the tired yield V of material balance timeeIt then can be by by above-mentioned substance equilibration time teIt substitutes into expression formula (1)
It is calculated.
For given oil/gas well type, Oil-gas Accumulation Types, inner boundary type, outer boundary type, fluid type, by upper
The double-log plate drafting function stated the relational model for the pressure and tired yield established in part and obtained is with tired yield
Abscissa can draw 4 curves, i.e., on log-log coordinate：Yield Regularization pressure curve, Regularization pressure derivative are bent
Line, yield Regularization pressure integral curve and yield Regularization pressure integral derivative curve, to obtain it is as shown in Figure 2 with
Tired yield is abscissa, using yield Regularization pressure as the double-log plate schematic diagram of coordinate.
It should be pointed out that although the method that the present embodiment is provided can draw yield Regularization pressure, Regularization
Pressure derivative, yield Regularization pressure integral and yield Regularization pressure 4 curves of integral derivative, but when practical application in order to
Convenient, this method can also be only regular using 2 curves of yield Regularization pressure and Regularization pressure derivative or yield simultaneously
Change pressure integral and yield Regularization pressure 2 curves of integral derivative.
When carrying out oil/gas well well test analysis using creation data, needs to carry out chart fitting and yield, the history of pressure are quasi-
It closes.As shown in Figure 1, the method yield Regularization pressure double-log that basis is got in step S104 that the present embodiment is provided
Measured data in plate is fitted yield Regularization double-log plate.
Specifically, for oil gas well yield Regularization pressure double-log plate, the theoretical mould on double-log plate is dragged
Type curve is allowed to be fitted with the double logarithmic curve of measured data.During chart fitting, the horizontal distance that plate drags is claimed
For tired yield match value, vertical range is known as yield Regularization pressure match value.
Wherein, in the present embodiment, tired production match value is defined as follows：
Yield Regularization pressure match value is defined as follows：
Wherein, VMAnd RNPMIt indicates to tire out yield match value and yield Regularization pressure match value, V respectivelyTBAnd VSJTable respectively
Show that the plate of match point in double logarithmic curve tires out yield data and yield data, RNP are tired out in actual measurementTBAnd RNPSJIndicate double right respectively
The plate yield Regularization pressure data of match point and actual measurement yield Regularization pressure data in number curve.
In the present embodiment, the fit procedure of yield Regularization pressure double-log plate be by adjusting tired yield match value and
Yield Regularization pressure match value makes fitting effect reach the process of better effects to adjust fitting effect.Work as fitting effect
When reaching preset requirement, the fit procedure of double-log plate namely stopping.
In order to enable the knot of the relevant parameter (such as reservoir, pit shaft and/or transformation parameter) of finally obtained well to be analyzed
Fruit is more accurate and reliable, as shown in Figure 1, the method that the present embodiment is provided can also be quasi- according to double-log plate in step S105
It closes result and carries out historical data fitting, to realize fitting synchronous with historical data to double-log plate, after obtaining optimization
Well to be analyzed relevant parameter.
Specifically, in the present embodiment, this method can obtain including that one group of tired yield is quasi- by this chart fitting process
Conjunction value VM, yield Regularization pressure match value RNPMAnd the reservoir of well to be analyzed, pit shaft, the chart fitting result that parameter is transformed.
For the creation data of any moment, yield V is tired out in the actual measurement that the corresponding moment can be calculated using expression formula (1)SJ, utilize
Expression formula (21) can tire out yield V according to above-mentioned actual measurementSJAnd tired yield match value VMIt is calculated on double-log plate corresponding
Tired yield plate value VTB。
Obtaining above-mentioned tired yield plate value VTBAfterwards, it can be found on double-log plate and tire out yield plate value with this
VTBThe plate value RNP of corresponding yield Regularization pressureTB(i.e. plate yield Regularization pressure data).Utilize expression formula
(22), according to the plate value RNP of above-mentioned yield Regularization pressureTBWith yield Regularization pressure match value RNPMIt can be calculated
The measured value RNP of corresponding yield Regularization pressureSJ(surveying the regular words pressure data of constant).
Since the daily output at corresponding moment is it is known that can be according to according to yield Regularization pressure using expression formula (14)
Measured value RNPSJThe fitting data of bottom pressure is calculated.Simultaneously as the actual measurement bottom pressure at corresponding moment is it is known that utilize
Expression formula (14) can be according to the shaft bottom producing pressure differential (i.e. reset pressure subtracts bottom pressure) and yield Regularization that actual measurement obtains
The measured value RNP of pressureSJDaily output fitting data is calculated.According to the daily output fitting data at the moment and last moment
Tired yield, the tired yield fitting data at the moment can also be calculated.
Based on same principle, the pressure that this method can obtain each moment from time 0 to production time most end point is quasi-
Close data, daily output fitting data and tired yield fitting data.After obtaining above-mentioned fitting data, this method can be by above-mentioned pressure
Power fitting data, daily output fitting data and tired yield fitting data and observed pressure data, actual measurement daily output data and
Actual measurement is tired out yield data and is fitted respectively.Wherein, if it is possible to the default fitting precision requirement of satisfaction (such as fitting data and reality
The absolute value of the difference of measured data is less than or equal to preset difference value threshold value), then then stopping being fitted, and this double-log plate is intended
Obtained relevant parameter (such as reservoir parameter, wellbore parameters and/or transformation parameter) is closed as final explanation results.And such as
Fruit cannot be satisfied default fitting precision requirement, then then adjust double-log plate fitting parameter and re-start chart fitting and
History matching process presets fitting precision requirement until meeting.
It should be pointed out that in the present embodiment, according to actual needs, double-log plate fit procedure used by this method
And historical data fit procedure can also be same as the prior art, therefore the phase no longer to step S104 and step S105 herein
Hold inside the Pass and is repeated.
In the present embodiment, using double-log chart fitting and history matching, the available well to be analyzed of this method
Relevant parameter preferably includes：Porosity, permeability, original formation pressure, formation temperature, reservoir thickness, fluid components, gas
Density, phase percolation curve (if selecting gas and water two-phase), crack series, fracture half-length, fracture interval, fracture condudtiviy and page
The parameters such as the corresponding gas supply area of rock gas well.Certainly, according to actual needs, the relevant parameter of the obtained well to be analyzed of this method
Can also only include a certain item or a few items in item listed above, the invention is not limited thereto.
Illustrate the application of the present invention by taking certain oil field a bite gas well as an example.The well is horizontal well, finishing drilling well depth 4168.0m, water
Flat segment length 1533m, 22 sections of pressure break, reservoir reset pressure are 36.7MPa, 81 DEG C of formation temperature, porosity 4%.The well produces
17184 hours, the daily output and pressure data were every 24 hours 1 points, totally 711 groups of data.
Using the method for the present invention, well test analysis is carried out to the creation data of the well, as shown in Figure 3 carries out creation data
It handles and is simultaneously fitted with theoretical plate, fitting result is that crack is averaged half a length of 77m, in-place permeability 0.024md, pit shaft
Skin factor 0.1, stratigraphic boundary are 1184m × 307m, tire out production match value VMIt is 9.46, pressure match value PMFor 0.054MPa-1。
In order to compare, using existing method, well test analysis is carried out to the creation data of the well, it as shown in Figure 5 will production
Data handle and be fitted with theoretical plate, and fitting result is that crack is averaged half a length of 67.6m, and in-place permeability is
0.028md, the wellbore skin factor 0.1, stratigraphic boundary are 1172m × 275m, time match value TMFor 14.19 (1/hr), pressure
Match value PMFor 0.0625MPa-1.
It is obtained according to chart fitting as a result, respectively apply the present invention method and existing method, carry out the daily output and
The history matching of pressure data, it is as a result as shown in Figure 4 and Figure 6 respectively.Wherein, if history matching obtained by the method for the present invention
As shown in figure 4, using the history matching of existing method, the results are shown in Figure 6.
By fig. 4 to fig. 6 as can be seen that since the daily output mutated the well at 4000 and 8000 hours, lead to it
The material balance time generates big fluctuation with actual time, destroys the monotonicity between material balance time and actual time,
When therefore being explained with existing method, go out on double-log theoretical curve of the yield Regularization pseudopressure with the material balance time
Apparent inflection point has been showed (at the corresponding about 4000 hour material balance time).Due to being different trend before and after the inflection point, so
Preferable fitting effect can only be all obtained before inflection point or after inflection point on Fig. 5 and Fig. 6, and is unable to get whole curve
Good fitting effect.Therefore the explanation results obtained by existing method are with a low credibility, and the multi-solution of result is stronger.
And the method for present invention the problem of avoiding the material balance time and not being monotonically changed with actual time, therefore in Fig. 3
Yield Regularization pseudopressure and fitting effect on the double-log plate of tired yield it is very good, also improve daily output shown in Fig. 4
Amount, tired yield and pressure history fitting effect.Therefore the credible result degree that method through the invention obtains is higher, and reduction finishes
The multi-solution of fruit.
By comparison as can be seen that the log-log graph proposed by the present invention based on yield Regularization (quasi-) pressure and tired yield
The creation data well test analysis method of version improves the problem of avoiding the material balance time and be not monotonically changed with actual time
Oil gas well production and pressure history fitting effect, substantially increase Explanation Accuracy, reduce multi-solution, to demonstrate the present invention
Correctness and practicability.The present invention has been improved particularly the production that creation data (the especially daily output) fluctuates larger oil/gas well
Data Explanation Accuracy, when creation data (the especially daily output) variation is smaller, this method is equally applicable, therefore the present invention has
Wide applicability.
As can be seen that the oil/gas well creation data analysis method that the present embodiment is provided is advised based on yield from foregoing description
The double-log plate of integralization (quasi-) pressure and tired yield passes through chart fitting, yield and number pressure using oil/gas well creation data
According to history matching carry out reservoir parameter (such as permeability, skin factor) and well and transformation parameter (such as fracturing fracture half be long, crack
Flow conductivity, SRV volumes) it explains, it can effectively avoid what the substance material balance time was not monotonically changed with actual time in this way
Problem is eliminated because creation data interruption, big fluctuation explain the negative effect brought to creation data, to improve production
Data (the especially daily output) interruption, the creation data Explanation Accuracy for fluctuating larger oil/gas well, improve yield and pressure are gone through
History fitting effect, substantially increases Explanation Accuracy, reduces multi-solution.
It should be understood that disclosed embodiment of this invention is not limited to specific structure disclosed herein or processing step
Suddenly, the equivalent substitute for these features that those of ordinary skill in the related art are understood should be extended to.It should also be understood that
It is that term as used herein is used only for the purpose of describing specific embodiments, and is not intended to limit.
" one embodiment " or " embodiment " mentioned in specification means the special characteristic described in conjunction with the embodiments, structure
Or characteristic includes at least one embodiment of the present invention.Therefore, the phrase " reality that specification various places throughout occurs
Apply example " or " embodiment " the same embodiment might not be referred both to.
Although above-mentioned example is used to illustrate principle of the present invention in one or more application, for the technology of this field
For personnel, without departing substantially from the principle of the present invention and thought, hence it is evident that can in form, the details of usage and implementation
It is upper that various modifications may be made and does not have to make the creative labor.Therefore, the present invention is defined by the appended claims.
Claims (9)
1. a kind of oil/gas well creation data analysis method, which is characterized in that the method includes：
Step 1: building well yield Regularization pressure function to be analyzed according to original formation pressure, flowing bottomhole pressure (FBHP) and the daily output
Or yield Regularization pseudopressure function；
Step 2: according to the yield Regularization pressure function or yield Regularization pseudopressure function, determined based on tired yield double
Logarithm plate drafting function draws yield Regularization pressure double-log plate or yield according to the double-log plate drafting function
Regularization pseudopressure double-log plate；
Step 3: according in the yield Regularization pressure double-log plate or yield Regularization pseudopressure double-log plate got
Measured data the yield Regularization pressure double-log plate or yield Regularization pseudopressure double-log plate are fitted,
And the relevant parameter of well to be analyzed is determined with pressure history fitting result according to double-log plate fitting result and yield.
2. the method as described in claim 1, which is characterized in that in the step 1,
If the well to be analyzed is oil well, the yield Regularization pressure function is built according to following expression：
If the well to be analyzed is gas well, the yield Regularization pseudopressure function is calculated according to following expression：
Wherein, RNP indicates yield Regularization pressure or yield Regularization pseudopressure, piIndicate original formation pressure, pwfIndicate well
Bottom pressure, t indicate that time, q indicate yield, ψiIndicate gas well prime stratum pseudopressure, ψwfIndicate gas well shaft bottom pseudopressure.
3. method as claimed in claim 1 or 2, which is characterized in that in the step 2,
Yield Regularization pressure derivative function is determined according to the yield Regularization pressure function, according to the yield Regularization pressure
Force function and yield Regularization pressure derivative function build the yield Regularization pressure double-log plate；Or,
Yield Regularization pressure integral function is determined according to the yield Regularization pressure function, and according to the yield Regularization
Pressure integral function determines yield Regularization pressure integral derivative function, according to the yield Regularization pressure integral function and production
Gauge integralization pressure integral derivative function builds the yield Regularization pressure double-log plate；Or,
Yield Regularization pressure derivative function and yield Regularization pressure integral are determined according to the yield Regularization pressure function
Function, and yield Regularization pressure integral derivative function is determined according to the yield Regularization pressure integral function, according to described
Yield Regularization pressure function, yield Regularization pressure derivative function, yield Regularization pressure integral function and yield Regularization
Pressure integral derivative function builds the yield Regularization pressure double-log plate.
4. method as claimed in claim 3, which is characterized in that determine that the yield Regularization pressure is led according to following expression
Number function：
Wherein, RNPdIndicate that yield Regularization pressure derivative, RNP indicate that yield Regularization pressure, V indicate to tire out yield.
5. method as described in claim 3 or 4, which is characterized in that determine the yield Regularization pressure according to following expression
Power integral function：
Wherein, RNPiIndicate yield Regularization pressure integral, VeIndicate that the tired yield in the material balance time, RNP indicate yield rule
Integralization pressure, V indicate to tire out yield.
6. the method as described in any one of claim 3~5, which is characterized in that determine the yield according to following expression
Regularization pressure integral derivative function：
Wherein, RNPiIndicate yield Regularization pressure integral, RNPidIndicate that yield Regularization pressure integral derivative, V indicate tired production
Amount.
7. such as method according to any one of claims 1 to 6, which is characterized in that in the step 3, by adjusting tired production
Match value and yield Regularization pressure match value adjust the fitting effect of yield Regularization pressure double-log plate.
8. method as claimed in claim 8, which is characterized in that determine the tired production match value according to following expression：
The yield Regularization pressure match value is determined according to following expression：
Wherein, VMAnd RNPMTired production match value and yield Regularization pressure match value, V are indicated respectivelyTBAnd VSJIndicate double right respectively
The plate of match point tires out yield data in number curve and yield data, RNP are tired out in actual measurementTBAnd RNPSJDouble logarithmic curve is indicated respectively
The plate yield Regularization pressure data and actual measurement yield Regularization pressure data of middle match point.
9. such as method according to any one of claims 1 to 8, which is characterized in that in the step 3, be also based on double-log
Chart fitting result carries out the history matching of yield and pressure data, by history matching come to the double-log chart fitting knot
Fruit optimizes, and realizes that plate data are synchronous with production history data and is fitted, to the phase of the well to be analyzed after being optimized
Related parameter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710212037.4A CN108664677B (en) | 2017-04-01 | 2017-04-01 | Oil and gas well production data analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710212037.4A CN108664677B (en) | 2017-04-01 | 2017-04-01 | Oil and gas well production data analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108664677A true CN108664677A (en) | 2018-10-16 |
| CN108664677B CN108664677B (en) | 2021-08-27 |
Family
ID=63784171
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710212037.4A Active CN108664677B (en) | 2017-04-01 | 2017-04-01 | Oil and gas well production data analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108664677B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110610288A (en) * | 2019-08-01 | 2019-12-24 | 裴雪皓 | Intelligent system analysis method for oil and gas well production data |
| CN111502652A (en) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | Yield decreasing and production dynamic prediction method for three-hole medium gas reservoir horizontal well |
| CN111550224A (en) * | 2020-04-27 | 2020-08-18 | 成都理工大学 | Hydraulic fracturing transformation system and method for unconventional reservoir oil and gas well |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0413873D0 (en) * | 2003-07-25 | 2004-07-21 | Schlumberger Holdings | Evaluation of fracture geometries in rock formations |
| CN105569646A (en) * | 2014-11-05 | 2016-05-11 | 中国石油化工股份有限公司 | Oil and gas well technical recoverable reserves prediction method |
| CN105760645A (en) * | 2014-12-18 | 2016-07-13 | 中国石油化工股份有限公司 | Well test interpretation method based on pressure recovery well test and production data well test |
-
2017
- 2017-04-01 CN CN201710212037.4A patent/CN108664677B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0413873D0 (en) * | 2003-07-25 | 2004-07-21 | Schlumberger Holdings | Evaluation of fracture geometries in rock formations |
| CN105569646A (en) * | 2014-11-05 | 2016-05-11 | 中国石油化工股份有限公司 | Oil and gas well technical recoverable reserves prediction method |
| CN105760645A (en) * | 2014-12-18 | 2016-07-13 | 中国石油化工股份有限公司 | Well test interpretation method based on pressure recovery well test and production data well test |
Non-Patent Citations (2)
| Title |
|---|
| BO SONG: "Rate-Normalized Pressure Analysis for Determination of Shale Gas Well Performance", 《THE NORTH AMERICAN UNCONVENTIONAL GAS CONFERENCE AND EXHIBITION》 * |
| WEI PANG: "Production Performance Modeling of Shale Gas Wells with Non-Uniform Fractures Based on Production Logging", 《SPE ANNUAL TECHNICAL CONFERENCE AND EXHIBITION》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110610288A (en) * | 2019-08-01 | 2019-12-24 | 裴雪皓 | Intelligent system analysis method for oil and gas well production data |
| CN111550224A (en) * | 2020-04-27 | 2020-08-18 | 成都理工大学 | Hydraulic fracturing transformation system and method for unconventional reservoir oil and gas well |
| CN111502652A (en) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | Yield decreasing and production dynamic prediction method for three-hole medium gas reservoir horizontal well |
| CN111502652B (en) * | 2020-07-01 | 2021-04-23 | 西南石油大学 | Yield decreasing and production dynamic prediction method for three-hole medium gas reservoir horizontal well |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108664677B (en) | 2021-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109594968B (en) | Shale gas multistage fracturing horizontal well post-fracturing fracture parameter evaluation method and system | |
| CN107506948B (en) | Shale oil gas comprehensive yield analysis method based on dynamic drainage volume | |
| Clarkson | Production data analysis of unconventional gas wells: Review of theory and best practices | |
| CN110321595B (en) | Fault sealing evaluation method for extracting static quality coefficient by logging | |
| CN106150477B (en) | A kind of method of the single well controlled reserves of determining fracture-pore reservoir | |
| CN108664679A (en) | A kind of oil/gas well creation data analysis method | |
| CN109033541B (en) | Method for evaluating heterogeneity of pressed shale gas reservoir based on EUR | |
| CN103912269B (en) | Method for determining formation fracture pressure gradient logging of shale gas reservoir | |
| CN107462936B (en) | Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law | |
| CN108518218B (en) | Unconventional oil and gas reservoir multi-section fracturing horizontal well single well dynamic reserve determination method | |
| CN107907910B (en) | Different lithologic reservoir transverse wave well logging determination method | |
| CN108614902B (en) | Shale gas well production data analysis method with changed production system | |
| CN104899411B (en) | A kind of reservoir productivity prediction model method for building up and system | |
| CN105930932B (en) | The acquisition methods of shale gas-bearing formation standardization open-flow capacity based on gassiness index | |
| CN104500016B (en) | Utilize the new method of shale gas reservoir pressing crack construction pressure drop segment data analysis reservoir properties | |
| CN108664677A (en) | A kind of oil/gas well creation data analysis method | |
| CN108133086A (en) | Water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir | |
| CN108319738A (en) | A kind of shale gas well yield prediction technique | |
| CN107503739B (en) | A kind of pressure monitoring method differentiating inflow direction for horizontal well | |
| CN109063232A (en) | One kind being used for shale gas well reservoir "sweet spot" evaluation method | |
| CN105822298A (en) | Method for acquiring absolute open flow of shale gas layer based on gas productivity index | |
| CN110043254B (en) | Method for obtaining stratum effective permeability based on cable stratum test data | |
| CN111927417A (en) | Shale gas staged fracturing horizontal well group reserve utilization condition evaluation method | |
| CN111764892B (en) | Shale gas well variable pressure and variable production stage production system adaptability judging method | |
| CN106096858B (en) | A kind of Shallow Heavy Oil Reservoir underground crude oil viscosity acquisition methods |
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 |