CN105117511B - A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter - Google Patents
A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter Download PDFInfo
- Publication number
- CN105117511B CN105117511B CN201510431255.8A CN201510431255A CN105117511B CN 105117511 B CN105117511 B CN 105117511B CN 201510431255 A CN201510431255 A CN 201510431255A CN 105117511 B CN105117511 B CN 105117511B
- Authority
- CN
- China
- Prior art keywords
- tracer
- well
- curve
- fracture hole
- flow
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to the Reservoir description field in oil field development, the characterizing method of especially a kind of fracture hole oil reservoir interwell communication passage and flow parameter comprises the following steps：(A) foundation describes the analytical expression of flowing law between fracture hole oil reservoir well；(B) tracer output concentration curve fit mathematics model is established；(C) cubic fitting model solves；According to principle described in above three step, fracture hole oil deposit tracking curve fitting software is worked out, theoretical tracer curve and actual measurement tracer curve are fitted by the method for the adjusting parameter on software interface, while fracture hole injection tracer quality m is calculatedik, flow velocity vik, area of passage Aik, flow tube length lik, so as to fulfill the characterization to communicating passage and flow parameter.The model for the description fracture hole oil deposit tracking agent flowing law that the present invention establishes is more in line with truth；Realize the quantitative computational representation to fracture-pore reservoir interwell communication passage and flow parameter；Fit procedure is combined by computer automatic Fitting with staff's tune ginseng.
Description
Technical field
The present invention relates to the Reservoir description field in oil field development, especially a kind of fracture hole oil reservoir interwell communication passage
And the characterizing method of flow parameter.
Background technology
At present both at home and abroad research inter well connectivity method mainly have wireline logging, Strata Comparison, geochemistry, well testing,
The modes such as the association analysis of injection-production well Production development, tracer flow, wherein first three methods belong to Static Assessment Method.For oil reservoir
For developer, what is focused more on is the dynamic connection situation during oil reservoir development between well.To oil reservoir continuity, it is particularly
The research of dynamic connectivity between injection-production well, is full appreciation reservoir geologic character, improves the element task content of oil recovery.
For the needs recognized oil reservoir, the research of inter well connectivity experienced qualitative analysis to the process for quantitatively seeking ginseng,
Related scholar studies the quantitative description of channel parameters well.Wang Mount Lushan et al. is directed to sandstone reservoirs high water-cut stage advantage
Seepage channel common problem, proposes that predominant pathway is identified using WATER FLOODING CHARACTERISTIC CURVE, and storage can be adopted using residue
Amount calculates predominant pathway volume, and then calculates predominant pathway permeability and pore throat radius.Once left a good name et al. and to use mathematics and oil reservoir
Engineering method establishes large pore path parameter quantitative calculation method.Zheng Qiang is basic by simultaneous mass-conservation equation, deliverability equation etc.
Equation derives the model that connected volume between injection-production well is calculated using Production performance data.Jiang Han bridges et al. establish base
In the Permeability Distribution probabilistic model of well log interpretation outcome table, it can be quantified by the model and calculate interwell communication predominant pathway
Parameter.Hu Shuyong et al. proposes the stochastic simulation mathematical model of macropore described in reservoir under the conditions of low speed laminar flow, and with should
The variation of large pore path parameter during modeling water drive.Wu Jun come et al. for conventional water-drive pool evaluation consider often
The problem of Statically Geologic Model, based on the simulation of streamline method, it is proposed that new water drive oil dynamic evaluation method realizes quantification
Relation between evaluation injection-production well.
One of the effective ways of inter-well tracer test technology as research interwell communication property, in grand celebration, huge port, the Liaohe River, Ji
The oil fields such as woods, North China, triumph, Xinjiang, Henan, Central Plains are used widely, and have been increasingly becoming the related sandstone oil reservoir well of description
Between the important means of anisotropism and flow behavior.Based on tracer flow data, domestic and foreign scholars are channel parameters well
Quantitatively characterizing method also studied.It is summed up, is broadly divided into two major class of analytic method and Method for Numerical.Brigham
The analytic method of the tracer quantitative interpretation communicating passage based on 5 well patterns is established earliest, for judging to test well group oil reservoir
Heterogeneous situation.Abbaszadeh and Brigham is on the basis of five-spot pattern flow of tracer characteristic prediction method, it is proposed that
With the analytic method of the inversion method in-place permeability of inter-well tracer test curve matching, pore throat radius and thickness.Zhang Yong establishes base
In the inter-well tracer test quantitative interpretation technology of streamline numerical simulation.Jia Jun mountains application enhancements streamline method numerical simulation technology describes flow field
Spread, characterize fluid neuron network rule.Computational methods of the inter-well tracer test technology in relation to interwell communication passage and flow parameter at present
Clastic rock oil reservoir is primarily directed to, theoretical foundation is seepage theory, and oil reservoir is regarded as in interpretation process is had not by multiple
Homogeneous substratum is formed in the layer of same permeability, pore throat radius and thickness, is fitted to obtain the corresponding ginseng of each layer by tracer curve
Number.
Earth's surface karst system is typical carbonate rock pipeline (crack)-hole-subterranean stream (lake) system, wherein
Fluid be mainly underground water, connection feature and fracture hole oil reservoir connectivity structure have certain similitude.Yang Lizheng, which is utilized, to be shown
Track curve is judged and is explained to karst waterpipe flow field structure, it is proposed that four kinds of karst connectivity types and its corresponding typical case
Tracer curve.R.G.Wilh and J.F.Quinlan is summarizing the analysis of karst area traces and is explaining in technical foundation, carries
8 kinds of tracer curve types and corresponding karst passage system feature are gone out.Mei Zhengxing is when summarizing domestic and international connection test method
It is also proposed that the tracer curve feature of 3 major class, 6 kinds of flow field structures.Zhang Zhenwu is theoretical according to Convention diffusion, by all kinds of karst waters
Pipe flow field geological conditions and the conclusion of tracer condition, give pipe flow field tracer mathematical model and analytic solutions.On connectivity types
Division be required for theoretic knowledge Karst System form, but considerably increase calculate Karst corridor parameter complexity,
Do not have operation possibility in actually calculating.
The development characteristics of fracture and vug carbonate reservoir reservoir have very strong anisotropism, cave, fracture development, reservoir
Thickness differs in size, cross direction profiles are unstable, and connective and floood conformance situation is increasingly complex, as shown in Figure 1.Fracture hole oil reservoir master
There is the three classes reservoir spaces such as seam, karst, solution cavity, this three space-like, which interconnects to form complicated oil gas and flow, to be led to
Road.
The content of the invention
The technical problems to be solved by the invention are to provide the table of a kind of fracture hole oil reservoir interwell communication passage and flow parameter
Sign method.
The technical solution adopted by the present invention to solve the technical problems is：A kind of fracture hole oil reservoir interwell communication passage and flowing
The characterizing method of parameter, comprises the following steps：
(A) foundation describes the analytical expression of flowing law between fracture hole oil reservoir well；
(B) tracer output concentration curve fit mathematics model is established；
(C) cubic fitting model solves；
(D) according to above three step, fracture hole oil deposit tracking curve fitting software is worked out, by being adjusted on software interface
The method of parameter is fitted theoretical tracer curve and actual measurement tracer curve, while fracture hole injection tracer injection is calculated
Quality, flow velocity, area of passage, flow tube length, so as to fulfill the characterization to communicating passage and flow parameter.
Preferably, the expression formula in the step (A) is：
In formula：C is that the concentration mg/L, D of tracer in flow field are dispersion coefficient m2/ d, v are for fluid flow velocity m/d, x
Length m, t are the sometime point h of tracer monitoring；
In formula：M is tracer quality mg, miFor tracer the quality mg, C of output in i-th of runnerbFor in body of groundwater
The concentration mg/L, C (t) of tracer are that the concentration mg/L, v of t moment withdrawal well tracer are fluid flow velocity m/d, and A is flowing
The equivalent cross-sectional area m of passage2, t is the sometime point h of tracer monitoring.
Preferably, the object function in the step (B) and constraints are respectively；
Purpose function：
In formula：Q is total number of days d of tracer monitoring, and t is the sometime point h of tracer monitoring, and p is to be monitored in a well group
To the number of the well of tracer output, k is containing p mouthfuls of withdrawal wells, C in a well groupk *(t) show for the well of serial number k in t moment
Field actual measurement obtains the concentration mg/L, C of tracerk(t) tracer is obtained by formula theoretical calculation in t moment for the well of serial number k
The concentration mg/L of agent；
In formula：CbkConcentration mg/L, n for tracer in k-th of well are that k-th of well surveys the wave crest occurred on tracer curve
Number, mikFor tracer the quality mg, v of output in i-th of runner of kth wellikIt is flowed for fluid in i-th of runner of k-th of well
Speed m/d, AikFor the equivalent cross-sectional area m of i-th of flow channel of kth well2, D is dispersion coefficient m2/ d, xikFor k-th of output
The length m, p of i-th of runner are the number for the well that tracer output is monitored in a well group between well and injection well；
Constraints：
In formula：M0For the initial total implantation quality mg of tracer of well group；
Lp≤xik≤10Lp
0 ＜ Aij≤5
In formula：LpThe straight line well spacing m, t between injection well and p-th of withdrawal wellTPFor after p-th of withdrawal well tracer output
The corresponding break through of one wave crest.
Preferably, the curve fitting algorithm calculation procedure of the step (C) is：
(1) an initial particle group containing 20 particles is generated；
(2) value of each particle is substituted into object function to calculate, more each target function value, searching meet local optimum and
The particle of global optimum；
(3) judge whether iteration reaches the condition of convergence, step (5) if the condition of convergence is reached, otherwise step (4)；
(4) calculate and determine iteration speed, generate the particle and Particle Swarm of a new generation, return to step (2)；
(5) optimizing terminates, and exports object function globally optimal solution and corresponding particle.
The beneficial effects of the invention are as follows：The model for the description fracture hole oil deposit tracking agent flowing law that the present invention establishes more accords with
Close truth；Realize the quantitative computational representation to fracture-pore reservoir interwell communication passage and flow parameter；Fit procedure by
Computer automatic Fitting is combined with staff's tune ginseng, calculates the combination for having taken into account theoretical calculation and practical experience.
Description of the drawings
Fig. 1 is TK221-TK214 measured curves schematic diagram in embodiment 1；
Fig. 2 is TK221-TK259 measured curves schematic diagram in embodiment 1；
Fig. 3 is TK221-TK250 measured curves schematic diagram in embodiment 1；
Fig. 4 is TK221-TK251 measured curves schematic diagram in embodiment 1；
Fig. 5 is TK634-TK747 measured curves form schematic diagram in embodiment 2；
Fig. 6 is TK634-TK744 measured curves form schematic diagram in embodiment 2；
Fig. 7 is TK634-TK711 measured curves form schematic diagram in embodiment 2；
Fig. 8 is TK634-TK625 measured curves form schematic diagram in embodiment 2；
Fig. 9 is TK634-T7-607 measured curves form schematic diagram in embodiment 2；
Figure 10 is TK634-S80 measured curves form schematic diagram in embodiment 2；
Figure 11 is TK634-T667 measured curves form schematic diagram in embodiment 2；
Figure 12 is tracer curve fitting software calculating parameter input interface schematic diagram；
Figure 13 is that tracer curve fitting software calculates main interface schematic diagram.
Specific embodiment
The present invention is described in further detail with reference to embodiment.
The characterizing method of a kind of fracture hole oil reservoir interwell communication passage and flow parameter, comprises the following steps：
(A) foundation describes the analytical expression of flowing law between fracture hole oil reservoir well；
Before to the flow of tracer rule description in the fracture hole reservoir system of underground, hypothesis below is carried out：
1. slight trace agent is ignored under reservoir conditions without chemical reaction, absorption and horizontal disperse, convection current is only considered
Influence of the diluting effect to tracer concentration；
2. assume that the slight trace agent of injection is all dissolved in water, and insoluble in oil；
3. it is migrated in seam_cavern type of the tracer aqueous solution injected between injection-production well in the form of pipe stream；
4. underground water is incompressible, flow velocity is the One-Dimensional flows of constant between injection-production well, tracer is along underground seam_cavern type
It migrates to withdrawal well (positive X-direction), and only flows to withdrawal well from injection well without external fluid replacement therebetween；
5. fracture hole oil reservoir interwell communication passage (fracture hole combining form) is reduced to the single channel and parallel port of circular tube shaped
Two big types characterize different communicating passage space sizes by different caliber and pipe range.
Based on assumed above, the fundamental differential of Tracer kinetics between injection-production well is obtained according to hydrodynamic dispersion equation
For：
In formula：Tracer concentration in C-flow field, mg/L；
D-dispersion coefficient, m2/d；
V-fluid flow velocity, m/d；
X-length, m；
T-moment point, h.
Section 1 represents the solute motion as caused by dispersion effect on the right of equation；Section 2 represents the convection current or total by liquid
Current caused by movement, last solute for being generated by chemical action of expression disappears and increases.If the tracer used
Agent is opposite nonabsorbable substance such as Cl-、NO3 -, chemical action does not occur with fracture hole medium when they migrate, then last can
To ignore, obtain
(1) single channel system
On the basis of foregoing 5 are assumed, then make the following assumptions：1. pipe flow path journey is semi-infinite long (0≤x ＜ ∞)；2. work as t
When=0, not tracer-containing between well；It is 3. instantaneous to the tracer that injection well input quality is M at x=0 in t=0.That
, mathematical model (3) as follows can be obtained
It is normal distribution that it, which is solved, i.e.,
Fixed in injection-production well position, each flow channel length x values are considered as fixed.T is variable, that is, is supervised
Survey time series.So formula (4) can be rewritten as with mathematic expectaion m=x, whose deviation markedForm, i.e.,
(2) parallel communication passage
Under certain geological conditions, underground flow tube be in net distribution, it is assumed that have n passage be parallel to injection well point A with
Between extraction well B.Each pipe stream length is respectively x1, x2…xn, flow velocity is respectively v1, v2…vn., for A points, quality is showing for M
Track agent is instantaneously injected, and the tracer quality that each passage is got between A-B is respectively m1, m2…mn, andObviously, every
Passage has the definite condition similar to single channel stream tracer field in tracer.Therefore, the concentration point of i-th passage output
Cloth can be write as
In withdrawal well B, the tracer concentration of parallel port flow field output exports the superposition of concentration for each runner.Then have
Formula (7) is a series of the sum of normal distributions formula.
Relatively simple and parallel port system tracer concentration output mathematic(al) representation understands that single channel is considered as
Situation of the parallel port expression formula in n=1.
When describing fracture hole oil reservoir connected space, the passage between injection-production well is equivalent to pipe-line system, concentration is led to
It can be represented with expression formula by formula (8,9).The flowing of fracture hole reservoir fluid is in the form of pipe stream, it is therefore assumed that water body seals in fracture hole
It closes, no outside fluid dilution.For a certain communicating passage between well, meet following expression in production end after tracer injection breakthrough
C0V0+CbVt=C ((t) (V0+Vt) (8)
Vt=vAt (9)
In formula：C0--is to launch initial time tracer and the mixed concentration of water, mg/L；
Cb--- for the concentration of tracer in body of groundwater, mg/L；
V0--- initial time reactor solution volume, m3；
C (t) --- the tracer concentration of t moment withdrawal well, mg/L；
Vt--- the volume that water body is supplemented after seam_cavern type is injected into for tracer, is approximately equal to fracture hole communicating passage body
Product, m3；
V --- fluid flow velocity, m/d；
A --- the equivalent cross-sectional area of flow channel, m2；
The sometime point of t --- tracer monitoring, h.
It can be obtained by (8)
C (t)=(C0V0+CbVt)/(V0+Vt) (10)
Due to V0<<Vt, therefore formula (10) can be reduced to
C0V0+CbVt=C (t) Vt (11)
Further obtain
Due to
M=C0V0 (13)
M --- the tracer quality that a certain communicating passage enters, mg.
Therefore, can be obtained by formula 9,12 and 13
For, there are n independent flow channels situation in parallel, production fluid concentration is each flow tube at a time concentration between well
Superposition
It can be obtained by formula (7) and formula (15), be in t moment extraction well production fluid concentration
Thus output concentration mathematic(al) representation of the fracture-pore reservoir inter-well tracer test curve in t moment is established.From formula 16,
The equation be there are the nonlinear equation of multiple variables, can not direct solution analysis solution.The tracer obtained by being fitted field monitoring
Agent output concentration curve can be with the relevant formation parameter of reverse.
(B) tracer output concentration curve fit mathematics model is established；
Fracture hole oil reservoir inter-well tracer test curve matching mathematical model is established, object function and constraints are as follows.
(1) object function
In formula：Total number of days of q-- tracer monitoring, unit d, the integer more than 0；
The sometime point of t-- tracer monitoring, h, the integer more than 0；
The number of the well of tracer output is monitored in p-- one well group, the integer more than 0；
Containing p mouthfuls of withdrawal wells, the wherein sequence number of some withdrawal well in k-- one well group, the integer more than 0；
The well of-- serial number k obtains the concentration of tracer, mg/L in t moment field measurement；
Ck(t) --- the well of serial number k obtains the concentration of tracer, mg/L in t moment by formula theoretical calculation.
In formula：CbkThe background concn of-- k-th well, mg/L；
The wave crest number (the flow channel number i.e. between well) occurred on the wells of n-- k-th actual measurement tracer curve, n >=1, integer；
mikThe tracer quality of output, mg in i-th of runner of-- kth well；
vikFluid-flow rate in-- i-th runner, m/d；
AikThe equivalent cross-sectional area of i-th of flow channel of-- kth well, m2, it is positive number；
D-- dispersion coefficient, m2/d；
xikThe length of i-th of runner between-- k-th withdrawal well well and injection well, m are positive number；
P --- the number of the well of tracer output is monitored in a well group；
(2) constraints
Merely from the point of view of mathematical model solution, mik、vik、Aik、xikThere is the possibility more solved, therefore need to be to its value model
It encloses and is defined.Parameter value scope to be asked is fitted spreadsheet dataset statistical result by substantial amounts of tracer curve and combines live expert
It is empirically determined, it is as follows：
1. quality constrains
M0The initial total implantation quality of the tracer of-- well group, mg；From m from the point of view of actual fitting resultik＜ ＜ M0, reflect
There is substantial amounts of be detained in underground in tracer.
2. flow velocity constrains
Lp--- straight line well spacing, m between injection well and p-th of withdrawal well；
tTP--- the corresponding break through of first wave crest after p-th of withdrawal well tracer output.
3. interwell communication passage length constrains
Lp≤xik≤10Lp
4. the equivalent area of passage constraint of interwell communication passage
0 ＜ Aij≤5
(C) cubic fitting model solves
From tracer curve fit mathematics model, the tracer quality m of each runner distribution between wellik, flow velocity vik, flowing
The equivalent cross-sectional area A of passageikWith length xikFor parameter to be asked.Object function is nonlinear equation, and 4 parameters to be asked are to connect
Continuous property variable, this mathematical model are the Optimization Solution problem of nonlinear function.Optimize for the solution generally use of this class function
Algorithm is realized.The present invention is used for the solution of tracer curve model of fit using particle swarm optimization.
There are the m particles without volume in D dimension spaces, a search group formed by several particles of the same generation,
The position coordinates of i-th of particle is Xi=(xi1,xi2,….xiD), component represents each variable respectively.
(1) generation of particle populations
The parameter to be asked of this paper is divided into four classes, is respectively quality mik, flow velocity vik, area of passage Aik, flow tube length lik, (k
=1 ..., p, p be withdrawal well number in well group), (i=1 ..., n, n are the wave crest number that a certain well needs are fitted).By parameter
Lower target variation understands that the quantity of parameter to be asked is variation, is determined by well number and wave crest number to be fitted.
The fracture hole parameter set combined by unknown variable in the case where meeting constraints is a particle Xi, wherein
Each variable corresponds to Xi one-components, and fracture hole parameter set just forms a particle populations.Initial particle group is random by system
Generate and meet 20 different particles of constraints, form a particle populations X_particles=X_particle (1),
X_particle (2) ..., X_particle (20) }, n_particle represent Particle Swarm in each particle number, value 1~
20, integer
(2) movement of Particle Swarm
Each particle is flown in search space with certain speed to desired positions (object function minimum value), residing for it
Position regard as when former generation part desired positions (scheme meets constraints, smaller than the preceding target function value once evolved)
It is denoted as Pi=(pi1,pi2,…piD), also referred to as pbest；Desired positions (the current evolutionary generation that all particles of group live through
Under, meet constraints works as the corresponding fracture hole parameter set of former generation minimum target function) it is when former generation globally optimal solution, it uses
Gbest is represented.
1. translational speed determines
The speed V of particle movementi=(vi1,vi2,..viD) represent, representation in components respectively ties up step-size in search.D dimensions point
Amount (1≤d≤D) changes according to equation below：
K represents the evolutionary generation that Particle Swarm is currently undergone in formula；ω is inertia weight, it makes particle keep motional inertia；
c1And c2For aceleration pulse, they make each particle be accelerated to the position of pbest and gbest, r1And r2For [0,1] scope
Interior variation is uniformly distributed pseudo random number.
ω's determines：
2. the generation of new particle
The generation formula of new particle
By taking the mass component in particle as an example
(3) fitness is evaluated
The part and global optimum for defining the iteration of kth _ iterate times are respectively y_p_best (k_iterate) and y_
G_best (k_iterate),
X_p_particle (k_iterate) be the corresponding particle of local best points, X_g_particle (k_iterate)
For the corresponding particle of globe optimum.
1. compare each particle target function value y_particle (n_particle) and part during kth _ iterate iteration
The size of optimal value, if using better than prior-generation local optimum there are particle object function when former generation and working as the optimal of former generation
Prior-generation local optimum is replaced, otherwise local optimum remains unchanged.
2. compare the local optimum under current iteration number and the global optimum under historical conditions.If when former generation local optimum
It is better than prior-generation global optimum, then replaces prior-generation global optimum with when former generation local optimum, otherwise global optimum remains unchanged.
(4) judge calculate whether meet the condition of convergence
Iterations searches for termination after reaching default maximum k_iterate_max in calculating.
To sum up, use particle swarm optimization solve curve fitting algorithm calculation procedure for：
(1) an initial particle group containing 20 particles is generated；
(2) value of each particle is substituted into object function to calculate, more each target function value, searching meet local optimum and
The particle of global optimum；
(3) judge whether iteration reaches the condition of convergence, step (5) if the condition of convergence is reached, otherwise step (4)
(4) calculate and determine iteration speed, generate the particle and Particle Swarm of a new generation, return to step (2)
(5) optimizing terminates, and exports object function globally optimal solution and corresponding particle.
(D) realization of fracture hole oil reservoir interwell communication passage and flow parameter characterization
According to principle described in foregoing three steps, establishment fracture hole oil deposit tracking curve fitting software such as Figure 12 and Figure 13 institutes
Show, for calculating interwell communication passage and flow parameter.
Theoretical tracer curve and actual measurement tracer curve are fitted by the method for the adjusting parameter on software interface, together
When fracture hole injection tracer implantation quality, flow velocity, area of passage, flow tube length is calculated, so as to fulfill to communicating passage and
The characterization of flow parameter.
Embodiment 1：TK221 well groups.
Injection well is TK221, and 4 mouthfuls of tracer withdrawal wells TK214, TK258, TK250CH, TK251CH are corresponded in supervised period.
From the point of view of tracer curve, TK214 wells monitor two independent wave crests, illustrate between TK221-TK214 wells that there are two independently to connect
Circulation passage, the tracer curve that remaining 3 mouthfuls of well each monitors is unimodal curve, and communicating passage is single between illustrating each injection-production well
Passage, as shown in measured curve in Fig. 1,2,3,4.Master data used in fitting (is respectively counted on actual measurement tracer curve as shown in table 1
According to slightly)
1 TK221 well groups tracer curve of table is fitted basic parameter
Well group fitting result is as shown in theoretical curve in table 2 and Fig. 1,2,3,4.
2 TK221 well groups of table connect parameter result of calculation table
The flow velocity of communicating passage, volume (overcurrent between each injection-production well of TK221 well groups can be obtained by tracer curve fitting
Area and passage length), the tracer quantum of output of each well and the injection fluid that is thus calculated be between each well of well group
Apportionment ratio.
Embodiment 2：TK634 well groups.
Injection well is TK634, and 7 mouthfuls of tracer withdrawal wells TK747, TK744, TK711, TK625, T7- are corresponded in supervised period
607、S80、TK667.Region is developed in the wellblock for weathering crust, and from the point of view of tracer curve is surveyed, each well tracer curve is by several
Mixed and disorderly random small wave crest composition, the flowing space really forms parallel connection by many small-sized communicating passages between reflecting each injection-production well
Channel design, as shown in measured curve form in Fig. 5 to 11.
Measured curve can need approximation to regard as with several main wave crests according to explanation.TK214 wells monitor two independences
Wave crest, illustrate that there are the tracer curves that two independent communication passages, remaining 3 mouthfuls of well each monitor between TK221-TK214 wells
It is unimodal curve, communicating passage is single channel between illustrating each injection-production well, as shown in measured curve in Fig. 1,2,3,4.Fitting
Master data used is as shown in table 1 (each point data is omited on actual measurement tracer curve)
It is as shown in table 3 to be fitted primary data
3 TK634 well groups tracer curve of table is fitted basic parameter
Well group fitting result is as shown in table 4 and Fig. 5,6,7,8,9,10,11.
4 TK634 well groups of table connect parameter result of calculation table
Claims (3)
1. the characterizing method of a kind of fracture hole oil reservoir interwell communication passage and flow parameter, which is characterized in that comprise the following steps：
(A) foundation describes the analytical expression of flowing law between fracture hole oil reservoir well；
(B) tracer output concentration curve fit mathematics model is established；
(C) cubic fitting model solves；
(D) according to above three step, fracture hole oil deposit tracking curve fitting software is worked out, passes through the adjusting parameter on software interface
Method to theoretical tracer curve and actual measurement tracer curve be fitted, while be calculated fracture hole injection tracer injection matter
Amount, flow velocity, area of passage, flow tube length, so as to fulfill the characterization to communicating passage and flow parameter；
Expression formula in the step (A) is：
In formula：C is that the concentration mg/L, D of tracer in flow field are dispersion coefficient m2/ d, v are length for fluid flow velocity m/d, x
M, t are the sometime point h of tracer monitoring；
In formula：M is tracer quality mg, miFor tracer the quality mg, C of output in i-th of runnerbFor tracer in body of groundwater
The concentration mg/L, C (t) of agent are that the concentration mg/L, v of t moment withdrawal well tracer are that fluid flow velocity m/d, A are flow channel
Equivalent cross-sectional area m2, t is the sometime point h of tracer monitoring.
2. the characterizing method of a kind of fracture hole oil reservoir interwell communication passage as described in claim 1 and flow parameter, feature exist
In object function and constraints in the step (B) are respectively；
Purpose function：
In formula：Q is total number of days d of tracer monitoring, and t is the sometime point h of tracer monitoring, and p is to monitor to show in a well group
The number of the well of track agent output, k are to contain p mouthfuls of withdrawal wells in a well group,It is real at t moment scene for the well of serial number k
Measure the concentration mg/L, C of tracerk(t) tracer is obtained by formula theoretical calculation in t moment for the well of serial number k
Concentration mg/L；
In formula：CbkConcentration mg/L, n for tracer in k-th of well are that k-th of well surveys the wave crest number occurred on tracer curve,
mikFor tracer the quality mg, v of output in i-th of runner of kth wellikFor fluid flowing speed in i-th of runner of k-th of well
Spend m/d, AikFor the equivalent cross-sectional area m of i-th of flow channel of kth well2, D is dispersion coefficient m2/ d, xikFor k-th of withdrawal well
The length m, p of i-th of runner are the number for the well that tracer output is monitored in a well group between injection well；
Constraints：
In formula：M0For the initial total implantation quality mg of tracer of well group；
Lp≤xik≤10Lp
0 ＜ Aik≤5
In formula：LpThe straight line well spacing m, t between injection well and p-th of withdrawal wellTPFor first after p-th of withdrawal well tracer output
The corresponding break through of wave crest.
3. the characterizing method of a kind of fracture hole oil reservoir interwell communication passage as claimed in claim 2 and flow parameter, feature exist
In the curve fitting algorithm calculation procedure of the step (C) is：
(1) an initial particle group containing 20 particles is generated；
(2) value of each particle is substituted into object function to calculate, more each target function value, searching meets local optimum and the overall situation
Optimal particle；
(3) judge whether iteration reaches the condition of convergence, step (5) if the condition of convergence is reached, otherwise step (4)；
(4) calculate and determine iteration speed, generate the particle and Particle Swarm of a new generation, return to step (2)；
(5) optimizing terminates, and exports object function globally optimal solution and corresponding particle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510431255.8A CN105117511B (en) | 2015-07-22 | 2015-07-22 | A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510431255.8A CN105117511B (en) | 2015-07-22 | 2015-07-22 | A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105117511A CN105117511A (en) | 2015-12-02 |
| CN105117511B true CN105117511B (en) | 2018-06-05 |
Family
ID=54665499
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510431255.8A Expired - Fee Related CN105117511B (en) | 2015-07-22 | 2015-07-22 | A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105117511B (en) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106761708B (en) * | 2017-01-20 | 2019-09-17 | 中国石油大学(北京) | A kind of water drive inter-well tracer test test interpretation method |
| CN109426657B (en) * | 2017-08-29 | 2022-09-30 | 中国石油化工股份有限公司 | Method and system for calculating inter-well communication coefficient of oil reservoir |
| CN108505981A (en) * | 2018-03-20 | 2018-09-07 | 西南石油大学 | A kind of velocity flow profile analysis method of Fractured reservoir injection profile agent |
| CN108999608A (en) * | 2018-06-14 | 2018-12-14 | 中国石油天然气股份有限公司 | Low-permeability conglomerate reservoir dominant channel identification method and system |
| CN108843300B (en) * | 2018-06-25 | 2022-03-01 | 中国石油天然气股份有限公司 | Method and device for determining type of main flow channel in complex porous medium |
| CN110965996A (en) * | 2018-09-28 | 2020-04-07 | 中国石油天然气股份有限公司 | Method for detecting dynamic communication strength between wells |
| CN109538193A (en) * | 2019-01-28 | 2019-03-29 | 西南石油大学 | A kind of runner phenetic analysis method applied to fracture-pore reservoir |
| CN109815600A (en) * | 2019-01-28 | 2019-05-28 | 西南石油大学 | A method of parsing fracture-pore reservoir tracer production curve feature |
| CN110162922A (en) * | 2019-06-03 | 2019-08-23 | 西南石油大学 | A kind of integrated recognition method of water-drive pool dominant flowing path |
| CN110259426B (en) * | 2019-07-02 | 2021-12-07 | 捷贝通石油技术集团股份有限公司 | Method for evaluating pressure channeling degree between unconventional platform wells |
| CN110348154B (en) * | 2019-07-17 | 2022-07-01 | 中国石油化工股份有限公司西北油田分公司 | Fracture-cavity type oil reservoir well group gas tracer interpretation and karst cave identification method |
| CN110593829B (en) * | 2019-08-14 | 2021-06-04 | 中国地质大学（武汉） | Automatic judgment method and device for interwell communication mode of fracture-cavity type oil reservoir |
| CN110644975B (en) * | 2019-09-27 | 2022-10-21 | 西安石油大学 | Fracture-cavity type oil reservoir tracer curve quantitative interpretation method |
| CN110671104B (en) * | 2019-11-01 | 2022-10-14 | 中国石油化工股份有限公司 | Interpretation method of interwell parameters of fracture-cavity type oil reservoir based on interference well testing interpretation |
| CN111859607B (en) * | 2020-05-09 | 2023-08-25 | 中国海洋石油集团有限公司 | Inter-well cross-flow channel quantification method and device considering in-channel heterogeneity |
| CN117076956B (en) * | 2023-10-16 | 2024-01-26 | 西安石油大学 | Fracture-cavity oil reservoir physical model similarity criterion optimization method and device |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102953726A (en) * | 2011-08-22 | 2013-03-06 | 中国石油大学(华东) | Method and device for water drive oilfield advantage channel recognition |
| CN103588801A (en) * | 2013-12-02 | 2014-02-19 | 天津大港油田圣达科技有限公司 | Trace element well-to-well tracer agent use method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6512371B2 (en) * | 1995-10-12 | 2003-01-28 | Halliburton Energy Services, Inc. | System and method for determining oil, water and gas saturations for low-field gradient NMR logging tools |
-
2015
- 2015-07-22 CN CN201510431255.8A patent/CN105117511B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102953726A (en) * | 2011-08-22 | 2013-03-06 | 中国石油大学(华东) | Method and device for water drive oilfield advantage channel recognition |
| CN103588801A (en) * | 2013-12-02 | 2014-02-19 | 天津大港油田圣达科技有限公司 | Trace element well-to-well tracer agent use method |
Non-Patent Citations (1)
| Title |
|---|
| "井组示踪产出曲线解释方法研究";查锋;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20110415(第4期);第14页至第51页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105117511A (en) | 2015-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105117511B (en) | A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter | |
| CN106837297B (en) | Method for identifying connectivity among wells and predicting oil-water dynamic state | |
| Eshraghi et al. | Optimization of miscible CO2 EOR and storage using heuristic methods combined with capacitance/resistance and Gentil fractional flow models | |
| CN103410502B (en) | A kind of acquisition methods of three-dimensional permeability fields of network-like fracture-pore reservoir | |
| Deli et al. | Data-driven optimization for fine water injection in a mature oil field | |
| MX2012009273A (en) | Method and system for predicting the effect of microbes injected into an oil-bearing reservoir. | |
| CN108603402A (en) | The variation of capillary pressure and relative permeability in porous media caused by mineral deposits and dissolving is modeled and predicted | |
| CN102041995A (en) | System for monitoring complicated oil deposit flooding conditions | |
| Li et al. | Modeling tracer flowback in tight oil reservoirs with complex fracture networks | |
| Lee et al. | Fracture-based modeling of complex flow and CO2 migration in three-dimensional fractured rocks | |
| CN104060985A (en) | Method and system for testing entering depth of stratified oil deposit profile control water plugging agent | |
| Li et al. | An analysis of tracer flowback profiles to reduce uncertainty in fracture-network geometries | |
| Zhai et al. | Production optimization in waterfloods with a new approach of inter-well connectivity modeling | |
| Xu et al. | Simulating seawater intrusion in a complex coastal karst aquifer using an improved variable-density flow and solute transport–conduit flow process model | |
| Liu et al. | A Discrete fracture–matrix model for pressure transient analysis in multistage fractured horizontal wells with discretely distributed natural fractures | |
| Liu et al. | A rapid waterflooding optimization method based on INSIM-FPT data-driven model and its application to three-dimensional reservoirs | |
| Gao et al. | Effect of density, trace length, aperture, and direction angle on permeability performance of fracture networks | |
| Chai et al. | Study on recovery factor and interlayer interference mechanism of multilayer co-production in tight gas reservoir with high heterogeneity and multi-pressure systems | |
| Du et al. | A deep learning framework using graph convolutional networks for adaptive correction of interwell connectivity and gated recurrent unit for performance prediction | |
| Chang et al. | Post-frac evaluation of deep shale gas wells based on a new geology-engineering integrated workflow | |
| Zeng et al. | Evaluating interwell connectivity in waterflooding reservoirs with graph-based cooperation-mission neural networks | |
| CN106481315A (en) | Land sandstone oil reservoir individual well recoverable reserves quickly determines model and method for building up | |
| Zhao et al. | A flownet-based method for history matching and production prediction of shale or tight reservoirs with fracturing treatment | |
| CN110644975A (en) | Fracture-cavity type oil reservoir tracer curve quantitative interpretation method | |
| US11501043B2 (en) | Graph network fluid flow modeling |
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: 20180605 Termination date: 20190722 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |