CN105626039B - A kind of method of production logging correlative flow data prediction - Google Patents
A kind of method of production logging correlative flow data prediction Download PDFInfo
- Publication number
- CN105626039B CN105626039B CN201511032493.8A CN201511032493A CN105626039B CN 105626039 B CN105626039 B CN 105626039B CN 201511032493 A CN201511032493 A CN 201511032493A CN 105626039 B CN105626039 B CN 105626039B
- Authority
- CN
- China
- Prior art keywords
- depth
- curve
- gamma curve
- time
- value
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 12
- 239000012530 fluid Substances 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000009795 derivation Methods 0.000 claims abstract description 5
- 230000009885 systemic effect Effects 0.000 claims abstract description 3
- 238000005259 measurement Methods 0.000 claims description 6
- 239000000523 sample Substances 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 239000000700 radioactive tracer Substances 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000011229 interlayer Substances 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000007921 spray Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011437 continuous method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Image Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of methods of production logging correlative flow data prediction, including to correlative flow data setting depth curve DEPTH, gamma curve GR and time graph TIME；FIR filtering carries out second order derivation to filtered gamma curve GR；When second dervative is less than 0, as gamma curve GR local peaking；Local peaking and systemic presupposition threshold value comparison, less than filtering out for threshold value, greater than the reservation of threshold value；Pseudo- depth, true depth and the time value that gamma curve GR is locally corresponded to depth are written in table together, calculate fluid velocity；The fluid velocity point value being calculated is depicted as curve.How its core by algorithm is accurately identifying GR curve local maxima, this processing simple and convenient operation flow, automatic identification accuracy is high, and treatment effect is preferable, has certain application value.
Description
Technical field
It is a kind of in the explanation of production logging injection profile the invention belongs to oil development and engineering field, correlative flow
The method of data prediction.
Background technique
With the development of production logging technology, correlative flow technology is more and more applied in the explanation of water injection profile.
Correlative flow be suitable in, Low-flow Wells, generally water injection well use, in producing well use since variations in flow patterns can make resolution ratio
Decline.In the well that turbine flowmeter cannot be used to measure, correlation flowmeters are generally used.Correlation flowmeters determine in the wellbore
The method of flow is tachometric method.Measurement installs instruments according to hole condition, and instrument, which is parked in spray between two perforation layers into pit shaft, to be shown
Then track agent measures the times required for tracer is transmitted in two points, refer generally between two detectors or injector is to detecting
Time between device thereby determines that the view flow velocity of each interpretation layer.For producing well, spray site should be close to the bottom of perforation layer interlayer
Portion then selects the top of interlayer for injection well.According to measurement method, tachometric method includes two methods, and one is static measurements
Method；Another kind is back tracking method.
It is accurate to determine injection tracer since the time of injection tracer is variation for gamma-ray detector
The time that agent reaches probe is more difficult, therefore generally uses continuous method for tracing.In this case, in an interlayer,
At least measured three times.Tracer bridge plug is also flowing while due to measurement, it is therefore necessary to which guarantee has sufficiently high flow velocity
Measure complete tracer slug.If it find that the displacement of first time is larger, then it should accelerate to test the speed, on the contrary then reduction is tested the speed.Flow velocity
Calculation method be
In formula, Δ H measures the distance (depth difference of peak value) of tracer slug displacement twice；
Δ t is the time needed for slug displacement；
What correlative flow pretreatment was done is exactly to find the local maximum of GR curve, that is, wave crest, then according to wave
The pseudo- depth at peak finds the corresponding true depth of Depth curve and Time curve corresponding time, then can calculate two o'clock
Flow velocity.
Summary of the invention
The purpose of the present invention is for currently without preferably by automatic program identification GR curve local maxima, craft
Interaction is picked up there is a problem of time-consuming and laborious again, and finding one kind by program and can set corresponding threshold value automatic identification part
The correlative flow preprocess method of maximum value and then calculating speed curve.The present invention is local most by automatic program identification GR curve
Big value obtains corresponding true depth, time according to the corresponding pseudo- depth of local maximum, so that flow velocity is calculated, operation of the present invention
Simplicity, automatic identification accuracy is high, has certain engineering application value.
The purpose of the present invention is what is realized by following technical proposals.
A kind of method of production logging correlative flow data prediction, comprising the following steps:
Step 1, to the production logging correlative flow data using single probe back tracking method measurement, setting depth curve DEPTH,
Gamma curve GR and time graph TIME；
Step 2, FIR filtering is carried out to gamma curve GR, removes flash removed and interference；
Step 3, second order derivation is carried out to filtered gamma curve GR:
GR "=GR [i+1]+GR [i-1] -2*GR [i]
In formula, i indicates current gamma curve data point, and i-1 indicates previous data point, and i+1 is the latter data point；
GR " is the second dervative of gamma curve；
Step 4, when second dervative is less than 0, the as position of gamma curve GR local maximum, that is, GR innings of gamma curve
Portion's peak value；
Step 5, gamma curve GR local peaking and systemic presupposition threshold value comparison are then not considered as gal less than threshold value
Horse curve GR local maximum, filters out；It is then required gamma curve GR local maximum greater than threshold value；
Step 6, required gamma curve GR local maximum step 5 obtained corresponds to the pseudo- depth, very deep of depth
Degree and time value are written in table together, calculate fluid velocity；
Step 7, the fluid velocity point value being calculated is depicted as curve.
Further, in step 2, there is limit for length's bandpass filter to carry out FIR filtering using 31 gamma curve GR, set GR
Discretization plot against time uses 4us, frequency filtering 1-8khz.
Further, in step 5, when program starts, gamma curve GR all values is read, therefrom count gamma curve GR most
/ 10th be worth greatly are used as default threshold.
Further, in step 6, pseudo- depth, true depth and the time value that all gamma curve GR are locally corresponded to depth are write
Enter table to fall into a trap fluid operator speed, be obtained by following formula:
V=60* (Depth [i-1]-Depth [i])/(Time [i-1]-Time [i])
Wherein, V indicates that the mean flow rate of Depth [i-1] to depth Depth [i], Depth [i-1] are right for (i-1)-th point
The true depth answered, Depth [i] are true depth corresponding to i-th point, and Time [i-1] is time corresponding to (i-1)-th point, Time
[i] is the time corresponding to i-th point.
Further, in step 7, the fluid velocity point value being calculated is depicted as curve, is accomplished in the following manner:
The fluid velocity of Depth1, Depth2 are respectively Vdepth1, Vdepth2, then fluid velocity VdepthAre as follows:
Vdepth=Vdepth2+(Depth-Depth2)*(Vdepth1-Vdepth2)/(Depth1-Depth2)。
Compared with the existing technology, the beneficial effects of the present invention are:
The present invention provides a kind of methods of production logging correlative flow data prediction.How its core is passing through algorithm
Accurately identifying GR curve local maxima, the present invention is using first FIR filtering is carried out to GR curve, then to the secondary derivation of curve,
By second derivative less than 0 be local maximum, local maximum is compared with given threshold, what is be less than filters out, then
According to the corresponding true depth of local maximum position acquisition, time, thus calculating speed curve.This process flow has innovation
Property.Execute-in-place is easy, and treatment effect is preferable, has certain application value.
Detailed description of the invention
Fig. 1 is processing flow schematic diagram.
Fig. 2 is parameter setting dialog box.
Fig. 3 is recognition effect figure, and red indicates the local maximum of identification.
Fig. 4 is the corresponding curve data of local maximum stored in a tabular form.
Fig. 5 is calculating speed curve dialog box.
Fig. 6 is the rate curve effect picture calculated.
Specific embodiment
Below with reference to particularly relevant flow rate log data, specific embodiment of the invention is illustrated.
As shown in Figure 1, a kind of method of production logging correlative flow data prediction of the present invention, includes the following steps:
Step 1, relevant parameter is set, as shown in Fig. 2, to the production logging correlative flow using single probe back tracking method measurement
The depth curve DEPTH, gamma curve GR, time graph TIME of input is arranged in data, and setting range coefficient 1 is the wave crest taken,
It is 603 that program calculates an input threshold size automatically.Setting direction is above to propose well logging, then program only counts calculating and above mentions direction
Maximum value.Depth is oriented to depth, when not instead of wave crest or trough that the amplitude of referring to takes, when wave crest is multiplied by range coefficient,
What value at this moment took is top half or lower half portion.
Step 2, it clicks and determines, program automatically processes.The process of program inter-process flow chart as shown in Figure 1, first
Filter GR using FIR, FIR filtering has limit for length's bandpass filter using 31, and the setting GR time uses 4us, filtering frequency
Rate is 1-8khz.
Step 3, to the GR discretization curve of filtering, second order derivation is carried out:
GR "=GR [i+1]+GR [i-1] -2*GR [i]
In formula, i indicates current gamma curve data point, and i-1 indicates previous data point, and i+1 is the latter data point；
GR " is the second dervative of gamma curve.
It step 4, is exactly the position of local maximum, that is, part GR peak value when second dervative is less than 0.
Step 5, local GR peak value and threshold value comparison are then not considered as local maximum less than threshold value, filtered out；
It is then required gamma curve GR local maximum greater than threshold value；
Wherein, when program starts, gamma curve GR all values is read, therefrom count gamma curve GR maximum value very
One of be used as default threshold.
Step 6, GR value, pseudo- depth, true depth, time value that all local maximums correspond to depth are written in table,
Calculate fluid velocity:
V=60* (Depth [i-1]-Depth [i])/(Time [i-1]-Time [i])
Wherein, V indicates that the mean flow rate of Depth [i-1] to depth Depth [i], Depth [i-1] are right for (i-1)-th point
The true depth answered, Depth [i] are true depth corresponding to i-th point, and Time [i-1] is time corresponding to (i-1)-th point, Time
[i] is the time corresponding to i-th point.
Step 7, the fluid velocity point value being calculated is depicted as curve, is accomplished in the following manner:
The fluid velocity of Depth1, Depth2 are respectively Vdepth1, Vdepth2, then fluid velocity VdepthAre as follows:
Vdepth=Vdepth2+(Depth-Depth2)*(Vdepth1-Vdepth2)/(Depth1-Depth2)。
The table of generation such as Fig. 4, the effect picture of peak value automatic Picking such as Fig. 3, it is Fig. 4 table that wave crest, which corresponds to horizontal line part,
Visualization display, corresponding is the position of GR peak value, and right side graph is track depth curve, and corresponding is true depth value,
What is calculated is the peak value for above mentioning direction, and corresponding track depth curve is the up increased part of depth.It can be seen that identification peak
It is very accurate to be worth position, does not omit.
According to the table of preservation, calculating speed curve and resultant curve；Fig. 5 is corresponding setting dialog box, and selection generates
TVAU01 table, click " transformation curve ", then generate corresponding curve, method is the speed table to generation, is carried out linear
Difference, formation speed curve.The speed effect figure of generation is as shown in Figure 6.
It is provided for the embodiments of the invention embodiment above to be described in detail, specific case used herein
The principle and embodiment of the embodiment of the present invention are expounded, the explanation of above embodiments is only applicable to help to understand this
The principle of inventive embodiments；At the same time, for those skilled in the art, according to an embodiment of the present invention, in specific embodiment party
There will be changes in formula and application range, in conclusion the contents of this specification are not to be construed as limiting the invention.
Claims (5)
1. a kind of method of production logging correlative flow data prediction, which comprises the following steps:
Step 1, to the production logging correlative flow data using single probe back tracking method measurement, depth curve DEPTH, gamma are set
Curve GR and time graph TIME；
Step 2, FIR filtering is carried out to gamma curve GR, removes flash removed and interference；
Step 3, second order derivation is carried out to filtered gamma curve GR:
GR "=GR [i+1]+GR [i-1] -2*GR [i]
In formula, i indicates current gamma curve data point, and i-1 indicates previous data point, and i+1 is the latter data point；GR " is
The second dervative of gamma curve；
Step 4, when second dervative is less than 0, the as position of gamma curve GR local maximum, that is, gamma curve GR local peaks
Value；
Step 5, by gamma curve GR local peaking and systemic presupposition threshold value comparison, less than threshold value, then it is not considered as gamma song
Line GR local maximum, filters out；It is then required gamma curve GR local maximum greater than threshold value；
Step 6, required gamma curve GR local maximum step 5 obtained correspond to the pseudo- depth of depth, true depth and
Time value is written in table together, calculates fluid velocity；
Step 7, the fluid velocity point value being calculated is depicted as curve.
2. the method according to claim 1, wherein in step 2, to gamma curve GR using 31 have limit for length with
Bandpass filter carries out FIR filtering, and setting GR discretization plot against time uses 4us, frequency filtering 1-8khz.
3. the method according to claim 1, wherein when program starts, reading gamma curve GR institute in step 5
There is value, therefrom count gamma curve GR maximum value 1/10th are used as default threshold.
4. the method according to claim 1, wherein all gamma curve GR are locally corresponded to depth in step 6
Pseudo- depth, true depth and time value write-in table fall into a trap fluid operator speed, obtained by following formula:
V=60* (Depth [i-1]-Depth [i])/(Time [i-1]-Time [i])
Wherein, V indicates that the mean flow rate of Depth [i-1] to depth Depth [i], Depth [i-1] are corresponding to (i-1)-th point
True depth, Depth [i] are true depth corresponding to i-th point, and Time [i-1] is the time corresponding to (i-1)-th point, and Time [i] is
Time corresponding to i-th point.
5. the method according to claim 1, wherein the fluid velocity point value being calculated is drawn in step 7
At curve, it is accomplished in the following manner:
The fluid velocity of Depth1, Depth2 are respectively Vdepth1, Vdepth2, then fluid velocity VdepthAre as follows:
Vdepth=Vdepth2+(Depth-Depth2)*(Vdepth1-Vdepth2)/(Depth1-Depth2)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511032493.8A CN105626039B (en) | 2015-12-31 | 2015-12-31 | A kind of method of production logging correlative flow data prediction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511032493.8A CN105626039B (en) | 2015-12-31 | 2015-12-31 | A kind of method of production logging correlative flow data prediction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105626039A CN105626039A (en) | 2016-06-01 |
| CN105626039B true CN105626039B (en) | 2019-01-18 |
Family
ID=56041330
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511032493.8A Active CN105626039B (en) | 2015-12-31 | 2015-12-31 | A kind of method of production logging correlative flow data prediction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105626039B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU365460A1 (en) * | 1970-12-11 | 1973-01-08 | Опытно конструкторское бюро геофизического приборостроени треста Бурмашремоит | METHOD OF MEASUREMENT OF OIL DEBIT IN WELLS WITH STAGENT WATER |
| CN2593215Y (en) * | 2003-01-13 | 2003-12-17 | 大庆油田有限责任公司 | Logging instrument relating to mixing and filling well tracer |
| CN203175536U (en) * | 2013-02-04 | 2013-09-04 | 西安思坦仪器股份有限公司 | Tracing correlation flowmeter |
| CN203499678U (en) * | 2013-10-16 | 2014-03-26 | 中国石油天然气股份有限公司 | Tracer related flow logging instrument based on double-channel releaser and single detection system |
| CN104265276A (en) * | 2014-09-12 | 2015-01-07 | 中国石油集团长城钻探工程有限公司测井公司 | Specific resistance tracer agent based flow measuring method and flowmeter |
| CN104481516A (en) * | 2014-11-26 | 2015-04-01 | 杭州泛太石油设备科技有限公司 | Continuous tracer logging method and logger thereof |
| WO2015097116A1 (en) * | 2013-12-23 | 2015-07-02 | Institutt For Energiteknikk | Particulate tracer materials |
-
2015
- 2015-12-31 CN CN201511032493.8A patent/CN105626039B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU365460A1 (en) * | 1970-12-11 | 1973-01-08 | Опытно конструкторское бюро геофизического приборостроени треста Бурмашремоит | METHOD OF MEASUREMENT OF OIL DEBIT IN WELLS WITH STAGENT WATER |
| CN2593215Y (en) * | 2003-01-13 | 2003-12-17 | 大庆油田有限责任公司 | Logging instrument relating to mixing and filling well tracer |
| CN203175536U (en) * | 2013-02-04 | 2013-09-04 | 西安思坦仪器股份有限公司 | Tracing correlation flowmeter |
| CN203499678U (en) * | 2013-10-16 | 2014-03-26 | 中国石油天然气股份有限公司 | Tracer related flow logging instrument based on double-channel releaser and single detection system |
| WO2015097116A1 (en) * | 2013-12-23 | 2015-07-02 | Institutt For Energiteknikk | Particulate tracer materials |
| CN104265276A (en) * | 2014-09-12 | 2015-01-07 | 中国石油集团长城钻探工程有限公司测井公司 | Specific resistance tracer agent based flow measuring method and flowmeter |
| CN104481516A (en) * | 2014-11-26 | 2015-04-01 | 杭州泛太石油设备科技有限公司 | Continuous tracer logging method and logger thereof |
Non-Patent Citations (2)
| Title |
|---|
| 基于ＭＡＴＬＡＢ相关流量测井处理软件的开发与应用;王现兵;《内蒙古石油化工》;20150630(第(2015年)12期);第28-30页 |
| 示踪流量测井技术在冀东油田的应用;董家东等;《石油钻采工艺》;20070731;第29卷(第(增刊)期);第106-109页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105626039A (en) | 2016-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20190316942A1 (en) | Methodologies and apparatus for the recognition of production tests stability | |
| CN104747152A (en) | Heavy oil deposit multi-element hot fluid throughput cross-well gas channeling degree detection method | |
| CN117077419B (en) | Novel formation pressure analysis method for fracture-cavity oil reservoir | |
| CN103352689B (en) | A kind of method utilizing radioactive tracer logging technology determination hole diameter | |
| CN111472760B (en) | A new method for extracting relative flow logging tracer peaks | |
| CN107842358B (en) | Method for extracting tracing logging information to form flow imaging graph | |
| CN105626039B (en) | A kind of method of production logging correlative flow data prediction | |
| US20210199000A1 (en) | Segmentation of time-frequency signatures for automated pipe defect discrimination | |
| CN102278109A (en) | Well leakage position ground and underground comprehensive detection system and method | |
| CN115906695B (en) | Multi-information fusion production profile logging interpretation optimization method | |
| CN112014881B (en) | Water flooding speed prediction method based on time lapse earthquake | |
| US11613986B1 (en) | Methods and systems for processing time-series well data using higher order channels to identify features therein and alter hydraulic fracturing operations based thereon | |
| CN104713631B (en) | The detection method of bulk sound velocity in a kind of oil well | |
| CN114483008A (en) | Method, system, device and storage medium for positioning rubber sealing ring of deep well tubing | |
| CN202596619U (en) | Rectangular tank type measuring device of drilling fluid returning flow | |
| CN103291279B (en) | The optimization method of the micro-abnormal signal of a kind of Gas Logging Value | |
| Huiqin et al. | Research on wavelet threshold denoising method of ultrasonic echo signal in ultrasonic drilling fluid leak detection system | |
| Gan et al. | Channeling analysis of wavelet threshold processing based on K-means clustering algorithm | |
| CN118818609B (en) | Method and device for detecting fluid mobility based on second-order synchronous extraction wavelet transformation | |
| Shi et al. | Determination of fluid properties and reservoir net pay cutoffs by production logging and conventional logs in exploration wells: a case study of the granite fractured reservoir in JZ oilfiled in Bohai sea | |
| CN119595063A (en) | Optical fiber flow measurement device, method and system based on peak heating | |
| CN113126165B (en) | Two-dimensional inclined shaft synthetic seismic record mosaic display method and device | |
| CN113534244B (en) | Method and device for determining the thickness of plumb bob encountered in drilling reservoir | |
| Jia et al. | Comprehensive model for multi-fracture localization based on water hammer signals: Evaluation and field application | |
| CN109992839B (en) | A method to avoid errors in inversion calculation caused by noise of water level data in water pumping experiments |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180608 Address after: 100007 Dongzhimen North Street, Dongcheng District, Dongcheng District, Beijing Applicant after: China National Petroleum Corporation Applicant after: China National Petroleum Group Logging Co., Ltd. Address before: 100007 Dongzhimen North Street, Dongcheng District, Dongcheng District, Beijing Applicant before: China National Petroleum Corporation Applicant before: Tubular Goods Research Center of China National Petroleum Corporation |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |