Multi-horizontal well fracturing design parameter optimization method based on optimal bridging joint network

Abstract

The invention provides a multi-horizontal well fracturing design parameter optimization method based on an optimal bridging fracture net, and the method comprises the steps: carrying out the fracture net form expansion of a target reservoir based on given fracturing design parameters, and obtaining the fracture point data of a plurality of fractured horizontal wells of the target reservoir; determining the fracture form of each fractured horizontal well according to the fracture point data of each fractured horizontal well; according to the fracture form of each fractured horizontal well, the initial bridging degree of the fracture network of the multiple fractured horizontal wells of the target reservoir is calculated, and the bridging degree of the fracture network of the multiple fractured horizontal wells is in negative correlation with the sum of the total unmodified volume and the total repeated modified volume of the multiple fractured horizontal wells; according to the initial bridging degree, a preset parameter optimization algorithm is used for optimizing the fracturing design parameters with the maximum bridging degree of the fracture network of the multiple fractured horizontal wells as the target, and the fracturing design parameters of the multiple fractured horizontal wells under the maximum bridging degree are obtained.

Description

Technical Field

[0001] The invention relates to the technical field of oil and gas reservoir production optimization, and in particular to a method for optimizing design parameters of multi-horizontal well fracturing based on an optimally bridged fracture network. Background Art

[0002] Fracturing plays an important role in the development of unconventional oil and gas. In recent years, deep coalbed methane exploration and development have made major breakthroughs. Now horizontal well fracturing has become an important technology for coalbed methane, especially deep coalbed methane development. Conventional fracturing has two main problems. First, it only targets a single well and cannot consider the interactive effects of the fracturing network between multiple wells. Second, there are a large number of unreconstructed areas within the well control range of multiple wells. Third, under multi-well conditions, there are often many repeated reconstructions between segments and between wells. Fourth, it is impossible to simultaneously optimize the segment cluster deployment design parameters (segment length, perforation position, fracture half-length) + fracturing construction parameters (displacement, liquid volume, etc.). Therefore, large-scale fracturing has problems of transitional reconstruction (resulting in cost waste) or local reconstruction blank areas (resulting in the inability to mobilize resources), which brings challenges to efficient development. Summary of the invention

[0003] In view of the problems in the prior art, an embodiment of the present invention provides a method for optimizing design parameters of multi-horizontal well fracturing based on an optimal bridged fracture network.

[0004] On the one hand, the present invention proposes a method for optimizing design parameters of multi-horizontal well fracturing based on an optimal bridged fracture network, comprising:

[0005] Based on the given fracturing design parameters, the fracture network of the target reservoir is expanded to obtain the fracture point data of multiple fracturing horizontal wells in the target reservoir;

[0006] Determine the fracture morphology of each fractured horizontal well according to the fracture point data of each fractured horizontal well;

[0007] Calculating the initial healing degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the fracture morphology of each fractured horizontal well, wherein the healing degree of the fracture network of multiple fractured horizontal wells is negatively correlated with the sum of the total unreformed volume and the total repeatedly reformed volume of the multiple fractured horizontal wells;

[0008] According to the initial bridging degree, the fracturing design parameters are optimized by using a preset parameter optimization algorithm with the goal of maximizing the bridging degree of the multiple fractured horizontal wells, so as to obtain the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree.

[0009] In some embodiments, the fracturing design parameters include construction parameters and segment cluster location parameters; wherein the construction parameters include construction displacement and / or construction fluid volume.

[0010] In some embodiments, the fracture patterns of the fracturing horizontal wells include the total well-controlled volume, the total effective stimulation volume, the total repeated stimulation volume, and / or the total un-stimulated volume; determining the fracture patterns of the fracturing horizontal wells according to the fracture point data of each fracturing horizontal well includes:

[0011] Obtaining the effective production range of each fracturing horizontal well through actual production dynamic simulation;

[0012] Determining the dynamic stimulation volume of a single-stage fracture of each fracturing horizontal well according to the effective production range and the fracture point data of each fracturing horizontal well;

[0013] Calculating the total effective stimulation volume, the total repeated stimulation volume, the total well-controlled volume, and / or the total un-stimulated volume of each fracturing horizontal well based on the background grid method according to the dynamic stimulation volume of a single-stage fracture of each fracturing horizontal well.

[0014] In some embodiments, calculating the total effective stimulation volume, the total repeated stimulation volume, the total well-controlled volume, and / or the total un-stimulated volume of each fracturing horizontal well based on the background grid method according to the dynamic stimulation volume of a single-stage fracture of each fracturing horizontal well includes:

[0015] Marking the index information of each fracture of each fracturing horizontal well based on the background grid method;

[0016] According to the index information, obtaining the total effective stimulation volume of each fracturing horizontal well by taking the union of the dynamic stimulation volumes of all single-stage fractures of each fracturing horizontal well; and / or

[0017] According to the index information, obtaining the total inter-stage repeated stimulation volume of each fracturing horizontal well by taking the intersection of the dynamic stimulation volumes of single-stage fractures in each fracturing horizontal well;

[0018] According to the index information, obtaining the total inter-well repeated stimulation volume of each fracturing horizontal well by taking the intersection of the fracture dynamic stimulation volumes of two-by-two fracturing horizontal wells;

[0019] Determining the total repeated stimulation volume of each fracturing horizontal well according to the total inter-stage repeated stimulation volume of each fracturing horizontal well and the total inter-well repeated stimulation volume of each fracturing horizontal well; and / or

[0020] According to the index information, obtaining the total well-controlled volume of each fracturing horizontal well by taking the union of the well-controlled volumes of each fracturing horizontal well;

[0021] According to the index information, obtaining the total un-stimulated volume of each fracturing horizontal well by taking the union of the volumes of un-stimulated grid units of each fracturing horizontal well.

[0022] In some embodiments, calculating the initial bridging degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture patterns of each fracturing horizontal well includes:

[0023] Calculate the initial bridging degree of the fracture networks of multiple fractured horizontal wells in the target reservoir according to the ratio of the sum of the total repeated reconstruction volume and the total unreconstructed volume of each fractured horizontal well to the total well control volume; or

[0024] Calculate the initial bridging degree of the fracture networks of multiple fractured horizontal wells in the target reservoir according to the ratio of the total effective reconstruction volume of each fractured horizontal well to the total well control volume.

[0025] In some embodiments, according to the initial bridging degree, using a preset parameter optimization algorithm, with the goal of maximizing the bridging degree of the fracture networks of the multiple fractured horizontal wells, optimize the fracture design parameters, and the fracture design parameters of the multiple fractured horizontal wells at the maximum bridging degree obtained include:

[0026] Optimize the fracture design parameters using a preset parameter optimization algorithm;

[0027] Based on the optimized fracture design parameters, expand the fracture network morphology again to obtain the fracture point data of each fractured horizontal well;

[0028] According to the fracture point data of each fractured horizontal well, determine the fracture morphology of each fractured horizontal well;

[0029] According to the fracture morphology of each fractured horizontal well, calculate the current bridging degree of the fracture networks of multiple fractured horizontal wells in the target reservoir;

[0030] If the current bridging degree is greater than the initial bridging degree, continue to use the preset parameter optimization algorithm to optimize the fracture design parameters again based on the optimized fracture design parameters;

[0031] Continue to iterate until a preset termination condition is reached, and obtain the maximum bridging degree and the fracture design parameters of the multiple fractured horizontal wells at the maximum bridging degree.

[0032] In some embodiments, the preset parameter optimization algorithm includes the Monte Carlo gradient approximation algorithm.

[0033] In some embodiments, the preset termination conditions include:

[0034] The bridging degree obtained from continuous X iterations no longer increases, where X is a positive integer; and / or

[0035] The iteration number threshold is reached.

[0036] On the other hand, the present invention proposes an optimization device for fracture design parameters of multiple horizontal wells based on the best bridging fracture network, including:

[0037] A morphology expansion module, which is used to expand the fracture network morphology of the target reservoir based on the given fracturing design parameters, and obtain the fracture point data of multiple fracturing horizontal wells in the target reservoir;

[0038] A fracture morphology determination module, which is used to determine the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well;

[0039] A bridging degree calculation module, which is used to calculate the initial bridging degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, wherein the bridging degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total untransformed volume and the total repeated transformation volume of the multiple fracturing horizontal wells;

[0040] A parameter optimization module, which is used to optimize the fracturing design parameters with the goal of maximizing the bridging degree of the fracture network of the multiple fracturing horizontal wells by using a preset parameter optimization algorithm according to the initial bridging degree, and obtain the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging degree.

[0041] In some embodiments, the fracturing design parameters include construction parameters and section cluster position parameters; wherein, the construction parameters include construction displacement and / or construction fluid volume.

[0042] In some embodiments, the fracture morphology of each fracturing horizontal well includes total well control volume, total effective transformation volume, total repeated transformation volume and / or total untransformed volume; the fracture morphology determination module is specifically used for:

[0043] Based on the actual production dynamic simulation, obtain the effective production range of each fracturing horizontal well;

[0044] According to the effective production range and the fracture point data of each fracturing horizontal well, determine the dynamic transformation volume of a single-section fracture of each fracturing horizontal well;

[0045] According to the dynamic transformation volume of a single-section fracture of each fracturing horizontal well, calculate the total effective transformation volume, total repeated transformation volume, total well control volume and / or total untransformed volume of each fracturing horizontal well based on the background grid method.

[0046] In some embodiments, the fracture morphology determination module calculates the total effective transformation volume, total repeated transformation volume, total well control volume and / or total untransformed volume of each fracturing horizontal well based on the background grid method according to the dynamic transformation volume of a single-section fracture of each fracturing horizontal well, including:

[0047] Mark the index information of each fracture of each fracturing horizontal well based on the background grid method;

[0048] According to the index information, by taking the union of the dynamic transformation volumes of each single-section fracture of all fracturing horizontal wells, obtain the total effective transformation volume of each fracturing horizontal well; and / or

[0049] According to the index information, by intersecting the dynamic transformation volumes of each single-stage fracture in each fractured horizontal well, the total inter-stage repeated transformation volume of each fractured horizontal well is obtained;

[0050] According to the index information, by intersecting the fracture dynamic transformation volumes of two fractured horizontal wells pairwise, the total inter-well repeated transformation volume of each fractured horizontal well is obtained;

[0051] According to the total inter-stage repeated transformation volume of each fractured horizontal well and the total inter-well repeated transformation volume of each fractured horizontal well, the total repeated transformation volume of each fractured horizontal well is determined; and / or

[0052] According to the index information, by taking the union of the well-controlled volumes of each fractured horizontal well, the total well-controlled volume of each fractured horizontal well is obtained;

[0053] According to the index information, by taking the union of the volumes of the untransformed grid cells of each fractured horizontal well, the total untransformed volume of each fractured horizontal well is obtained.

[0054] In some embodiments, the bridging degree calculation module is specifically configured to:

[0055] Calculate the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the ratio of the sum of the total repeated transformation volume and the total untransformed volume of each fractured horizontal well to the total well-controlled volume; or

[0056] Calculate the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the ratio of the total effective transformation volume of each fractured horizontal well to the total well-controlled volume.

[0057] In some embodiments, the parameter optimization module is specifically configured to:

[0058] Optimize the fracturing design parameters by using a preset parameter optimization algorithm;

[0059] Again, based on the optimized fracturing design parameters, perform fracture network morphology expansion to obtain the fracture point data of each fractured horizontal well;

[0060] According to the fracture point data of each fractured horizontal well, determine the fracture morphology of each fractured horizontal well;

[0061] According to the fracture morphology of each fractured horizontal well, calculate the current bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir;

[0062] If the current bridging degree is greater than the initial bridging degree, continue to use the preset parameter optimization algorithm to optimize the fracturing design parameters again on the basis of the optimized fracturing design parameters;

[0063] Continue to iterate until a preset termination condition is reached, and obtain the maximum bridging degree and the fracturing design parameters of the multi-well fracturing horizontal wells under the maximum bridging degree.

[0064] In some embodiments, the preset parameter optimization algorithm includes a Monte Carlo gradient approximation algorithm.

[0065] In some embodiments, the preset termination conditions include:

[0066] The bridging degree obtained from continuous X iterations no longer increases, where X is a positive integer; and / or

[0067] The iteration number threshold is reached.

[0068] On the other hand, the present invention provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the steps of the multi-horizontal well fracturing design parameter optimization method based on the best bridging fracture network described in any of the above embodiments are implemented.

[0069] On another hand, the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the multi-horizontal well fracturing design parameter optimization method based on the best bridging fracture network described in any of the above embodiments are implemented.

[0070] The multi-horizontal well fracturing design parameter optimization method and device based on the best bridging fracture network provided by the embodiments of the present invention expand the fracture network morphology of the target reservoir based on the given fracturing design parameters to obtain the fracture point data of multiple fracturing horizontal wells in the target reservoir; determine the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; calculate the initial bridging degree of the fracture network of the multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, where the bridging degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fracturing horizontal wells; optimize the fracturing design parameters with the maximum bridging degree of the fracture network of the multiple fracturing horizontal wells as the goal by using a preset parameter optimization algorithm according to the initial bridging degree, and obtain the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging degree. In this way, optimizing the fracturing design parameters with the maximization of the bridging degree of the fracture network of multiple fracturing horizontal wells as the goal, that is, optimizing the fracturing design parameters with the minimization of the total unmodified volume and the total repeated modification volume of multiple fracturing horizontal wells as the goal, realizes the full utilization of resources under the best benefit. Description of the Drawings

[0071] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can be obtained based on these drawings. In the accompanying drawings:

[0072] Figure 1 It is a schematic diagram of the optimized variables of the bridging fracture network of a fractured horizontal well obtained in an embodiment of the present invention.

[0073] Figure 2 It is a schematic flow chart of a method for optimizing the fracturing design parameters of multiple horizontal wells based on the optimal bridging fracture network provided by an embodiment of the present invention.

[0074] Figure 3 It is a schematic diagram of the fracture propagation mechanism provided by an embodiment of the present invention.

[0075] Figure 4 It is a schematic diagram of determining the fracture distribution according to the grid background method provided by an embodiment of the present invention.

[0076] Figure 5 It is a partial schematic flow chart of a method for optimizing the fracturing design parameters of multiple horizontal wells based on the optimal bridging fracture network provided by an embodiment of the present invention.

[0077] Figure 6 It is a schematic diagram of the fracture morphology and ESRV of each section of a fractured horizontal well provided by an embodiment of the present invention.

[0078] Figure 7 It is a schematic diagram of the effective pressure propagation range around the fracture provided by an embodiment of the present invention.

[0079] Figure 8 It is a schematic diagram of various volumes provided by an embodiment of the present invention.

[0080] Figure 9 It is a schematic flow chart of optimizing the fracturing design parameters by using a preset parameter optimization algorithm provided by an embodiment of the present invention.

[0081] Figure 10a and Figure 10b are respectively schematic diagrams of the permeability and porosity distribution fields of a coalbed methane reservoir.

[0082] Figure 11a and Figure 11b are respectively schematic diagrams of the perforation positions and fracturing fracture morphologies of two horizontal wells before and after optimization.

[0083] Figure 12 It is a schematic diagram of the change in the bridging degree during the optimization process of the fracturing design parameters in an embodiment of the present invention.

[0084] Figure 13 、 Figure 14 、 Figure 15 are the change data of the displacement, liquid volume of each section of two horizontal wells and the horizontal position of a certain cluster in each section during the iterative process, respectively.

[0085] Figure 16 is a schematic structural diagram of an optimization device for multi-horizontal well fracturing design parameters based on the best bridging fracture network provided by an embodiment of the present invention.

[0086] Figure 17 is a schematic physical structure diagram of an electronic device provided by an embodiment of the present invention. Specific Embodiments

[0087] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer and more understandable, the following further elaborates on the embodiments of the present invention with reference to the accompanying drawings. Herein, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but not to limit the present invention. It should be noted that, without conflict, the embodiments in the present invention and the features in the embodiments can be arbitrarily combined with each other.

[0088] To better understand the present invention, the following first details the research background of the present invention.

[0089] Facing two main problems in conventional multi-well fracturing, increasing the inter-stage repeated transformation volume and the inter-well repeated transformation volume has little impact on increasing the coalbed methane production, while reducing the volume of the untransformed area can significantly increase the coalbed methane production. Therefore, to achieve the full utilization of resources under the best benefit, with the goal of minimizing the untransformed volume and the repeated transformation volume, the maximum bridging degree of the multi-well fracture network can be ensured. Therefore, a method that can cover all possible situations as much as possible can be adopted in the setting of fracturing design parameters. Specifically, according to microseismic data, a fracture propagation simulator can be used to simulate different fracture network morphologies according to different construction parameters (such as construction displacement and construction liquid volume), and the effect of non-uniform fracture distribution in multiple wells can be achieved by changing the segment cluster position. Then, the fracturing design parameters are optimized integrally according to the intelligent optimization algorithm to achieve the purpose of the best bridging fracture network.

[0090] To achieve the goal of the best bridging fracture network, the well parameters and fracture parameters can be considered separately, that is, the construction displacement (dis), construction liquid volume (vol), segment length (spac), and cluster (perforation) position (loc) are used as independent optimization variables. In this way, it can ensure that these parameters are adjusted flexibly to meet different geological conditions and engineering requirements, thereby maximizing the oil and gas production benefit and avoiding unnecessary resource waste. Figure 1Shows a schematic diagram of optimizing variables for bridging fracture networks in fractured horizontal wells, where the section length refers to the length of each section, the cluster (perforation) position refers to the position of each fracture, and the displacement and liquid volume are variables that control the fracture network morphology.

[0091] In multi - well horizontal well fracturing, in order to achieve the ultimate utilization of oil and gas resources, while avoiding excessive inter - stage and inter - well repeated stimulation areas and un - stimulated areas caused by uniform conceptual design, the present invention proposes a method aiming at maximizing the bridging degree of fracture networks in multiple fractured horizontal wells (or minimizing the un - stimulated volume and repeated stimulation volume), which is achieved by constructing an objective function for the bridging degree of the fracture network. This objective function can reduce costs and increase efficiency in multi - well horizontal well fracturing development, reduce the area of un - stimulated areas within the well - controlled area, increase the oil (gas) drainage area of the reservoir, and at the same time reduce the inter - stage and inter - well repeated stimulation areas, thus achieving the ultimate utilization of oil and gas resources. This method takes into account the changes in geological conditions, the irregularity of the fracture network, and the influence of the fracture network morphology of individual horizontal wells and between wells on the bridging degree.

[0092] The basic idea of the Monte Carlo gradient approximation algorithm (MCGA) is as follows: First, use the Monte Carlo method to generate several realizations of random variables around the current optimal variables, then calculate the values of the performance index functions corresponding to each variable realization respectively, and finally use each realization and its performance index function value to estimate the gradient of the performance index function (such as the objective function of the bridging degree in the present invention).

[0093] For reservoir production optimization, let the optimal control variable obtained in the l - th iteration step be u l , and its corresponding performance index function be F(u l , y). According to the principle of the MCGA method, N l perturbation variables are generated around u r , and their expression is:

[0094] u l,j = u l + γ × r j , j = 1, 2, …, N r (1)

[0095] In the formula, u l,j represents the j - th control variable realization; γ is the perturbation constant; r j is the j - th perturbation vector, and its component r i j (i = 1, 2, …, N u ) is generally a variable conforming to the standard normal distribution, where N u represents the total number of control variables, that is, r ij ~N(0,1).

[0096] Substitute each implementation u l,j into the reservoir simulator, and the corresponding performance index F(u l,j , y) can be obtained through calculation. Let b j be the difference between F(u l,j , y) and the current optimal performance index F(u l,j , y), that is

[0097] b j = F(u l,j , y) - F(u l , y), j = 1, 2, …, N r (2)

[0098] Then, the gradient estimate l of the index function F(u l at u is obtained by applying the MCGA method, and the expression is:

[0099]

[0100] In the formula, represents the value of the i-th component. Next, the properties of will be analyzed.

[0101] Consider the first-order Taylor expansion of F(u l,j , y) at u l . Since the value of γ is usually small, then

[0102] F(u l,j , y) = F(u l , y) + γ(r j ) T g(u l ) (4)

[0103] In the formula, g(u l ) is the true gradient of F(u l , y) at u l . Substitute this formula into b j , then b j becomes

[0104]

[0105] where g s (u l ) represents the s-th component of the true gradient, s = 1, 2, …, N u . Substitute the above formula into the expression, and we have

[0106]

[0107] Taking the expectation value on both sides of the above equation, we can obtain

[0108]

[0109] Since the components in r j are random variables that conform to the standard normal distribution, they satisfy the following conditions:

[0110]

[0111] According to the above equation, the expectation value of

[0112]

[0113] Obviously, the expectation value of is the i-th component of the true gradient, so the expectation value of the gradient estimated based on the MCGA method is the true gradient

[0114]

[0115] where u l+1 is the control variable after iterative optimization; represents the infinity norm of

[0116] Figure 2 is the flow chart of the multi-horizontal well fracturing design parameter optimization method based on the best bridging fracture network provided by an embodiment of the present invention. As Figure 2 shown, the multi-horizontal well fracturing design parameter optimization method based on the best bridging fracture network provided by the embodiment of the present invention includes:

[0117] S101. Expand the fracture network morphology of the target reservoir based on the given fracturing design parameters to obtain the fracture point data of multiple fracturing horizontal wells in the target reservoir;

[0118] In step S101, based on the coalbed methane fracture propagation method, the expansion of the fracture network morphology with the given fracturing design parameters is realized; specifically, the coalbed methane fracture propagation simulation process is as follows:

[0119] The reservoir geological parameters, fracturing construction parameters, etc. are equivalent to potential parameters, and the mechanical problems in the process of fracture propagation are equivalent to the electrical problems in the process of lightning breakdown. The simulation idea is adopted to carry out fracture propagation simulation. The modified maximum circumferential tensile stress is used to judge the initiation of fracture: fractures can initiate at positions where the circumferential tensile stress is greater than the critical stress for initiation (i.e., fractures can propagate in multiple directions); the fractal index is introduced to calculate the probability distribution of fracture at each node at the fracture tip, and a random number is used to judge the fracture propagation direction. For the fracture propagation mechanism, see Figure 3 。

[0120] S102. Determine the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well;

[0121] In step S102, the fracture morphology can be characterized by a grid system, that is, the fracture point data in the meshless system is characterized by the background grid to represent its fracture morphology. For the determination of fracture distribution by the grid background method, see Figure 4 。

[0122] S103. Calculate the initial closure degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, where the closure degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fracturing horizontal wells;

[0123] In step S103, the closure degree of multiple fracturing horizontal wells is related to the fracture morphology of each fracturing horizontal well. Specifically, the closure degree of multiple fracturing horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fracturing horizontal wells, and is positively correlated with the total effective modification volume. Therefore, the initial closure degree of the fracture network of the multiple fracturing horizontal wells under the given fracturing design parameters can be calculated according to the fracture morphology of each fracturing horizontal well.

[0124] S104. According to the initial closure degree, use the preset parameter optimization algorithm to optimize the fracturing design parameters with the maximum closure degree of the fracture network of the multiple fracturing horizontal wells as the goal, and obtain the fracturing design parameters of the multiple fracturing horizontal wells under the maximum closure degree.

[0125] In step S104, the preset parameter optimization algorithm is used to iteratively optimize the fracturing design parameters. Each time the fracturing design parameters are optimized, calculate the closure degree of the fracture network of the multiple fracturing horizontal wells under the fracturing design parameters, and compare this closure degree with the closure degree of the fracture network of the multiple fracturing horizontal wells under the fracturing design parameters of the previous iteration. If the closure degree calculated this time is greater than the closure degree calculated in the previous iteration, continue to use the preset parameter optimization algorithm to optimize the fracturing design parameters until the iteration termination condition is reached, and obtain the maximum closure degree and the fracturing design parameters of the multiple fracturing horizontal wells under the maximum fitting degree.

[0126] The optimization method for fracturing design parameters of multi - horizontal wells based on the best - fitting fracture network provided by the present invention expands the fracture network morphology of the target reservoir based on the given fracturing design parameters to obtain the fracture point data of multiple fracturing horizontal wells in the target reservoir; determines the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; calculates the initial fitting degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, wherein the fitting degree of the fracture network of multiple fracturing horizontal wells is negatively correlated with the sum of the total untransformed volume and the total repeated transformation volume of multiple fracturing horizontal wells; optimizes the fracturing design parameters with the maximum fitting degree of the fracture network of multiple fracturing horizontal wells as the goal by using a preset parameter optimization algorithm, and obtains the fracturing design parameters of multiple fracturing horizontal wells under the maximum fitting degree. In this way, optimizing the fracturing design parameters with the maximum fitting degree of the fracture network of multiple fracturing horizontal wells as the goal, that is, optimizing the fracturing design parameters with the minimum sum of the total untransformed volume and the total repeated transformation volume of multiple fracturing horizontal wells as the goal, realizes the full utilization of resources under the best benefit.

[0127] In some embodiments, the fracturing design parameters include construction parameters and section - cluster position parameters; wherein, the construction parameters include construction displacement and / or construction fluid volume. Among them, according to the finally optimized section - cluster position, the section length, section spacing, cluster spacing, etc. can be correspondingly obtained to guide the on - site fracturing construction.

[0128] As Figure 5 shown, in some embodiments, the fracture morphology of each fracturing horizontal well includes total well - controlled volume, total effective transformation volume, total repeated transformation volume and / or total untransformed volume; the determining of the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well includes:

[0129] S1021. Obtaining the effective production range of each fracturing horizontal well through actual production dynamic simulation;

[0130] In step S1021, the stimulated reservoir volume (SRV) greatly enhances the original formation permeability and has a decisive impact on the production of production wells. SRV can be calculated through microseismic imaging, secondary fracture monitoring technology and mathematical models based on fracture propagation laws. However, there are significant differences between the SRV calculated by the above - mentioned methods and the effective stimulated volume (ESRV) estimated through dynamic production curves or actual oil and gas field production. Usually, SRV is much larger than ESRV. Compared with SRV, ESRV is a more important decisive condition for determining the success or failure of hydraulic fracturing and predicting the production improvement effect of fracturing. As Figure 6As shown, the effective stimulation volume (ESRV) is the volume of the reservoir fluid that can be effectively mobilized within the range affected by pressure after pressure propagation during pressure drawdown production, based on the fracture patterns in each stage. By establishing a matrix-fracture model, the effective pressure propagation range around the fractures during pressure drawdown production is simulated. As Figure 7 shown, the effective stimulation volume (ESRV) can be calculated by combining the fracture patterns. Specifically, the effective mobilization range of each fractured horizontal well can be obtained through dynamic simulation of actual production.

[0131] S1022. Determine the dynamic stimulation volume of the single-stage fractures of each fractured horizontal well according to the effective mobilization range and the fracture point data of each fractured horizontal well;

[0132] In step S1022, the irregular boundaries of each well are calculated based on the fracture point data, and the dynamic stimulation volume of each single-stage fracture is calculated based on the effective mobilization range.

[0133] S1023. Calculate the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume, and / or total non-stimulated volume of each fractured horizontal well based on the background grid method according to the dynamic stimulation volume of the single-stage fractures of each fractured horizontal well.

[0134] In step S1023, the index information of each fracture can be marked based on the background grid method. Specifically, in the background grid, 0 can be used to represent the non-stimulated volume, 1 can be used to represent the natural fracture-controlled volume and the artificial fracture-controlled volume, and 2 can be used to represent the inter-stage and inter-well repeated stimulation volume. By substituting the statistical position information of the index marks into the respective calculation formulas of the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume, and total non-stimulated volume, the above four volumes can be calculated respectively.

[0135] In some embodiments, the calculating the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume, and / or total non-stimulated volume of each fractured horizontal well based on the background grid method according to the dynamic stimulation volume of the single-stage fractures of each fractured horizontal well includes:

[0136] Mark the index information of each fracture of each fractured horizontal well based on the background grid method;

[0137] According to the index information, obtain the total effective stimulation volume of each fractured horizontal well by taking the union of the dynamic stimulation volumes of the single-stage fractures of all fractured horizontal wells; and / or

[0138] According to the index information, obtain the total inter-stage repeated stimulation volume of each fractured horizontal well by taking the intersection of the dynamic stimulation volumes of the single-stage fractures in each fractured horizontal well;

[0139] According to the index information, obtain the total inter-well repeated stimulation volume of each fractured horizontal well by taking the intersection of the fracture dynamic stimulation volumes of pairwise fractured horizontal wells;

[0140] Determine the total repeated stimulation volume of each fractured horizontal well according to the total inter-stage repeated stimulation volume and the total inter-well repeated stimulation volume of each fractured horizontal well; and / or

[0141] According to the index information, obtain the total well-controlled volume of each fractured horizontal well by taking the union of the well-controlled volumes of each fractured horizontal well;

[0142] According to the index information, obtain the total non-stimulated volume of each fractured horizontal well by taking the union of the non-stimulated grid cell volumes of each fractured horizontal well.

[0143] Specifically, according to the simulation results of the propagation of fracturing fractures, obtain the fracture point data of each fractured horizontal well, and then based on the actual production dynamic simulation, obtain the effective production range of each fractured horizontal well. Finally, determine the fracture distribution under the grid by the grid background method to calculate the effective stimulation volume, repeated stimulation volume (the sum of the inter-stage repeated stimulation volume and the inter-well repeated stimulation volume), well-controlled volume, and non-stimulated volume. Schematic diagrams of various volumes are shown in Figure 8 as shown. Among them, for the well-controlled volume V c it can be obtained by summing the well-controlled volumes of each fractured horizontal well:

[0144]

[0145] In the formula, V c is the total well-controlled volume; N is the number of well-controlled grids in the block; is the volume of the i-th unit grid; dx i is the grid step size of the i-th unit grid in the x direction; dy i is the grid step size of the i-th unit grid in the y direction; dz i is the grid step size of the i-th unit grid in the z direction.

[0146] By summing the control volumes with index information of 0 in each unit grid, the total non-stimulated volume V 0 can be obtained:

[0147]

[0148] In the formula, V 0 is the total non-stimulated volume; N w is the number of fractured horizontal wells; N u | λ＝0 (w i ) is the number of non-stimulated grid cells of the w i fractured horizontal well; is the volume of the i-th unit grid; dx i is the grid step size of the i-th unit grid in the x direction; dy i is the grid step size of the i-th unit grid in the y direction; dzi is the grid step size in the z direction of the i-th unit grid.

[0149] Based on the fracture propagation of each fractured horizontal well and the effective utilization range of the actual fractures, the effective control volume and the total effective control volume of each fractured horizontal well are calculated as follows:

[0150]

[0151]

[0152] In the formula, is for w i The control volume with the index information of 1 for the fractured horizontal well, that is, the effective control volume of a single well; V 1 is the control volume with the index information of 1 for all fractured horizontal wells, that is, the total effective control volume; uni(·) is the unique value function; N ut | λ＝1 is the total number of (effective) transformed grid cells; N u | λ＝1 (w i ) is the number of (effective) transformed grid cells of w i for the fractured horizontal well. N w is the number of fractured horizontal wells; dx i is the grid step size in the x direction of the i-th unit grid; dy i is the grid step size in the y direction of the i-th unit grid; dz i is the grid step size in the z direction of the i-th unit grid.

[0153] Furthermore, the total repeated transformation volume can be calculated based on the effective control volume of each fractured horizontal well:

[0154]

[0155] In the formula, V 2 is the total repeated transformation volume; N r | λ＝2 is the number of repeated transformation grid cells. However, when calculating the actual bridging degree, the total repeated transformation volume V can be calculated by taking the intersection of the effective transformation volumes of each fractured horizontal well 2 .

[0156] In some embodiments, calculating the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the fracture morphology of each fractured horizontal well includes:

[0157] Calculating the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the ratio of the sum of the total repeated transformation volume and the total untransformed volume of each fractured horizontal well to the total well control volume; or

[0158] Calculate the initial bridging degree of the fracture networks of multiple fractured horizontal wells in the target reservoir according to the ratio of the total effective stimulation volume to the total well control volume of each fractured horizontal well.

[0159] Specifically, the definition of the bridging degree of the fracture networks of multiple fractured horizontal wells can be as follows:

[0160]

[0161] In the formula, M is the bridging degree of the fracture network, dimensionless; N u | λ＝0 (w i ) is the total number of stage repeated stimulations and inter-well repeated stimulations of the fractured horizontal well; N i | ut is the total number of (effective) stimulated grid cells; uni(·) is the unique value function; N λ＝1 | rt is the total number of repeated stimulated grid cells; N is the number of well control grids in the block; for the definitions of other parameters, see above. λ＝2

[0162] It can be seen that the results of the total well control volume calculated by formula (11), the total non-stimulated volume calculated by formula (12), and the total repeated stimulation volume calculated by formula (15) can be substituted into formula (16), or the results of the total well control volume calculated by formula (11) and the total effective control volume calculated by formula (14) can be substituted into formula (16) to calculate the bridging degree of the fracture network. Both methods can calculate the bridging degree M under the current number of wells and fracture patterns.

[0163] In some embodiments, based on the initial bridging degree, using a preset parameter optimization algorithm, with the goal of maximizing the bridging degree of the fracture networks of the multiple fractured horizontal wells, optimizing the fracture design parameters, the fracture design parameters of the multiple fractured horizontal wells at the maximum bridging degree include:

[0164] Optimize the fracture design parameters using a preset parameter optimization algorithm;

[0165] Based on the optimized fracture design parameters, expand the fracture network morphology again to obtain the fracture point data of each fractured horizontal well;

[0166] According to the fracture point data of each fractured horizontal well, determine the fracture morphology of each fractured horizontal well;

[0167] According to the fracture morphology of each fractured horizontal well, calculate the current bridging degree of the fracture networks of the multiple fractured horizontal wells in the target reservoir;

[0168] If the current closing degree is greater than the initial closing degree, continue to use the preset parameter optimization algorithm to optimize the fracturing design parameters again based on the optimized fracturing design parameters;

[0169] Continue to iterate until a preset termination condition is reached, and obtain the maximum closing degree and the fracturing design parameters of the multiple fracturing horizontal wells at the maximum closing degree.

[0170] Specifically, as Figure 9 shown, the process of optimizing the fracturing design parameters using the preset parameter optimization algorithm is as follows:

[0171] Step 1: Input the initial displacement, liquid volume, and perforation positions of each fracturing horizontal well.

[0172] Step 2: Based on the fracture propagation inversion simulator, obtain the fracture morphology of each fracturing horizontal well, and determine the effective production range through actual production simulation. Next, calculate the effective stimulation volume, repeated stimulation volume, and well control volume, and then use Equation (16) to calculate the initial closing degree M 0 .

[0173] Step 3: Use the preset parameter optimization algorithm (such as the MCGA algorithm) to iteratively update the optimization variables.

[0174] Step 4: Call the fracture propagation simulator again to obtain the updated fracture morphology, and combine it with the effective production range to calculate the new closing degree M 1 .

[0175] Step 5: Judge that when M 1 > M 0 , then perform the next operation; otherwise, execute Step 3.

[0176] Step 6: Judge that when the termination condition is reached, output the optimal displacement, liquid volume, each perforation position, and the optimal closing degree M* that meet the conditions; otherwise, use the optimized fracture morphology as the initial fracture morphology, and assign the closing degree M 1 to M 0 , and execute Step 3.

[0177] In some embodiments, the preset termination conditions include: the closing degree obtained from continuous X iterations no longer increases, where X is a positive integer; and / or reaching the iteration number threshold.

[0178] To better understand the present invention, the following uses a specific embodiment to illustrate the multi-horizontal well fracturing design parameter optimization method based on the best closed fracture network provided by the present invention.

[0179] Taking two horizontal wells (JS8-8P01 and JS8-8P02) in a heterogeneous coalbed methane reservoir as the actual model, the fracture pattern obtained from the current fracturing scheme is used as the initial scheme. The MCGA algorithm is used to optimize the displacement, liquid volume, and cluster position. The fracture is simulated using an extended numerical simulator for the fracture pattern, and the bridging degree of multiple horizontal wells is calculated by combining the effective utilization range and the background grid method. In addition, the displacement, liquid volume, and cluster position are continuously iteratively updated through the MCGA algorithm to continuously find the optimal bridging degree.

[0180] The permeability and porosity distribution fields of this coalbed methane reservoir are as Figure 10a and Figure 10b shown. The maximum permeability value of this reservoir is 0.1 mD, the minimum permeability value is 0.01 mD, and the average permeability value is 0.05 mD. The maximum porosity is 0.05, the minimum porosity is 0.01, and the average porosity is 0.036. The reservoir parameters and simulation information in the actual example are shown in Table 1.

[0181] Table 1 Reservoir parameters and simulation information table in the conceptual example

[0182]

[0183] Taking the maximum bridging degree of multi-well fractured horizontal wells as the optimization objective function, the MCGA algorithm is used to perturb and iteratively update each displacement, liquid volume, and cluster position to achieve the integrated optimization of the optimization parameters and obtain the fracture deployment scheme with the current optimal bridging degree. Figure 11a and Figure 11b show the comparison of the fracture network morphology before and after optimization during the optimization process of displacement, liquid volume, and cluster position. From Figure 11a and Figure 11b it can be clearly seen that the horizontal positions of each fracture change after optimization and the fracture lengths generally show an increasing trend. From Figure 12 it can be seen that the convergence rate of the multi-well fracture network optimization based on MCGA is relatively fast, starting to converge around the 15th step, and the final bridging degree is 0.5638, with an increase in the bridging degree of 40.95%. Figure 13 、 Figure 14 、 Figure 15 respectively show the changes in the displacement, liquid volume, and cluster position of each section of the two horizontal wells during the iterative process. Taking the optimization result of the first cluster of each section as an example, when the fracture position is optimized, the section length, section spacing, and cluster spacing between each cluster can be obtained.

[0184] Figure 16 is the structural schematic diagram of the multi-horizontal well fracturing design parameter optimization device based on the optimal bridging fracture network provided by an embodiment of the present invention, as Figure 16As shown in the figure, the multi-horizontal well fracturing design parameter optimization device based on the best bridging fracture network provided by the embodiment of the present invention includes:

[0185] A morphology expansion module 21, configured to perform fracture network morphology expansion on a target reservoir based on given fracturing design parameters to obtain fracture point data of multiple fracturing horizontal wells in the target reservoir;

[0186] A fracture morphology determination module 22, configured to determine the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well;

[0187] A bridging degree calculation module 23, configured to calculate the initial bridging degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, wherein the bridging degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fracturing horizontal wells;

[0188] A parameter optimization module 24, configured to optimize the fracturing design parameters with the maximum bridging degree of the fracture network of the multiple fracturing horizontal wells as the target by using a preset parameter optimization algorithm according to the initial bridging degree, so as to obtain the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging degree.

[0189] The multi-horizontal well fracturing design parameter optimization device based on the best bridging fracture network provided by the embodiment of the present invention performs fracture network morphology expansion on a target reservoir based on given fracturing design parameters to obtain fracture point data of multiple fracturing horizontal wells in the target reservoir; determines the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; calculates the initial bridging degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, wherein the bridging degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fracturing horizontal wells; and optimizes the fracturing design parameters with the maximum bridging degree of the fracture network of the multiple fracturing horizontal wells as the target by using a preset parameter optimization algorithm according to the initial bridging degree, so as to obtain the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging degree. In this way, the fracturing design parameters are optimized with the maximum bridging degree of the fracture network of multiple fracturing horizontal wells as the target, that is, the fracturing design parameters are optimized with the minimum total unmodified volume and total repeated modification volume of multiple fracturing horizontal wells as the target, realizing the full utilization of resources under the best benefit.

[0190] In some embodiments, the fracturing design parameters include construction parameters and stage cluster position parameters; wherein, the construction parameters include construction displacement and / or construction fluid volume.

[0191] In some embodiments, the fracture morphology of each fracturing horizontal well includes the total well control volume, the total effective stimulation volume, the total repeated stimulation volume, and / or the total un-stimulated volume; specifically, the fracture morphology determination module is configured to:

[0192] Obtain the effective production range of each fracturing horizontal well based on actual production dynamic simulation;

[0193] Determine the dynamic stimulation volume of a single-stage fracture of each fracturing horizontal well according to the effective production range and the fracture point data of each fracturing horizontal well;

[0194] Calculate the total effective stimulation volume, the total repeated stimulation volume, the total well control volume, and / or the total un-stimulated volume of each fracturing horizontal well based on the background grid method according to the dynamic stimulation volume of a single-stage fracture of each fracturing horizontal well.

[0195] In some embodiments, the fracture morphology determination module calculates the total effective stimulation volume, the total repeated stimulation volume, the total well control volume, and / or the total un-stimulated volume of each fracturing horizontal well based on the background grid method according to the dynamic stimulation volume of a single-stage fracture of each fracturing horizontal well, including:

[0196] Mark the index information of each fracture of each fracturing horizontal well based on the background grid method;

[0197] According to the index information, obtain the total effective stimulation volume of each fracturing horizontal well by taking the union of the dynamic stimulation volumes of single-stage fractures of all fracturing horizontal wells; and / or

[0198] According to the index information, obtain the total inter-stage repeated stimulation volume of each fracturing horizontal well by taking the intersection of the dynamic stimulation volumes of single-stage fractures in each fracturing horizontal well;

[0199] According to the index information, obtain the total inter-well repeated stimulation volume of each fracturing horizontal well by taking the intersection of the fracture dynamic stimulation volumes of two-by-two fracturing horizontal wells;

[0200] Determine the total repeated stimulation volume of each fracturing horizontal well according to the total inter-stage repeated stimulation volume of each fracturing horizontal well and the total inter-well repeated stimulation volume of each fracturing horizontal well; and / or

[0201] According to the index information, obtain the total well control volume of each fracturing horizontal well by taking the union of the well control volumes of each fracturing horizontal well;

[0202] According to the index information, obtain the total un-stimulated volume of each fracturing horizontal well by taking the union of the un-stimulated grid cell volumes of each fracturing horizontal well.

[0203] In some embodiments, the bridging degree calculation module is specifically configured to:

[0204] Calculate the initial bridging degree of the fracture networks of multiple fractured horizontal wells in the target reservoir according to the ratio of the sum of the total repeated stimulation volume and the total non-stimulated volume of each fractured horizontal well to the total well control volume; or

[0205] Calculate the initial bridging degree of the fracture networks of multiple fractured horizontal wells in the target reservoir according to the ratio of the total effective stimulation volume of each fractured horizontal well to the total well control volume.

[0206] In some embodiments, the parameter optimization module is specifically configured to:[[]]

[0207] Optimize the fracturing design parameters by using a preset parameter optimization algorithm;

[0208] Perform fracture network morphology expansion again based on the optimized fracturing design parameters to obtain the fracture point data of each fractured horizontal well;

[0209] Determine the fracture morphology of each fractured horizontal well according to the fracture point data of each fractured horizontal well;

[0210] Calculate the current bridging degree of the fracture networks of multiple fractured horizontal wells in the target reservoir according to the fracture morphologies of each fractured horizontal well;

[0211] If the current bridging degree is greater than the initial bridging degree, continue to optimize the fracturing design parameters again based on the optimized fracturing design parameters by using the preset parameter optimization algorithm;

[0212] Continue the iteration until a preset termination condition is reached, and obtain the maximum bridging degree and the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree.

[0213] In some embodiments, the preset parameter optimization algorithm includes the Monte Carlo gradient approximation algorithm.

[0214] In some embodiments, the preset termination conditions include:

[0215] The bridging degree obtained from X consecutive iterations no longer increases, where X is a positive integer; and / or

[0216] The iteration number threshold is reached.

[0217] The embodiments of the device provided by the embodiments of the present invention can specifically be used to execute the processing procedures of the above method embodiments, and their functions will not be elaborated here. Reference can be made to the detailed descriptions of the above method embodiments.

[0218] Figure 17 It is a schematic physical structure diagram of an electronic device provided by an embodiment of the present invention, as Figure 17As shown in the figure, the electronic device 600 may include: a processor 100 and a memory 140. The memory 140 is coupled to the processor 100. The processor 100 may call the logical instructions in the memory 140 to execute the following method: by expanding the fracture network morphology of the target reservoir based on given fracturing design parameters, obtaining fracture point data of multiple fracturing horizontal wells in the target reservoir; determining the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; calculating the initial bridging degree of the fracture network of the multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, wherein the bridging degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fracturing horizontal wells; according to the initial bridging degree, using a preset parameter optimization algorithm to optimize the fracturing design parameters with the goal of maximizing the bridging degree of the fracture network of the multiple fracturing horizontal wells, and obtaining the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging degree.

[0219] This embodiment discloses a computer program product. The computer program product includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the methods provided in the above method embodiments, for example, including: by expanding the fracture network morphology of the target reservoir based on given fracturing design parameters, obtaining fracture point data of multiple fracturing horizontal wells in the target reservoir; determining the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; calculating the initial bridging degree of the fracture network of the multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, wherein the bridging degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fracturing horizontal wells; according to the initial bridging degree, using a preset parameter optimization algorithm to optimize the fracturing design parameters with the goal of maximizing the bridging degree of the fracture network of the multiple fracturing horizontal wells, and obtaining the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging degree.

[0220] This embodiment provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and the computer program causes the computer to execute the methods provided in the above method embodiments. For example, it includes: expanding the fracture network morphology of a target reservoir based on given fracturing design parameters to obtain fracture point data of multiple fracturing horizontal wells in the target reservoir; determining the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; calculating the initial bridging degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, where the bridging degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total untransformed volume and the total repeated transformation volume of the multiple fracturing horizontal wells; optimizing the fracturing design parameters with the maximum bridging degree of the fracture network of the multiple fracturing horizontal wells as the goal by using a preset parameter optimization algorithm according to the initial bridging degree, so as to obtain the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging degree.

[0221] As Figure 17 shown, the electronic device 600 may further include: a communication module 110, an input unit 120, an audio processor 130, a display 160, and a power supply 170. It should be noted that the electronic device 600 does not necessarily have to include Figure 17 all the components shown in Figure 17 ; in addition, the electronic device 600 may further include Figure 17 components not shown in Figure 17 . Reference may be made to the prior art. It should be noted that this figure is exemplary; other types of structures may also be used to supplement or replace this structure to achieve telecommunication functions or other functions.

[0222] As Figure 17 shown, the processor 100 is sometimes also referred to as a controller or an operation control unit. It may include a microprocessor or other processor devices and / or logic devices. The processor 100 receives inputs and controls the operations of the various components of the electronic device 600.

[0223] Among them, the memory 140 may be, for example, one or more of a buffer, a flash memory, a hard drive, a removable medium, a volatile memory, a non-volatile memory, or other suitable devices. It can store the above information related to failures, and can also store programs for executing relevant information. And the processor 100 can execute the program stored in the memory 140 to achieve information storage or processing, etc.

[0224] The input unit 120 provides inputs to the processor 100. The input unit 120 is, for example, a key or a touch input device. The power supply 170 is used to supply power to the electronic device 600. The display 160 is used for displaying display objects such as images and texts. The display 160 may be, for example, an LCD display, but is not limited thereto.

[0225] The memory 140 can be a solid-state memory, for example, a read-only memory (ROM), a random access memory (RAM), a SIM card, etc. It can also be a memory that stores information even when powered off, can be selectively erased and has more data. Examples of the memory 140 are sometimes referred to as EPROM, etc. The memory 140 can also be some other type of device. The memory 140 includes a buffer 141 (sometimes referred to as a buffer memory). The memory 140 can include an application / function storage unit 142, which is used to store application programs and function programs or the processes for operating the electronic device 600 through the processor 100.

[0226] The memory 140 can also include a data storage unit 143, which is used to store data, such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. The driver storage unit 144 of the memory 140 can include various drivers for the communication functions of the electronic device and / or for performing other functions of the electronic device (such as a messaging application, an address book application, etc.).

[0227] The communication module 110 includes a transmitter / receiver that transmits and receives signals via the antenna 111. The communication module 110 is coupled to the processor 100 to provide input signals and receive output signals, which can be the same as in the case of a conventional mobile communication terminal.

[0228] Based on different communication technologies, multiple communication modules 110 can be provided in the same electronic device, such as a cellular network module, a Bluetooth module, and / or a wireless local area network module, etc. The communication module 110 is also coupled to the speaker 131 and the microphone 132 via the audio processor 130 to provide an audio output via the speaker 131 and receive an audio input from the microphone 132, so as to achieve normal telecommunication functions. The audio processor 130 can include any suitable buffer, decoder, amplifier, etc. In addition, the audio processor 130 is also coupled to the processor 100, so that recording can be performed on the local machine through the microphone 132, and the sound stored on the local machine can be played through the speaker 131.

[0229] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk memories, CD-ROMs, optical memories, etc.) that contain computer-usable program codes.

[0230] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram, and the combination of flows and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing device produce means for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or means for implementing the functions specified in one block or multiple blocks.

[0231] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufactured article including instruction means for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or means for implementing the functions specified in one block or multiple blocks.

[0232] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or means for implementing the functions specified in one block or multiple blocks.

[0233] In the description of this specification, the description with reference to terms such as "one embodiment", "a specific embodiment", "some embodiments", "for example", "example", "specific example", or "some examples", etc. means that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0234] The above specific embodiments have further elaborated on the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above are only specific embodiments of the present invention and are not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. An optimization method for fracturing design parameters of multi - horizontal wells based on the best - fitting fracture network, characterized in that, it includes: Expanding the fracture network morphology of the target reservoir based on the given fracturing design parameters to obtain the fracture point data of multiple fracturing horizontal wells in the target reservoir; Determining the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; Calculating the initial fitting degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well, wherein the fitting degree of the fracture network of the multiple fracturing horizontal wells is negatively correlated with the sum of the total untransformed volume and the total repeated transformation volume of the multiple fracturing horizontal wells; According to the initial fitting degree, using a preset parameter optimization algorithm to optimize the fracturing design parameters with the goal of maximizing the fitting degree of the fracture network of the multiple fracturing horizontal wells, and obtaining the fracturing design parameters of the multiple fracturing horizontal wells under the maximum fitting degree.

2. The method according to claim 1, characterized in that, the fracturing design parameters include construction parameters and section - cluster position parameters; wherein, the construction parameters include construction displacement and / or construction fluid volume.

3. The method according to claim 2, characterized in that, the fracture morphology of each fracturing horizontal well includes total well - controlled volume, total effective transformation volume, total repeated transformation volume and / or total untransformed volume; The determining the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well includes: Based on actual production dynamic simulation, obtaining the effective production range of each fracturing horizontal well; According to the effective production range and the fracture point data of each fracturing horizontal well, determining the dynamic transformation volume of a single - stage fracture of each fracturing horizontal well; Based on the dynamic transformation volume of a single - stage fracture of each fracturing horizontal well, calculating the total effective transformation volume, total repeated transformation volume, total well - controlled volume and / or total untransformed volume of each fracturing horizontal well based on the background grid method.

4. The method according to claim 3, characterized in that, the calculating the total effective transformation volume, total repeated transformation volume, total well - controlled volume and / or total untransformed volume of each fracturing horizontal well based on the background grid method according to the dynamic transformation volume of a single - stage fracture of each fracturing horizontal well includes: Marking the index information of each fracture of each fracturing horizontal well based on the background grid method; According to the index information, obtaining the total effective transformation volume of each fracturing horizontal well by taking the union of the dynamic transformation volumes of single - stage fractures of all fracturing horizontal wells; and / or According to the index information, obtaining the total inter - stage repeated transformation volume of each fracturing horizontal well by taking the intersection of the dynamic transformation volumes of single - stage fractures in each fracturing horizontal well; According to the index information, obtaining the total inter - well repeated transformation volume of each fracturing horizontal well by taking the intersection of the fracture dynamic transformation volumes of two - by - two fracturing horizontal wells; Determining the total repeated transformation volume of each fracturing horizontal well according to the total inter - stage repeated transformation volume of each fracturing horizontal well and the total inter - well repeated transformation volume of each fracturing horizontal well; and / or According to the index information, obtaining the total well - controlled volume of each fracturing horizontal well by taking the union of the well - controlled volumes of each fracturing horizontal well. According to the index information, the total unmodified volume of each fractured horizontal well is obtained by taking the union of the unmodified grid cell volumes of each fractured horizontal well.

5. The method according to claim 1, wherein, the calculating of the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the fracture morphology of each fractured horizontal well includes: calculating the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the ratio of the sum of the total repeated modification volume and the total unmodified volume of each fractured horizontal well to the total well-controlled volume; or calculating the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the ratio of the total effective modification volume of each fractured horizontal well to the total well-controlled volume.

6. The method according to claim 1, wherein, the optimizing of the fracturing design parameters with the maximum bridging degree of the fracture network of the multiple fractured horizontal wells as the target according to the initial bridging degree by using a preset parameter optimization algorithm includes: optimizing the fracturing design parameters by using a preset parameter optimization algorithm; performing fracture network morphology expansion again based on the optimized fracturing design parameters to obtain the fracture point data of each fractured horizontal well; determining the fracture morphology of each fractured horizontal well according to the fracture point data of each fractured horizontal well; calculating the current bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the fracture morphology of each fractured horizontal well; if the current bridging degree is greater than the initial bridging degree, continue to optimize the fracturing design parameters again based on the optimized fracturing design parameters by using a preset parameter optimization algorithm; continue to iterate until a preset termination condition is reached, to obtain the maximum bridging degree and the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree.

7. The method according to claim 6, wherein, the preset parameter optimization algorithm includes the Monte Carlo gradient approximation algorithm.

8. The method according to claim 6, wherein, the preset termination condition includes: the bridging degree obtained by continuous X iterations no longer increases, where X is a positive integer; and / or reaching the iteration number threshold.

9. A device for optimizing the fracturing design parameters of multiple horizontal wells based on the best-bridged fracture network, wherein, it includes: a morphology expansion module, configured to perform fracture network morphology expansion on the target reservoir based on given fracturing design parameters to obtain the fracture point data of multiple fractured horizontal wells in the target reservoir; a fracture morphology determination module, configured to determine the fracture morphology of each fractured horizontal well according to the fracture point data of each fractured horizontal well; a bridging degree calculation module, configured to calculate the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the fracture morphology of each fractured horizontal well, wherein the bridging degree of the fracture network of the multiple fractured horizontal wells is negatively correlated with the sum of the total unmodified volume and the total repeated modification volume of the multiple fractured horizontal wells. A parameter optimization module, configured to optimize the fracturing design parameters according to the initial bridging degree by using a preset parameter optimization algorithm with the goal of maximizing the bridging degree of the fracture network of the multi-well fractured horizontal well, so as to obtain the fracturing design parameters of the multi-well fractured horizontal well under the maximum bridging degree.

10. The apparatus according to claim 9, wherein, the fracturing design parameters include construction parameters and section cluster position parameters; wherein, the construction parameters include construction displacement and / or construction fluid volume.

11. The apparatus according to claim 10, wherein, the fracture patterns of the fractured horizontal wells include total well-controlled volume, total effective stimulation volume, total repeated stimulation volume and / or total un-stimulated volume; the fracture pattern determination module is specifically configured to: obtain the effective production range of each fractured horizontal well based on actual production dynamic simulation; determine the dynamic stimulation volume of the single-section fracture of each fractured horizontal well according to the effective production range and the fracture point data of each fractured horizontal well; calculate the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume and / or total un-stimulated volume of each fractured horizontal well based on the background grid method according to the dynamic stimulation volume of the single-section fracture of each fractured horizontal well.

12. The apparatus according to claim 11, wherein, the fracture pattern determination module calculates the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume and / or total un-stimulated volume of each fractured horizontal well based on the background grid method according to the dynamic stimulation volume of the single-section fracture of each fractured horizontal well, including: marking the index information of each fracture of each fractured horizontal well based on the background grid method; obtaining the total effective stimulation volume of each fractured horizontal well by taking the union of the dynamic stimulation volumes of the single-section fractures of all fractured horizontal wells according to the index information; and / or obtaining the total inter-section repeated stimulation volume of each fractured horizontal well by taking the intersection of the dynamic stimulation volumes of the single-section fractures in each fractured horizontal well according to the index information; obtaining the total inter-well repeated stimulation volume of each fractured horizontal well by taking the intersection of the fracture dynamic stimulation volumes of two-by-two fractured horizontal wells according to the index information; determining the total repeated stimulation volume of each fractured horizontal well according to the total inter-section repeated stimulation volume of each fractured horizontal well and the total inter-well repeated stimulation volume of each fractured horizontal well; and / or obtaining the total well-controlled volume of each fractured horizontal well by taking the union of the well-controlled volumes of each fractured horizontal well according to the index information; obtaining the total un-stimulated volume of each fractured horizontal well by taking the union of the un-stimulated grid unit volumes of each fractured horizontal well according to the index information.

13. The apparatus according to claim 9, wherein, the bridging degree calculation module is specifically configured to: calculate the initial bridging degree of the fracture network of the multi-well fractured horizontal wells in the target reservoir according to the ratio of the sum of the total repeated stimulation volume and the total un-stimulated volume of each fractured horizontal well to the total well-controlled volume; or calculate the initial bridging degree of the fracture network of the multi-well fractured horizontal wells in the target reservoir according to the ratio of the total effective stimulation volume of each fractured horizontal well to the total well-controlled volume.

14. The apparatus according to claim 9, wherein, the parameter optimization module is specifically configured to: Optimize the fracturing design parameters by using a preset parameter optimization algorithm; Based on the optimized fracturing design parameters, expand the fracture network morphology again to obtain the fracture point data of each fracturing horizontal well; Determine the fracture morphology of each fracturing horizontal well according to the fracture point data of each fracturing horizontal well; Calculate the current closing degree of the fracture network of multiple fracturing horizontal wells in the target reservoir according to the fracture morphology of each fracturing horizontal well; If the current closing degree is greater than the initial closing degree, continue to use the preset parameter optimization algorithm to optimize the fracturing design parameters again based on the optimized fracturing design parameters; Continue to iterate until a preset termination condition is reached, and obtain the maximum closing degree and the fracturing design parameters of the multiple fracturing horizontal wells under the maximum closing degree.

15. The device according to claim 14, wherein, the preset parameter optimization algorithm includes a Monte Carlo gradient approximation algorithm.

16. The device according to claim 14, wherein, the preset termination condition includes: the closing degree obtained by continuous X iterations no longer increases, where X is a positive integer; and / or reaching the iteration number threshold.

17. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein, when the processor executes the computer program, the steps of the method according to any one of claims 1 to 8 are implemented.

18. A computer-readable storage medium, on which a computer program is stored, wherein, when the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 8 are implemented.

Images