CN105183929B - A kind of computational methods of producing well instantaneous flow - Google Patents
A kind of computational methods of producing well instantaneous flow Download PDFInfo
- Publication number
- CN105183929B CN105183929B CN201510323955.5A CN201510323955A CN105183929B CN 105183929 B CN105183929 B CN 105183929B CN 201510323955 A CN201510323955 A CN 201510323955A CN 105183929 B CN105183929 B CN 105183929B
- Authority
- CN
- China
- Prior art keywords
- pressure
- derivative
- flow
- amplitude
- producing well
- 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.)
- Active
Links
Abstract
An embodiment of the present invention provides a kind of computational methods of producing well instantaneous flow, the amplitude of pressure first derivative is determined by calculating the first derivative of pressure data, and the time that the amplitude of each pressure first derivative occurs corresponds to a flow changed time, so according to the amplitude of pressure first derivative can be divided flowing period, and the amplitude of pressure first derivative and the average proportions coefficient of flow volume change values are calculated, finally calculate producing well instantaneous flow.By using above-mentioned technical proposal the error brought using reservoir model is avoided without establishing reservoir model；This method can be used for the arbitrary period of oil reservoir flowing simultaneously, overcome the drawbacks of production index method can be used only in oil reservoir stable state or quasi-stable state；This method has a wide range of application, and the calculated instantaneous flow result precision of institute is high.
Description
Technical field
The present embodiments relate to oil reservoir technical field more particularly to a kind of computational methods of producing well instantaneous flow.
Background technology
The data on flows of producing well is the important parameter of oil field development, is joined to yield monitoring, the reservoir characteristics of producing well
Number is explained and reservoir model history matching is most important.In Oil Field, every mouthful of producing well is flowed using flowmeter
Gauge amount is unpractical, because flowmeter is expensive, gross investment is larger, and precision is not high, inconvenient maintenance.Therefore, it produces
The instantaneous flow data of well often lack seriously.
Existing Oil Field, the main total output that several mouthfuls or tens mouthfuls wells are measured using measuring station, by adjusting flow,
The flow of different producing wells can be monitored.Specific practice is：The total flow of several mouthfuls or tens mouthfuls wells is measured with measuring station, then is carried out
Flow is split point.Such method of calculating flux has prodigious error.What main cause was that this method surveyed is that a bite producing well exists
Average flow rate in a period of time.Moreover, there are many quantity of producing well, each well is required for carrying out flow rate test, therefore to same
For a bite producing well, the time interval of flow rate test is very long twice, and producing well capacity meeting in the meantime, great changes will take place,
The changes in flow rate that cannot reflect producing well in real time, causes the error of flow monitoring larger.
In addition, the Oil Field also having calculates the flow of producing well using the production index PI of each well.The production index
(PI) the ratio between production Q and pressure drop Δ p when oil reservoir reaches stable state or quasi-stable state, PI=Q/ Δ p, wherein Δ p=P are defined ase-
Pwf, PeIt is oil reservoir average pressure, PwfIt is flowing bottomhole pressure (FBHP).According to known oil reservoir average pressure PeWith flowing bottomhole pressure (FBHP) PwfIt can calculate
The pressure drop Δ p of different moments, according to the known production well flow index, flow can be calculated by following formula：QΔ
p×PI.It is also larger that flow error is calculated using the production index.The production index can be used only in stable state or quasi-stable state period, Bu Nengji
Calculate the transient change of flow.And oil reservoir average pressure can change at any time, need constantly to update.It counts with this method
When calculating flow, the production index must be known, but produce well capacity and can change at any time, so the production index is not yet
It is constant.The variation of the production index can equally bring prodigious error to flow rate calculation.
Invention content
The purpose of the present invention is to propose to a kind of computational methods of producing well instantaneous flow, to solve to use existing producing well
Instantaneous flow calculates the instantaneous flow value that method is calculated and there are problems that large error.
For this purpose, the present invention uses following technical scheme：
A kind of computational methods of producing well instantaneous flow, including step：
S1, original pressure data is obtained；
S2, the original pressure data is pre-processed, obtains analyzable pressure data；
S3, the first derivative for calculating the pressure data, and determine the amplitude of pressure first derivative；
S4, flowing period is divided according to the amplitude of the pressure first derivative, wherein the amplitude of each pressure first derivative
The time of appearance corresponds to a flow changed time；
S5, calculate the pressure first derivative amplitude and flow volume change values average proportions coefficient；
S6, producing well instantaneous flow is calculated.
Further, the first derivative of pressure is found out by following formula：
Wherein, pderivativeFor the first derivative of pressure, p1And p2For two adjacent pressure, t1For p1The corresponding time,
t2For p2The corresponding time.
Further, the average proportions coefficient is found out by following formula：
Wherein, b is average proportions coefficient, and Q is total output, q1For initial flow, Δ tiIt is q for instantaneous flowiWhen it is corresponding
Production time interval, AjFor the amplitude of pressure first derivative, n is the number of different instantaneous flows.
Further, the producing well instantaneous flow is found out by following formula：
Wherein, qiFor producing well instantaneous flow, q1For initial flow, AjFor the amplitude of pressure first derivative, b is average specific
Example coefficient.
Further, which is characterized in that the original pressure data is pre-processed in S2, including：
Remove noise in the original pressure data and singular point and amount of compressed data.
The computational methods of a kind of producing well instantaneous flow provided in an embodiment of the present invention, by the single order for calculating pressure data
Derivative determines the amplitude of pressure first derivative, and the time that the amplitude of each pressure first derivative occurs corresponds to a flow hair
The time for changing so according to the amplitude of pressure first derivative can be divided flowing period, and calculates the width of pressure first derivative
The average proportions coefficient of value and flow volume change values, finally calculates producing well instantaneous flow.By using above-mentioned technical proposal, nothing
Reservoir model need to be established, the error brought using reservoir model is avoided；This method can be used for the arbitrary of oil reservoir flowing simultaneously
Period overcomes the drawbacks of production index method can be used only in oil reservoir stable state or quasi-stable state；This method has a wide range of application, and is calculated
The instantaneous flow result precision gone out is high.
Description of the drawings
Fig. 1 is a kind of flow chart of the computational methods for producing well instantaneous flow that the embodiment of the present invention one provides；
Fig. 2 is the image that the pressure that the embodiment of the present invention one provides changes over time；
Fig. 3 be the embodiment of the present invention one provide first derivative is sought to pressure in Fig. 2 after obtained pressure single order lead
The image that number changes over time；
Fig. 4 is a kind of flow chart of the computational methods of producing well instantaneous flow provided by Embodiment 2 of the present invention；
Fig. 5 is the image that flowing bottomhole pressure (FBHP) provided by Embodiment 2 of the present invention changes over time；
Fig. 6 be it is provided by Embodiment 2 of the present invention first derivative is sought to pressure in Fig. 4 after obtained pressure single order lead
The image that number changes over time.
Specific implementation mode
Technical solution to further illustrate the present invention below with reference to the accompanying drawings and specific embodiments.
Embodiment one
Fig. 1 is a kind of flow chart of the computational methods for producing well instantaneous flow that the embodiment of the present invention one provides, such as Fig. 1 institutes
Show, this method comprises the following steps：
S1, original pressure data is obtained.
Illustratively, original flowing bottomhole pressure (FBHP) data can be obtained in this step.
S2, original pressure data is pre-processed, obtains analyzable pressure data.
Illustratively, the original pressure data is pre-processed, it may include：Remove the noise in original pressure data
And singular point, amount of compressed data.It is under the premise of keeping the integrality of reset pressure information, obtain light to carry out pretreatment
Sliding analyzable pressure data.
S3, the first derivative for calculating pressure data, and determine the amplitude of pressure first derivative.
Specifically, the first derivative of pressure can be found out by following formula：
Wherein, pderivativeFor the first derivative of pressure, p1And p2For two adjacent pressure, t1For p1The corresponding time,
t2For p2The corresponding time.
S4, flowing period is divided according to the amplitude of pressure first derivative, wherein the amplitude of each pressure first derivative occurs
Time correspond to a flow changed time.
This step is illustrated with specific data instance below.Fig. 2 be the embodiment of the present invention one provide pressure at any time
Between the image that changes, Fig. 3 be the embodiment of the present invention one provide obtained pressure one after first derivative is sought to pressure in Fig. 2
The image that order derivative changes over time.Comparison diagram 2 and Fig. 3 can be seen that at the time of transient change occurs for pressure, and the one of pressure
Order derivative will appear an amplitude.When the pressure drops, the amplitude of corresponding pressure first derivative is negative value；When pressure rises
When, the amplitude of corresponding pressure first derivative is positive value.On this basis, different flowing periods is divided, that is, marks off difference
Production time section.
S5, calculate pressure first derivative amplitude and flow volume change values average proportions coefficient.
Specifically, average proportions coefficient can be found out by following formula：
Wherein, b is average proportions coefficient, and Q is total output, q1For initial flow, Δ tiIt is q for instantaneous flowiWhen it is corresponding
Production time interval, AjFor the amplitude of pressure first derivative, n is the number of different instantaneous flows.
S6, producing well instantaneous flow is calculated.
Specifically, producing well instantaneous flow can be found out by following formula：
Wherein, qiFor producing well instantaneous flow, q1For initial flow, AjFor the amplitude of pressure first derivative, b is average specific
Example coefficient.
The operation principle involved in the embodiment of the present invention is described in detail below：
The embodiment of the present invention find for the first time pressure single order inverse amplitude and flow variation it is proportional, it was demonstrated that such as
Under.According to the principle of stacking of oil reservoir, flowing bottomhole pressure (FBHP) p (t) is equal to the pleat of flow q (τ) and reservoir system Instantaneous Phase function g (t)
Product indicates as follows：
In formula：p0For original formation pressure, q (τ) is the sandface flow rate of producing well, and p (t) is the flowing bottomhole pressure (FBHP) of producing well, g
(t) it is reservoir system instantaneous function.Wherein, the value of g (t) decays rapidly at any time.In sandface flow rate changed wink
Between, flowing bottomhole pressure (FBHP) also changes therewith.Δ t, instantaneous pressure change, elta p in the changed very little time interval of flow
Mainly caused by changes in flow rate Δ q：
When time interval Δ t is identical,It is a definite value, therefore, instantaneous pressure change, elta p and flow
Changes delta q is directly proportional：
Δp∝Δq (3)
Comparison diagram 2 and Fig. 3 can be seen that at the time of transient change occurs for pressure, and the first derivative of pressure will appear one
A amplitude.When the pressure drops, the amplitude of corresponding pressure first derivative is negative value；When pressure rises, corresponding pressure one
The amplitude of order derivative is positive value.Instantaneous pressure variation is bigger, and pressure derivative is bigger, the amplitude A of pressure single order inverse and instantaneous pressure
The size of power changes delta p is directly proportional：
A∝Δp (4)
In conjunction with formula (3) and formula (4), it is directly proportional to changes in flow rate Δ q that first derivative amplitude A can be obtained：
A∝Δq (5)
Assuming that proportionality coefficient is b, then have：
A=Δs q × b (6)
Assuming that flow initial value is q1, production time section is Δ t1, changes in flow rate q2, production time section is Δ t2, according to
Formula (6), when flow is from q1Variation is q2When, the amplitude A of pressure first derivative1It can be expressed as：
A1=(q2-q1)×b (7)
Then flow q2Q can be used1It is indicated with b：
Similarly, when flow is from q2Variation is q3When, the amplitude of pressure first derivative is A2It can be expressed as：
A2=(q3-q2)×b (9)
Then flow q3Q can be used2It is indicated with b：
In conjunction with formula (8), formula (10) can be expressed as：
Similarly, flow qiQ can be used1It is indicated with b：
Assuming that flow is respectively q1, q2... qnWhen, the production time is respectively Δ t1, Δ t2... Δ tn, total flow Q can be with
It is expressed as：
Formula (12) is substituted into formula (13), formula (13) becomes：
According to formula (14), average proportions coefficient b can be in the hope of：
The average proportions coefficient b acquired, which is substituted into formula (12), can acquire producing well any time uninterrupted.
The computational methods for a kind of producing well instantaneous flow that the embodiment of the present invention one provides, by calculate pressure data one
Order derivative determines the amplitude of pressure first derivative, and the time that the amplitude of each pressure first derivative occurs corresponds to a flow
The changed time so according to the amplitude of pressure first derivative can be divided flowing period, and calculates pressure first derivative
The average proportions coefficient of amplitude and flow volume change values, finally calculates producing well instantaneous flow.By using above-mentioned technical proposal,
Without establishing reservoir model, the error brought using reservoir model is avoided；This method can be used for appointing for oil reservoir flowing simultaneously
Meaning period overcomes the drawbacks of production index method can be used only in oil reservoir stable state or quasi-stable state；This method has a wide range of application, and is counted
The instantaneous flow result precision of calculating is high.
Embodiment two
In the present embodiment, the specific implementation mode of the embodiment of the present invention is illustrated with specific data instance.This reality
The production total time for applying producing well in example is 192 hours, and total oil production is 987.87 cubic metres, initial flow 0.Specifically,
Fig. 4 is a kind of flow chart of the computational methods of producing well instantaneous flow provided by Embodiment 2 of the present invention, as shown in figure 4, the party
Method includes the following steps：
S10, original flowing bottomhole pressure (FBHP) data are obtained；
S20, by removing noise and the pretreatments such as singular point and amount of compressed data in original flowing bottomhole pressure (FBHP) data, obtain
To analyzable flowing bottomhole pressure (FBHP) data.
Fig. 5 is the image that flowing bottomhole pressure (FBHP) provided by Embodiment 2 of the present invention changes over time, as shown, by pretreatment
The curve smoother that flowing bottomhole pressure (FBHP) afterwards changes over time is conducive to subsequent calculating and analyzes.
S30, the first derivative for calculating flowing bottomhole pressure (FBHP) data, and determine the amplitude of pressure first derivative.
Fig. 6 be it is provided by Embodiment 2 of the present invention first derivative is sought to pressure in Fig. 5 after obtained pressure single order lead
The image that number changes over time.
S40, flowing period is divided according to the amplitude of pressure first derivative.
Comparison diagram 5 and Fig. 6 are it is found that time one flow of correspondence that the amplitude of each pressure first derivative occurs changes
Time, the correspondence between changes in flow rate time and the amplitude of pressure first derivative is had recorded in table 1, wherein each when
Between be spaced a corresponding flowing period.
Table 1
S50, calculate pressure first derivative amplitude and flow volume change values average proportions coefficient.
Known total flow Q=987.87m3, initial flow q1=0, it can be acquired using formula (15)：
S60, producing well instantaneous flow is calculated.
Utilize formulaProducing well any time instantaneous flow, result of calculation can be calculated
As shown in table 2.
Table 2
| Time (hour) | Time interval (hour) | Real traffic (m3) | Calculate flow (m3) | Error (%) |
| 0 | 0 | 0 | 0 | |
| 12.024 | 12.024 | 47.696 | 47.57 | -0.26 |
| 14.424 | 2.4 | 79.492 | 79.348 | -0.18 |
| 16.824 | 2.4 | 158.98 | 158.85 | -0.08 |
| 48.024 | 31.2 | 0 | 0 | 0 |
| 60.024 | 12 | 79.494 | 79.11 | -0.48 |
| 60.264 | 0.24 | 95.393 | 95.356 | -0.04 |
| 60.504 | 0.24 | 127.19 | 127.38 | 0.15 |
| 62.424 | 1.92 | 158.98 | 159.25 | 0.17 |
| 120.024 | 57.6 | 95.368 | 95.609 | 0.25 |
| 122.424 | 2.4 | 0 | 0 | 0 |
| 124.824 | 2.4 | 158.99 | 158.75 | -0.15 |
| 180.024 | 55.2 | 0 | 0 | 0 |
As can be seen from Table 2, the calculating flow results obtained using computational methods provided in an embodiment of the present invention with
Real traffic is smaller compared to error, and worst error is within 0.5%.
The foregoing is merely the preferred embodiment of the present invention, are not intended to restrict the invention, for those skilled in the art
For, the present invention can have various modifications and changes.Each embodiment of the present invention can phase on the basis of not violating logic
Mutually combination.All any modification, equivalent replacement, improvement and so within spirit and principles of the present invention should be included in this
Within the protection domain of invention.
Claims (3)
1. a kind of computational methods of producing well instantaneous flow, which is characterized in that including step：
S1, original pressure data is obtained；
S2, the original pressure data is pre-processed, obtains analyzable pressure data；
S3, the first derivative for calculating the pressure data, and determine the amplitude of pressure first derivative；
S4, flowing period is divided according to the amplitude of the pressure first derivative, wherein the amplitude of each pressure first derivative occurs
Time correspond to a flow changed time；
S5, calculate the pressure first derivative amplitude and flow volume change values average proportions coefficient；
S6, producing well instantaneous flow is calculated；
Wherein, it includes removing noise in the original pressure data and unusual to carry out pretreatment to the original pressure data
Point and amount of compressed data；
The first derivative of the pressure is found out by following formula：
Wherein, pderivativeFor the first derivative of pressure, p1And p2For two adjacent pressure, t1For p1Corresponding time, t2For p2
The corresponding time.
2. according to the method described in claim 1, it is characterized in that, the average proportions coefficient is found out by following formula：
Wherein, b is average proportions coefficient, and Q is total output, q1For initial flow, Δ tiIt is q for instantaneous flowiWhen corresponding life
Produce time interval, AjFor the amplitude of pressure first derivative, n is the number of different instantaneous flows.
3. according to the method described in claim 1, it is characterized in that, the producing well instantaneous flow is found out by following formula：
Wherein, qiFor producing well instantaneous flow, q1For initial flow, AjFor the amplitude of pressure first derivative, b is average proportions system
Number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510323955.5A CN105183929B (en) | 2015-06-12 | 2015-06-12 | A kind of computational methods of producing well instantaneous flow |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510323955.5A CN105183929B (en) | 2015-06-12 | 2015-06-12 | A kind of computational methods of producing well instantaneous flow |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105183929A CN105183929A (en) | 2015-12-23 |
| CN105183929B true CN105183929B (en) | 2018-08-21 |
Family
ID=54906009
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510323955.5A Active CN105183929B (en) | 2015-06-12 | 2015-06-12 | A kind of computational methods of producing well instantaneous flow |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105183929B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202381071U (en) * | 2011-11-25 | 2012-08-15 | 中国石油天然气股份有限公司 | Wellhead check valve type digital wireless transmission flowmeter of low-production oil producing well |
| CN103362486A (en) * | 2012-04-05 | 2013-10-23 | 中国石油天然气股份有限公司 | Digital underground layering water injection testing and adjusting process pipe column and operation method thereof |
| WO2014031223A1 (en) * | 2012-08-23 | 2014-02-27 | Exxonmobil Upstream Research Company | Systems and methods for re-completing multi-zone wells |
| CN104594860A (en) * | 2014-11-28 | 2015-05-06 | 航天科工惯性技术有限公司 | Flow monitoring controller for water injecting pipeline |
-
2015
- 2015-06-12 CN CN201510323955.5A patent/CN105183929B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202381071U (en) * | 2011-11-25 | 2012-08-15 | 中国石油天然气股份有限公司 | Wellhead check valve type digital wireless transmission flowmeter of low-production oil producing well |
| CN103362486A (en) * | 2012-04-05 | 2013-10-23 | 中国石油天然气股份有限公司 | Digital underground layering water injection testing and adjusting process pipe column and operation method thereof |
| WO2014031223A1 (en) * | 2012-08-23 | 2014-02-27 | Exxonmobil Upstream Research Company | Systems and methods for re-completing multi-zone wells |
| CN104594860A (en) * | 2014-11-28 | 2015-05-06 | 航天科工惯性技术有限公司 | Flow monitoring controller for water injecting pipeline |
Non-Patent Citations (1)
| Title |
|---|
| Unknown Rate History Calculation from Down-hole Transient Pressure Data UsingWavelet Transform;FuyongWang 等;《Transport in Porous Media》;20121218;第96卷;第547-566页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105183929A (en) | 2015-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105464652A (en) | Dynamic reserve calculation method and system for fracture and cave type carbonatite oil reservoir fracture and cave unit | |
| CA2665116A1 (en) | System and method for performing oilfield simulation operations | |
| CN103335682B (en) | Measuring method for gas flow of natural gas | |
| CN105900099A (en) | Optimizing flow control device properties for a liquid injection well using a coupled wellbore-reservoir model | |
| Haouche et al. | Virtual flow meter pilot: Based on data validation and reconciliation approach | |
| CN105183929B (en) | A kind of computational methods of producing well instantaneous flow | |
| CN106481315B (en) | Land sandstone oil reservoir individual well recoverable reserves quickly determines model and method for building up | |
| CN112554864A (en) | Method for calculating single-well control reserve of water-producing gas well | |
| CN103674541A (en) | Method for testing performance of differential pressure valve | |
| CN109063258A (en) | Satellitosis telemetry parameter storage method and system | |
| CN107340714A (en) | A kind of adaptive inverse control of nanometer of servo-drive system | |
| CN204027618U (en) | A kind of flow measurement device that is installed in parallel areal survey wide range ratio | |
| CN104237101B (en) | Bridge imitation and flow measurement avoidance type cross rock core permeability measurement method and device | |
| CN108562944A (en) | A kind of parameter of pore structure computational methods based on Electrical imaging porosity spectrum | |
| CN204902902U (en) | Low velocity of flow to high flow rate wide range clings to flowmeter than pitot | |
| CN110765660B (en) | Method for quickly determining movable reserves of low-permeability gas under different gas production speeds | |
| CN204461518U (en) | A kind of heating type turbine integral type flowmeter | |
| CN107358541A (en) | A kind of method for detecting abnormality of electric analog data | |
| CN103164530B (en) | Modeling based on measuring point | |
| CN202420763U (en) | Remote calibration intelligent pressure transducer system | |
| CN107330207B (en) | Sliding bearing flow parameter calculation method corrected by multi-factor coupling test | |
| CN105302169A (en) | Flow control method | |
| US20210073320A1 (en) | Video, audio, and historical trend data interpolation | |
| Rafiee et al. | Subsurface Back Allocation: Calculating Production and Injection Allocation by Layer in a Multilayered Waterflood Using a Combination of Machine Learning and Reservoir Physics | |
| JP2006190264A (en) | Method for simulating semiconductor integrated circuit |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |