Method and apparatus for optimizing fracture design parameters of multiple horizontal wells on basis of closed fracture network

Abstract

A method and apparatus for optimizing fracture design parameters of multiple horizontal wells on the basis of a closed fracture network. The method comprises: performing fracture network form expansion on a target reservoir on the basis of given fracture design parameters, so as to obtain fracture point data of multiple fracture horizontal wells of the target reservoir (S101); on the basis of the fracture point data of the fracture horizontal wells, determining fracture forms of the fracture horizontal wells (S102); on the basis of the fracture forms of the fracture horizontal wells, calculating an initial degree of closure of a fracture network of the multiple fracture horizontal wells of the target reservoir, wherein the degree of closure of the fracture network of the multiple fracture horizontal wells is negatively correlated with the sum of a total unmodified volume and a total repeatedly modified volume of the multiple fracture horizontal wells (S103); and on the basis of the initial degree of closure, optimizing the fracture design parameters by using a preset parameter optimization algorithm and taking the maximization of the degree of closure of the fracture network of the multiple fracture horizontal wells as an objective, so as to obtain fracture design parameters of the multiple fracture horizontal wells at the maximum degree of closure (S104).

Description

Multi-horizontal well fracturing design parameter optimization method and device based on bridging fracture network

[0001] Related applications

[0002] This application claims priority to Chinese patent application No. 202311620079.3 filed on November 30, 2023, and cites the contents disclosed in the above patent application as part of this application. Technical Field

[0003] The present application relates to the technical field of oil and gas reservoir production optimization, and in particular to a method and device for optimizing design parameters of multi-horizontal well fracturing based on a bridging fracture network. Background Art

[0004] Hydraulic fracturing plays a vital role in the development of unconventional oil and gas. In recent years, significant breakthroughs have been made in the exploration and development of deep coalbed methane. Horizontal well fracturing has become a crucial technology for coalbed methane, particularly deep coalbed methane development. Conventional fracturing, however, has two major issues: first, it targets only a single well and fails to consider the interactive effects of the fracturing network across multiple wells. Second, there are large unreformed areas within the well-control range of multiple wells. Third, under multi-well conditions, there is often a high incidence of repeated re-reforms between segments and between wells. Fourth, it is impossible to simultaneously optimize segment cluster deployment design parameters (segment length, perforation location, fracture half-length) and fracturing operation parameters (displacement rate, fluid volume, etc.). Therefore, large-scale fracturing is subject to issues of either over-reform (resulting in cost waste) or localized reformation blanks (resulting in unused resources), posing challenges to efficient development.

[0005] Summary of the Invention

[0006] In response to the problems in the prior art, an embodiment of the present application provides a method for optimizing design parameters of multi-horizontal well fracturing based on bridging fracture networks.

[0007] On the one hand, the present application proposes a method for optimizing design parameters of multi-horizontal well fracturing based on bridging fracture networks, comprising:

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

[0009] Determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0010] Based on the fracture morphology of each fractured horizontal well, the initial healing degree of the fracture network of multiple fractured horizontal wells in the target reservoir is calculated. The healing degree of the fracture network of multiple fractured horizontal wells is negatively correlated with the sum of the total un-fractured volume and the total repeatedly-fractured volume of the multiple fractured horizontal wells.

[0011] According to the initial bridging degree, the preset parameter optimization algorithm is used to optimize the fracturing design parameters with the goal of maximizing the bridging degree of the fracture network of multiple fractured horizontal wells, and the fracturing design parameters of multiple fractured horizontal wells under the maximum bridging degree are obtained.

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

[0013] In some embodiments, the fracture morphology of each fractured horizontal well includes the total well-controlled volume, the total effective stimulated volume, the total repeatedly stimulated volume, and / or the total unstimulated volume. Determining the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well includes:

[0014] Based on actual production dynamic simulation, the effective producing range of each fractured horizontal well is obtained;

[0015] Determine the dynamic stimulation volume of a single fracture segment in each fractured horizontal well based on the effective producing range and fracture point data of each fractured horizontal well;

[0016] According to the dynamic stimulation volume of a single fracture segment of each fractured horizontal well, the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume and / or total unstimulated volume of each fractured horizontal well are calculated based on the background grid method.

[0017] In some embodiments, according to the dynamic stimulation volume of a single fracture segment of each fractured horizontal well, calculating the total effective stimulation volume, the total repeated stimulation volume, the total well-controlled volume and / or the total unstimulated volume of each fractured horizontal well based on the background grid method includes:

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

[0019] According to the index information, the total effective stimulation volume of each fractured horizontal well is obtained by combining the dynamic stimulation volumes of each single fracture segment of all fractured horizontal wells; and / or

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

[0021] According to the index information, the total inter-well repeated stimulation volume of each fractured horizontal well is obtained by intersecting the fracture dynamic stimulation volumes of each fractured horizontal well.

[0022] Determine the total repetitive stimulation volume of each fractured horizontal well based on the total inter-stage repetitive stimulation volume of each fractured horizontal well and the total inter-well repetitive stimulation volume of each fractured horizontal well; and / or

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

[0024] According to the index information, the total unstimulated volume of each fractured horizontal well is obtained by combining the unstimulated grid cell volumes of each fractured horizontal well.

[0025] In some embodiments, calculating the initial healing degree of a fracture network of multiple fractured horizontal wells in a target reservoir based on the fracture morphology of each fractured horizontal well includes:

[0026] Calculate the initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir based on the ratio of the sum of the total re-stimulated volume and the total unstimulated volume of each fractured horizontal well to the total well-controlled volume; or

[0027] The initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir is calculated based on the ratio of the total effective stimulation volume of each fractured horizontal well to the total well-controlled volume.

[0028] In some embodiments, based on the initial bridging degree, a preset parameter optimization algorithm is used to optimize the fracturing design parameters with the goal of maximizing the bridging degree of the fracture network of multiple fractured horizontal wells. The fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree include:

[0029] Optimize fracturing design parameters using preset parameter optimization algorithms;

[0030] The fracture network morphology is expanded again based on the optimized fracturing design parameters to obtain the fracture point data of each fracturing horizontal well;

[0031] Determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0032] Calculate the current degree of healing of the fracture network of multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well;

[0033] If the current degree of bridging is greater than the initial degree of bridging, the preset parameter optimization algorithm is continued to be used to optimize the fracturing design parameters again based on the optimized fracturing design parameters;

[0034] The iteration is continued until the preset termination condition is reached, and the maximum bridging degree and the fracturing design parameters of multiple fractured horizontal wells under the maximum bridging degree are obtained.

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

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

[0037] The degree of bridging obtained after X consecutive iterations does not increase, where X is a positive integer; and / or

[0038] The iteration threshold is reached.

[0039] On the other hand, the present application proposes a multi-horizontal well fracturing design parameter optimization device based on bridging fracture networks, comprising:

[0040] The morphology expansion module is used to expand the fracture network morphology of the target reservoir based on the given fracturing design parameters to obtain the fracture point data of multiple fractured horizontal wells in the target reservoir;

[0041] A fracture morphology determination module is used to determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0042] The healing degree calculation module is used to calculate the initial healing degree of the fracture network of multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well. The healing degree of the fracture network of multiple fractured horizontal wells is negatively correlated with the sum of the total un-reformed volume and the total repeatedly re-reformed volume of the multiple fractured horizontal wells.

[0043] The parameter optimization module is used to optimize the fracturing design parameters based on the initial bridging degree and the preset parameter optimization algorithm with the goal of maximizing the bridging degree of the fracture network of multiple fractured horizontal wells, and obtain the fracturing design parameters of multiple fractured horizontal wells under the maximum bridging degree.

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

[0045] In some embodiments, the fracture morphology of each fractured horizontal well includes the total well-controlled volume, the total effective stimulation volume, the total repeated stimulation volume and / or the total unstimulated volume; the fracture morphology determination module is specifically configured to:

[0046] Based on actual production dynamic simulation, the effective producing range of each fractured horizontal well is obtained;

[0047] Determine the dynamic stimulation volume of a single fracture segment in each fractured horizontal well based on the effective producing range and fracture point data of each fractured horizontal well;

[0048] According to the dynamic stimulation volume of a single fracture segment of each fractured horizontal well, the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume and / or total unstimulated volume of each fractured horizontal well are calculated based on the background grid method.

[0049] In some embodiments, the fracture morphology determination module calculates the total effective stimulated volume, the total repeated stimulated volume, the total well-controlled volume, and / or the total unstimulated volume of each fractured horizontal well based on the background grid method according to the dynamic stimulated volume of a single fracture segment of each fractured horizontal well, including:

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

[0051] According to the index information, the total effective stimulation volume of each fractured horizontal well is obtained by combining the dynamic stimulation volumes of each single fracture segment of all fractured horizontal wells; and / or

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

[0053] According to the index information, the total inter-well repeated stimulation volume of each fractured horizontal well is obtained by intersecting the fracture dynamic stimulation volumes of each fractured horizontal well.

[0054] Determine the total repetitive stimulation volume of each fractured horizontal well based on the total inter-stage repetitive stimulation volume of each fractured horizontal well and the total inter-well repetitive stimulation volume of each fractured horizontal well; and / or

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

[0056] According to the index information, the total unstimulated volume of each fractured horizontal well is obtained by combining the unstimulated grid cell volumes of each fractured horizontal well.

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

[0058] Calculate the initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir based on the ratio of the sum of the total re-stimulated volume and the total unstimulated volume of each fractured horizontal well to the total well-controlled volume; or

[0059] The initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir is calculated based on the ratio of the total effective stimulation volume of each fractured horizontal well to the total well-controlled volume.

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

[0061] Optimize fracturing design parameters using preset parameter optimization algorithms;

[0062] The fracture network morphology is expanded again based on the optimized fracturing design parameters to obtain the fracture point data of each fracturing horizontal well;

[0063] Determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0064] Calculate the current degree of healing of the fracture network of multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well;

[0065] If the current degree of bridging is greater than the initial degree of bridging, the preset parameter optimization algorithm is continued to be used to optimize the fracturing design parameters again based on the optimized fracturing design parameters;

[0066] The iteration is continued until the preset termination condition is reached, and the maximum bridging degree and the fracturing design parameters of multiple fractured horizontal wells under the maximum bridging degree are obtained.

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

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

[0069] The degree of bridging obtained after X consecutive iterations does not increase, where X is a positive integer; and / or

[0070] The iteration threshold is reached.

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

[0072] On the other hand, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the multi-horizontal well fracturing design parameter optimization method based on bridging fracture networks described in any of the above embodiments.

[0073] The embodiments of the present application provide a method and device for optimizing fracturing design parameters of multiple horizontal wells based on a bridged fracture network. The method and device expand the fracture network morphology of a target reservoir based on given fracturing design parameters to obtain fracture point data of multiple fractured horizontal wells in the target reservoir; determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well; calculate the initial bridge degree of the fracture network of multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well, wherein the bridge 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; based on the initial bridge degree, optimize the fracturing design parameters using a preset parameter optimization algorithm with the goal of maximizing the bridge degree of the fracture network of multiple fractured horizontal wells, and obtain the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridge degree. In this way, the fracturing design parameters are optimized with the goal of maximizing the degree of bridging the fracture network of multiple fractured horizontal wells, that is, the fracturing design parameters are optimized with the goal of minimizing the total unreformed volume and the total repeated reformed volume of multiple fractured horizontal wells, thereby achieving full utilization of resources with optimal benefits. BRIEF DESCRIPTION OF THE DRAWINGS

[0074] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following briefly introduces the drawings required for the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present application. For those skilled in the art, other drawings can be obtained based on these drawings without creative work. In the drawings:

[0075] FIG1 is a schematic diagram of optimized variables for bridging the fracture network of a fractured horizontal well obtained in one embodiment of the present application.

[0076] FIG2 is a flow chart of a method for optimizing design parameters of multi-horizontal well fracturing based on a bridged fracture network according to an embodiment of the present application.

[0077] FIG3 is a schematic diagram of a crack propagation mechanism provided in an embodiment of the present application.

[0078] FIG4 is a schematic diagram of determining crack distribution based on a grid background method according to an embodiment of the present application.

[0079] FIG5 is a partial flow diagram of a method for optimizing design parameters of multi-horizontal well fracturing based on a bridged fracture network provided in one embodiment of the present application.

[0080] FIG6 is a schematic diagram of the fracture morphology and ESRV of each section of a fractured horizontal well provided in one embodiment of the present application.

[0081] FIG7 is a schematic diagram of the effective pressure sweep range around a crack provided by an embodiment of the present application.

[0082] FIG8 is a schematic diagram of various volumes provided in an embodiment of the present application.

[0083] FIG9 is a flow chart of optimizing fracturing design parameters using a preset parameter optimization algorithm according to an embodiment of the present application.

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

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

[0086] FIG12 is a schematic diagram showing changes in the degree of bridging during optimization of fracturing design parameters according to an embodiment of the present application.

[0087] Figures 13, 14, and 15 show the change data of displacement, liquid volume, and horizontal position of a cluster in each section of the two horizontal wells during the iteration process.

[0088] FIG16 is a schematic structural diagram of a multi-horizontal well fracturing design parameter optimization device based on a bridged fracture network provided in one embodiment of the present application.

[0089] FIG17 is a schematic diagram of the physical structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0090] To make the purpose, technical solutions, and advantages of the embodiments of the present application more clearly understood, the embodiments of the present application are further described in detail below with reference to the accompanying drawings. The illustrative embodiments of the present application and their descriptions are used to explain the present application but are not intended to limit the present application. It should be noted that, unless there is a conflict, the embodiments and features in the embodiments of the present application may be combined with each other in any manner.

[0091] In order to better understand this application, the research background of this application is first introduced in detail below.

[0092] Conventional multi-well fracturing presents two major challenges. Increasing the volume of repeated stimulation between stages and between wells has little impact on increasing CBM production, while reducing the volume of unstimulated areas can significantly increase CBM production. Therefore, to achieve optimal resource utilization, minimizing the unstimulated and repeated stimulation volumes can ensure maximum bridging of the multi-well fracture network. Therefore, a comprehensive approach can be adopted in setting fracturing design parameters. Specifically, based on microseismic data, a fracture propagation simulator can be used to simulate different fracture network morphologies based on different operation parameters (such as operation flow rate and operation fluid volume). By varying the location of the stage clusters, a non-uniform fracture distribution effect can be achieved across multiple wells. Then, an intelligent optimization algorithm can be used to comprehensively optimize the fracturing design parameters to achieve optimal fracture bridging.

[0093] To achieve the goal of optimally bridging the fracture network, well parameters and fracture parameters can be considered separately, with displacement (dis), fluid volume (vol), stage length (spacing of stages), clusters (perforations), and location (loc) being considered independent optimization variables. This ensures flexible adjustment of these parameters to meet diverse geological conditions and engineering requirements, maximizing oil and gas production efficiency while avoiding unnecessary resource waste. Figure 1 shows a schematic diagram of the optimization variables for bridging the fracture network in fractured horizontal wells. Stage length refers to the length of each stage, cluster (perforation) location refers to the location of each fracture, and displacement and fluid volume are variables controlling the fracture network morphology.

[0094] In multi-horizontal well fracturing, in order to achieve the ultimate utilization of oil and gas resources while avoiding excessive inter-stage repeated stimulation, inter-well repeated stimulation areas, and unstimulated areas caused by uniform conceptual design, this application proposes a method with the goal of maximizing the degree of bridging of the fracture network of multiple fractured horizontal wells (or minimizing the unstimulated volume and the repeatedly stimulated volume), which is achieved by constructing an objective function for the degree of fracture network bridging. This objective function can reduce the cost and increase the efficiency of multi-horizontal well fracturing development, reduce the area of ​​unstimulated areas within the well control area, increase the oil (gas) leakage area of ​​the reservoir, and simultaneously reduce the inter-stage repeated stimulation and inter-well repeated stimulation areas to achieve the ultimate utilization of oil and gas resources. This method takes into account the changes in geological conditions, the irregularities of the fracture network, and the impact of the fracture network morphology of individual horizontal wells and between wells on the degree of bridging.

[0095] The basic idea of ​​the Monte Carlo gradient approximation (MCGA) algorithm is to first use the Monte Carlo method to generate several realizations of random variables around the current optimal variable, then calculate the performance indicator function value corresponding to each variable realization, and finally use each realization and its performance indicator function value to estimate the gradient of the performance indicator function (for example, the objective function of the degree of bridging in this application).

[0096] For oil and gas reservoir production optimization, let the optimal control variable obtained in the lth iteration be u l , and its corresponding performance index function is F(u l ,y). According to the principle of MCGA method, in u l Generate N around r The disturbance variable is expressed as: l,j =u l +γ×r j ,j=1,2,…,N r (1)

[0097] Where u l,j represents the jth control variable realization; γ is the disturbance constant; r j is the jth perturbation vector, where the component r i j (i=1,2,…,N u ) is generally a variable that conforms to the standard normal distribution, where N u Represents the total number of control variables, i.e. r i j ~N(0,1).

[0098] Each implementation u l,j Substituted into the oil and gas reservoir simulator and calculated, the corresponding performance index F(u l,j ,y), let bj F(u l,j ,y) and the current optimal performance index F(u l,j ,y), that is, b j =F(u l,j ,y)-F(u l ,y),j=1,2,…,N r (2)

[0099] Then the MCGA method is used to obtain the index function F(u l ,y) in u l Gradient estimate at The expression is:

[0100] Where, express The i-th component value. Analyze the nature of .

[0101] Consider F(u l,j ,y) in u l Perform a first-order Taylor expansion. Since the value of γ is usually small, we have F(u l,j ,y)=F(u l ,y)+γ(r j ) T g(u l ) (4)

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

[0103] Among them, g s (u l ) represents the sth component of the true gradient s=1,2…,N u Substitute the above formula into In the expression, there is

[0104] Taking the expected value on both sides of the above equation, we can get

[0105] Due to r j The components in are random variables that conform to the standard normal distribution, then they satisfy the following conditions:

[0106] According to the above formula, The expected value of becomes,

[0107] Obviously, The expected value of the true gradient is the i-th component, then the expected value of the gradient estimate based on the MCGA method is the true gradient Therefore, the gradient of the MCGA method can be used to iteratively optimize the performance index function, and its expression is:

[0108] Where u l+1 is the control variable after iterative optimization; express The infinite norm of ; α is the search step size. During the optimization process, the step size α is determined by a simple inexact line search method. If the performance index value of the current step does not strictly increase, α will be halved until the performance index value increases.

[0109] FIG2 is a flow chart of a method for optimizing design parameters of multi-horizontal well fracturing based on a bridged fracture network according to an embodiment of the present application. As shown in FIG2 , the method for optimizing design parameters of multi-horizontal well fracturing based on a bridged fracture network according to an embodiment of the present application includes:

[0110] S101, expanding the fracture network of a target reservoir based on given fracturing design parameters to obtain fracture point data of multiple fractured horizontal wells in the target reservoir;

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

[0112] Reservoir geological parameters, fracturing operation parameters, and other parameters are equivalently converted into electrical potential parameters. The mechanical problems in the crack expansion process are equivalently converted into electrical problems during lightning breakdown. This simulation approach is used to simulate crack propagation. A modified maximum circumferential tensile stress is used to determine crack initiation: cracks can initiate at locations where the circumferential tensile stress exceeds the critical stress for crack initiation (i.e., cracks can propagate in multiple directions). A fractal index is introduced to calculate the probability distribution of rupture at each node at the crack tip, and random numbers are used to determine the direction of crack propagation. The crack propagation mechanism is shown in Figure 3.

[0113] S102, determining the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0114] In step S102, a grid system can be used to characterize the fracture morphology. That is, the fracture point data in the gridless system can be characterized by a background grid. The grid background method determines the fracture distribution, as shown in Figure 4.

[0115] S103. Calculating the initial healing degree of the fracture network of multiple fractured horizontal wells in the target reservoir based on 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 un-stimulated volume and the total repeatedly stimulated volume of the multiple fractured horizontal wells;

[0116] In step S103, the degree of healing of a multi-fractured horizontal well is related to the fracture morphology of each fractured horizontal well. Specifically, the degree of healing is negatively correlated with the sum of the total unstimulated volume and the total repeatedly stimulated volume of the multi-fractured horizontal well, and positively correlated with the total effective stimulated volume. Therefore, the initial degree of healing of the fracture network of a multi-fractured horizontal well can be calculated based on the fracture morphology of each fractured horizontal well under given fracturing design parameters.

[0117] S104. Based on the initial bridging degree, a preset parameter optimization algorithm is used to optimize the fracturing design parameters with the goal of maximizing the bridging degree of the fracture network of the multiple fractured horizontal wells, thereby obtaining the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree.

[0118] 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, the degree of bridging of the multiple fractured horizontal well network under the fracturing design parameters is calculated, and the bridging degree is compared with the bridging degree of the multiple fractured horizontal well network under the fracturing design parameters of the previous iteration. If the bridging degree calculated this time is greater than the bridging degree calculated in the previous iteration, the preset parameter optimization algorithm is continued to be used to optimize the fracturing design parameters until the iteration termination condition is reached, and the fracturing design parameters of the multiple fractured horizontal wells under the maximum bridging degree and the maximum fitting degree are obtained.

[0119] The present application provides a method for optimizing fracturing design parameters of multiple horizontal wells based on bridged fracture networks, which obtains fracture point data of multiple fractured horizontal wells in the target reservoir by expanding the fracture network morphology of the target reservoir based on given fracturing design parameters; determines the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well; calculates the initial bridge degree of the fracture network of multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well, wherein the bridge 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; based on the initial bridge degree, optimizes the fracturing design parameters using a preset parameter optimization algorithm with the goal of maximizing the bridge degree of the fracture network of multiple fractured horizontal wells, and obtains the fracturing design parameters of multiple fractured horizontal wells at the maximum bridge degree. In this way, the fracturing design parameters are optimized with the goal of maximizing the degree of bridging the fracture network of multiple fractured horizontal wells, that is, the fracturing design parameters are optimized with the goal of minimizing the total unreformed volume and the total repeated reformed volume of multiple fractured horizontal wells, thereby achieving full utilization of resources with optimal benefits.

[0120] In some embodiments, fracturing design parameters include operation parameters and segment cluster location parameters; wherein the operation parameters include operation displacement and / or operation fluid volume. Based on the optimized segment cluster location, segment length, segment spacing, and cluster spacing can be calculated to guide on-site fracturing operation.

[0121] As shown in FIG5 , in some embodiments, the fracture morphology of each fractured horizontal well includes the total well-controlled volume, the total effective stimulated volume, the total repeatedly stimulated volume, and / or the total unstimulated volume. Determining the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well includes:

[0122] S1021. Determine the effective producing range of each fractured horizontal well based on actual production dynamic simulation;

[0123] In step S1021, the stimulated reservoir volume (SRV) is calculated, which greatly enhances the original formation permeability and has a decisive impact on the production of production wells. SRV can be calculated using microseismic imaging, secondary fracture monitoring technology, and mathematical models based on fracture propagation laws. However, the SRV calculated using the above methods differs significantly from the effective stimulated reservoir volume (ESRV) estimated through dynamic production curves or actual oil and gas field production. SRV is typically much larger than ESRV. Compared to SRV, ESRV is a more important determinant of the success of hydraulic fracturing and the prediction of fracturing-enhanced production effects. As shown in Figure 6, the effective stimulated reservoir volume (ESRV) is the volume of the effective range of reservoir fluid after pressure sweep during depressurization production, based on the morphology of each fracture segment. By establishing a matrix-fracture model, the effective pressure sweep range around the fracture during depressurization production is simulated, as shown in Figure 7. The effective stimulated reservoir volume (ESRV) can be calculated based on the fracture morphology. Specifically, the effective range of each fractured horizontal well can be obtained based on actual production dynamic simulation.

[0124] S1022. Determine the dynamic stimulation volume of a single fracture segment in each fractured horizontal well based on the effective producing range and fracture point data of each fractured horizontal well;

[0125] 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 fracture segment is calculated based on the effective producing range.

[0126] S1023. Calculate the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume and / or total unstimulated volume of each fractured horizontal well based on the background grid method according to the dynamic stimulation volume of a single fracture segment of each fractured horizontal well.

[0127] In step S1023, the index information of each fracture can be marked based on the background grid method. Specifically, in the background grid, 0 represents the unstimulated volume, 1 represents the volume controlled by natural fractures and artificial fractures, and 2 represents the volume repeatedly stimulated between sections and wells. The position information of the index marks is statistically incorporated into the calculation formulas for the total effective stimulated volume, total repeatedly stimulated volume, total well-controlled volume, and total unstimulated volume to calculate the aforementioned four volumes.

[0128] In some embodiments, according to the dynamic stimulation volume of a single fracture segment of each fractured horizontal well, calculating the total effective stimulation volume, the total repeated stimulation volume, the total well-controlled volume and / or the total unstimulated volume of each fractured horizontal well based on the background grid method includes:

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

[0130] According to the index information, the total effective stimulation volume of each fractured horizontal well is obtained by combining the dynamic stimulation volumes of each single fracture segment of all fractured horizontal wells; and / or

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

[0132] According to the index information, the total inter-well repeated stimulation volume of each fractured horizontal well is obtained by intersecting the fracture dynamic stimulation volumes of each fractured horizontal well.

[0133] Determine the total repetitive stimulation volume of each fractured horizontal well based on the total inter-stage repetitive stimulation volume of each fractured horizontal well and the total inter-well repetitive stimulation volume of each fractured horizontal well; and / or

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

[0135] According to the index information, the total unstimulated volume of each fractured horizontal well is obtained by combining the unstimulated grid cell volumes of each fractured horizontal well.

[0136] Specifically, based on the results of the fracture propagation simulation, the fracture point data of each fractured horizontal well is obtained, and then the effective utilization range of each fractured horizontal well is obtained based on the actual production dynamic simulation. Finally, the fracture distribution under the grid is determined by the grid background method to calculate the effective stimulation volume, repeated stimulation volume (the sum of the repeated stimulation volume between sections and the repeated stimulation volume between wells), well-controlled volume and unstimulated volume. The schematic diagram of each type of volume is shown in Figure 8. c It can be obtained by summing the well-controlled volumes of each fractured horizontal well:

[0137] Where 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 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.

[0138] The total unmodified volume V can be obtained by summing the control volumes with index information 0 in each unit grid. 0 :

[0139] Where V 0 is the total unmodified volume; N w is the number of fractured horizontal wells; N u | λ＝0 (w i ) is w i The number of unstimulated grid cells in the fractured horizontal well; is the volume of the i-th unit grid; 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.

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

[0141] Where, w i The control volume of the fractured horizontal well with index information of 1, that is, the effective control volume of a single well; V 1 is the control volume of all fractured horizontal wells with index information of 1, that is, the total effective control volume; uni(·) is a unique value function; N ut | λ＝1 is the total number of (effective) transformed grid cells; N u | λ＝1 (w i ) is w i The number of (effective) stimulation grid cells for fractured horizontal wells. 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 iis the grid step size in the z direction of the i-th unit grid.

[0142] The total repeated stimulation volume can then be calculated based on the effective control volume of each fractured horizontal well:

[0143] Where V 2 is the total repeated transformation volume; N r | λ＝2 The number of grid cells to be transformed repeatedly is the number of grid cells to be transformed repeatedly. However, when calculating the actual degree of bridging, the effective transformation volume of each fractured horizontal well can be obtained by intersection. Calculate the total repeated transformation volume V 2 .

[0144] In some embodiments, calculating the initial healing degree of a fracture network of multiple fractured horizontal wells in a target reservoir based on the fracture morphology of each fractured horizontal well includes:

[0145] Calculate the initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir based on the ratio of the sum of the total re-stimulated volume and the total unstimulated volume of each fractured horizontal well to the total well-controlled volume; or

[0146] The initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir is calculated based on the ratio of the total effective stimulation volume of each fractured horizontal well to the total well-controlled volume.

[0147] Specifically, the degree of fracture network bridging in multi-fractured horizontal wells can be defined as follows:

[0148] Where, M is the degree of seam network bridging, dimensionless; N u | λ＝0 (w i ) is w i Total number of repeated stimulations between sections and wells in fractured horizontal wells; N ut | λ＝1 is the total number of (effective) transformed grid cells; uni(·) is a unique value function; N rt | λ＝2 is the total number of grid cells that are repeatedly transformed; N is the number of well-controlled grids in the block; the definitions of other parameters refer to the above.

[0149] It can be seen from this that the total well-controlled volume calculated by equation (11), the total unreformed volume calculated by equation (12), and the total re-reformed volume calculated by equation (15) can be substituted into equation (16), or the total well-controlled volume calculated by equation (11) and the total effective control volume calculated by equation (14) can be substituted into equation (16) to calculate the degree of fracture network bridging. Both methods can calculate the degree of bridging M under the current number of wells and fracture morphology.

[0150] In some embodiments, based on the initial bridging degree, a preset parameter optimization algorithm is used to optimize the fracturing design parameters with the goal of maximizing the bridging degree of the fracture network of multiple fractured horizontal wells. The fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree include:

[0151] Optimize fracturing design parameters using preset parameter optimization algorithms;

[0152] The fracture network morphology is expanded again based on the optimized fracturing design parameters to obtain the fracture point data of each fracturing horizontal well;

[0153] Determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0154] Calculate the current degree of healing of the fracture network of multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well;

[0155] If the current degree of bridging is greater than the initial degree of bridging, the preset parameter optimization algorithm is continued to be used to optimize the fracturing design parameters again based on the optimized fracturing design parameters;

[0156] The iteration is continued until the preset termination condition is reached, and the maximum bridging degree and the fracturing design parameters of multiple fractured horizontal wells under the maximum bridging degree are obtained.

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

[0158] Step 1: Input the initial displacement, fluid volume, and perforation positions of each fractured horizontal well.

[0159] Step 2: Based on the fracture propagation inversion simulator, the fracture morphology of each fractured horizontal well was obtained, and the effective production range was determined through actual production simulation. Next, the effective stimulation volume, repeated stimulation volume, and well control volume were calculated, and the initial bridging degree M0 was calculated using Equation (16).

[0160] Step 3: Use a preset parameter optimization algorithm (such as MCGA algorithm - mapping crossover genetic algorithm) to iteratively update the optimization variables.

[0161] Step 4: Call the crack extension simulator again to obtain the updated crack morphology, and calculate the new bridging degree M1 based on the effective mobilization range.

[0162] Step 5: If M1>M0, proceed to the next step; otherwise, proceed to step 3.

[0163] Step 6: When the termination condition is met, the optimal displacement, liquid volume, perforation positions, and optimal healing degree M* are output. Otherwise, the optimized fracture morphology is used as the initial fracture morphology, and the healing degree M1 is assigned to M0, and step 3 is executed.

[0164] In some embodiments, the preset termination conditions include: the degree of bridging obtained after X consecutive iterations no longer increases, where X is a positive integer; and / or a threshold number of iterations is reached.

[0165] For a better understanding of the present application, the multi-horizontal well fracturing design parameter optimization method based on bridging fracture networks provided by the present application is described below through a specific embodiment.

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

[0167] The permeability and porosity distributions of this coalbed methane reservoir are shown in Figures 10a and 10b. The maximum permeability is 0.1 mD, the minimum is 0.01 mD, and the average is 0.05 mD. The maximum porosity is 0.05, the minimum is 0.01, and the average is 0.036. The reservoir parameters and simulation information for this example are shown in Table 1.

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

[0169] With maximizing the degree of bridging in a multi-well fractured horizontal well as the optimization objective, the MCGA algorithm was used to perturb and iteratively update the flow rates, fluid volumes, and cluster locations, achieving integrated optimization of the optimization parameters and obtaining the fracture deployment plan with the optimal degree of bridging. Figures 11a and 11b compare the fracture network morphology before and after optimization of the flow rates, fluid volumes, and cluster locations. Figures 11a and 11b clearly show that the horizontal positions of the fractures change after optimization, and the fracture lengths generally increase. Figure 12 shows that the MCGA-based multi-well fracture network optimization converges rapidly, starting around step 15. The final degree of bridging is 0.5638, an increase of 40.95%. Figures 13, 14, and 15 respectively illustrate the changes in flow rates, fluid volumes, and cluster locations for each segment of two horizontal wells during the iterative process. Taking the optimization results for the first cluster of each segment as an example, once the fracture locations are determined, the segment lengths, segment spacing, and cluster spacing between clusters can be determined.

[0170] FIG16 is a schematic diagram of the structure of a multi-horizontal well fracturing design parameter optimization device based on a bridged fracture network according to an embodiment of the present application. As shown in FIG16 , the multi-horizontal well fracturing design parameter optimization device based on a bridged fracture network according to an embodiment of the present application includes:

[0171] The morphology expansion module 21 is used to expand the fracture network morphology of the target reservoir based on the given fracturing design parameters to obtain the fracture point data of multiple fractured horizontal wells in the target reservoir;

[0172] The fracture morphology determination module 22 is used to determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0173] The healing degree calculation module 23 is used to calculate the initial healing degree of the fracture network of multiple fractured horizontal wells in the target reservoir based on 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 un-reformed volume and the total repeatedly re-reformed volume of the multiple fractured horizontal wells;

[0174] The parameter optimization module 24 is used to optimize the fracturing design parameters based on the initial bridging degree using a preset parameter optimization algorithm with the goal of maximizing the bridging degree of the fracture network of multiple fractured horizontal wells, and obtain the fracturing design parameters of the multiple fractured horizontal wells under the maximum bridging degree.

[0175] The embodiment of the present application provides a multi-horizontal well fracturing design parameter optimization device based on bridged fracture networks, which obtains fracture point data of multiple fractured horizontal wells in the target reservoir by expanding the fracture network morphology of the target reservoir based on given fracturing design parameters; determines the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well; calculates the initial bridge degree of the fracture network of the multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well, wherein the bridge degree of the fracture network of the 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; based on the initial bridge degree, optimizes the fracturing design parameters using a preset parameter optimization algorithm with the goal of maximizing the bridge degree of the fracture network of the multiple fractured horizontal wells, and obtains the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridge degree. In this way, the fracturing design parameters are optimized with the goal of maximizing the degree of bridging the fracture network of multiple fractured horizontal wells, that is, the fracturing design parameters are optimized with the goal of minimizing the total unreformed volume and the total repeated reformed volume of multiple fractured horizontal wells, thereby achieving full utilization of resources with optimal benefits.

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

[0177] In some embodiments, the fracture morphology of each fractured horizontal well includes the total well-controlled volume, the total effective stimulation volume, the total repeated stimulation volume and / or the total unstimulated volume; the fracture morphology determination module is specifically configured to:

[0178] Based on actual production dynamic simulation, the effective producing range of each fractured horizontal well is obtained;

[0179] Determine the dynamic stimulation volume of a single fracture segment in each fractured horizontal well based on the effective producing range and fracture point data of each fractured horizontal well;

[0180] According to the dynamic stimulation volume of a single fracture segment of each fractured horizontal well, the total effective stimulation volume, total repeated stimulation volume, total well-controlled volume and / or total unstimulated volume of each fractured horizontal well are calculated based on the background grid method.

[0181] In some embodiments, the fracture morphology determination module calculates the total effective stimulated volume, the total repeated stimulated volume, the total well-controlled volume, and / or the total unstimulated volume of each fractured horizontal well based on the background grid method according to the dynamic stimulated volume of a single fracture segment of each fractured horizontal well, including:

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

[0183] According to the index information, the total effective stimulation volume of each fractured horizontal well is obtained by combining the dynamic stimulation volumes of each single fracture segment of all fractured horizontal wells; and / or

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

[0185] According to the index information, the total inter-well repeated stimulation volume of each fractured horizontal well is obtained by intersecting the fracture dynamic stimulation volumes of each fractured horizontal well.

[0186] Determine the total repetitive stimulation volume of each fractured horizontal well based on the total inter-stage repetitive stimulation volume of each fractured horizontal well and the total inter-well repetitive stimulation volume of each fractured horizontal well; and / or

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

[0188] According to the index information, the total unstimulated volume of each fractured horizontal well is obtained by combining the unstimulated grid cell volumes of each fractured horizontal well.

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

[0190] Calculate the initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir based on the ratio of the sum of the total re-stimulated volume and the total unstimulated volume of each fractured horizontal well to the total well-controlled volume; or

[0191] The initial degree of bridging of the fracture network of multiple fractured horizontal wells in the target reservoir is calculated based on the ratio of the total effective stimulation volume of each fractured horizontal well to the total well-controlled volume.

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

[0193] Optimize fracturing design parameters using preset parameter optimization algorithms;

[0194] The fracture network morphology is expanded again based on the optimized fracturing design parameters to obtain the fracture point data of each fracturing horizontal well;

[0195] Determine the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well;

[0196] Calculate the current degree of healing of the fracture network of multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well;

[0197] If the current degree of bridging is greater than the initial degree of bridging, the preset parameter optimization algorithm is continued to be used to optimize the fracturing design parameters again based on the optimized fracturing design parameters;

[0198] The iteration is continued until the preset termination condition is reached, and the maximum bridging degree and the fracturing design parameters of multiple fractured horizontal wells under the maximum bridging degree are obtained.

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

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

[0201] The degree of bridging obtained after X consecutive iterations does not increase, where X is a positive integer; and / or

[0202] The iteration threshold is reached.

[0203] The embodiments of the apparatus provided in the embodiments of the present application can be specifically used to execute the processing flow of the above-mentioned method embodiments. Its functions will not be described in detail here, and reference can be made to the detailed description of the above-mentioned method embodiments.

[0204] FIG17 is a schematic diagram of the physical structure of an electronic device provided in an embodiment of the present application. As shown in FIG17 , 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 invoke logic instructions in the memory 140 to execute the following method: obtaining fracture point data for multiple fractured horizontal wells in the target reservoir by expanding the fracture network morphology of the target reservoir based on given fracturing design parameters; determining the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well; calculating the initial degree of bridging of the fracture network of the multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well, wherein the degree of bridging of the fracture network of the multiple fractured horizontal wells is negatively correlated with the sum of the total unstimulated volume and the total repeatedly stimulated volume of the multiple fractured horizontal wells; and optimizing the fracturing design parameters based on the initial degree of bridging, using a preset parameter optimization algorithm, with the goal of maximizing the degree of bridging of the fracture network of the multiple fractured horizontal wells, to obtain the fracturing design parameters for the multiple fractured horizontal wells at the maximum degree of bridging.

[0205] The present embodiment discloses a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the methods provided by the above-mentioned method embodiments, for example, including: obtaining fracture point data of multiple fractured horizontal wells in the target reservoir by expanding the fracture network morphology of the target reservoir based on given fracturing design parameters; determining the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well; calculating the initial bridging degree of the fracture network of the multiple fractured horizontal wells in the target reservoir based on 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 unreformed volume and the total repeatedly reformed volume of the multiple fractured horizontal wells; based on the initial bridging degree, optimizing the fracturing design parameters using a preset parameter optimization algorithm with the goal of maximizing the bridging degree of the fracture network of the multiple fractured horizontal wells, and obtaining the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree.

[0206] This embodiment provides a computer-readable storage medium, which stores a computer program. The computer program enables a computer to execute the methods provided by the above-mentioned method embodiments, for example, including: obtaining fracture point data of multiple fractured horizontal wells in the target reservoir by expanding the fracture network morphology of the target reservoir based on given fracturing design parameters; determining the fracture morphology of each fractured horizontal well based on the fracture point data of each fractured horizontal well; calculating the initial healing degree of the fracture network of the multiple fractured horizontal wells in the target reservoir based on the fracture morphology of each fractured horizontal well, wherein the healing degree of the fracture network of the 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; based on the initial healing degree, optimizing the fracturing design parameters using a preset parameter optimization algorithm with the goal of maximizing the healing degree of the fracture network of the multiple fractured horizontal wells, and obtaining the fracturing design parameters of the multiple fractured horizontal wells at the maximum healing degree.

[0207] As shown in Figure 17 , 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 is worth noting that electronic device 600 does not necessarily include all of the components shown in Figure 17 ; furthermore, electronic device 600 may also include components not shown in Figure 17 , for which reference may be made to existing technologies. It is worth noting that this figure is illustrative only; other types of structures may be used to supplement or replace this structure to implement telecommunication or other functions.

[0208] As shown in FIG. 17 , the processor 100 is sometimes also referred to as a controller or an operation control, and may include a microprocessor or other processor device and / or logic device. The processor 100 receives input and controls the operation of various components of the electronic device 600 .

[0209] Memory 140 may be, for example, one or more of a cache, flash memory, a hard drive, removable media, volatile memory, non-volatile memory, or other suitable devices. It may store the aforementioned failure-related information and may also store a program for executing the relevant information. Processor 100 may execute the program stored in memory 140 to implement information storage or processing.

[0210] Input unit 120 provides input to processor 100. Input unit 120 is, for example, a keypad or touch input device. Power supply 170 is used to provide power to electronic device 600. Display 160 is used to display objects such as images and text. Display 160 can be, for example, an LCD display, but is not limited thereto.

[0211] The memory 140 may be a solid-state memory, such as a read-only memory (ROM), a random access memory (RAM), or a SIM card. Alternatively, it may be a memory that retains information even when power is off, can be selectively erased, and is provided with more data. Examples of memory 140 are sometimes referred to as EPROMs. The memory 140 may also be some other type of device. The memory 140 includes a buffer 141 (sometimes referred to as a buffer memory). The memory 140 may include an application / function storage unit 142 for storing application programs and function programs or processes for executing the operation of the electronic device 600 via the processor 100.

[0212] The memory 140 may also include a data storage unit 143 for storing 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 may include various driver programs for communication functions of the electronic device and / or for executing other functions of the electronic device (such as messaging applications, address book applications, etc.).

[0213] The communication module 110 includes a transmitter / receiver for transmitting and receiving signals via an antenna 111. The communication module 110 is coupled to the processor 100 to provide input signals and receive output signals, which may be the same as in the case of a conventional mobile communication terminal.

[0214] Based on different communication technologies, multiple communication modules 110 may be provided in the same electronic device, such as a cellular network module, a Bluetooth module, and / or a wireless local area network module. The communication module 110 is also coupled to a speaker 131 and a microphone 132 via an audio processor 130 to provide audio output via the speaker 131 and receive audio input from the microphone 132, thereby implementing common telecommunication functions. The audio processor 130 may include any suitable buffer, decoder, amplifier, etc. Furthermore, the audio processor 130 is coupled to the processor 100, enabling local recording via the microphone 132 and playback of stored audio via the speaker 131.

[0215] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0216] The present application is described with reference to the flow chart and / or block diagram of the method, device (system), and computer program product according to the embodiment of the present application. It should be understood that each flow process and / or box in the flow chart and / or block diagram and the combination of the flow process and / or box in the flow chart and / or block diagram can be realized by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processing machine or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for realizing the function specified in one flow chart flow or multiple flows and / or one box or multiple boxes of the block diagram.

[0217] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a product including an instruction device that implements the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

[0218] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce 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 or more processes in the flowchart and / or one or more boxes in the block diagram.

[0219] In the description of this specification, reference to the terms "one embodiment," "a specific embodiment," "some embodiments," "for example," "example," "specific example," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. 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 may be combined in any appropriate manner in any one or more embodiments or examples.

[0220] The specific embodiments described above further illustrate the purpose, technical solutions and beneficial effects of the present application. It should be understood that the above description is only a specific embodiment of the present application and is not intended to limit the scope of protection of the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application should be included in the scope of protection of the present application.

Claims

1. A method for optimizing fracturing design parameters based on bridging fracture networks for multi - horizontal wells, characterized in that, it includes: Expanding the fracture network morphology of the target reservoir based on the 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 networks 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 networks of the multiple fracturing horizontal wells is negatively correlated with the sum of the total un - transformed volume and the total repeated transformation 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 networks of the multiple fracturing horizontal wells, and obtaining the fracturing design parameters of the multiple fracturing horizontal wells under the maximum bridging 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 un - transformed 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 - section fracture of each fracturing horizontal well; Based on the dynamic transformation volume of a single - section fracture of each fracturing horizontal well, calculating the total effective transformation volume, total repeated transformation volume, total well - controlled volume and / or total un - transformed 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 un - transformed 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 includes: Marking the index information of each fracture of each fracturing horizontal well based on the background grid method; According to the index information, by taking the union of the dynamic transformation volumes of single - section fractures of all fracturing horizontal wells, obtaining the total effective transformation volume of each fracturing horizontal well; and / or According to the index information, by taking the intersection of the dynamic transformation volumes of single - section fractures in each fracturing horizontal well, obtaining the total inter - section repeated transformation volume of each fracturing horizontal well; According to the index information, by taking the intersection of the fracture dynamic transformation volumes of two fracturing horizontal wells pairwise, obtaining the total inter - well repeated transformation volume of each fracturing horizontal well; Determining the total repeated transformation volume of each fracturing horizontal well according to the total inter - section 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, by taking the union of the well - controlled volumes of each fracturing horizontal well, obtaining the total well - controlled volume 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 volumes of the unmodified grid cells of each fractured horizontal well.

5. The method according to claim 1, wherein, the calculating the initial bridging degree of the fracture network of multiple fractured horizontal wells in the target reservoir according to the fracture patterns 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 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 multiple fractured horizontal wells to obtain the fracturing design parameters of the multiple fractured horizontal wells at the maximum bridging degree 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 pattern 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 patterns 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 on the basis of the optimized fracturing design parameters by using a preset parameter optimization algorithm; continue to iterate 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.

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 from continuous X iterations no longer increases, where X is a positive integer; and / or reaching the iteration number threshold.

9. A device for optimizing fracturing design parameters of multiple horizontal wells based on a bridged fracture network, wherein, it includes: a morphology expansion module, configured to perform fracture network morphology expansion on a target reservoir based on given fracturing design parameters to obtain fracture point data of multiple fractured horizontal wells in the target reservoir; a fracture pattern determination module, configured to determine the fracture pattern 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 patterns 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 device 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 device according to claim 10, wherein, the fracture patterns of the fractured horizontal wells include total well-controlled volume, total effective stimulated volume, total repeated stimulated 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 stimulated 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 stimulated volume, total repeated stimulated 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 stimulated volume of the single-section fracture of each fractured horizontal well.

12. The device according to claim 11, wherein, the fracture pattern determination module calculates the total effective stimulated volume, total repeated stimulated 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 stimulated 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 stimulated volume of each fractured horizontal well by taking the union of the dynamic stimulated volumes of the single-section fractures of all fractured horizontal wells according to the index information; and / or obtaining the total inter-section repeated stimulated volume of each fractured horizontal well by taking the intersection of the dynamic stimulated volumes of the single-section fractures in each fractured horizontal well according to the index information; obtaining the total inter-well repeated stimulated volume of each fractured horizontal well by taking the intersection of the fracture dynamic stimulated volumes of two fractured horizontal wells pairwise according to the index information; determining the total repeated stimulated volume of each fractured horizontal well according to the total inter-section repeated stimulated volume of each fractured horizontal well and the total inter-well repeated stimulated 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 cell volumes of each fractured horizontal well according to the index information.

13. The device 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 stimulated 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 stimulated volume of each fractured horizontal well to the total well-controlled volume.

14. The device 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 on the basis of 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 the iteration number threshold is reached.

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