CN106980758A - A kind of quick calculation method of flooding pattern flow field velocity - Google Patents

A kind of quick calculation method of flooding pattern flow field velocity Download PDF

Info

Publication number
CN106980758A
CN106980758A CN201710159951.7A CN201710159951A CN106980758A CN 106980758 A CN106980758 A CN 106980758A CN 201710159951 A CN201710159951 A CN 201710159951A CN 106980758 A CN106980758 A CN 106980758A
Authority
CN
China
Prior art keywords
computational methods
flow
speed
flow field
well
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710159951.7A
Other languages
Chinese (zh)
Other versions
CN106980758B (en
Inventor
郭文敏
吕爱华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changzhou University
Original Assignee
Changzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changzhou University filed Critical Changzhou University
Priority to CN201710159951.7A priority Critical patent/CN106980758B/en
Publication of CN106980758A publication Critical patent/CN106980758A/en
Application granted granted Critical
Publication of CN106980758B publication Critical patent/CN106980758B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Abstract

The present invention provides a kind of quick calculation method of flooding pattern flow field velocity, including note adopts streamline approximate calculation method, streamline and adopts correction factor M computational methods, the computational methods of flooding pattern flow field velocity of the angle to speed to the correction factor C computational methods of speed, note around tomography computational methods, flow curvature.According to it is described note adopt streamline approximate calculation method, streamline around tomography computational methods set up injection-production well between Flow Field Distribution, correction factor of the angle to the correction factor M computational methods computational flow speed of speed is adopted to the correction factor C computational methods of speed, note using the flow curvature, flow field velocity calculating is carried out using the computational methods of the flooding pattern flow field velocity.Permeation fluid mechanics method can not embody oil reservoir Heterogeneous Characteristics and method for numerical simulation process complexity, length research cycle, to personnel requirement limitation during the present invention solves Study on Flow Field so that the quick calculating of flooding pattern fluid velocity is achieved.

Description

A kind of quick calculation method of flooding pattern flow field velocity
Technical field
The present invention relates to the oil reservoir Flow Field Distribution and flow field velocity research field in oil reservoir development adjustment, in particular it relates to A kind of quick calculation method field of flooding pattern flow field velocity.
Background technology
Especially no matter harvesting is improved in oilfield development process using conventional hydrodynamics method in ultra-high water cut stage Rate, or using the tertiary oil recovery mode of injection chemical agent, top priority is to get oil reservoir injection clear to adopt flow field and flow velocity size, bright True water breakthrough direction determines development contradiction, according to flow field contradiction, carries out targetedly Well pattern edjustment conceptual design and liquid measure optimization is set Meter.It can be seen that on the basis of flooding pattern, a kind of quick calculation method of flooding pattern flow field velocity is set up, for enriching hydrodynamic force Method, perfect scheme design theory, improve conceptual design efficiency and effect there is very useful engineering theory and valency is applied Value.
The research method in flooding pattern flow field mainly include Interwell tracer, well-logging method, reservoir engineering method, based on most Inter well connectivity evaluation method, method for numerical simulation and the streamline simulation method of the principle of optimality, wherein Interwell tracer, well testing side Method construction cost is high, research cycle is long;Reservoir engineering method requires high to researcher's professional qualities, and analysis result is only limited to In qualitatively;Optimal method is a kind of fuzzy Judgment method completely based on mathematical theory, and as a result ambiguity can use by force Property is low;Numerical simulation and streamlined impeller analogy method are stressed the influence of field computation precision, Method And Principle science but resultant error Larger actual operation is poor.The limitation of various methods is directed to as can be seen here, it is necessary to which setting up a set of can cover geology, well pattern Well spacing, note adopt dynamic flooding pattern fluid velocity field computation method that is quick, simple, quantifying to be had for reservoir engineering research Very important meaning.
For tracer testing method the high periodically long, reservoir engineering method of cost qualitatively quantify by force it is weak, optimal The features such as strong periodicity length with a low credibility and numerical simulation streamlined impeller method of the ambiguity of change method is error big, to ooze Flow based on theory of mechanics, reservoir heterogeneity is for seepage flow approximate processing, consideration well are carried out streamline track injection-production well Influence, the change of consideration streamline Macroscopic Curvature and influence of the angle to flow stream velocity between injection-production well of speed, are built using principle of stacking It is vertical to consider the quick calculation method that reservoir properties, Well Pattern And Spacing, note adopt dynamic flooding pattern fluid velocity.
The content of the invention
The main purpose of the embodiment of the present invention is to provide a kind of quick calculation method of flooding pattern flow field velocity, to solve Either certainly existing flow field velocity computational methods method is complicated, cost is high, the cycle is long, or it is strong qualitatively to quantify weak, ambiguity by force Defect with a low credibility.
To achieve these goals, the embodiment of the present invention provides a kind of quick calculation method of flooding pattern flow field velocity, Including:
Note adopts streamline approximate calculation method, streamline and the correction factor C of speed is calculated around tomography computational methods, flow curvature Method, note adopt correction factor M computational methods, the computational methods of flooding pattern flow field velocity of the angle to speed.
According to it is described note adopt streamline approximate calculation method, streamline around tomography computational methods set up injection-production well between Flow Field Distribution;
Respectively using the flow curvature correction factor C computational methods of speed, note are adopted angle to the correction of speed because Sub- M computational methods, obtain flow field velocity correction factor C and correction factor M;
Flow field velocity calculating is carried out using the computational methods of the flooding pattern flow field velocity.
The angle being association of activity and inertia by means of above-mentioned technical proposal, the present invention from oil reservoir, is carried out to flooding pattern streamline track The change of track Macroscopic Curvature, note adopt influence of the angle for flow field velocity between anisotropism, well between approximate processing, consideration well, profit With principle of stacking, set up and consider that reservoir properties, Well Pattern And Spacing, note adopt a kind of dynamic quick calculating of flooding pattern flow field velocity Method.Compared to prior art feature, the present invention solves method for numerical simulation Complicated Periodic length, tracer testing technical costs The high cycle is long, and reservoir engineering method qualitatively quantifies by force contradiction weak, that optimal method ambiguity is strong so that flooding pattern flow field The calculating of speed is simple, quick, quantify to be achieved, and further enriches hydrodynamics research method, perfect conceptual design It is theoretical, improve conceptual design efficiency and effect, the invention has very useful engineering theory application value.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, embodiment will be described below In required for the accompanying drawing that uses be briefly described, it should be apparent that, drawings in the following description are only some of the present invention Embodiment, for those of ordinary skill in the art, without having to pay creative labor, can also be according to these Accompanying drawing obtains other accompanying drawings.
Fig. 1 is the schematic flow sheet that the flooding pattern flow field velocity that the embodiment of the present invention one is provided quickly is calculated;
Fig. 2 is the rhombus model trajectory calculation principle schematic that the embodiment of the present invention one is provided;
Fig. 3 is the straight line model trajectory calculation principle schematic that the embodiment of the present invention one is provided;
Fig. 4 is the passing fault streamline schematic diagram that the embodiment of the present invention one is provided;
Fig. 5 is the Macroscopic Curvature correction principle schematic diagram that the embodiment of the present invention one is provided;
Fig. 6 is the circular arc stream tube model decomposing schematic representation that the embodiment of the present invention one is provided;
Fig. 7 is the circular expanding angle calculation of arc length model that the embodiment of the present invention one is provided;
Fig. 8 is the expanding angle calculation of straight line model that the embodiment of the present invention one is provided;
Fig. 9 is the expanding angle calculation of rhombus model that the embodiment of the present invention one is provided;
Figure 10 is that the note that the embodiment of the present invention one is provided adopts angle calcu-lation schematic diagram;
Figure 11 is the heterogeneous linear model schematic diagram of transverse direction that the embodiment of the present invention one is provided;
Figure 12 be the embodiment of the present invention one provide anisotropism 4.:3.=400:400、④:3.=550:250、④:③ =700:100 flow field velocity figures.
Embodiment
Below in conjunction with 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 Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made Embodiment, belongs to the scope of protection of the invention.
Embodiment one
The present embodiment provides a kind of quick calculation method of flooding pattern flow field velocity, as shown in figure 1, this method includes:
Step 1:Note adopts streamline approximate calculation method
If P coordinates in any point are (X between injection-production wellp, Yp), any point P to water injection well Inj (X is calculated respectively1, Y1) Apart from d1And to producing well Pro (X2, Y2) apart from d2
(1) rhombus model
For straight line model as shown in Fig. 2 straightway d1Slope K1For:
Straightway d1Equation be:Y=K1(X-X1)+Y1
Injection-production well line midpoint O coordinates:
The vertical line slope K of injection-production well line2For:
The vertical line equation for obtaining midpoint O and injection-production well line is:Y=K2(X-X0)+Y0
Order:Straightway d1It is (X with crossing midpoint O with injection-production well line vertical line intersection point T coordinatesT,YT), thus obtain:
YT=K1(XT-X1)+Y1Or YT=K2(XT-X0)+Y0
It is (X by T coordinatesT,YT), calculating is obtained:
(2) straight line model
For straight line model as shown in figure 3, directly calculating any point P (X between injection-production wellp, Yp) arrive water injection well Inj (X1, Y1) apart from d1With to producing well Pro (X2, Y2) apart from d2, this method is compared with rhombus model calculating speed faster.
Step 2:Streamline is around tomography computational methods
Micro-small fault is there may be for fault block oil reservoir, between water injection well and producing well to block, when tomography extension is not enough to When playing packing effect for injection-production well flow field, injection water will form injection water around tomography around passing fault endpoint stream to producing well Phenomenon is specific as shown in figure 4, now, for arbitrfary point P in oil reservoir, can be divided into two kinds of situations progress streamlines approximate.
(1) P1 points are in the position between water injection well and tomography
Streamline can flow to tomography end points, i.e. Inj-M1-A and I-M1-B by two kinds of approach, because ∠ I-M1-A are small It is turn back line in 90 °, can not possibly exists among reality, it is seen that the simplification path of the position M1 points between water injection well and tomography Mode can carry out approximately, thereby determining that water injection well Inj by fault barrier effect and producing well using Inj-M1-B-Pro Pro streamline should approximately meet following condition:
①∠Inj-M1-B>90°
②∠M1-B-Pro>90°
3. Dis (Inj-M1-B-Pro)=d01+d02+d1<Dmax
Wherein:Dis (Inj-M1-B-Pro) is the length in Inj-M1-B-Pro paths;Dmax is well testing injection-production well control half 2 times of footpath.
(2) P2 points are in the position between producing well and tomography
For producing well side, streamline can flow to producing well Pro points by two kinds of approach, i.e. Inj-A-M2-Pro and Inj-B-M2-Pro, is turn back line because ∠ Pro-M2-B are less than 90 °, can not possibly exist among reality, it is seen that in life Production well can be carried out approximate with the simplification path fashion of the position M2 points between tomography using Inj-A-M2-Pro.Thereby determine that Producing well P should be met following condition by fault barrier effect is approximate with water injection well I streamline:
①∠Pro-M2-A>90°
②∠M2-A-Inj>90°
3. Dis (Inj-A-M2-Pro)=d0+d11+d12<Dmax
Step 3:Correction factor C computational methods of the flow curvature to speed
(1) Macroscopic Curvature correction principle
The gradation type bend pipe flowing shown in Fig. 6 can be expressed as flow tube unit model during plane seepage flow Fig. 5, For the flowing of gradation type bend pipe, it can be decomposed into:Equal-diameter elbow flow process+flaring type flow process+gradual shrinkage flow process.
1. for equal-diameter elbow flow process
Flowed in bend pipe, flow direction changes, and speed is constant, the waterhead fall that bend pipe is produced depends on In the angle, θ and the ratio of radius of curvature size and pipe diameter size of bend pipe, bend pipe often uses resistance coefficient expression formula:
In formula:ξ is resistance coefficient;R, R are respectively elbow radius and radius of curvature;α is bend pipe angle.
For the calculating of 180 ° of bend pipe resistance coefficient, 2 times of right-angle elbow pipe resistance coefficient approximate calculation can be used.
2. for flaring type flow process
Expression formula for increaser local loosening ξ is calculated as:
In formula:λ is friction drag coefficient, for laminar flow,V is flowing velocity, m/s;D is interior Footpath, m;θ is expanding angle, °;A1、A2For inlet/outlet sectional area;
3. for gradual shrinkage flow process
For gradual shrinkage flow process, resistance coefficient change is little, typically takes ξ=0.04.
(2) calculating of expanding angle, θ
For the calculating of expanding angle, θ, circular arc length model, straight line model, three kinds of models and rhombus model can be passed through Calculated, comprised the following steps that.
1. circular arc length model
Calculating for the circular expanding angle, θ of arc length model is as shown in fig. 7, make two arc distances from L3It is set as that grid is maximum wide Degree, expanding angle, θ computational methods:
Because:L4=R1-R1cosθ1-L3
Because:So:θ=θ12
2. straight line model
For the expanding angle, θ of straight line model calculating as shown in figure 8, calculating process is as follows.
3. rhombus model
For the expanding angle, θ of rhombus model calculating as shown in figure 9, calculating process is as follows.
It can be seen that, the extra coefficient of resistance caused by flow tube is bent with flow tube gradual change is:
(3) curvature correction coefficient C is determined
Frictional resistance change ratio is that curvature correction coefficient C is:
Step 4:Note adopts correction factor M computational methods of the angle to speed
In the case where flooding pattern is controlled (A, B, C be producing well, D, E be water injection well) as shown in Figure 10, on reservoir Plane Fluid Field The size of any point O flow velocitys by static parameters such as distance, reservoir properties in addition to being influenceed, also by the point and water injection well, production The angle ∠ AOD of well line formation, ∠ BOD, ∠ COD, ∠ AOE, ∠ BOE, ∠ COE sizes are relevant, it can be seen that:
(1) ∠ BOE, ∠ COE are respectively less than 90 °, it is seen that O points are influenceed very little by control between water injection well E and producing well B, C well, It can ignore;
(2) ∠ AOE are more than 90 °, from injection-production relation, and O points are had one by control between water injection well E and producing well A wells Fixed possibility;
(3) ∠ AOD, ∠ BOD, ∠ COD are all higher than 90 °, wherein ∠ BOD>∠COD>∠ AOD, from injection-production relation, angle Degree is bigger, and control ability is stronger;
As can be seen here, any point O in reservoir plane, when the angle of itself and injection-production well line formation is less than 90 °, can be neglected For the control ability of the point slightly between injection-production well, when angle is more than 90 °, with the increase of angle, control ability is incrementally increased, When angle is equal to 180 °, O points are on main stream line, and most strong by controlling extent between injection-production well well, percolation flow velocity is bigger.
Consider that water filling that the difference on injection-production well line perpendicular direction is caused is advanced by leaps and bounds and injected water horizontal proliferation, is obtained Note adopts correction factor M of the angle to speed:
M=(cos δ)n
In formula:M is angle modifying factor between injection-production well;δ is the point and the angle of water injection well producing well line formation;N is Injection-production well line permeability in plane with being perpendicularly to the direction the ratio between permeability.
Step 5:The computational methods of flooding pattern flow field velocity
Know that any point O corresponds to i-th in oil reservoir by the principle of stacking of the Isoquant multi-source remittance infinity plane radial direction force and velocity of a current Mouthful well is in the velocity that the point is produced:
It is preferred that, in above-mentioned steps 5, it particularly may be divided into following Four processes:
Process 1:Non-homogeneous model computing permeability
For linear seepage flow, it is assumed that streamwise is heterogeneous, length is respectively L1、L2、L3, core intersection is h, is oozed Saturating rate is respectively K1、K2、K3, then:
Stratum total length L=L1+L2+L3, Δ P=Δs P1+ΔP2+ΔP3, Q=Q1=Q2=Q3
According to Darcy formula:
Obtaining linear seepage flow mean permeability K is:
Process 2:The calculating of injection production pressure difference
(1) injection production pressure difference when giving bottom pressure is calculated
ΔpMij=pInjMij-pproMij
(2) injection production pressure difference given under note Discharge rate is calculated
In formula:Subscript Inji、projRespectively water injection well i and producing well j ambient parameters, a is drainage radius.
(3) injection production pressure difference given under rice note Discharge rate is calculated
(4) drainage radius a determination
The calculating of oil-water well drainage radius is using infinitely great formation plane radial flow Darcy formula:
By oil and gas well testing, in specific permeability k' and a', calculate specific pressure drop Δ p', then it is pressed than Example is corrected to total pressure drop corresponding during a → ∞, and particular value such as 5000m can be set in this ∞;It is actual according to each mouth well afterwards Situation k, h carries out pressure drop correction, and Δ p, this pressure drop drop in the actual pressure obtained under the actual note Discharge rate of each individual well at present As Δ pMij
Process 3:O point flooding patterns flow field velocity is calculated
O points corresponding i, j well formation vector velocity be:
In formula:I is the corresponding water injection well numbering of O points;J is the corresponding producing well numbering of O points;CijFor the corresponding macroscopic view of O points Curvature correction coefficient;For the corresponding track mean permeability of O points;aijFor distance between water injection well i and producing well j;For Respectively O points are to water injection well i, producing well j distances.
The flow field velocity of O points is vector superposed to be:
Process 4:Flooding pattern flow field velocity calculation procedure
Specific flooding pattern flow field velocity calculation procedure is as follows:
1. according to numerical simulation mesh generation result and well pattern control form, appoint and take a grid element center point O;
2. corresponding water injection well i and producing well j is searched for;
3. trajectory parameters model is utilized, point O and water injection well i track and the track to producing well j is calculated;
4. the average reservoir physical parameter thickness in track, permeability, saturation degree etc. are calculated;
5. Macroscopic Curvature correction coefficient C is calculatedij
6. calculate note and adopt angle influence correction coefficient (cos δMij)n
7. all water injection well is corresponding with O points, producing well j strength of water drive are calculated repeatedly
8. vector calculating is carried out to the vector flow field velocity of all calculating, obtains the flow field velocity of the point;
9. 1.~8. process is repeated, all grid nodes are calculated, obtains considering that static parameter, well pattern form, note are adopted The flooding pattern flow field velocity of state.
Embodiment two
In order to which the application effect of the quick calculation method of the flooding pattern flow field velocity of the offer of embodiment one has more intuitively Understand, the embodiment of the present invention is now illustrated by taking the calculating process using the above method as an example.
Triangle well pattern is adopted for a note two, the heterogeneous design of model is interpolated using the different ratios of permeability between injection-production well After form non-homogeneous model, design parameter Proportionality design is as shown in table 1
The note two of table 1 one adopts the design of triangle well pattern anisotropism
Anisotropism Model 1 Model 2 Model 3
④:③ 400:400 550:250 700:100
Three kinds of models of design are adopted for a note two, calculate respectively remaining oil saturation when moisture content reaches 98%, Strength of water drive, water drive degree dividing elements result are as shown in figure 12, it can be seen that when water injection well and production interwell permeability bigger than normal, After this mouthful of well moisture content reaches 98%, obvious water breakthrough passage is formed between this two mouthfuls of wells, flow field velocity is larger, and remaining oil is in water Cross-hole region is altered substantially without potentiality, REGION OF WATER INJECTION OILFIELD is concentrated mainly on flow field velocity between producing well and between non-water breakthrough well compared with cell Domain.
Particular embodiments described above, has been carried out further in detail to the purpose of the present invention, technical scheme and beneficial effect Describe in detail it is bright, should be understood that the foregoing is only the present invention specific embodiment, the guarantor being not intended to limit the present invention Scope is protected, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc. should be included in this Within the protection domain of invention.

Claims (6)

1. a kind of quick calculation method of flooding pattern flow field velocity, it is characterised in that including:Note adopts streamline approximate calculation side Method, streamline adopt correction of the angle to speed around tomography computational methods, flow curvature to the correction factor C computational methods of speed, note The computational methods of factor M computational methods, flooding pattern flow field velocity.
According to it is described note adopt streamline approximate calculation method, streamline around tomography computational methods set up injection-production well between Flow Field Distribution, utilize The flow curvature adopts angle and the correction factor M computational methods of speed are calculated to the correction factor C computational methods of speed, note to flow The correction factor of field speed degree, flow field velocity calculating is carried out using the computational methods of the flooding pattern flow field velocity.
2. according to the method described in claim 1, it is characterised in that the note adopts streamline approximate calculation method, it is specially:Homogeneous The uniform thickness infinity stratum Isoquant source remittance streamline that any point P is produced in the earth formation is one group of circular arc, by residing for arc track Grid search process calculating speed influences, the rhombus mould that arc model is reduced to be constituted with water injection well -3 points of point P- producing wells Type or straight line model, to shorten the calculation of streamlines time.
3. according to the method described in claim 1, it is characterised in that the streamline is around tomography computational methods, specially:Work as tomography When curtailment adopts packing effect to form note, form injection water and flow into producing well phenomenon around tomography, for any between injection-production well Point P, can be divided into point P and be in tomography homonymy and point P with producing well in two kinds of situations progress streamlines of tomography homonymy with water injection well Calculate, Computing Principle is that the track of discounting of water injection well-point P (tomography end points)-tomography end points (point P)-producing well formation comes approximate Simulate streamline.
4. according to the method described in claim 1, it is characterised in that correction factor C calculating side of the flow curvature to speed Method, be specially:The flow tube unit model of streamline formation can be expressed as the flowing of gradation type bend pipe, by the flowing of gradation type bend pipe again It is secondary to be decomposed into equal-diameter elbow flow process+flaring type flow process+gradual shrinkage flow process, three flow process are calculated respectively to rub Resistance, calculates percolation flow velocity correction factor C.
5. according to the method described in claim 1, it is characterised in that the note adopts correction factor M calculating side of the angle to speed Method, be specially:The size of any point P percolation flow velocities by static parameters such as distance, reservoir properties except being influenceed in reservoir plane Outside, go back receptor site P relevant with the corner dimension of water injection well, the formation of producing well line, when line angle is less than 90 °, can ignore The percolation flow velocity between injection-production well, when angle is more than 90 °, with the increase of angle, speed influence increase, when angle is equal to At 180 °, the point is on main stream line, and speed influence is maximum, and the power side of the angle cosine value thus is adopted into angle as note Speed correction factor M.
6. according to the method described in claim 1, it is characterised in that the computational methods of the flooding pattern flow field velocity, specifically For:By the principle of stacking of the Isoquant multi-source remittance infinity plane radial direction force and velocity of a current, any point P exists corresponding to i-th mouthful of well in oil reservoir The velocity that the point is produced, obtains final flooding pattern flow field speed after the speed calibrated factor C and correction factor M correction Degree, wherein the calculating for pressure difference is divided into three kinds of situations, that is, under the injection production pressure difference, given note Discharge rate when giving bottom pressure Injection production pressure difference computational methods under injection production pressure difference computational methods and given rice note Discharge rate.
CN201710159951.7A 2017-03-17 2017-03-17 Rapid calculation method for flow field velocity of injection-production well pattern Active CN106980758B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710159951.7A CN106980758B (en) 2017-03-17 2017-03-17 Rapid calculation method for flow field velocity of injection-production well pattern

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710159951.7A CN106980758B (en) 2017-03-17 2017-03-17 Rapid calculation method for flow field velocity of injection-production well pattern

Publications (2)

Publication Number Publication Date
CN106980758A true CN106980758A (en) 2017-07-25
CN106980758B CN106980758B (en) 2020-06-26

Family

ID=59338274

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710159951.7A Active CN106980758B (en) 2017-03-17 2017-03-17 Rapid calculation method for flow field velocity of injection-production well pattern

Country Status (1)

Country Link
CN (1) CN106980758B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108049849A (en) * 2017-09-07 2018-05-18 中国石油化工股份有限公司 Water-drive pool Plane Fluid Field regulates and controls design method
CN108868745A (en) * 2018-07-09 2018-11-23 中国石油大学(华东) A kind of oil reservoir flow field matching evaluation method
CN109933899A (en) * 2019-03-13 2019-06-25 中国石油大学(华东) A kind of water drive old filed cyclic waterflooding development cycle calculates plate and determines method
CN110826171A (en) * 2019-09-18 2020-02-21 深圳先进技术研究院 Method and device for simulating external flow field of fan and computer readable storage medium
CN116341406A (en) * 2023-03-03 2023-06-27 常州大学 Calculation method for determining injection and production interference time of inclined oil reservoir

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1224848A (en) * 1998-08-03 1999-08-04 辽河石油勘探局钻采工艺研究院 Inter-well tracing determination technology using chemical tracer
CN101661514A (en) * 2008-05-21 2010-03-03 中国石化股份胜利油田分公司地质科学研究院 Oil deposit black oil model numerical simulation system
CN103970989A (en) * 2014-04-15 2014-08-06 昆明理工大学 Immersing boundary flow field calculation method based on fluid/solid interface consistency
CN105822271A (en) * 2015-09-09 2016-08-03 中国石油化工股份有限公司 Adjustment method of variable flow line vector of thick-stratum oil deposit

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1224848A (en) * 1998-08-03 1999-08-04 辽河石油勘探局钻采工艺研究院 Inter-well tracing determination technology using chemical tracer
CN101661514A (en) * 2008-05-21 2010-03-03 中国石化股份胜利油田分公司地质科学研究院 Oil deposit black oil model numerical simulation system
CN103970989A (en) * 2014-04-15 2014-08-06 昆明理工大学 Immersing boundary flow field calculation method based on fluid/solid interface consistency
CN105822271A (en) * 2015-09-09 2016-08-03 中国石油化工股份有限公司 Adjustment method of variable flow line vector of thick-stratum oil deposit

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
POLLOCK D W. ET AL: "Semi-analytical computation of path lines for finite-difference models", 《GROUND WATER》 *
THIELE M.R.ET AL: "Water Injection Optimization Using a Streamline-Based Workflow", 《SOCIETY OF PETROLEUM ENGINEERS》 *
喻鹏: "静观2区块高凝油油藏流动单元预测表征及剩余油分布研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *
郭文敏等: "二维可视化水驱物理模拟井间参数解释方法", 《石油天然气学报(江汉石油学院学报)》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108049849A (en) * 2017-09-07 2018-05-18 中国石油化工股份有限公司 Water-drive pool Plane Fluid Field regulates and controls design method
CN108049849B (en) * 2017-09-07 2019-11-29 中国石油化工股份有限公司 Water-drive pool Plane Fluid Field regulates and controls design method
CN108868745A (en) * 2018-07-09 2018-11-23 中国石油大学(华东) A kind of oil reservoir flow field matching evaluation method
CN108868745B (en) * 2018-07-09 2021-10-01 中国石油大学(华东) Oil reservoir flow field matching evaluation method
CN109933899A (en) * 2019-03-13 2019-06-25 中国石油大学(华东) A kind of water drive old filed cyclic waterflooding development cycle calculates plate and determines method
CN109933899B (en) * 2019-03-13 2023-04-25 中国石油大学(华东) Method for determining water flooding old oilfield period water flooding development period calculation chart
CN110826171A (en) * 2019-09-18 2020-02-21 深圳先进技术研究院 Method and device for simulating external flow field of fan and computer readable storage medium
CN116341406A (en) * 2023-03-03 2023-06-27 常州大学 Calculation method for determining injection and production interference time of inclined oil reservoir
CN116341406B (en) * 2023-03-03 2024-04-23 常州大学 Calculation method for determining injection and production interference time of inclined oil reservoir

Also Published As

Publication number Publication date
CN106980758B (en) 2020-06-26

Similar Documents

Publication Publication Date Title
CN106980758A (en) A kind of quick calculation method of flooding pattern flow field velocity
CN104594872B (en) A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy
CN109948272A (en) Tune based on inter well connectivity blocks up dynamic prediction method and system
CN103886187B (en) A kind of urban river water sand real-time predicting method based on data assimilation
Alazba et al. Field assessment of friction head loss and friction correction factor equations
CN109284531A (en) An a kind of two-dimentional hydrodynamics coupling process based on space topological
CN102087130B (en) Computational fluid dynamics (CFD) technology-based method for optimizing acoustic path of multi-path ultrasonic flowmeter assembled in elbow pipe
WO2021253902A1 (en) Method for calculating head loss coefficient and branch pipe flow distribution of multi-point water intake port
CN104915530B (en) The method for building up of oil reservoir interwell communication relation
CN109271672A (en) A kind of river-lake-pumping plant influences each other the river water surface curve calculation method under acting on
CN102630426B (en) Structure optimization method for Venturi fertilizer injector
CN106446438B (en) A kind of nearly far field approach of coupled numerical simulation of power plant&#39;s warm water discharge deep water discharge
Josiah et al. An experimental investigation of head loss through trash racks in conveyance systems
CN104091065A (en) Intermittent flow numerical simulation method for solving shallow water problem
Ghare et al. Cylindrical central baffle flume for flow measurements in open channels
Li et al. Physical and numerical modeling of the hydraulic characteristics of type-A piano key weirs
Khanahmadi et al. Hydraulic of curved type-B piano key weirs characteristics under free flow conditions
CN201707087U (en) Vortex street mass flow meter
Rai et al. Analysis of water distribution network using EPANET
Vatankhah et al. Discussion of “New method for modeling thin-walled orifice flow under partially submerged conditions” by David Brandes and William T. Barlow
CN105444828B (en) The design method of the rectifier of ultrasonic flowmeter
CN108534846A (en) A kind of drawing type portable weir body and flow-measuring method
CN108506622A (en) A kind of lower resistance threeway component based on arc flow deflector
CN110991084B (en) Reservoir permeability calculation method based on streamline numerical value well test
Rihan et al. Physical simulation for the flow in straight and rectangular loop manifolds

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant