CN103246778B - Curve fitting processing method - Google Patents
Curve fitting processing method Download PDFInfo
- Publication number
- CN103246778B CN103246778B CN201310177697.5A CN201310177697A CN103246778B CN 103246778 B CN103246778 B CN 103246778B CN 201310177697 A CN201310177697 A CN 201310177697A CN 103246778 B CN103246778 B CN 103246778B
- Authority
- CN
- China
- Prior art keywords
- curve
- actual measurement
- peak value
- peak
- derivative curve
- 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.)
- Expired - Fee Related
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000005259 measurement Methods 0.000 claims description 42
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims 1
- 238000007906 compression Methods 0.000 claims 1
- 238000013507 mapping Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 101000878457 Macrocallista nimbosa FMRFamide Proteins 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a curve fitting processing method. By means of the method, the peak position of an actually measured derivative curve can be obtained according to a coordinate mapping relationship by utilizing a mouse position on the basis of uniform oil deposit plate characteristics, and the corresponding CDe2S parameters can be obtained accordingly; and a theoretical curve and a corresponding theoretical derivative curve are calculated and drawn according to the obtained CDe2S parameters. The aim of fast curve fitting can be achieved by controlling cursor for once, and the working efficiency is greatly improved.
Description
Technical field
The present invention relates to computer technology and petroleum works field, more particularly to a kind of processing method of curve matching.
Background technology
Existing well test analysis method mainly includes：Straightway analysis method, curve-fitting method and non-linear regression method.
Wherein, the essence of existing curve-fitting method is to realize curve matching by trial and error procedure；That is, parameter is constantly changed,
After often changing once, fitting result is checked, until the profile of curve is approximate with radian；Again moving curve makes it overlap.The method
Defect be mainly：1）It is that method is gathered in examination, to explaining that engineer's level has high requirements；2）Need repetitious change ginseng
Number, and accordingly calculated, work efficiency is relatively low；3）When parameter area is larger（The scope 10 of relevant parameter-6～10200）More
The number of times for changing is very more, it is impossible to be rapidly performed by curve matching.
The content of the invention
It is an object of the invention to provide a kind of processing method of curve matching, improves work efficiency.
The purpose of the present invention is achieved through the following technical solutions：
A kind of processing method of curve matching, the method includes：
Measured curve and corresponding actual measurement derivative curve are plotted on log-log graph so that the measured curve and actual measurement
The initial part of derivative curve, and through origin and slope overlaps for the oblique line of predetermined value s；And the ending of the actual measurement derivative curve
Horizontal line of the part with numerical value for predetermined value n overlaps；
Cursor control instruction is received, cursor is controlled and moved according to the instruction, when cursor is moved to the actual measurement derivative curve
Peak value at when, calculate the corresponding zero dimension combination parameter C of the actual measurement derivative curve peak valueDe2S；
According to the C for calculatingDe2SParameter calculates and draws a theoretical curve and corresponding theoretical derivative curve, completes curve
Fitting.
As seen from the above technical solution provided by the invention, based on coordinate mapping relations, so as to utilize mouse position
The peak of actual measurement derivative curve is acquired, and corresponding C is calculated according to peakDe2SParameter；According to the C for obtainingDe2S
Parameter calculates and draws a theoretical curve and corresponding theoretical derivative curve；Reach control cursor and be once capable of achieving curve
The purpose of Fast Fitting, greatly improves the efficiency of work.Compared with the curve-fitting method that examination is gathered, this method is according to homogenizing
Oil reservoir plate feature, it is artificial to obtain derivative curve peak value, relevant parameter is calculated, reach the successful effect of once fitting.
Description of the drawings
In order to be illustrated more clearly that the technical scheme of the embodiment of the present invention, below will be to use needed for embodiment description
Accompanying drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the present invention, for this
For the those of ordinary skill in field, on the premise of not paying creative work, can be obtaining other according to these accompanying drawings
Accompanying drawing.
Fig. 1 is a kind of flow chart of the processing method of curve matching provided in an embodiment of the present invention；
Fig. 2 is the schematic diagram that boost line is drawn in log-log graph provided in an embodiment of the present invention；
Fig. 3 is the schematic diagram that measured curve and actual measurement derivative curve are drawn in log-log graph provided in an embodiment of the present invention；
Fig. 4 is the schematic diagram of curve fitting process provided in an embodiment of the present invention；
Fig. 5 is the schematic diagram completed after curve matching provided in an embodiment of the present invention.
Specific embodiment
With reference to the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Ground description, it is clear that described embodiment is only a part of embodiment of the invention, rather than the embodiment of whole.Based on this
Inventive embodiment, the every other enforcement that those of ordinary skill in the art are obtained under the premise of creative work is not made
Example, belongs to protection scope of the present invention.
Fig. 1 is a kind of flow chart of the processing method of curve matching provided in an embodiment of the present invention.As shown in figure 1, the party
Method mainly comprises the steps：
Step 11, measured curve and corresponding actual measurement derivative curve are plotted on log-log graph so that the actual measurement is bent
The initial part of line and actual measurement derivative curve, and through origin and slope overlaps for the oblique line of predetermined value s；And the actual measurement derivative is bent
Horizontal line of the ending of line with numerical value for predetermined value n overlaps.
As shown in Fig. 2 in log-log graph（For example, abscissa is the time, in the range from 10-2～104, vertical coordinate is pressure,
In the range from 10-2～102）Two boost lines of upper drafting：One slope is predetermined value s（For example, 1）Oblique line, another numerical value
For predetermined value n（For example, 0.5）Horizontal line.
The measured curve can be bottom pressure data and curves, and corresponding actual measurement derivative curve is the pressure data derivation
Curve afterwards.As shown in figure 3, by two curve plottings on log-log graph, and cause the measured curve and actual measurement derivative bent
The initial part of line, and through origin and slope overlaps for the oblique line of s；And the ending of the actual measurement derivative curve is with numerical value
The horizontal line of n overlaps；The position of above-mentioned two curves is demarcated by boost line, be easy to the curve matching work in later stage
Make.
Step 12, reception cursor control instruction, control and move cursor according to the instruction, lead when cursor is moved to the actual measurement
When at the peak value of number curve, the corresponding zero dimension combination parameter C of the actual measurement derivative curve peak value is calculatedDe2S。
Because the curve matching of the present embodiment is based on zero dimension combination parameter CDe2S, and derivative curve is relative to CDe2SParameter
It is more sensitive, it is therefore desirable to calculate the C of actual measurement derivative curveDe2SParameter；There is corresponding relation with the peak value of curve in the parameter,
After the peak value of curve is known, the corresponding C of acquisition can be calculatedDe2SParameter.In the same manner, when learn actual measurement derivative curve peak value it is corresponding
CDe2SAfter parameter, then the theoretical curve that is completely superposed with measured curve and corresponding actual measurement derivative curve and corresponding can be calculated
Theoretical derivative curve.
Below for the peak C of actual measurement derivative curveDe2SThe solution procedure of parameter is described in detail：
Exemplary, if measured curve is the bottom pressure data and curves on homogenizing infinity stratum, dimensionless shaft bottom pressure
Power expression formula is：
Wherein, J0(u) and J1U () is that 0 rank surveys class Bessel function with 1 rank（Bessel function）；Y0(u) and Y1U () is
0 rank and 1 rank theory class Bessel function；For dimensionless expression formula；For dimensionless
Time；For dimensionless well hole storage constant；S is skin factor；pwfFor bottom hole flowing pressure（Mpa）；piFor
Original formation pressure（Mpa）；H is Effective thickness of formation（m）；μ is fluid viscosity（Mpa.s）；φ is formation porosity；CtFor comprehensive
Close coefficient of compressibility（1/Mpa）.
Solve the dimensionless bottom pressure expression formula and obtain the corresponding time t of actual measurement Derivative peakD：
Wherein, Δ=[CDuJ0(u)-(1-CDu2S)J1(u)]2+[CDuY0(u)-(1-CDu2S)Y1(u)]2。
It is right（1）-（2）Carry out to be drawn after induction and conclusion：
Lg(CDe2S)=2.351187hD+4.73227×10-2；（3）
tD/CD=0.8968505Lg (CDe2S)+3.604613.（4）
Wherein, hDFor the peak value for surveying derivative curve after nondimensionalization；（4）Formula represents the expression of time after nondimensionalization
Formula.
By（3）Formula understands, as the h for learning actual measurement derivative curveDAfterwards, then C can be calculatedDe2SParameter.The embodiment of the present invention
Based on coordinate mapping relations, using mouse position the h of actual measurement derivative curve is calculatedD, and then calculate corresponding CDe2SParameter.Tool
Body method is as follows：First, cursor control instruction is received, cursor is controlled and moved according to the instruction（I.e. user clicks on cursor and protects
Click state movement light timestamp is held, computer is operated accordingly according to the instruction for receiving）, when cursor is moved to the actual measurement
When at the peak value of derivative curve, the coordinate at peak value is calculated using coordinate mapping relations, its formula is：
Wherein, y is the vertical coordinate that light is marked in screen, y0For the corresponding screen vertical coordinate of origin vertical coordinate of log-log graph, lyFor screen
Height.Peak value h is calculated further according to peak coordinateD, its formula is：Finally combine（3）Formula calculates CDe2SGinseng
Number.
On the other hand, for the ease of the result of Real Time Observation curve matching, in cursor actual measurement derivative curve peak value is moved to
Before, it is also possible to calculate theoretical curve and theoretical derivative curve in real time, as shown in Figure 4.So, user just can be according to theoretical bent
Whether line and theoretical derivative curve judge whether light target position is real with the fitting degree of measured curve and actual measurement derivative curve
Survey at the peak value of derivative curve.
Step 13, according to the C for calculatingDe2SParameter calculates and draws a theoretical curve and corresponding theoretical derivative curve,
Complete curve matching.
As shown in figure 5, when cursor is moved at the peak value of the actual measurement derivative curve, according to the C for calculatingDe2SParameter meter
Calculate and draw a theoretical curve and corresponding theoretical derivative curve；Now, four curves are completely superposed, and complete curve matching, can
To terminate to light target clicking operation, the purpose that mobile cursor once completes curve matching is realized by the way, reduce
The difficulty of work and time, improve the efficiency of work.
The embodiment of the present invention is based on coordinate mapping relations, so as to acquire the peak value of actual measurement derivative curve using mouse position
Position, and the corresponding C of acquisition is calculated with thisDe2SParameter；According to the C for obtainingDe2SParameter calculates and draws a theoretical curve and right
The theoretical derivative curve answered；The purpose that control cursor is once capable of achieving the Fast Fitting of curve is reached, has been greatly improved
The efficiency of work.
Through the above description of the embodiments, those skilled in the art can be understood that above-described embodiment can
To be realized by software, it is also possible to realize by the mode of software plus necessary general hardware platform.Based on such understanding,
The technical scheme of above-described embodiment can be embodied in the form of software product, the software product can be stored in one it is non-easily
The property lost storage medium（Can be CD-ROM, USB flash disk, portable hard drive etc.）In, including some instructions use is so that a computer sets
It is standby（Can be personal computer, server, or network equipment etc.）Perform the method described in each embodiment of the invention.
The above, the only present invention preferably specific embodiment, but protection scope of the present invention is not limited thereto,
Any those familiar with the art in the technical scope of present disclosure, the change or replacement that can be readily occurred in,
All should be included within the scope of the present invention.Therefore, protection scope of the present invention should be with the protection model of claims
Enclose and be defined.
Claims (1)
1. a kind of processing method of curve matching, it is characterised in that the method includes：
Measured curve and corresponding actual measurement derivative curve are plotted on log-log graph so that the measured curve and actual measurement derivative
The initial part of curve, and through origin and slope overlaps for the oblique line of predetermined value s；And the ending of the actual measurement derivative curve
Overlap with horizontal line of the numerical value for predetermined value n；
Cursor control instruction is received, cursor is controlled and moved according to the instruction, when cursor is moved to the peak of the actual measurement derivative curve
When at value, the corresponding zero dimension combination parameter C of the actual measurement derivative curve peak value is calculatedDe2S；
According to the C for calculatingDe2SParameter calculates and draws a theoretical curve and corresponding theoretical derivative curve, completes curve plan
Close；
The measured curve includes：
The bottom pressure data and curves on homogenizing infinity stratum, its dimensionless bottom pressure expression formula is：
Wherein, J0(u) and J1U () is 0 rank and 1 rank the primal Bessel function；Y0(u) and Y1U () is 0 rank and 1 rank Equations of The Second Kind
Bessel functions；tDFor nondimensional time；CDFor dimensionless well hole storage constant；S is skin factor；pwfFor flowing bottom hole pressure
Power；piFor original formation pressure；H is Effective thickness of formation；μ is fluid viscosity；φ is formation porosity；CtFor comprehensive compression system
Number；
The corresponding time t of actual measurement Derivative peakD：
Wherein, Δ=[CDuJ0(u)-(1-CDu2S)J1(u)]2+[CDuY0(u)-(1-CDu2S)Y1(u)]2；
It is described when cursor is moved at the peak value of the actual measurement derivative curve, calculate the actual measurement derivative curve peak value it is corresponding without because
Secondary combination parameter CDe2SIncluding：
When cursor is moved at the peak value of the actual measurement derivative curve, actual measurement described in the position calculation in screen is marked on according to light and is led
The peak coordinate of number curve, and the actual measurement derivative curve peak value and correspondence are calculated according to the peak coordinate of the actual measurement derivative curve
CDe2SParameter；
The peak coordinate of derivative curve is surveyed described in the position calculation being marked on according to light in screen, and according to the actual measurement derivative
The peak coordinate of curve calculates the actual measurement derivative curve peak value and corresponding CDe2SParameter includes：
The peak value vertical coordinate that derivative curve is surveyed described in the position calculation in screen is marked on according to light, its formula is：Wherein, y is the vertical coordinate that light is marked in screen, y0Origin vertical coordinate for log-log graph is corresponding
Screen vertical coordinate, lyFor the height of screen；
Peak value h is calculated according to peak value vertical coordinateD, its formula is：
By peak value hDCalculate corresponding CDe2SParameter, its formula is：lg(CDe2S)=2.351187hD+4.73227×10-2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310177697.5A CN103246778B (en) | 2013-05-14 | 2013-05-14 | Curve fitting processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310177697.5A CN103246778B (en) | 2013-05-14 | 2013-05-14 | Curve fitting processing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103246778A CN103246778A (en) | 2013-08-14 |
| CN103246778B true CN103246778B (en) | 2017-04-19 |
Family
ID=48926296
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310177697.5A Expired - Fee Related CN103246778B (en) | 2013-05-14 | 2013-05-14 | Curve fitting processing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103246778B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104392068B (en) * | 2014-12-10 | 2017-09-29 | 中国科学院合肥物质科学研究院 | A kind of processing method of the oil reservoir of closed boundary containing straight line curve matching |
| CN112731161B (en) * | 2021-02-08 | 2021-10-26 | 中南大学 | Nonlinear data feature extraction and classification prediction method based on small amount of data mixed insertion |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4607524A (en) * | 1985-04-09 | 1986-08-26 | Scientific Software-Intercomp, Inc. | Method for obtaining a dimensionless representation of well pressure data without the use of type-curves |
| CN101216564A (en) * | 2008-01-16 | 2008-07-09 | 中国科学技术大学 | Early well testing analytical method based on small signal extraction technology |
| CN101446196A (en) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | Well test analysis method and device of treble medium oil pool branch horizontal well |
| CN102305062A (en) * | 2011-07-25 | 2012-01-04 | 中国科学技术大学 | Method for interpreting saturation degree of multiphase liquid |
-
2013
- 2013-05-14 CN CN201310177697.5A patent/CN103246778B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4607524A (en) * | 1985-04-09 | 1986-08-26 | Scientific Software-Intercomp, Inc. | Method for obtaining a dimensionless representation of well pressure data without the use of type-curves |
| CN101216564A (en) * | 2008-01-16 | 2008-07-09 | 中国科学技术大学 | Early well testing analytical method based on small signal extraction technology |
| CN101446196A (en) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | Well test analysis method and device of treble medium oil pool branch horizontal well |
| CN102305062A (en) * | 2011-07-25 | 2012-01-04 | 中国科学技术大学 | Method for interpreting saturation degree of multiphase liquid |
Non-Patent Citations (3)
| Title |
|---|
| 特征点拟合法在非牛顿流试井解释中的应用;梁光跃 等;《西南石油大学学报（自然科学版）》;20111231;第33卷(第6期);摘要，引言，第3-4节，附图1-3,6-7 * |
| 考虑井筒热效应的气体井底压力计算;卢德唐 等;《中国科学技术大学学报》;20100331;第40卷(第3期);全文 * |
| 考虑热效应的CO2气井试井压力模型和计算方法;侯彬彬;《中国科学技术大学学报》;20110331;第41卷(第3期);参见第1,3部分 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103246778A (en) | 2013-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106599449B (en) | Well test interpretation method for karst cave volume calculation | |
| CN103065356B (en) | A kind of method for building up of offshore oil facility three-dimensional model | |
| CN106886046B (en) | Determine that fractured-cavernous reservoir is not gone into operation the methods of block available reserves | |
| CN105631935A (en) | BIM-based three-dimensional geological automatic modeling method | |
| CN104236893B (en) | Performance parameter test system and performance parameter test method of hydraulic damper | |
| CN103699634A (en) | Database creation method and database creation system for one-stop pipeline detection | |
| CN103866769A (en) | Dam grouting construction analysis and control method based on 3D (three-dimensional) geological model and real-time monitoring | |
| CN108875273A (en) | Base pit engineering Top-Down Design and operation management system based on BIM | |
| CN104181548A (en) | Method for monitoring water and soil conservation condition based on three-dimensional laser scanning technology | |
| CN104405365A (en) | Pumping unit indicator diagram liquid production capacity measurement technology | |
| WO2016141685A1 (en) | Telecommunications engineering drawing generation method and drawing design client | |
| CN104764435A (en) | Intelligent surveying and mapping system and method | |
| CN103246778B (en) | Curve fitting processing method | |
| CN108845663A (en) | Design and construction system | |
| CN105067050A (en) | Method for improving accuracy of natural gas measurement | |
| CN103047975B (en) | Method of accurate and rapid measurement lofting system | |
| CN108665172A (en) | A kind of interactive building engineering construction quality control method | |
| CN103156646B (en) | Method and device of dual-screen ultraphonic medical treatment measurement | |
| CN205138544U (en) | Intelligence is measured and is managed integration system based on internet | |
| CN105045958A (en) | Implementation system and method of GPS (Global Positioning System) elevation fitting on the basis of BP (Back Propagation) neural network | |
| CN104408218A (en) | Visual examination method of oil gas pipeline design achievement | |
| CN104392068B (en) | A kind of processing method of the oil reservoir of closed boundary containing straight line curve matching | |
| CN108664678B (en) | Yield prediction method | |
| CN105302979A (en) | Modeling method and system of valve groups in two-phase fluid network model | |
| CN106021571A (en) | Method and device for generating oil production engineering scheme |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170419 |