CN107939366B

CN107939366B - A kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method and device

Info

Publication number: CN107939366B
Application number: CN201711182968.0A
Authority: CN
Inventors: 田冷; 高星星; 顾岱鸿; 董俊林; 朱贤峰
Original assignee: China University of Petroleum Beijing
Current assignee: China University of Petroleum Beijing
Priority date: 2017-11-23
Filing date: 2017-11-23
Publication date: 2019-10-11
Anticipated expiration: 2037-11-23
Also published as: CN107939366A

Abstract

The application provides a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method and device.The described method includes: constructing the corresponding double medium model of well section to be analyzed using fractal theory；According to the double medium model, the Pressure Dynamic Model of the well section to be analyzed is established using fractional calculus method；According to the Pressure Dynamic Model, pressure dynamic analysis is carried out to the well section to be analyzed.Utilize embodiment each in the application, the Pressure Dynamic Model of well section to be analyzed is established using fractal theory, fractional calculus method, and it can use semi analytical method and solve Pressure Dynamic Model, the solution of Pressure Dynamic Model can be quickly and accurately got, provides accurate data basis for the subsequent pressure dynamic analysis for carrying out well section to be analyzed.Precision of analysis is improved while reducing the calculation amount of Dynamic Analysis of Horizontal Well Pressure.

Description

A kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method and device

Technical field

The application belongs to oil/gas well Well Test Technology field more particularly to a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure side Method and device.

Background technique

Tight gas is as one of three big Unconventional gas, and stock number is abundant, and potentiality to be exploited is big, will become China's natural gas The valuable source of industrial fast and stable development.Tight gas is generally found in the densification under the specific conditions such as low hole, hypotonic and low pressure In reservoir, gas well natural production ability is low, does not often have commercial mining only with the exploitation of the fracturing technique of horizontal well and routine Value, it is therefore necessary to which research uses novel massive hydraulic fracture renovation technique.Volume fracturing technology is formed because of its pressure break with major fracture For " the netted seam " in length and breadth of trunk, seepage area and reservoir permeability are increased, improve tight gas reservoir initil output and is finally adopted Yield, it has also become the necessary reservoir reconstruction means of compact reservoir Efficient Development.

In the prior art, more and more using volume fracturing technology exploitation dense gas, using volume fracturing technology In horizontal well, Pressure behaviour variation of the research level well in recovery process is needed to improve water to ensure the safety of recovery process Flat well yield.The method for carrying out Dynamic Analysis of Horizontal Well Pressure in the prior art is mainly: in conjunction with numerical simulation technology to body The overstocked Production development for splitting horizontal well is simulated.But this method is computationally intensive, and needs deep theoretical grounding in basic skills It is realized substantially with the ability of a large amount of practical experience.And this method needs to obtain the crack of accurate volume fracturing horizontal well Distributed intelligence, still, in the prior art, the micro-seismic monitoring information of generally up to resolution ratio is also difficult to provide high-precision crack System distributed intelligence will lead to the pressure dynamic analysis inaccuracy of volume fracturing horizontal well in this way.Therefore, volume fracturing water is reduced The calculation amount of horizontal well Pressure behaviour mutation analysis further increases the accuracy of pressure dynamic analysis, skill referred to as urgently to be resolved Art problem.

Summary of the invention

The application is designed to provide a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method and device, is managed using a point shape By the Pressure Dynamic Model for establishing volume fracturing horizontal well with fractional calculus method, volume fracturing horizontal well pressure dynamic is reduced While the calculation amount of model foundation, the accuracy of volume fracturing horizontal well pressure Dynamic Model is improved, is further increased The accuracy of volume fracturing Dynamic Analysis of Horizontal Well Pressure.

On the one hand this application provides a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure methods, comprising:

The corresponding double medium model of well section to be analyzed is constructed using fractal theory；

According to the double medium model, the Pressure behaviour mould of the well section to be analyzed is established using fractional calculus method Type；

According to the Pressure Dynamic Model, pressure dynamic analysis is carried out to the well section to be analyzed.

Further, described to construct the well section to be analyzed using fractal theory in another embodiment of the method Corresponding double medium model, comprising:

It is using fractal theory using the corresponding microcrack network of the well section to be analyzed as fractal, the fractal is embedding Enter into Euclid's matrix；

Euclid's matrix comprising the fractal is divided into inner region and outskirt, constructs the double medium model；

The inner region includes waterpower major fracture, Euclid's matrix, point shape microcrack, and the outskirt includes Euclid's base Matter and divide shape intrinsic fracture, described point of shape microcrack includes the intrinsic fracture of induced fractures and activation, described to divide shape intrinsic fracture Including original intrinsic fracture.

Further, in another embodiment of the method, the inner region is round radial compound region, the inner region In waterpower major fracture be vertical lap seam.

Further, described according to the double medium model in another embodiment of the method, utilize fractional order The method fluxions establish the Pressure Dynamic Model of the well section to be analyzed, comprising:

According to the double medium model, inner region fractal number rank is established respectively using the fractional calculus method and is spread Model and outskirt fractal number rank diffusion model；

Using the combination of the inner region fractal number rank diffusion model and the outskirt fractal number rank diffusion model as institute State Pressure Dynamic Model.

Further, in another embodiment of the method, the inner region fractal number rank diffusion model includes:

In above formula, r_DIndicate dimensionless radius, θ₁Indicate unusual diffusion index, ψ_Df1Indicate the Fracture System of the inner region without Dimensionless pressure, β₁=d_mf1-θ₁- 1, d_mf1Indicate the mass-fractal dimension of the fracture network of the inner region,M₁₂Table Show the crack mobility ratio between the inner region and outskirt；W₁₂Indicate that ratio, t are held in the storage between the inner region and outskirt_DIndicate it is no because The secondary time,φ_fr1Indicate porosity of the crack of the inner region under reference length, C_tf1Table Show the Fracture System system compressibility of the inner region, φ_m1Indicate the porosity of the matrix system of the inner region, C_tm1Indicate institute State the matrix system system compressibility of inner region, γ₁=2/ (2+ θ₁), ψ_Dm1Indicate the matrix system zero dimension pressure of the inner region Power, λ₁Indicate the interporosity flow coefficient of the inner region；

The outskirt fractal number rank diffusion model includes:

In above formula, r_DIndicate dimensionless radius, θ₂Indicate unusual diffusion index, ψ_Df2Indicate the Fracture System of the outskirt without Dimensionless pressure, β₂=d_mf2-θ₂- 1, d_mf2Indicate the mass-fractal dimension of the fracture network of the outskirt, t_DWhen indicating zero dimension Between,φ_fr2Indicate porosity of the crack of the outskirt under reference length, C_tf2It indicates The Fracture System system compressibility of the outskirt, φ_m2Indicate the porosity of the matrix system of the outskirt, C_tm2Described in expression The matrix system system compressibility of outskirt, γ₂=2/ (2+ θ₂), ψ_Dm2Indicate the matrix system zero dimension pressure of the outskirt, λ₂Indicate the interporosity flow coefficient of the outskirt.

Further, described according to the Pressure Dynamic Model in another embodiment of the method, to described wait divide It analyses well section and carries out pressure dynamic analysis, comprising:

Solve the Pressure Dynamic Model using semi analytical method, according to the solution of the Pressure Dynamic Model, draw it is described to Analyze the Pressure behaviour curve of well section；

According to the Pressure behaviour curve, pressure dynamic analysis is carried out to the well section to be analyzed.

Further, described to solve the Pressure behaviour mould using semi analytical method in another embodiment of the method Type, comprising:

Using Laplace transformation, the Pressure Dynamic Model is transformed to Laplace space Pressure Dynamic Model；

The Laplace space pressure is solved using line source function, discrete method, principle of stacking and Stehfest Numerical inverse method Power dynamic model obtains the solution of the Pressure Dynamic Model.

Further, described according to the Pressure behaviour curve in another embodiment of the method, to described wait divide It analyses well section and carries out pressure dynamic analysis, comprising:

By the Pressure behaviour curve carry out flowing stage division, and according to the Pressure behaviour curve to it is described to point It analyses well section pressure and carries out parameters sensitivity analysis.

On the other hand, this application provides a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure device, described device includes:

Dielectric model establishes module, for constructing the corresponding double medium model of well section to be analyzed using fractal theory；

Pressure Dynamic Model establishes module, for being established using fractional calculus method according to the double medium model The Pressure Dynamic Model of the well section to be analyzed；

Pressure dynamic analysis module, for it is dynamic to carry out pressure to the well section to be analyzed according to the Pressure Dynamic Model State analysis.

In another aspect, present invention also provides a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure devices, comprising: processor And the memory for storage processor executable instruction, the processor realize above-mentioned volume fracturing when executing described instruction Dynamic Analysis of Horizontal Well Pressure method.

Volume fracturing Dynamic Analysis of Horizontal Well Pressure method and device provided by the present application is established using fractal theory wait divide The corresponding double medium model of well section is analysed, the structure of well section microcrack network to be analyzed and answering for distribution can be accurately described Polygamy, scrambling and heterogeneity provide accurately theoretical base for the subsequent pressure dynamic analysis for carrying out well section to be analyzed Plinth improves the accuracy of volume fracturing Dynamic Analysis of Horizontal Well Pressure.Meanwhile it can be with part wait divide using fractal theory The feature for analysing the microcrack network of well section, states the feature of the microcrack network of full well section, reduces and subsequent establishes well to be analyzed The calculation amount of section Pressure Dynamic Model, further reduces the calculation amount of volume fracturing Dynamic Analysis of Horizontal Well Pressure.In addition, this Application can establish the Pressure Dynamic Model of volume fracturing horizontal well using fractional calculus method, and can use semi analytical method Solve Pressure Dynamic Model.The accuracy of Pressure Dynamic Model building can be improved in fractional calculus method, and semi analytical method has The result of analytic method is accurately and numerical method handles the more flexible advantage in boundary.Therefore, volume fracturing water provided by the present application Horizontal well pressure dynamic analysis method and device can improve analysis while reducing the calculation amount of Dynamic Analysis of Horizontal Well Pressure As a result accuracy.

Detailed description of the invention

In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property Under, it is also possible to obtain other drawings based on these drawings.

Fig. 1 is that a kind of method flow of volume fracturing Dynamic Analysis of Horizontal Well Pressure embodiment of the method provided by the present application shows It is intended to；

Fig. 2 is the structural schematic diagram of volume fracturing horizontal well double medium model in a kind of embodiment of the application；

Fig. 3 is the structural schematic diagram of volume fracturing horizontal well double medium model in the application another kind embodiment；

Fig. 4 is the line source top view of volume fracturing horizontal well double medium model in the application one embodiment；

Fig. 5 is the line source front view of volume fracturing horizontal well double medium model in the application one embodiment；

Fig. 6 is the contrast schematic diagram of Pressure Dynamic Model in one embodiment provided by the present application；

Fig. 7 is Pressure behaviour curve synoptic diagram in one embodiment provided by the present application；

Fig. 8 is that the mass-fractal dimension of the fracture network of the inner region in the application one embodiment influences Pressure behaviour Schematic diagram；

Fig. 9 is that the mass-fractal dimension of the fracture network of the outskirt in the application one embodiment influences Pressure behaviour Schematic diagram；

Figure 10 is the schematic diagram that the inner region abnormality diffusion index in the application one embodiment influences Pressure behaviour；

Figure 11 is the schematic diagram that the outskirt abnormality diffusion index in the application one embodiment influences Pressure behaviour；

Figure 12 is the schematic diagram that the inside and outside compound radius in the application one embodiment influences Pressure behaviour；

Figure 13 is the influence schematic diagram that interior outskirt crack mobility ratio changes Pressure behaviour in the application one embodiment；

Figure 14 is that the schematic diagram for comparing Pressure behaviour influence is held in interior outskirt storage in the application one embodiment；

Figure 15 is the module of the volume fracturing Dynamic Analysis of Horizontal Well Pressure device in a kind of embodiment provided by the present application Structural schematic diagram；

Figure 16 is the module knot of another volume fracturing Dynamic Analysis of Horizontal Well Pressure Installation practice provided by the present application Structure schematic diagram.

Specific embodiment

In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common The application protection all should belong in technical staff's every other embodiment obtained without creative efforts Range.

Formula in each embodiment provided herein is only merely a kind of embodiment, can also be as needed to each A formula is deformed or is converted, and the application is not especially limited.

Fig. 1 is a kind of a kind of method stream of embodiment of volume fracturing Dynamic Analysis of Horizontal Well Pressure method provided by the present application Journey schematic diagram, as shown in Figure 1, volume fracturing Dynamic Analysis of Horizontal Well Pressure method provided by the present application includes:

S1, the corresponding double medium model of well section to be analyzed is constructed using fractal theory.

Fractal theory can refer to that part is known as a point shape with whole body similar in some way, and introducing fractal theory can be with The structure of the microcrack network of well section to be analyzed and the complexity of distribution, irregular shape and expense mean value etc. are described.Specifically, exist In the application one embodiment, it can use the structure progress Fractal process that fractal theory is analysed to well section, construct to be analyzed The corresponding double medium model of well section.Such as: the portion with representing characteristic can be chosen according to the practical structures of well section to be analyzed The architectural characteristic of the entire well section to be analyzed of the representation of point well section, can also utilize geometry similar with well section structure to be analyzed The structure of shape representation well section to be analyzed constructs the corresponding double medium model of well section to be analyzed.

It is described to construct the corresponding dual media of the well section to be analyzed using fractal theory in the application one embodiment Model may include:

Specifically, Euclid's matrix includes the matrix of conventional geometric dimension, and dimension includes one-dimensional, two dimension, three-dimensional, this Shen Please in Euclid's matrix can be two-dimensional matrix.In the application one embodiment, can use fractal theory will be wait divide The corresponding microcrack network of well section is analysed as fractal, the corresponding fractal of microcrack network (is referred to as a point shape fine fisssure Seam network) it is embedded into Euclid's matrix.Euclid's matrix comprising point shape microcrack is divided into inner region (and to be properly termed as The area SRV, stimulated reservoir volume, reservoir reconstruction volume area) and outskirt (and it is properly termed as the area USRV, Unstimulated reservoir volume, volume area is not transformed), construct double medium model.Inner region includes waterpower master Crack, Euclid's matrix and divide shape microcrack, outskirt includes Euclid's matrix and point shape intrinsic fracture.Dual media includes Water-bearing media with crack and hole two media property, double medium model can accurately describe fluid in well to be analyzed Flow behavior in the corresponding matrix of section.

In the application one embodiment, fractal theory is quoted, simulates the knot of microcrack network in volume fracturing horizontal well Structure, the microcrack network being analysed in well section are divided into waterpower major fracture (major fracture generated when volume fracturing), point shape fine fisssure It stitches, point shape intrinsic fracture.Point shape microcrack includes the day of induced fractures (crack generated is induced during volume fracturing), activation Right crack (intrinsic fracture started during volume fracturing), a point shape intrinsic fracture includes original intrinsic fracture.Fig. 2 is this Shen Please in a kind of embodiment volume fracturing horizontal well double medium model structural schematic diagram, as shown in Fig. 2, double medium model Inner region include: waterpower major fracture, induced fractures, activation intrinsic fracture, remaining blank space is Euclid's matrix, and outskirt includes Euclid's matrix of original intrinsic fracture and blank space.

In another embodiment of the application, the inner region is round radial compound region, the waterpower in the inner region Major fracture is vertical lap seam, can simplify double medium model in this way, reduces subsequent calculation amount.Fig. 3 is the application another kind The structural schematic diagram of volume fracturing horizontal well double medium model in embodiment, as shown in figure 3, by inner region in double medium model It is reduced to a round radial compound region.And it is possible to the induced fractures in inner region and the intrinsic fracture of activation be merged, system One, as shape microcrack is divided, can also be reduced to the waterpower major fracture in inner region the vertical lap seam of rule compositor.Fig. 4 is the application The line source top view of volume fracturing horizontal well double medium model in one embodiment, Fig. 5 are bodies in the application one embodiment Overstock the line source front view for splitting horizontal well double medium model, as shown in Figure 4 and Figure 5, r_mIt can indicate interior outskirt compound radius, r_eIt can indicate outer boundary radius, R_sourceIt can indicate that the inner boundary radius of line source, ε can indicate that mathematical limit symbol, h can To indicate the effective thickness on stratum.Assuming that r_eLevel off to infinity, R_sourceLevel off to 0.

S2, according to the double medium model, the pressure for establishing the well section to be analyzed using fractional calculus method is dynamic States model.

Specifically, fractional calculus is a branch of calculus, including carries out fractional order differential integral to function, can To be used to describe the dynamic behavior of physical system and material.Pressure Dynamic Model can be used to describe well section to be analyzed working Pressure behaviour variation in the process.It, can be according to the double of foundation after the completion of the corresponding double medium model of well section to be analyzed is established Dense media model establishes the Pressure Dynamic Model of well section to be analyzed using fractional calculus method.In the application one embodiment In, the Pressure Dynamic Model of well section to be analyzed can be using fractional calculus method, according to fluid in well section to be analyzed The diffusion equation that diffusion property is established.The application can establish point of well section to be analyzed based on fractal theory and fractional calculus Shape Fractional Diffusion Equation, using the fractal number rank diffusion equation of foundation as the Pressure Dynamic Model of well section to be analyzed.

S3, according to the Pressure Dynamic Model, pressure dynamic analysis is carried out to the well section to be analyzed.

Specifically, after establishing the Pressure Dynamic Model of well section to be analyzed, according to the Pressure Dynamic Model of foundation, to be analyzed Well section carries out dynamic analysis.Such as: the solution of Pressure Dynamic Model can be obtained by solving the Pressure Dynamic Model established.Root According to the solution of Pressure Dynamic Model, the Pressure behaviour curve of well section to be analyzed is drawn, the pressure dynamic analysis of well section to be analyzed is completed.

It is described according to the Pressure Dynamic Model in the application one embodiment, the well section to be analyzed is pressed Power dynamic analysis may include:

Specifically, semi analytical method can be a kind of method for solving including a variety of analytic methods and numerical method, have analytic method Result accurately and numerical method handles the more flexible advantage in boundary.Therefore, while semi analytical method can reduce calculation amount Improve the precision of calculated result.After the application obtains the Pressure Dynamic Model of well section to be analyzed, pressure is solved using semi analytical method Dynamic model obtains the solution of Pressure Dynamic Model.The pressure that well section to be analyzed can be drawn according to the solution of Pressure Dynamic Model is dynamic State curve can carry out the pressure dynamic analysis of well section to be analyzed according to the Pressure behaviour curve drawn out.

It is described according to the Pressure behaviour curve in one embodiment of the application, the well section to be analyzed is carried out Pressure dynamic analysis may include: the division that the Pressure behaviour curve is carried out to flowing stage, and according to the Pressure behaviour Curve carries out parameters sensitivity analysis to the well section pressure to be analyzed.The Pressure behaviour curve of acquisition can be subjected to flowing rank The division of section analyzes the Pressure behaviour variation in the well section course of work to be analyzed according to the flowing stage of division.It is also possible that root Parameters sensitivity analysis is carried out to well section to be analyzed according to the Pressure behaviour curve of well section to be analyzed, analyzes parameters to be analyzed The influence of pressure change in well section, to ensure work safety, is improved the standard to adjust parameters during oil-gas mining Well yield.

It is described to solve the Pressure Dynamic Model using semi analytical method in the application one embodiment, may include:

Specifically, after the Pressure Dynamic Model for getting well section to be analyzed, Laplace transformation can be used, pressure is moved States model is transformed to Laplace space Pressure Dynamic Model (Pressure Dynamic Model i.e. in Laplace space).Obtain Pressure Dynamic Model Laplace space Pressure Dynamic Model after, can use line source function, discrete method, principle of stacking and Stehfest numerical inversion Method solves Laplace space Pressure Dynamic Model, obtains the solution of Pressure Dynamic Model.

Volume fracturing Dynamic Analysis of Horizontal Well Pressure method provided by the present application establishes well section to be analyzed using fractal theory Corresponding double medium model, can accurately describe the structure of well section microcrack network to be analyzed and the complexity of distribution, Scrambling and heterogeneity.Accurate theoretical basis is provided for the subsequent pressure dynamic analysis for carrying out well section to be analyzed, is mentioned The high accuracy of volume fracturing Dynamic Analysis of Horizontal Well Pressure.Meanwhile the well to be analyzed of part can be used using fractal theory The feature of the microcrack network of section, states the feature of the microcrack network of full well section.Reduce and subsequent establishes well section pressure to be analyzed The calculation amount of power dynamic model further reduces the calculation amount of volume fracturing Dynamic Analysis of Horizontal Well Pressure.In addition, the application The Pressure Dynamic Model of volume fracturing horizontal well can be established using fractional calculus method, and can be solved using semi analytical method Pressure Dynamic Model can also improve the accurate of analysis result while reducing the calculation amount of Dynamic Analysis of Horizontal Well Pressure Property.

On the basis of the above embodiments, described according to the double medium model, it is established using fractional calculus method The Pressure Dynamic Model of the well section to be analyzed may include:

Specifically, in one embodiment of the application, according to the double medium model of foundation, fractional calculus is utilized Method establishes outskirt in the corresponding inner region fractal number rank diffusion model of inner region in double medium model and double medium model respectively Corresponding outskirt fractal number rank diffusion model.It can be corresponding outer by the inner region fractal number rank diffusion model and outskirt of inner region Distinguish Pressure Dynamic Model of the combination of shape fractional order diffusion model as well section to be analyzed.

In the application one embodiment, inner region fractal number rank diffusion model may include:

Outskirt fractal number rank diffusion model may include:

The primary condition of inner region fractal number rank diffusion model and outskirt fractal number rank diffusion model, boundary condition and Interface condition can be expressed as:

Zero dimension primary condition: ψ_Df1(r, 0)=ψ_Df2(r, 0)=ψ_Dm1(r, 0)=ψ_Dm2(r, 0)=0

Zero dimension internal boundary condition:

Zero dimension Outer Boundary Conditions:

Zero dimension interface condition: In the above formulas, r_DIt can indicate dimensionless radius；θ₁It can indicate the unusual diffusion index of inner region, θ₂It can indicate outskirt Unusual diffusion index, unusual diffusion index are also referred to as abnormal conduction index, flow road along complicated Disordered Media to describe fluid The feature of diameter diffusion, represents the unusual degree of diffusion process deviation from the norm diffusion；ψ_Df1Can indicate the Fracture System of inner region without Dimensionless pressure；t_DIndicate non dimensional time；λ₁It can indicate the interporosity flow coefficient of inner region；ψ_Dm1Can indicate the matrix system of inner region without Dimensionless pressure；ψ_Dm2It can indicate the matrix system zero dimension pressure of outskirt；ψ_Df2It can indicate the Fracture System zero dimension of outskirt Pressure；λ₂It can indicate the interporosity flow coefficient of outskirt；For the yield of zero dimension line source, r_mDIt can indicate outskirt in zero dimension Compound radius, i.e. the outer boundary radius of zero dimension inner region.

In the above formulas, each parameter can be calculated using following formula:

In above formula, r can indicate the radial coordinate of arbitrary point in double medium model；L_refIt can indicate reference length；ψ_D It can indicate zero dimension pseudopressure；H can indicate Effective thickness of formation, and unit can be m；T_scIt can indicate under the status of criterion Formation temperature, unit can be K；q_scIt can indicate that the yield under the status of criterion, unit can be m³/s；p_scIt can indicate standard Strata pressure under situation, unit can be MPa；ψ_iIt can indicate zero dimension prime stratum pseudopressure；ψ can indicate quasi- pressure Power；d_mf1It can indicate the mass-fractal dimension of the fracture network of inner region；d_mf2It can indicate the quality point of the fracture network of outskirt Shape dimension；M₁₂It can indicate the crack mobility ratio between inner region and outskirt；W₁₂It can indicate that the storage between inner region and outskirt is held Than；k_fr1It can indicate permeability of the crack of inner region under reference length；k_fr2The crack of outskirt can be indicated in reference length Under permeability；φ_fr1It can indicate porosity of the crack of inner region under reference length；φ_fr2It can indicate the crack of outskirt Porosity under reference length；C_tf1It can indicate the Fracture System system compressibility of inner region；C_tf2It can indicate outskirt Fracture System system compressibility；φ_m1It can indicate the porosity of the matrix system of inner region；φ_m2It can indicate the matrix of outskirt The porosity of system；C_tm1It can indicate the matrix system system compressibility of inner region；C_tm2It can indicate the matrix system of outskirt System compressibility；μ can indicate gas viscosity, and unit can be mPas；α can indicate that form factor, unit can be m^-2。

In above formula, undefined parameter such as: β₁、β₂、η₁₂、ω₁、ω₂、γ₁、γ₂Deng being arranged for simplified model Reduced parameter, without specific meaning.

In the application one embodiment, above-mentioned inner region fractal number rank diffusion model and outskirt fractal number rank are being established Before diffusion model, point shape permeability and the porosity of microcrack net in double medium model can be first calculated, can specifically be adopted It is calculated and is obtained with following formula:

In above formula, l=1,2, when can respectively indicate inner region and outskirt, i.e. l=1, it can indicate inner region, it, can when l=2 To indicate outskirt；k_frlAnd φ_frlIt can indicate the permeability and porosity under reference length；L_refIt can indicate reference length；d It can indicate Euclidean dimension, for cylindrical system, value is equal to 2；d_mf1It can indicate the mass fractal dimension of fracture network Number, represents a point complexity for shape fracture network structure and distribution；θ_lIt can indicate unusual diffusion index, it is also referred to as abnormal to pass Index is led, to describe fluid along the feature of complicated Disordered Media flow path diffusion, represents the expansion of diffusion process deviation from the norm Scattered unusual degree；R is the radial coordinate of arbitrary point in double medium model.

Dynamic Analysis of Horizontal Well Pressure method provided by the present application utilizes parting theory and fractional calculus method, building The Pressure Dynamic Model of volume fracturing horizontal well can accurately describe the structure and distribution of well section microcrack network to be analyzed Complexity, scrambling and heterogeneity.Analysis knot is improved while reducing the calculation amount of Dynamic Analysis of Horizontal Well Pressure The accuracy of fruit.

Obtain above-mentioned inner region fractal number rank diffusion model and outskirt fractal number rank diffusion model well section i.e. to be analyzed After Pressure Dynamic Model, using semi analytical method, solves above-mentioned inner region fractal number rank diffusion model and outskirt fractal number rank is expanded Dissipate model.Specifically (Laplace transform) can be converted by Laplace, Pressure Dynamic Model is transformed to Laplace space pressure Dynamic model obtains the fractal number rank diffusion model of Laplace space.In the application one embodiment, pass through Laplace The fractal number rank diffusion model, that is, Laplace space Pressure Dynamic Model for converting the Laplace space obtained may include inner region Laplace Spatial Fractal fractional order diffusion model and outskirt Laplace space fractal number rank diffusion model, are specifically as follows following expression sides Formula:

Inner region Laplace space fractal number rank diffusion model:

Outskirt Laplace space fractal number rank diffusion model:

The boundary condition and interface condition of Laplace space Pressure Dynamic Model are as follows:

Internal boundary condition:

Outer Boundary Conditions:

Interface condition:

In the above formulas, r_DIt can indicate dimensionless radius；θ₁It can indicate the unusual diffusion index of inner region, θ₂It can be with table Show that the unusual diffusion index of outskirt, unusual diffusion index are also referred to as abnormal conduction index, to describe fluid along complicated unordered Jie The feature of matter flow path diffusion, represents the unusual degree of diffusion process deviation from the norm diffusion；It can indicate Laplace sky Between lower inner region Fracture System zero dimension pressure；t_DIndicate non dimensional time；λ₁It can indicate the interporosity flow coefficient of inner region；μ can be with table Show gas viscosity, unit can be mPas；It can indicate the matrix system zero dimension pressure of inner region under Laplace space；It can indicate the matrix system zero dimension pressure of outskirt；It can indicate the Fracture System zero dimension pressure of outskirt； λ₂It can indicate the interporosity flow coefficient of outskirt；For the yield of zero dimension line source, r_mDIt can indicate compound half of outskirt in zero dimension Diameter, i.e. the outer boundary radius of zero dimension inner region.

In the above formulas, undefined parameter is such as: β₁、β₂、η₁₂、ω₁、ω₂、γ₁、γ₂Deng, be for simplified model and The reduced parameter of setting, without specific meaning.The circular of parameters can be consistent with the various embodiments described above, this Place repeats no more.

Obtain above-mentioned Laplace space Pressure Dynamic Model i.e.: inner region Laplace space fractal number rank diffusion model and outskirt are drawn It, can be anti-by line source function, discrete method, superposition principle and Stehfest numerical value after family name's Spatial Fractal fractional order diffusion model It drills method etc. to be solved, the solution of the mystery gas reservoir volume fracturing horizontal well pressure dynamic model of acquisition can be expressed as shape Formula:

The general solution of inner region fractal number rank diffusion model:

The general solution of outskirt fractal number rank diffusion model:

Line-source solution:

In above formula, each parameter can be calculated using following formula:

In the above formulas, r_DIt can indicate dimensionless radius；θ₁It can indicate the unusual diffusion index of inner region, θ₂It can be with table Show that the unusual diffusion index of outskirt, unusual diffusion index are also referred to as abnormal conduction index, to describe fluid along complicated unordered Jie The feature of matter flow path diffusion, represents the unusual degree of diffusion process deviation from the norm diffusion；It can indicate inner region Fracture System zero dimension pressure；t_DIndicate non dimensional time；λ₁It can indicate the interporosity flow coefficient of inner region；μ can indicate that gas is viscous Degree, unit can be mPas；It can indicate the matrix system zero dimension pressure of inner region；It can indicate outskirt Matrix system zero dimension pressure；It can indicate the Fracture System zero dimension pressure of outskirt；λ₂It can indicate altering for outskirt Flow coefficient；For the yield of zero dimension line source, r_mDIt can indicate the compound radius of outskirt in zero dimension, i.e. zero dimension inner region Outer boundary radius；I_v1、K_v1、I_v2、K_v2Belong to Bessel function special symbol, v1, v2 can indicate the order of Bessel function.

In the above formulas, undefined parameter is such as: β₁、β₂、η₁₂、ω₁、ω₂、γ₁、γ₂、A₁、B₁、f₁(u)、g₁(r_D, u)、A₂、B₂、f₂(u)、g₂(r_D, u) etc., it is the reduced parameter being arranged for simplified model, without specific meaning.

By above-mentioned each embodiment, the pressure for establishing well section to be analyzed using fractal theory, fractional calculus method is dynamic States model, and can use semi analytical method and solve Pressure Dynamic Model, it can quickly and accurately get Pressure behaviour mould The solution of type provides accurate data basis for the subsequent pressure dynamic analysis for carrying out well section to be analyzed.

Fig. 6 is the contrast schematic diagram of Pressure Dynamic Model in one embodiment provided by the present application, as shown in fig. 6, at this Apply in one embodiment, the Pressure Dynamic Model for the well section to be analyzed that the application establishes is moved with a kind of pressure in the prior art States model (discuss Subcommittee-to and study carefully in article by the multiple dimensioned unstable Mongolian stream reason of the shale gas reservoir fractured well based on complicated percolation mechanism of Zhao Yulong Pressure Dynamic Model) compare.As shown in fig. 6, the solid line of L1 and the solid line of L2 are the application d_mf1=d_mf2=2 and θ₁= θ₂Fractal number rank diffusion model (Pressure Dynamic Model) when=0 as a result, result with Zhao model in the prior art (as shown in circular lines in L1 and L2) are completely coincident.Therefore, what the application established divides shape tight gas reservoir volume fracturing horizontal well The corresponding Pressure Dynamic Model of compound dual media is correct.

In Fig. 6, parameters and the value of parameter when generating figure may include: C_D=10 can indicate zero dimension well Cylinder reserve coefficient, S=0.01 can indicate skin factor, r_m=1500m can indicate interior outskirt compound radius, and L=1000m can To indicate horizontal well length, L_F=80m can indicate hydraulic fracture length, and δ=90 ° can indicate crack and pit shaft angle, M= 2 can indicate that man-made fracture item number, N=5 can indicate that crack is divided into 2N equal portions, L_ref=0.1m can indicate reference length, λ₁ =10^-5It can indicate the interporosity flow coefficient of inner region, λ₂=10^-4It can indicate the interporosity flow coefficient of outskirt, M₁₂=5 can indicate inner region Crack mobility ratio between outskirt, W₁₂=0.8 can indicate that ratio, d are held in the storage between inner region and outskirt_mf1It can indicate inner region Fracture network mass fractal dimension, d_mf2It can indicate the mass-fractal dimension of the fracture network of outskirt, θ₁It can indicate inner region Unusual diffusion index, θ₂It can indicate the unusual diffusion index of outskirt.In Fig. 6, curve L1 and L3 are overlapped, curve L2 and L4 weight It closes, ψ_D'It can indicate pressure curve, ψ '_WDt_D/C_D'It can indicate derivative curve.L1 is indicated using the present processes in d_mf1= d_mf2=2, θ₁=θ₂Pressure curve when=0, the pressure curve that the method that L3 indicates the prior art obtains.L2 indicates to use this Shen Method please is in d_mf1=d_mf2=2, θ₁=θ₂Derivative curve when=0, L4 indicate the derivative curve that art methods obtain.

Fig. 7 is Pressure behaviour curve synoptic diagram in one embodiment provided by the present application, as shown in fig. 7, in the application Pressure behaviour curve can be divided into pressure curve and derivative curve, in order to observe complete flowing stage, take to relevant parameter Theoretical value.Specifically, dotted line L1 represents d_mf1=d_mf2=2 and θ₁=θ₂The result of fractal number rank diffusion model when=0；Point Scribing line L2 represents d_mf1=d_mf2=1.95 and θ₁=θ₂The result of fractal number rank diffusion model when=0.05；Dotted line L3 is represented d_mf1=1.95, d_mf2=2, θ₁=0.05, θ₂The result of fractal number rank diffusion model when=0；Solid line L4 represents d_mf1= 1.9, d_mf2=1.95, θ₁=0.1, θ₂The result of fractal number rank diffusion model when=0.05.d_mf1It can indicate splitting for inner region Stitch the mass fractal dimension of network, d_mf2It can indicate the mass-fractal dimension of the fracture network of outskirt, θ₁It can indicate the anti-of inner region Normal diffusion index, θ₂It can indicate the unusual diffusion index of outskirt.The meaning of parameters is identical in Fig. 6 in Fig. 7, herein not It repeats again.

As shown in fig. 7, the fractal number rank of point shape tight gas reservoir volume fracturing horizontal well can be flowed can be divided into I-Ⅸ 9 flowing stage: pure bottom-hole storage stage, transition flow stage early stage, inner region early stage pseudo-linear stream stage, inner region quasi-stable state Channelling stage, inner region mid-term pseudoradial flow stage, inner region mid-term pseudo-linear stream stage, inner region advanced stage pseudoradial flow stage, advanced stage Transition flow stage, composite system advanced stage pseudoradial flow stage.Pressure dynamic analysis, tool are carried out to each flowing stage according to Fig. 7 Body analyzes result can be with reference to following statement:

Flowing stage I: in the pure bottom-hole storage stage, mainly being controlled by well bore storage effect, pressure curve and derivative curve weight It is combined into the straight line that slope is 1.

Flowing stage II: there is a hump in transition flow stage early stage, derivative curve, and the height and width of hump depend on The size of bottom-hole storage coefficient and skin factor.

Flowing stage III: pseudo-linear stream stage early stage, derivative curve show as straight line.In this stage, interior differentiation Gas flowing in shape microcrack network is approximately the linear flow perpendicular to each single waterpower major fracture, can be referred to as pseudo-linear Stream.It is worth noting that, the derivative curve in this stage is no longer the straight line of fixed slope, slope value is with inner region microcrack Fractal dimension d_mf1With inner region abnormality diffusion index θ₁Change and change, and slope value is typically larger than 1/2, only works as d_mf1= d_mf2=2 and θ₁=θ₂When=0, the value of straight slope is just equal to 1/2.d_mf2It can indicate the mass fractal of the fracture network of outskirt Dimension.

Flowing stage IV: the inner region quasi-stable state channelling stage, at this time in territorial matrix with divided pressure difference in shape microcrack network Establish, the gas in matrix in a manner of quasi-stable state to shape microcrack system channelling is divided because being hypothesized matrix and point shape fine fisssure Stitching is quasi-stable state channelling between network, is occurred apparent " caviton " on derivative curve, and the position of " caviton " on derivative curve It sets high than " caviton " position of traditional Euclid's model.

Flowing stage V: the inner region mid-term pseudoradial flow stage divides the gas in shape intrinsic fracture network in each single waterpower It is flowed in a manner of pseudoradial flow to each single waterpower major fracture around major fracture, interferes between waterpower major fracture not yet formed at this time. It it is worth noting that, derivative curve here is inclined straight line, and is no longer value for 1/ (2MM₁₂) horizontal line.Derivative curve The position of straightway with inner region microcrack fractal dimension d_mf1With inner region abnormality diffusion index θ₁Change and change, usually compares 1/(2MM₁₂) horizontal line wants high.Particularly, work as d_mf1=d_mf2=2 and θ₁=θ₂When=0, it is 1/ (2MM that derivative curve, which is degenerated,₁₂) Horizontal line.d_mf2It can indicate the mass-fractal dimension of the fracture network of outskirt, M₁₂It can indicate the crack between inner region and outskirt Mobility ratio, M can indicate man-made fracture item number.

Flowing stage VI: inner region mid-term pseudo-linear stream stage, derivative curve also show as straight line section.When pressure wave passes It is multicast to except waterpower major fracture region, starts generation between adjacent hydraulic fracture at this time and interfere with each other, divide in shape microcrack network Gas flowed in a manner of being approximately parallel to the linear flow in hydraulic fracture face to entire hydraulic fracturing area, referred to as pseudo-linear stream. Due to considering point shape microcrack and unusual diffusion, the slope of derivative curve straightway is typically larger than the 1/2 of traditional Euclidean model, And work as d_mf1=d_mf2=2 and θ₁=θ₂When=0, the slope of straight line is equal to 1/2.Since this phase duration is shorter, so The feature of curved section is less obvious.d_mf1It can indicate the mass-fractal dimension of the fracture network of inner region, d_mf2It can indicate outskirt Fracture network mass-fractal dimension.

Flowing stage VII: inner region advanced stage pseudoradial flow stage, derivative curve show as inclined straight line, and are no longer that value is 1/(2M₁₂) horizontal line.Similarly, the position of derivative curve straightway is usually above 1/ (2M₁₂) horizontal line, and d_mf1= d_mf2=2 and θ₁=θ₂When=0, it is 1/ (2M that straight line, which is degenerated,₁₂) horizontal line.In fact, the duration in this stage is very It is short, only when interior outskirt compound radius is observed much larger than waterpower major fracture half is Shi Caineng long.In order to make typical curve feature Obviously, interior outskirt compound radius can be taken as 1500m, a length of 40m of waterpower major fracture half by the application.d_mf1It can indicate inner region Fracture network mass-fractal dimension, d_mf2It can indicate the mass-fractal dimension of the fracture network of outskirt, M₁₂In can indicating Crack mobility ratio between area and outskirt, M can indicate man-made fracture item number.

Flowing stage VIII: in the transition flow stage in advanced stage, when pressure wave is diffused into interior outskirt compound radius, the flowing of outskirt is opened Begin, the transition flow regime of appearance.

Flowing stage Ⅸ: the composite system advanced stage pseudoradial flow stage works as pressure since the lateral outer boundary of gas reservoir is infinity Wave, which has propagated to, remotely divides the gas of shape tight gas reservoir in a manner of approximate pseudoradial flow to gas reservoir center flow, whole system (i.e. double medium model) reaches quasi- stable state.Since inner region and outskirt all consider point shape fracture network and unusual diffusion, Derivative curve is no longer the horizontal line that fixed value is 1/2, but inclined straight line.Similarly, the position of derivative curve straightway Horizontal line usually above 1/2, and work as d_mf1=d_mf2=2 and θ₁=θ₂When=0, straight line degenerate be 1/2 horizontal line.d_mf1It can To indicate the mass-fractal dimension of the fracture network of inner region, d_mf2It can indicate the mass-fractal dimension of the fracture network of outskirt.

After the Pressure behaviour curve for obtaining well section to be analyzed, sensitivity to parameter point can also be carried out according to Pressure behaviour curve Analysis, influence of the analysis parameters to the Pressure behaviour variation of volume fracturing horizontal well.Fig. 8 is in the application one embodiment The schematic diagram that the mass-fractal dimension of the fracture network of inner region influences Pressure behaviour, Fig. 8 shows θ₁=1.5, θ₂=0, d_mf2=2 When, the mass-fractal dimension d of different inner region microcrack networks_mf1Corresponding Pressure behaviour curve.Curve L1-L4 is respectively indicated d_mf1=2, d_mf1=1.95, d_mf1=1.9, d_mf1=1.85 corresponding pressure curves, curve L5-L8 respectively indicate d_mf1=2, d_mf1 =1.95, d_mf1=1.9, d_mf1=1.85 corresponding derivative curves.

As shown in figure 8, the mass-fractal dimension d of inner region microcrack network_mf1To in addition to pure bottom-hole storage stage and early stage mistake The almost entire flowing stage for crossing the stream stage suffers from very big influence.With the mass-fractal dimension of inner region microcrack network d_mf1Reduction, the value of zero dimension pseudopressure and pseudopressure derivative also increases with it, show as the corresponding stage of typical curve, The position of derivative curve is higher or the slope of straightway is bigger.The meaning of parameters is identical in Fig. 6 in Fig. 8, herein not It repeats again.

Fig. 9 is that the mass-fractal dimension of the fracture network of the outskirt in the application one embodiment influences Pressure behaviour Schematic diagram, Fig. 9 indicate d_mf1=1.95, θ₁=0.05, θ₂When=0, the mass-fractal dimension d of different outskirt intrinsic fractures_mf2It is right The Pressure behaviour curve answered.Curve L1-L4 respectively indicates d_mf2=2, d_mf2=1.95, d_mf2=1.9, d_mf2=1.85 corresponding pressures Force curve, curve L5-L8 respectively indicate d_mf2=2, d_mf2=1.95, d_mf2=1.9, d_mf2=1.85 corresponding derivative curves.Such as Shown in Fig. 9, the mass-fractal dimension d of outskirt intrinsic fracture_mf2The pressure for only influencing the pseudoradial flow stage in composite system advanced stage is dynamic State.In the composite system advanced stage pseudoradial flow stage, derivative curve is no longer the horizontal line that fixed value is 1/2, but inclined straight Line, the slope of straight line with intrinsic fracture mass-fractal dimension d_mf2Reduce and increases.Same Fig. 6 of the meaning of parameters in Fig. 9 In it is identical, details are not described herein again.

Figure 10 is the schematic diagram that the inner region abnormality diffusion index in the application one embodiment influences Pressure behaviour, Figure 10 Indicate d_mf1=1.95, d_mf2=2, θ₂When=0, different inner region abnormality diffusion index θ₁Corresponding Pressure behaviour curve.Curve L1-L4 respectively indicates θ₁=0, θ₁=0.05, θ₁=0.1, θ₁=0.2 corresponding pressure curve, curve L5-L8 respectively indicate θ₁= 0、θ₁=0.05, θ₁=0.1, θ₁=0.2 corresponding derivative curve.As shown in Figure 10, inner region abnormality diffusion index θ₁To pressure song The influence of line and derivative curve was mainly reflected in after the inner region early stage transition flow stage, inner region abnormality diffusion index θ₁Not only change The value for having become zero dimension pseudopressure and pseudopressure derivative has also postponed pressure response time, the position of pressure curve and derivative curve It sets and also moves up and move to right therewith.Inner region abnormality diffusion index θ₁Value it is bigger, pressure response spot is more late.Each ginseng in Figure 10 Several meanings is identical in Fig. 6, and details are not described herein again.

Figure 11 is the schematic diagram that the outskirt abnormality diffusion index in the application one embodiment influences Pressure behaviour, Figure 11 Indicate d_mf1=1.9, d_mf2=1.95, θ₁When=0.1, different outskirt abnormality diffusion index θ₂Corresponding Pressure behaviour curve.It is bent Line L1-L4 respectively indicates θ₂=0, θ₂=0.05, θ₂=0.1, θ₂=0.2 corresponding pressure curve, curve L5-L8 respectively indicate θ₂ =0, θ₂=0.05, θ₂=0.1, θ₂=0.2 corresponding derivative curve.As shown in figure 11, outskirt abnormality diffusion index θ₂Only shadow Ring the Pressure behaviour in advanced stage transition flow stage and composite system advanced stage pseudoradial flow stage.Outskirt abnormality diffusion index θ₂Value is got over Greatly, the duration in advanced stage transition flow stage is longer, and the pressure response time of origin in composite system advanced stage pseudoradial flow stage is got over Evening.Meanwhile outskirt abnormality diffusion index θ₂It is worth bigger, the pressure curve and derivative curve in composite system advanced stage pseudoradial flow stage Position it is higher.The meaning of parameters is identical in Fig. 6 in Figure 11, and details are not described herein again.

Figure 12 is the schematic diagram that the inside and outside compound radius in the application one embodiment influences Pressure behaviour, and Figure 12 is indicated Different interior outskirt compound radius r_mCorresponding Pressure behaviour curve.Curve L1-L4 respectively indicates r_m=500m, r_m=1000m, r_m=1500m, r_mThe corresponding pressure curve of=2000m, curve L5-L8 respectively indicate r_m=500m, r_m=1000m, r_m= 1500m、r_mThe corresponding derivative curve of=2000m.As shown in figure 12, interior outskirt compound radius r_mIt is main to influence inner region advanced stage quasi- diameter To the duration in stream stage, interior outskirt compound radius r_mBigger, the duration in inner region advanced stage pseudoradial flow stage is longer, evening Transition flow stage phase is also more late with regard to spot.The meaning of parameters is identical in Fig. 6 in Figure 12, and details are not described herein again.

Figure 13 is the influence schematic diagram that interior outskirt crack mobility ratio changes Pressure behaviour in the application one embodiment, figure 13 indicate different interior outskirt crack mobility ratio M₁₂Corresponding Pressure behaviour curve.Curve L1-L4 respectively indicates M₁₂=1, M₁₂= 2、M₁₂=5, M₁₂=10 corresponding pressure curves, curve L5-L8 respectively indicate M₁₂=1, M₁₂=2, M₁₂=5, M₁₂=10 is corresponding Derivative curve.As shown in figure 13, interior outskirt crack mobility ratio M₁₂Intend in addition to pure bottom-hole storage stage and composite system advanced stage The Pressure behaviour of nearly all flowing stage except the radial flow stage suffers from significant impact.Specific manifestation are as follows: inside and outside area Crack mobility ratio M₁₂Bigger, zero dimension pseudopressure and pseudopressure derivative are smaller, and the position of pressure curve and derivative curve is lower.This It is because dimensionless variable as defined herein is all based on outskirt parameter definition, interior outskirt crack mobility ratio M₁₂More Greatly, the permeability of inner region microcrack is higher.Again since horizontal well is fixed output quota production, the permeability of inner region microcrack is higher, seeps Flow resistance power is smaller, and pressure drop is smaller, zero dimension pseudopressure and pseudopressure derivative it is smaller.Same Fig. 6 of the meaning of parameters in Figure 13 In it is identical, details are not described herein again.

Figure 14 is that the schematic diagram for comparing Pressure behaviour influence is held in interior outskirt storage in the application one embodiment, and Figure 14 is indicated not Same Pressure behaviour is all held by the storage of interior outskirt and compares W₁₂Corresponding Pressure behaviour curve.Curve L1-L4 respectively indicates W₁₂=0.5, W₁₂ =0.75, W₁₂=1, W₁₂=1.75 corresponding pressure curves, curve L5-L8 respectively indicate W₁₂=0.5, W₁₂=0.75, W₁₂= 1、W₁₂=1.75 corresponding derivative curves.As shown in figure 14, in addition to pure bottom-hole storage stage and composite system advanced stage pseudoradial flow Stage, the Pressure behaviour of almost all of flowing stage are all held by the storage of interior outskirt and compare W₁₂Influence.From transition flow stage early stage To the entire flowing period between the transition flow stage in advanced stage, holds with the storage of interior outskirt and compare W₁₂Increase, pressure curve and derivative Curve generates one and significantly moves to right, i.e., pressure response time is postponed.Equally, this is because dimensionless variable as defined herein It is all based on outskirt parameter definition, so the storage of interior outskirt, which is held, compares W₁₂Bigger, the porosity compressed coefficient product of inner region is bigger.This Sample results in the pressure diffusion speed of inner region just slower, therefore pressure response time is postponed, and pressure curve is moved to right with derivative curve. The meaning of parameters is identical in Fig. 6 in Figure 14, and details are not described herein again.

Volume fracturing Dynamic Analysis of Horizontal Well Pressure method provided by the present application, utilizes fractal theory and fractional calculus Method establishes the Pressure Dynamic Model of volume fracturing horizontal well, using and can use semi analytical method solve Pressure Dynamic Model, The solution of Pressure Dynamic Model can quickly and accurately be got.According to the solution of the Pressure Dynamic Model got, can obtain The fractal number rank flowing of horizontal well can be divided into 9 streams according to dynamic pressure profile is obtained by horizontal well pressure performance graph Dynamic stage, further progress parameters sensitivity analysis.While reducing the calculation amount of Dynamic Analysis of Horizontal Well Pressure, improve The accuracy of volume fracturing Dynamic Analysis of Horizontal Well Pressure.

Based on volume fracturing Dynamic Analysis of Horizontal Well Pressure method described above, this specification one or more embodiment A kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure device is also provided.The device may include having used this specification implementation The example system (including distributed system) of the method, software (application), module, component, server, client etc. and in conjunction with must The device for the implementation hardware wanted.Based on same innovation thinking, in one or more embodiments that this specification embodiment provides Device is as described in the following examples.Since the implementation that device solves the problems, such as is similar to method, this specification is implemented The implementation of the specific device of example may refer to the implementation of preceding method, and overlaps will not be repeated.It is used below, term The combination of the software and/or hardware of predetermined function may be implemented in " unit " or " module ".Although described in following embodiment Device is preferably realized with software, but the realization of the combination of hardware or software and hardware is also that may and be contemplated.

Specifically, Figure 15 is the volume fracturing Dynamic Analysis of Horizontal Well Pressure device in a kind of embodiment provided by the present application Modular structure schematic diagram, as shown in figure 15, volume fracturing Dynamic Analysis of Horizontal Well Pressure device provided by the present application include: be situated between Matter model building module 151, Pressure Dynamic Model establish module 152, pressure dynamic analysis module 153.

Dielectric model establishes module 151, can be used for constructing the corresponding dual media of well section to be analyzed using fractal theory Model；

Pressure Dynamic Model establishes module 152, can be used for utilizing fractional calculus according to the double medium model Method establishes the Pressure Dynamic Model of the well section to be analyzed；

Pressure dynamic analysis module 153 can be used for carrying out the well section to be analyzed according to the Pressure Dynamic Model Pressure dynamic analysis.

Volume fracturing Dynamic Analysis of Horizontal Well Pressure device provided by the present application establishes well section to be analyzed using fractal theory Corresponding double medium model, can accurately describe the structure of well section microcrack network to be analyzed and the complexity of distribution, Scrambling and heterogeneity provide accurate theoretical basis for the subsequent pressure dynamic analysis for carrying out well section to be analyzed, mention The high accuracy of volume fracturing Dynamic Analysis of Horizontal Well Pressure.Meanwhile the well to be analyzed of part can be used using fractal theory The feature of the microcrack network of section, states the feature of the microcrack network of full well section, reduces and subsequent establishes well section pressure to be analyzed The calculation amount of power dynamic model further reduces the calculation amount of volume fracturing Dynamic Analysis of Horizontal Well Pressure.The application can be with The Pressure Dynamic Model of volume fracturing horizontal well is established using fractional calculus method, reduces volume fracturing horizontal well pressure dynamic While the calculation amount of model foundation, the accuracy of volume fracturing horizontal well pressure Dynamic Model is improved, is further increased The accuracy of volume fracturing Dynamic Analysis of Horizontal Well Pressure.

The above-mentioned volume fracturing Dynamic Analysis of Horizontal Well Pressure method or apparatus that this specification embodiment provides can counted Corresponding program instruction is executed to realize, such as using the c++ language of windows operating system in the end PC reality by processor in calculation machine Existing, Linux system is realized or other are for example realized using android, iOS system programming language in intelligent terminal, and Processing logic realization based on quantum computer etc..A kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure dress that this specification provides In another embodiment set, Figure 16 is that another volume fracturing Dynamic Analysis of Horizontal Well Pressure device provided by the present application is implemented The modular structure schematic diagram of example, as shown in figure 16, the volume fracturing horizontal well pressure dynamic point that another embodiment of the application provides Analysis apparatus may include processor 161 and the memory 162 for storage processor executable instruction,

Processor 161 and memory 162 pass through bus 163 and complete mutual communication；

The processor 161 is used to call the program instruction in the memory 162, to execute above-mentioned each volume fracturing water Method provided by horizontal well pressure dynamic analysis embodiment of the method, for example, construct well section pair to be analyzed using fractal theory The double medium model answered；According to the double medium model, the well section to be analyzed is established using fractional calculus method Pressure Dynamic Model；According to the Pressure Dynamic Model, pressure dynamic analysis is carried out to the well section to be analyzed.

It should be noted that specification device described above can also include it according to the description of related method embodiment His embodiment, concrete implementation mode are referred to the description of embodiment of the method, do not repeat one by one herein.In the application Various embodiments are described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, often What a embodiment stressed is the difference from other embodiments.For hardware+program class embodiment, Since it is substantially similar to the method embodiment, so being described relatively simple, related place is said referring to the part of embodiment of the method It is bright.

This specification embodiment is not limited to meet industry communication standard, standard computer data processing sum number According to situation described in storage rule or this specification one or more embodiment.The right way of conduct is made in certain professional standards or use by oneself In formula or the practice processes of embodiment description embodiment modified slightly also may be implemented above-described embodiment it is identical, it is equivalent or The implementation result being anticipated that after close or deformation.Using these modifications or deformed data acquisition, storage, judgement, processing side The embodiment of the acquisitions such as formula still may belong within the scope of the optional embodiment of this specification embodiment.

In the 1990s, the improvement of a technology can be distinguished clearly be on hardware improvement (for example, Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So And with the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit. Designer nearly all obtains corresponding hardware circuit by the way that improved method flow to be programmed into hardware circuit.Cause This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, programmable logic device (Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate Array, FPGA)) it is exactly such a integrated circuit, logic function determines device programming by user.By designer Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, designs and makes without asking chip maker Dedicated IC chip.Moreover, nowadays, substitution manually makes IC chip, this programming is also used instead mostly " is patrolled Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development, And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language (Hardware Description Language, HDL), and HDL is also not only a kind of, but there are many kind, such as ABEL (Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL (Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language) etc., VHDL (Very-High-Speed is most generally used at present Integrated Circuit Hardware Description Language) and Verilog.Those skilled in the art also answer This understands, it is only necessary to method flow slightly programming in logic and is programmed into integrated circuit with above-mentioned several hardware description languages, The hardware circuit for realizing the logical method process can be readily available.

Controller can be implemented in any suitable manner, for example, controller can take such as microprocessor or processing The computer for the computer readable program code (such as software or firmware) that device and storage can be executed by (micro-) processor can Read medium, logic gate, switch, specific integrated circuit (Application Specific Integrated Circuit, ASIC), the form of programmable logic controller (PLC) and insertion microcontroller, the example of controller includes but is not limited to following microcontroller Device: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320 are deposited Memory controller is also implemented as a part of the control logic of memory.It is also known in the art that in addition to Pure computer readable program code mode is realized other than controller, can be made completely by the way that method and step is carried out programming in logic Controller is obtained to come in fact in the form of logic gate, switch, specific integrated circuit, programmable logic controller (PLC) and insertion microcontroller etc. Existing identical function.Therefore this controller is considered a kind of hardware component, and to including for realizing various in it The device of function can also be considered as the structure in hardware component.Or even, it can will be regarded for realizing the device of various functions For either the software module of implementation method can be the structure in hardware component again.

System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity, Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer.Specifically, computer for example may be used Think personal computer, laptop computer, vehicle-mounted human-computer interaction device, cellular phone, camera phone, smart phone, individual Digital assistants, media player, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or The combination of any equipment in these equipment of person.

Although this specification one or more embodiment provides the method operating procedure as described in embodiment or flow chart, It but may include more or less operating procedure based on conventional or without creativeness means.The step of being enumerated in embodiment Sequence is only one of numerous step execution sequence mode, does not represent and unique executes sequence.Device in practice or When end product executes, can be executed according to embodiment or the execution of method shown in the drawings sequence or parallel (such as it is parallel The environment of processor or multiple threads, even distributed data processing environment).The terms "include", "comprise" or its Any other variant is intended to non-exclusive inclusion so that include the process, methods of a series of elements, product or Equipment not only includes those elements, but also including other elements that are not explicitly listed, or further include for this process, Method, product or the intrinsic element of equipment.In the absence of more restrictions, being not precluded is including the element There is also other identical or equivalent elements in process, method, product or equipment.The first, the second equal words are used to indicate name Claim, and does not indicate any particular order.

For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this The function of each module can be realized in the same or multiple software and or hardware when specification one or more, it can also be with The module for realizing same function is realized by the combination of multiple submodule or subelement etc..Installation practice described above is only It is only illustrative, for example, in addition the division of the unit, only a kind of logical function partition can have in actual implementation Division mode, such as multiple units or components can be combined or can be integrated into another system or some features can be with Ignore, or does not execute.Another point, shown or discussed mutual coupling, direct-coupling or communication connection can be logical Some interfaces are crossed, the indirect coupling or communication connection of device or unit can be electrical property, mechanical or other forms.

The present invention be referring to according to the method for the embodiment of the present invention, the process of device (system) and computer program product Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.

These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.

These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.

In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net Network interface and memory.

Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/or The forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable medium Example.

Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data. The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM), Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage, graphene stores or other Magnetic storage device or any other non-transmission medium, can be used for storage can be accessed by a computing device information.According to herein In define, computer-readable medium does not include temporary computer readable media (transitory media), such as the data of modulation Signal and carrier wave.

It will be understood by those skilled in the art that this specification one or more embodiment can provide as method, system or calculating Machine program product.Therefore, this specification one or more embodiment can be used complete hardware embodiment, complete software embodiment or The form of embodiment combining software and hardware aspects.Moreover, this specification one or more embodiment can be used at one or It is multiple wherein include computer usable program code computer-usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) on the form of computer program product implemented.

This specification one or more embodiment can computer executable instructions it is general on It hereinafter describes, such as program module.Generally, program module includes executing particular task or realization particular abstract data type Routine, programs, objects, component, data structure etc..This this specification one can also be practiced in a distributed computing environment Or multiple embodiments, in these distributed computing environments, by being held by the connected remote processing devices of communication network Row task.In a distributed computing environment, program module can be located at the local and remote computer including storage equipment In storage medium.

All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method Part explanation.In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", The description of " specific example " or " some examples " etc. means specific features described in conjunction with this embodiment or example, structure, material Or feature is contained at least one embodiment or example of this specification.In the present specification, to the signal of above-mentioned term Property statement be necessarily directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described It may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other, this The technical staff in field can be by the spy of different embodiments or examples described in this specification and different embodiments or examples Sign is combined.

The foregoing is merely the embodiments of this specification one or more embodiment, are not limited to book explanation Book one or more embodiment.To those skilled in the art, this specification one or more embodiment can have various Change and variation.All any modification, equivalent replacement, improvement and so within the spirit and principles of the present application should all include Within scope of the claims.

Claims

1. a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method, which is characterized in that the described method includes:

According to the double medium model, the Pressure Dynamic Model of the well section to be analyzed is established using fractional calculus method；

According to the Pressure Dynamic Model, pressure dynamic analysis is carried out to the well section to be analyzed；

It is wherein, described to construct the corresponding double medium model of the well section to be analyzed using fractal theory, comprising:

Using fractal theory using the corresponding microcrack network of the well section to be analyzed as fractal, the fractal is embedded into In Euclid's matrix；

The inner region includes waterpower major fracture, Euclid's matrix, point shape microcrack, the outskirt include Euclid's matrix and Divide shape intrinsic fracture, described point of shape microcrack includes the intrinsic fracture of induced fractures and activation, and described point of shape intrinsic fracture includes Original intrinsic fracture；

It is described according to the double medium model, the Pressure behaviour mould of the well section to be analyzed is established using fractional calculus method Type, comprising:

According to the double medium model, inner region fractal number rank diffusion model is established respectively using the fractional calculus method With outskirt fractal number rank diffusion model；

Using the combination of the inner region fractal number rank diffusion model and the outskirt fractal number rank diffusion model as the pressure Power dynamic model；

The inner region fractal number rank diffusion model includes:

In above formula, r_DIndicate dimensionless radius, θ₁Indicate unusual diffusion index, ψ_Df1Indicate the Fracture System zero dimension of the inner region Pressure, β₁=d_mf1-θ₁- 1, d_mf1Indicate the mass-fractal dimension of the fracture network of the inner region,M₁₂Indicate institute State the crack mobility ratio between inner region and outskirt；W₁₂Indicate that ratio, t are held in the storage between the inner region and outskirt_DWhen indicating zero dimension Between,φ_fr1Indicate porosity of the crack of the inner region under reference length, C_tf1It indicates The Fracture System system compressibility of the inner region, φ_m1Indicate the porosity of the matrix system of the inner region, C_tm1Described in expression The matrix system system compressibility of inner region, γ₁=2/ (2+ θ₁), ψ_Dm1Indicate the matrix system zero dimension pressure of the inner region, λ₁Indicate the interporosity flow coefficient of the inner region；

The outskirt fractal number rank diffusion model includes:

In above formula, r_DIndicate dimensionless radius, θ₂Indicate unusual diffusion index, ψ_Df2Indicate the Fracture System zero dimension of the outskirt Pressure, β₂=d_mf2-θ₂- 1, d_mf2Indicate the mass-fractal dimension of the fracture network of the outskirt, t_DIndicate non dimensional time,φ_fr2Indicate porosity of the crack of the outskirt under reference length, C_tf2Described in expression The Fracture System system compressibility of outskirt, φ_m2Indicate the porosity of the matrix system of the outskirt, C_tm2Indicate the outskirt Matrix system system compressibility, γ₂=2/ (2+ θ₂), ψ_Dm2Indicate the matrix system zero dimension pressure of the outskirt, λ₂Table Show the interporosity flow coefficient of the outskirt.

2. a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method as described in claim 1, which is characterized in that the inner region For round radial compound region, the waterpower major fracture in the inner region is vertical lap seam.

3. such as a kind of described in any item volume fracturing Dynamic Analysis of Horizontal Well Pressure methods of claim 1-2, which is characterized in that It is described according to the Pressure Dynamic Model, pressure dynamic analysis is carried out to the well section to be analyzed, comprising:

The Pressure Dynamic Model is solved using semi analytical method, according to the solution of the Pressure Dynamic Model, is drawn described to be analyzed The Pressure behaviour curve of well section；

4. a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method as claimed in claim 3, which is characterized in that the utilization Semi analytical method solves the Pressure Dynamic Model, comprising:

It is dynamic that the Laplace space pressure is solved using line source function, discrete method, principle of stacking and Stehfest Numerical inverse method States model obtains the solution of the Pressure Dynamic Model.

5. a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method as claimed in claim 4, which is characterized in that the basis The Pressure behaviour curve carries out pressure dynamic analysis to the well section to be analyzed, comprising:

The Pressure behaviour curve is carried out to the division of flowing stage, and according to the Pressure behaviour curve to the well to be analyzed Section pressure carries out parameters sensitivity analysis.

6. a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure device, which is characterized in that described device includes:

Pressure Dynamic Model establishes module, is used for according to the double medium model, using described in the foundation of fractional calculus method The Pressure Dynamic Model of well section to be analyzed；

Pressure dynamic analysis module, for carrying out Pressure behaviour point to the well section to be analyzed according to the Pressure Dynamic Model Analysis；

The inner region fractal number rank diffusion model includes:

The outskirt fractal number rank diffusion model includes:

7. a kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure device, which is characterized in that including processor and at storage The memory of device executable instruction is managed, the processor is realized when executing described instruction such as any one of claim 1 to 5 institute The step of stating method.