Method and device for evaluating high crossing layer of lamination crack of stratified reservoir

Abstract

The invention provides a stratified reservoir lamination fracture high-crossing evaluation method and device, and relates to the technical field of oil-gas exploration, and the method comprises the steps: carrying out the combined analysis of a geomechanical model based on a target position of a research area, a first fracture propagation analysis result and a second fracture propagation analysis result, and obtaining a first fracture propagation analysis result and a second fracture propagation analysis result; determining conditions of crack induced bedding opening and crossing expansion under different reservoir and interlayer stress differences and interlayer thicknesses, and drawing a crack height crossing evaluation chart corresponding to the target position by taking the reservoir and interlayer stress differences and the interlayer thicknesses as two-dimensional coordinate axes; and performing layer crossing evaluation before fracturing construction of the to-be-built production block by using the seam height layer crossing evaluation chart. According to the method, the fracture height crossing evaluation chart is determined by combining the formation mechanical characterization, the fracture propagation physical simulation experiment and the multi-cluster fracturing fracture propagation numerical simulation experiment in the horizontal well, then the crossing evaluation of the to-be-built production block is carried out, and the post-fracturing production capacity, the overall production effect of the reservoir and the overall development benefit of the stratified reservoir are improved.

Description

Technical Field

[0001] This invention belongs to the field of oil and gas exploration technology, and specifically relates to a method for evaluating high-penetration fractures in laminar reservoirs, a device for evaluating high-penetration fractures in laminar reservoirs, a machine-readable storage medium, and a computer device. Background Technology

[0002] Laminated shale oil reservoirs suffer from complex hydraulic fracture propagation patterns, unclear controlling factors, and a lack of quantitative evaluation methods, resulting in low overall development efficiency for laminated shale oil.

[0003] Research has revealed that many scholars have conducted studies on the propagation law of hydraulic fractures in layered shale oil. For example, Chinese patent application No. 202210627720.5 proposes a simulation test method for directional cross-layer fracturing in horizontal wells with soft coal seam roof strata. This method, targeting coalbed methane reservoirs, uses a layered loading true triaxial hydraulic fracturing simulation device to conduct horizontal well cross-layer fracturing simulation experiments. After the experiment, the experimental data and fracture morphology data are integrated to analyze the law of directional cross-layer fracturing. This method solves the problems of not being able to obtain the rock mechanical parameters of the soft coal seam under experimental conditions and the difficulty in processing and manufacturing the "roof-coal seam" cementation physical model of the soft coal seam. The manufactured "roof-coal seam" cementation physical model based on rock mechanical strength can more realistically reflect the actual formation and reservoir characteristics. It has the characteristics of high success rate and high efficiency in directional cross-layer fracturing physical simulation, making the experimental parameters used more consistent with the actual formation parameters, accurately reflecting the actual fracturing situation of the formation, and providing stronger guidance for the experimental results. However, this method does not quantify the specific implementation process of geological and engineering parameters. Chinese patent application No. 202010973323.4 proposes a method for cross-layer fracturing in thin interbedded reservoirs and its application. This method addresses the problems of large fluid loss along the bedding planes, ineffective fracturing of multi-lithological reservoirs, inadequate proppant placement, and weak propping effect during fracturing of thin interbedded oil and gas reservoirs in coal-bearing strata. However, this method does not provide a quantitative evaluation or implementation method for cross-layer fracturing through experiments, numerical simulations, or field applications. Chinese patent application No. 202010494993.8 proposes a feasibility evaluation method for horizontal well cross-layer fracturing. This method targets tight, low-permeability oil and gas reservoirs and can evaluate whether fracturing the target oil and gas layer can simultaneously consider the permeability of adjacent oil and gas layers and the flow conductivity of proppant-carrying fluids in artificial fractures within mudstone interlayers. However, this method cannot account for the propagation patterns of artificial fractures in shale laminae.

[0004] Among the aforementioned existing technologies, some provide evaluation or implementation methods for cross-layer fracturing through experiments, numerical simulations, or field applications. However, none of them provide specific implementation steps that quantify geological and engineering parameters, nor can they obtain artificial fracture propagation patterns that take into account shale laminae. Therefore, they cannot guide the hydraulic fracturing construction process, making it difficult to improve the overall fracturing and stimulation effect of the reservoir.

[0005] In summary, there is an urgent need to propose a feasible scheme to quantitatively characterize the propagation law of artificial fractures in laminar reservoirs, so as to carry out evaluation of fracturing through laminar flow. Summary of the Invention

[0006] In view of this, the purpose of this invention is to provide a method for evaluating high-penetration fractures in laminar reservoirs, a device for evaluating high-penetration fractures in laminar reservoirs, a machine-readable storage medium, and a computer device, in order to overcome the technical problem in the prior art that it is impossible to conduct fracturing penetration evaluation on laminar reservoirs such as laminar shale oil.

[0007] To achieve the above objectives, a first aspect of the present invention provides a method for evaluating high-penetration fractures in laminar reservoirs, the method comprising:

[0008] By combining the geomechanical model based on the target location of the study area, the results of the first fracture propagation analysis, and the results of the second fracture propagation analysis, the conditions for fracture-induced bedding opening and cross-layer propagation under different reservoir stress differences and interlayer thicknesses were determined. The reservoir stress difference and interlayer thickness were used as two-dimensional coordinate axes to draw the fracture height cross-layer evaluation chart corresponding to the target location.

[0009] The aforementioned fracture height and layer-penetration evaluation chart is used to evaluate the layer-penetration before fracturing construction in the production block to be built.

[0010] Specifically, the first fracture propagation analysis result is determined based on the monitoring data from the physical simulation experiment of fracture propagation under different lamination characteristics, and the second fracture propagation analysis result is determined based on the output data from the numerical simulation experiment of multi-cluster fracturing fracture propagation in horizontal wells under different lamination characteristics.

[0011] Further improved, the process of constructing the geomechanical model includes:

[0012] Quantitative characterization of the laminar mechanical features at the target location in the study area;

[0013] A geomechanical model was constructed based on the quantitatively characterized laminar mechanical features.

[0014] Further improvements include the mechanical characteristics of the laminations, including the mechanical characteristics of the interlayer interfaces, the mechanical characteristics of the rock matrix, and the characteristics of in-situ stress.

[0015] Further improved, the process for determining the results of the first crack propagation analysis includes:

[0016] Based on the rock samples obtained from the target location in the study area, fracturing specimens with different lamination characteristics were constructed. Physical simulation experiments on crack propagation were carried out on the fracturing specimens to obtain monitoring data of the experimental process.

[0017] The first crack propagation analysis result is generated based on the monitoring data;

[0018] Among them, the results of the first crack propagation analysis reflect the laws and main controlling factors of artificial crack competition initiation, penetration, and induced bedding opening under different lamination characteristics.

[0019] Further improved, the process for determining the results of the second crack propagation analysis includes:

[0020] After characterizing the bedding features of the target location in the study area, a three-dimensional model of crack propagation was constructed.

[0021] Numerical simulation experiments on the propagation of multiple clusters of hydraulic fracturing fractures in horizontal wells under different lamination characteristics were carried out using the aforementioned three-dimensional model of fracture propagation, and the output data of the experimental process were obtained.

[0022] A second crack propagation analysis result is generated based on the output data;

[0023] The second crack propagation analysis results reflect the crack morphology and the law affecting crack height propagation under different laminar flow characteristics.

[0024] Further improvements include constructing a three-dimensional model of crack propagation after characterizing the bedding feature parameters of the target location in the study area, which includes:

[0025] Based on discrete element mechanics theory, a discrete element mechanics model considering the bedding of the target location in the study area is established;

[0026] After the injected fracturing fluid passes through the wellbore and perforation, it enters each cluster of fractures. The fluid inflow into each cluster of fractures is controlled by the parameters of the wellbore, perforation, and fractures to implement multi-cluster fracturing.

[0027] A fluid-structure interaction model for crack propagation is established by coupling rock deformation and fluid flow within cracks.

[0028] Characterization of bedding feature parameters for strata with bedding development at different scales.

[0029] Further improvements include bedding characteristic parameters such as bedding density, bedding width, bedding spacing, bedding permeability, and bedding mechanical properties.

[0030] Further improvements include characterizing the bedding characteristic parameters of strata with bedding development at different scales, including:

[0031] The total permeability of multiple bedding planes in the target location of the study area is equivalently represented as the permeability of a single bedding plane, thus obtaining the equivalent bedding plane permeability.

[0032] The bedding density of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding permeability, and the equivalent bedding density is obtained.

[0033] The equivalent bedding spacing of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding density.

[0034] The equivalent bedding width is obtained by equivalent characterizing the bedding width of the strata at the target location in the study area based on the equivalent bedding permeability.

[0035] Further improvements include the minimum displacement during fracturing operations as a condition for cross-layer propagation.

[0036] In a further improvement, the fracture height penetration evaluation chart is divided into a direct penetration zone, a T-shaped penetration zone, and an impenetrable suppression zone based on the stress difference and thickness of the reservoir layer.

[0037] A second aspect of the present invention provides a device for evaluating high-penetration fractures in lamellar reservoirs, the device comprising:

[0038] The fracture height and cross-layer evaluation chart construction module is used to combine and analyze the geomechanical model based on the target location of the study area, the results of the first fracture propagation analysis and the results of the second fracture propagation analysis, to determine the conditions for fracture-induced bedding and cross-layer propagation under different reservoir stress differences and interlayer thicknesses, and to draw the fracture height and cross-layer evaluation chart corresponding to the target location using the reservoir stress difference and interlayer thickness as two-dimensional coordinate axes.

[0039] The layer-penetration evaluation module is used to evaluate the layer-penetration of the production block before fracturing construction using the fracture height layer-penetration evaluation chart.

[0040] Specifically, the first fracture propagation analysis result is determined based on the monitoring data from the physical simulation experiment of fracture propagation under different lamination characteristics, and the second fracture propagation analysis result is determined based on the output data from the numerical simulation experiment of multi-cluster fracturing fracture propagation in horizontal wells under different lamination characteristics.

[0041] A third aspect of the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the high-penetration evaluation method for lamellar reservoir pressure fractures described in the first aspect of the present invention.

[0042] A fourth aspect of the present invention provides a machine-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method for evaluating high-penetration fractures in laminar reservoirs as described in the first aspect of the present invention.

[0043] In the aforementioned technical solution, by comprehensively applying methods for characterizing the mechanical features of laminae, physical simulation experiments of fracture propagation under different laminae characteristics, and numerical simulation experiments of multi-cluster hydraulic fracturing fracture propagation in horizontal wells, a high-penetration evaluation chart for hydraulically fractured laminae reservoirs is constructed. This chart uses the reservoir stress difference and interlayer thickness as two-dimensional coordinate axes and the conditions for fracture-induced bedding opening and cross-layer propagation as graphical information. This enables a quantitative evaluation of the propagation law and main controlling factors of artificial fractures under different lithological sequences, mechanical parameters, and interlayer interface characteristics. Before fracturing operations are carried out in the production area to be developed, the cross-layer propagation conditions can be obtained based on the high-penetration evaluation chart and the determined reservoir stress difference and interlayer thickness. These cross-layer propagation conditions can be used to establish fracturing operation parameters or optimize existing fracturing operation parameters, thereby improving the post-fracturing production capacity of horizontal wells and the overall utilization effect of the reservoir. The differentiated fracturing scheme based on the high-penetration evaluation chart can controllably reduce single-well investment and improve the overall development benefits of laminae reservoirs.

[0044] Meanwhile, the above-mentioned technical solution makes cross-layer evaluation convenient and quick. For example, in the early stages of development, an initial fracture height cross-layer evaluation chart applicable to each sweet spot location in the target production area can be quickly established. At the same time, during field practice, by continuously improving the understanding of the lithological sequence, mechanical parameters and interlayer interface characteristics of the target production area, the fracture height cross-layer evaluation chart can be updated and iterated, thereby guiding the rapid establishment of differentiated fracturing schemes and targeted optimization of fracturing processes and key parameters in the production blocks to be developed.

[0045] Other features and advantages of the embodiments of the present invention will be described in detail in the following detailed description section. Attached Figure Description

[0046] The accompanying drawings are provided to further illustrate embodiments of the present invention and form part of the specification. They are used together with the following detailed description to explain the embodiments of the present invention, but do not constitute a limitation thereof. In the drawings:

[0047] Figure 1 A flowchart illustrating a method for evaluating high-penetration fractures in lamellar reservoirs according to an embodiment of the present invention is shown.

[0048] Figure 2 This schematically illustrates another flowchart of a method for evaluating high-penetration fractures in lamellar reservoirs according to an embodiment of the present invention;

[0049] Figure 3A schematic diagram of the geostress profile is shown.

[0050] Figure 4 A schematic diagram of the geomechanical model is shown.

[0051] Figure 5 The diagram illustrates a comparison of crack propagation morphology under different stress differences.

[0052] Figure 6 The diagram illustrates the construction process of the three-dimensional model of crack propagation.

[0053] Figure 7 A schematic diagram of the evaluation chart for seam height penetration is shown;

[0054] Figure 8 A schematic diagram illustrating the components of a high-penetration evaluation device for pressure fractures in lamellar reservoirs according to an embodiment of the present invention is shown.

[0055] Figure 9 A schematic block diagram of a computer device according to an embodiment of the present invention is shown. Detailed Implementation

[0056] The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are for illustration and explanation only and are not intended to limit the scope of the present invention.

[0057] Example 1

[0058] See Figure 1 This invention provides a method for evaluating hydraulic fracturing penetration in laminar reservoirs, used to evaluate hydraulic fracturing penetration in, for example, laminar shale oil reservoirs. The method includes the following steps:

[0059] Step S100 involves combining the geomechanical model of the target location in the study area, the first fracture propagation analysis results determined by monitoring data from physical simulation experiments of fracture propagation under different laminar flow characteristics, and the second fracture propagation analysis results determined by output data from numerical simulation experiments of multi-cluster fracturing fracture propagation in horizontal wells under different laminar flow characteristics. This analysis determines the conditions for fracture-induced bedding opening and cross-layer propagation under different reservoir stress differences and interlayer thicknesses. The reservoir stress difference and interlayer thickness are then used as two-dimensional coordinate axes to plot the fracture height and cross-layer evaluation chart corresponding to the target location. The combined analysis refers to the sensitivity analysis of horizontal well multi-cluster fracturing construction parameters to induced bedding opening and cross-layer propagation.

[0060] During hydraulic fracturing, the target location is usually selected as the sweet spot. Therefore, in step S100, a fracture height and cross-layer evaluation chart corresponding to each sweet spot location can be generated, that is, a fracture height and cross-layer evaluation chart under different lithological sequences, mechanical parameters and interlayer interface characteristics.

[0061] As can be seen, the monitoring results of physical simulation experiments on fracture propagation can provide information such as the length, width, and height of artificially created fractures. Correspondingly, the first fracture propagation analysis results can reflect the laws and main controlling factors of artificial fracture competition initiation, penetration, and induced bedding opening under different laminar flow characteristics. The output results of numerical simulation experiments on multi-cluster fracturing fracture propagation in horizontal wells can provide information such as fracture morphology. Correspondingly, the second fracture propagation analysis results can reflect the influence of fracture morphology and fracture height propagation under different laminar flow characteristics.

[0062] It is important to understand that, in order to construct a geomechanical model of the target location in the study area, it is necessary to quantitatively characterize the laminar mechanical features within the laminar reservoir at the target location. These laminar mechanical features include not only in-situ stress characteristics and rock matrix mechanical characteristics, but also the interlayer interface mechanical characteristics. For example, the process of constructing a geomechanical model of the target location in the study area may include the following implementation steps:

[0063] By conducting rock mechanics tests at target locations in the study area, quantitative characterization of microfractures in rock cores, and dynamic rock mechanics calculations, the mechanical characteristics of interlayer interfaces, rock matrix mechanical characteristics, and geostress characteristics under different lithologies and bedding were quantitatively characterized.

[0064] Based on the interlayer interface mechanical characteristics, rock matrix mechanical characteristics, and geostress characteristics obtained from quantitative characterization, a geomechanical model corresponding to the target location is constructed.

[0065] For example, in one specific embodiment, the process for determining the result of the first crack propagation analysis is as follows:

[0066] Based on the rock samples obtained from the target location in the study area, fracturing samples with different laminar characteristics were constructed. Physical simulation experiments of multi-cluster fracturing fracture propagation in horizontal wells were carried out on the fracturing samples to obtain monitoring data of the experimental process.

[0067] The first crack propagation analysis results were generated based on the monitoring data.

[0068] As is known, rock samples from the target location in the study area are obtained through rock sampling, typically from downhole cores and outcrops. Based on the lithological sequence, mechanical parameters, and interlayer interface characteristics of the rock samples, systematic testing and benchmarking can be performed to construct fracturing specimens with different laminar characteristics. Monitoring data during the experiment can be obtained using acoustic emission technology and high-precision CT scanning technology, thereby enabling quantitative characterization of the length, width, and height of the formed artificial fractures.

[0069] For example, in one specific embodiment, the process for determining the second crack propagation analysis result is as follows:

[0070] After characterizing the bedding features of the target location in the study area, a three-dimensional model of crack propagation was constructed.

[0071] Numerical simulation experiments on the propagation of multiple clusters of hydraulic fracturing fractures in horizontal wells under different lamination characteristics were carried out using a three-dimensional model of fracture propagation, and the output data of the experimental process were obtained.

[0072] The second crack propagation analysis results are generated based on the output data.

[0073] It is important to understand that bedding characteristic parameters include bedding density, bedding width, bedding spacing, bedding permeability, and bedding mechanical properties. When characterizing these bedding characteristic parameters, since the three-dimensional model of fracture propagation cannot be set according to the actual number of bedding lines in the strata, multiple bedding lines in the actual strata need to be equivalently represented as a single bedding line in the model. Therefore, in a specific embodiment, the bedding characteristic parameters are equivalently characterized, specifically including the following implementation steps:

[0074] The total permeability of multiple bedding planes in the target location of the study area is equivalently represented as the permeability of a single bedding plane, thus obtaining the equivalent bedding plane permeability.

[0075] The bedding density of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding permeability, and the equivalent bedding density is obtained.

[0076] The equivalent bedding spacing of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding density.

[0077] The equivalent bedding width is obtained by equivalent characterizing the bedding width of the strata at the target location in the study area based on the equivalent bedding permeability.

[0078] For example, in a specific embodiment, the equivalent characterization process of the bedding characteristic parameters, combined with the specific characteristics of the layered shale oil reservoir, is as follows:

[0079] 1) The equivalent bedding density is calculated based on the relationship between bedding density and bedding permeability, as shown in the following formula:

[0080] ρ′ BP =ρ BP ×k BP / k′ BP (Formula 1);

[0081] In Formula 1, ρ′ BP ρ represents the equivalent bedding density in the crack propagation 3D model. BP k represents the bedding density of the actual strata at the target location in the study area. BPk′ represents the bedding permeability of the actual strata at the target location in the study area. BP Equivalent bedding permeability in the three-dimensional model of crack propagation;

[0082] 2) The equivalent bedding spacing is calculated based on the relationship between bedding spacing and bedding density, as shown in the following formula:

[0083] d BP =1 / ρ′ BP (Formula 2);

[0084] In Formula 2, d BP The equivalent bedding spacing in the three-dimensional model of crack propagation;

[0085] 3) The equivalent bedding width is calculated based on the relationship between bedding width and bedding permeability, as shown in the following formula:

[0086]

[0087] In Formula 3, w′ BP The equivalent bedding width in the crack propagation 3D model.

[0088] As an improvement to the above embodiment, a three-dimensional model of crack propagation is constructed after characterizing the bedding feature parameters of the target location in the study area. This specifically includes the following implementation steps:

[0089] Step S01: Based on the discrete element mechanics theory, establish a discrete element mechanics model that considers the stratification of the target location in the study area.

[0090] As is known, in order to implement numerical simulation of hydraulic fracturing fracture propagation, rock deformation and rock fracture criteria are usually characterized by governing equations.

[0091] For example, in reservoirs, the block motion of the rock matrix obeys Newton's second law. Considering the deformability of the blocks, it can be represented by a small-displacement linear elastic dynamic equation, as shown in the following formula:

[0092] σ ij,j +b i -ρu i,tt -αu i,t =0 (Formula 4);

[0093] In Formula 4, σ ij =D ij,st ×ε st , σ ij,j Represents the Cauchy tensor; b i The unit of force is N; ρ represents the density of rock, with units of kg / m³. 3 ; α represents the damping coefficient; ui,t u i,tt u i,j and u j,i ε represents displacement, in mm. ij and ε st D represents strain, in mm. ij,st This represents the Hooke tensor.

[0094] The rock fracture criterion can be expressed as follows:

[0095]

[0096] In Formula 5, F n F represents the normal force. s This represents tangential force, measured in N; K n Indicates the normal spring stiffness; K s Indicates the tangential spring stiffness; Δu n Δu represents the normal relative displacement between adjacent nodes, in mm. s This represents the tangential relative displacement between adjacent nodes, in mm.

[0097] Step S02: For the multi-cluster fracturing problem, the injected fracturing fluid passes through the wellbore and perforation before entering each cluster of fractures. The fluid inflow into each cluster of fractures is controlled by the "wellbore-perforation-fracture" system.

[0098] In one specific embodiment, since the flow friction in the horizontal wellbore is relatively small, the frictional resistance along the wellbore is temporarily ignored. Accordingly, the inlet pressure of each cluster of fractures satisfies the following formula:

[0099] p w =p p,k +p in，k k = 1, 2, ..., n f (Formula 6);

[0100] In Formula 6, p p,k p represents the perforation friction of the k-crack, in MPa. in,k p represents the inlet pressure of crack k, in MPa. w This represents the pressure at the crack inlet, in MPa; n f Indicates the number of cracks.

[0101] If the fracturing fluid within the fracture is considered as an incompressible Newtonian fluid flowing laminarly within a flat plate, then it satisfies the continuity equation shown in Equation 7 and the global mass conservation equation shown in Equation 8:

[0102]

[0103]

[0104] In Formulas 7 and 8, p represents fluid pressure in MPa; w represents dynamic crack width in meters; t represents time in seconds; s represents the coordinates of any point within the crack; q l This indicates the fracturing fluid filtrate loss, in units of m³. 3 / s, because the matrix has ultra-low permeability, the effect of fracturing fluid loss is ignored in the examples, i.e., q l =0; Q0 represents the construction displacement, in meters. 3 / s; μ represents fluid viscosity, with units of Pa·s.

[0105] Step S03: Couple rock deformation and fluid flow within the crack to establish a fluid-structure interaction model for crack propagation.

[0106] Since rock deformation and fluid flow within the fracture interact, it is necessary to couple rock deformation with fluid flow within complex fractures to establish a three-dimensional model of fracture propagation considering fluid-structure interaction. The relationship between fracturing fluid flow and matrix block deformation within the fracture is as follows: 1) Changes in fluid pressure at the matrix block boundary will affect the deformation of the matrix block, thus causing changes in fracture width; 2) Changes in fracture width cause changes in flow rate, ultimately affecting the distribution of fluid pressure within the fracture. A weak coupling method is used to achieve the above iterative process. Before rock fracturing, the permeability (k) of the initial fracture (bedding, natural fracture) unit is... f ) should be related to the matrix permeability (k m Since the permeability and width of the fracture are equal, the relationship between the fracture permeability and width, k, can be used. f =w 2 / 12 gives the initial width of the crack. To provide an initial flow path for the fracturing fluid, the fracture width is constant within a single time step of calling the fluid flow solver. The pressure distribution in the fracture is solved using equations and converted into external loads on the surface nodes of the matrix block fracture. Then, the stress, strain, displacement, and normal relative displacement Δu between adjacent blocks are calculated for the block system under fluid pressure and initial stress field. n Tangential relative displacement Δu s The dynamic crack width w and shear displacement of the hydraulic crack are used as known conditions for the next fluid calculation time step. Then, the coupled iteration of the fluid pressure field and matrix block deformation is calculated, and the pressure convergence of all nodes is used as the judgment condition.

[0107] In one specific embodiment, the Picard iterative method can be used to solve the fluid-structure interaction problem, as shown in the following equation:

[0108]

[0109] In Formula 9, p represents the fluid pressure, in MPa; w represents the slit width, in mm; pk′+1 This represents the pressure at the next iteration step, numbered k′+1, in MPa; u k′+1 This represents the displacement of the next iteration step, numbered k′+1, in mm; w k′+1 The slit width at the next iteration step (k′+1) is represented in meters (m); t represents time (s); A represents the global stiffness matrix; α′ represents the empirical coefficient; w k′ This represents the slit width in meters (m) for the current iteration step, k′. k′ The duration of the current iteration step, k′, is represented in seconds; F represents the frictional force, in MPa.

[0110] Furthermore, the aforementioned governing equations can be solved using the dynamic relaxation method, and the aforementioned continuity equations can be calculated using the finite element method. Both are solved using the aforementioned Picard iteration method.

[0111] Step S04: Characterize the bedding characteristic parameters of strata with bedding development at different scales.

[0112] In the above technical solution, the discrete element method is used to overcome the problem of medium discontinuity caused by stratification, improve the accuracy of the three-dimensional model of crack propagation, and thus improve the accuracy of the cross-layer evaluation results based on the crack height cross-layer evaluation chart.

[0113] In one specific embodiment, the condition for cross-layer propagation includes a minimum fracturing displacement. For example... Figure 7 As shown, the evaluation chart for fracture height and interlayer penetration is plotted with the stress difference between the reservoir and interlayer as the horizontal axis and the interlayer thickness as the vertical axis. Figure 7 The condition shown for cross-layer propagation is the minimum hydraulic fracturing displacement required to allow the fracture to penetrate the layer. Figure 7 Based on the differences in reservoir stress and interlayer thickness, the fracture height penetration evaluation chart is divided into a direct penetration zone (penetration zone), a cross-shaped penetration zone (transition zone), and an inaccessible suppression zone.

[0114] Step S200: Use the fracture height and layer-penetration evaluation chart to conduct a layer-penetration evaluation before fracturing construction in the production block to be built.

[0115] As can be seen, after conducting cross-layer evaluation on the block to be developed, the fracturing construction parameters can be designed or the existing fracturing construction parameters can be optimized based on the cross-layer evaluation results. After fracturing based on the designed or optimized fracturing construction parameters, the post-fracturing production capacity can be evaluated, thereby determining the effectiveness of the fracture height cross-layer evaluation method implemented in this embodiment for production capacity optimization.

[0116] See Figure 2 The following describes the application of the above embodiments to the lamellar shale oil reservoir in the first specific region, including the following implementation process:

[0117] Step A1, construct a high-precision geomechanical model:

[0118] Based on core observations, CT scans, and well logging curve interpretations at various sweet spot locations, refined core descriptions and micro- and macroscopic testing analyses were conducted to clarify the interfacial mechanical characteristics, rock matrix mechanical characteristics, and geostress characteristics under different laminae and lithologies. A high-precision geomechanical model with multiple vertically continuous sublayers was constructed for different sweet spot locations in the region.

[0119] Step A2: Conduct a physical simulation experiment on crack propagation.

[0120] Based on the systematic testing and benchmarking of lithological sequence, mechanical parameters, and interlayer interface characteristics of downhole cores and outcrops, fracturing samples with different lamination characteristics were constructed. Physical simulation experiments of multi-cluster fracture propagation in horizontal well sections were carried out using a one-step segmented fracturing string and a true triaxial layered pressurization device. Based on acoustic emission monitoring and high-precision CT scanning technology, the propagation morphology of artificial fractures was quantitatively characterized, and the laws and main controlling factors of artificial fracture competition initiation, penetration, and induced lamination opening under different lamination characteristics were clarified.

[0121] Step A3: Conduct numerical simulation experiments on the propagation of multiple clusters of hydraulic fracturing fractures in horizontal wells.

[0122] By comprehensively applying outcrop profiles and core CT scan images, an equivalent characterization method for bedding feature parameters is established. The discrete element method is used to solve the problem of medium discontinuity caused by bedding, thereby constructing a three-dimensional model of fracture propagation. Using the constructed three-dimensional model of fracture propagation, numerical simulation experiments of multi-cluster fracture propagation in horizontal well sections under different bedding characteristics are carried out to quantitatively evaluate fracture morphology and understand the influence of fracture height propagation.

[0123] Step A4, create a seam height and layer penetration evaluation chart:

[0124] By combining and analyzing the geomechanical models and hydraulic fracture propagation laws obtained in steps A1 to A3, the conditions for hydraulic fracture-induced bedding opening and cross-layer propagation under different reservoir stress differences and interlayer thicknesses are clarified. Curves are plotted with reservoir stress difference as the abscissa and interlayer thickness as the ordinate to obtain shale oil fracture high cross-layer evaluation charts under different lithological sequences, mechanical characteristics, and interlayer interface characteristics.

[0125] It is important to understand that the development of reservoir-interstitial layers varies significantly at different sweet spot locations in lamellar shale oil reservoirs, and the distribution of microseismic event points differs markedly across different wellbore trajectories or perforation locations. This directly impacts the morphology of artificial fracture propagation and the effectiveness of the stimulation. Therefore, different sweet spot locations correspond to different fracture height and cross-layer evaluation charts. Based on this, a fracture height and cross-layer evaluation chart for a specific sweet spot location encompasses the lithological sequence, mechanical parameters, interlayer interface characteristics, reservoir-interstitial layer stress difference, interstitial layer thickness, and conditions for induced bedding opening and cross-layer propagation at that location.

[0126] For example, in this application example, the geological and engineering parameters covered by the evaluation chart of the fracture height of a certain dessert location are shown in Table 1.

[0127] Table 1

[0128]

[0129]

[0130] Step A5 guides the design or optimization of fracturing schemes for the blocks to be developed:

[0131] After obtaining the reservoir-interstitial stress difference and interstitial thickness of the proposed production area, a corresponding fracture height cross-layer evaluation chart is selected based on the lithological sequence, mechanical parameters, and interlayer interface characteristics of the proposed production block for cross-layer evaluation. The cross-layer evaluation results are then used to design or optimize fracturing parameters before fracturing operations. For example, based on the corresponding fracture height cross-layer evaluation chart, after knowing the reservoir-interstitial stress difference and interstitial thickness of the proposed production block, it is determined whether cross-layer fracturing is feasible. If cross-layer fracturing is feasible, the minimum fracturing displacement required for cross-layer fracturing is determined from the fracture height cross-layer evaluation chart.

[0132] Taking well T1 in the region as an example, the comparison between the optimized fracturing parameters and production results is shown in Table 2:

[0133] Table 2

[0134]

[0135] As shown in Table 1, calculations show that the stress difference in the reservoir layer where well T1 is located is 8 MPa, and the layer thickness is 4 m. Based on the fracture height and layer crossing evaluation chart corresponding to well T1, the single-fracture drilling displacement needs to be maintained at 3.5 m³ / s. 3 Only hydraulic fractures with a speed of 8-12 m / min or higher can achieve cross-layer fracturing. Before optimizing the hydraulic fracturing construction parameters, well D1 used 8-12 m / min. 3 Fracturing operations were performed using a displacement of / min, with a single-stage design of 5 clusters. After optimization, well T1 adopted a displacement of 15-18m. 3 The fracturing operation was carried out at a displacement of / min, and the design still adopted a single-stage 5-cluster configuration. After optimizing the fracturing parameters, the fracture height in a typical section of well T1 increased by 21.4% compared to well D1 without optimized parameters. Compared to well D1, with the same cumulative production days, well T1's average daily oil production increased by 30.1. Therefore, it is evident that the optimized fracturing parameters based on the above embodiments significantly improved fracturing production efficiency.

[0136] Combine Figures 2 to 7 As shown, the following content describes the application of the above embodiments in the lamellar shale oil reservoir in the second specific region, including the following implementation process:

[0137] Step B1, construct a high-precision geomechanical model:

[0138] Based on core observations, CT scans, and well logging interpretations at various sweet spot locations, refined core descriptions and micro-to-macroscopic testing analyses were conducted to clarify the interfacial mechanical characteristics, matrix mechanical characteristics, and in-situ stress characteristics under different laminae and lithologies. A high-precision geomechanical model of multiple sublayers with continuous longitudinal stratigraphy was constructed for different sweet spot locations in this area. (See details at [link to relevant documentation]). Figure 3 and Figure 4 As shown.

[0139] Step B2: Conduct a physical simulation experiment on crack propagation.

[0140] Based on systematic testing and benchmarking of lithological sequence, mechanical parameters, and interlayer interface characteristics of downhole cores and outcrops, fracturing samples with different lamination characteristics were constructed. Physical simulation experiments of multi-cluster fracture propagation in horizontal well sections were conducted using a one-step segmented fracturing string and a true triaxial layered pressurization device. The morphology of artificial fracture propagation was quantitatively characterized using acoustic emission monitoring and high-precision CT scanning technology. The laws and main controlling factors of artificial fracture competition initiation, trans-layer movement, and induced lamination opening under different lamination characteristics were clarified. For details, please refer to [link to relevant documentation]. Figure 5 As shown.

[0141] Step B3: Conduct numerical simulation experiments on the propagation of multiple clusters of hydraulic fracturing fractures in a horizontal well.

[0142] By comprehensively utilizing outcrop profiles and core CT scan images, an equivalent characterization method for bedding feature parameters was established. The discrete element method (DEM) was employed to address the medium discontinuity problem caused by bedding, thereby constructing a three-dimensional model of fracture propagation. Using this model, numerical simulation experiments were conducted on multi-cluster fracture propagation within horizontal well sections under different bedding characteristics to quantitatively evaluate fracture morphology and understand the influence of fracture height propagation. For details, please refer to [link to relevant documentation]. Figure 6 As shown.

[0143] Step B4, establish a seam height and layer-penetration evaluation chart:

[0144] Based on the geomechanical model and hydraulic fracture propagation law obtained in steps B1 to B3, the conditions for hydraulic fracture-induced bedding opening and cross-layer propagation under different reservoir stress differences and interlayer thicknesses were clarified. Curves were plotted with reservoir stress difference as the x-axis and interlayer thickness as the y-axis to obtain evaluation charts for high-penetration shale oil fractures under different lithological sequences, mechanical parameters, and interlayer interface characteristics. For details, please refer to [reference needed]. Figure 7 As shown.

[0145] Step B5 guides the design or optimization of fracturing schemes for the blocks to be developed:

[0146] After obtaining the reservoir stress difference and interlayer thickness of the production area to be developed, the corresponding fracture height cross-layer evaluation chart is selected based on the lithological sequence, mechanical parameters and interlayer interface characteristics of the production area to be developed for cross-layer evaluation. The cross-layer evaluation results are obtained, and the fracturing construction parameters are designed or optimized before fracturing construction based on the cross-layer evaluation results.

[0147] For example, if a reservoir in a block to be developed is located in the lower sweet spot, the corresponding fracture height cross-layer evaluation chart is as follows: Figure 7 As shown, based on the well network and spacing deployment and the distribution of interlayers, and after knowing the stress difference and thickness of the reservoir interlayers, it was determined that the artificial fractures in this reservoir are mainly for localized penetration. Under the premise of overall utilization, it is necessary to ensure that the artificial fractures can penetrate the layers and guarantee the height and conductivity of the artificial fractures. Therefore, the horizontal wells in this block should appropriately increase the drilling flow rate and the proportion of pre-fracturing fluid. Based on this principle, the existing fracturing operation parameters should be optimized.

[0148] Taking wells V1 and J1 in this reservoir as examples, the comparison of reservoir physical parameters, reservoir thickness and mechanical properties between wells V1 and J1 is shown in Table 3.

[0149] Table 3

[0150]

[0151] Taking well V1 as an example, the comparison between the optimized fracturing construction parameters and production results is shown in Table 4.

[0152] Table 4

[0153]

[0154]

[0155] As shown in Table 3, calculations show that the stress difference in the reservoir layer where well V1 is located is 9 MPa, and the layer thickness is 4.5 m. Based on the fracture height and layer crossing evaluation chart corresponding to well V1, the single-fracture drilling displacement needs to be maintained at 4.0 m³ / s. 3 Only hydraulic fractures with a speed of over 10 m / min can achieve cross-layer penetration. Before optimizing the hydraulic fracturing parameters, well J1 used 10-12 m / min. 3 Fracturing operations were performed using a displacement of / min, with a single-stage 4-cluster design. After optimization, well V1 adopted a displacement of 16-18m. 3 The fracturing operation was carried out at a displacement of / min, and the design still adopted a single-stage 4-cluster configuration. After optimizing the fracturing parameters, the fracture height in a typical section of well V1 increased by 50% compared to well J1 without optimized parameters. Compared to well J1, with the same cumulative production days, well V1's average daily oil production increased by 29.2%. Therefore, it is evident that the optimized fracturing parameters based on the above embodiments significantly improved fracturing production efficiency.

[0156] Example 2

[0157] See Figure 8 This invention provides a device 400 for evaluating high-penetration fractures in laminar reservoirs, comprising a fracture high-penetration evaluation map construction module 410 and a penetration evaluation module 420 connected to each other, wherein:

[0158] The fracture height and cross-layer evaluation chart construction module 410 is used to combine and analyze the geomechanical model based on the target location of the study area, the results of the first fracture propagation analysis and the results of the second fracture propagation analysis, to determine the conditions for fracture-induced bedding opening and cross-layer propagation under different reservoir stress differences and interlayer thicknesses, and to draw the fracture height and cross-layer evaluation chart corresponding to the target location using the reservoir stress difference and interlayer thickness as two-dimensional coordinate axes.

[0159] The cross-layer evaluation module 420 is used to evaluate the cross-layers before fracturing construction in the production block using the fracture height cross-layer evaluation chart.

[0160] The first fracture propagation analysis result was determined based on the monitoring data from the physical simulation experiment of fracture propagation under different lamination characteristics, and the second fracture propagation analysis result was determined based on the output data from the numerical simulation experiment of multi-cluster fracturing fracture propagation in horizontal wells under different lamination characteristics.

[0161] In one specific embodiment, the process of constructing the geomechanical model includes:

[0162] Quantitative characterization of the laminar mechanical features at the target location in the study area;

[0163] A geomechanical model was constructed based on the quantitatively characterized laminar mechanical features.

[0164] In one specific embodiment, the laminar mechanical characteristics include interlayer interface mechanical characteristics, rock matrix mechanical characteristics, and geostress characteristics.

[0165] In one specific embodiment, the process of determining the result of the first crack propagation analysis includes:

[0166] Based on the rock samples obtained from the target location in the study area, fracturing specimens with different lamination characteristics were constructed. Physical simulation experiments on crack propagation were carried out on the fracturing specimens to obtain monitoring data of the experimental process.

[0167] The first crack propagation analysis results were generated based on the monitoring data;

[0168] Among them, the results of the first crack propagation analysis reflect the laws and main controlling factors of artificial crack competition initiation, cross-layering, and induced bedding opening under different laminar characteristics.

[0169] In one specific embodiment, the process of determining the results of the second crack propagation analysis includes:

[0170] After characterizing the bedding features of the target location in the study area, a three-dimensional model of crack propagation was constructed.

[0171] Numerical simulation experiments on the propagation of multiple clusters of hydraulic fracturing fractures in horizontal wells under different lamination characteristics were carried out using a three-dimensional model of fracture propagation, and the output data of the experimental process were obtained.

[0172] The second crack propagation analysis results are generated based on the output data;

[0173] Among them, the results of the second crack propagation analysis reflect the crack morphology and the law affecting crack height propagation under different laminar features.

[0174] In one specific embodiment, after characterizing the bedding feature parameters of the target location in the study area, a three-dimensional model of crack propagation is constructed, including:

[0175] Based on discrete element mechanics theory, a discrete element mechanics model considering the bedding of the target location in the study area is established;

[0176] After the injected fracturing fluid passes through the wellbore and perforation, it enters each cluster of fractures. The fluid inflow into each cluster of fractures is controlled by the parameters of the wellbore, perforation, and fractures to implement multi-cluster fracturing.

[0177] A fluid-structure interaction model for crack propagation is established by coupling rock deformation and fluid flow within cracks.

[0178] Characterization of bedding feature parameters for strata with bedding development at different scales.

[0179] In one specific embodiment, the bedding characteristic parameters include bedding density, bedding width, bedding spacing, bedding permeability, and bedding mechanical properties.

[0180] In one specific embodiment, the bedding characteristic parameters of strata with bedding development at different scales are characterized, including:

[0181] The total permeability of multiple bedding planes in the target location of the study area is equivalently represented as the permeability of a single bedding plane, thus obtaining the equivalent bedding plane permeability.

[0182] The bedding density of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding permeability, and the equivalent bedding density is obtained.

[0183] The equivalent bedding spacing of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding density.

[0184] The equivalent bedding width is obtained by equivalent characterizing the bedding width of the strata at the target location in the study area based on the equivalent bedding permeability.

[0185] In one specific embodiment, the condition for cross-layer propagation includes a minimum fracturing operation displacement.

[0186] In one specific embodiment, the fracture height permeation evaluation chart is divided into a direct permeation zone, a T-shaped permeation zone, and an impenetrable suppression zone based on the stress difference and thickness of the reservoir layer.

[0187] On the other hand, embodiments of the present invention also provide a machine-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the above-described methods for evaluating high-penetration fractures in lamellar reservoirs.

[0188] In another aspect, embodiments of the present invention also provide a processor for running a program, wherein the program executes the above-described methods for evaluating high-penetration fractures in lamellar reservoirs.

[0189] In another aspect, embodiments of the present invention also provide a computer device, which may be a terminal, and its internal structure diagram may be as follows: Figure 9 As shown in the figure, the computer device includes a processor A01, a network interface A02, a display screen A04, an input device A05, and a memory (not shown) connected via a system bus. The processor A01 provides computing and control capabilities. The memory includes internal memory A03 and a non-volatile storage medium A06. The non-volatile storage medium A06 stores an operating system B01 and a computer program B02. The internal memory A03 provides an environment for the operation of the operating system B01 and the computer program B02 stored in the non-volatile storage medium A06. The network interface A02 is used for communication with external terminals via a network connection. When the computer program is executed by the processor A01, it implements a method for evaluating high-penetration fractures in layered reservoirs. The display screen A04 can be a liquid crystal display or an e-ink display. The input device A05 can be a touch layer covering the display screen, buttons, a trackball, or a touchpad on the computer device casing, or an external keyboard, touchpad, or mouse.

[0190] In one embodiment, the high-penetration evaluation device 400 for laminated reservoir pressure fractures provided by the present invention can be implemented as a computer program, which can be implemented in, for example... Figure 9 The computer device shown runs on this device. The computer device's memory can store the various program modules that make up the high-penetration evaluation device 400 for laminated reservoir pressure fractures, for example, Figure 8 The diagram shows a high-penetration fracture evaluation template construction module 410 and a penetration evaluation module 420. The computer program, comprised of these modules, causes the processor to execute the steps in the high-penetration fracture evaluation method for laminar reservoir pressure fractures described in this specification.

[0191] Figure 9 The computer device shown can be used as follows Figure 8 The fracture high-penetration evaluation chart construction module 410 in the layered reservoir pressure fracture high-penetration evaluation device 400 shown executes step S100, and the computer equipment can execute step S200 through the penetration evaluation module 420.

[0192] This invention also provides a computer program product that, when executed on a data processing device, is suitable for executing a program that initializes the following method steps:

[0193] By combining the geomechanical model based on the target location of the study area, the results of the first fracture propagation analysis, and the results of the second fracture propagation analysis, the conditions for fracture-induced bedding opening and cross-layer propagation under different reservoir stress differences and interlayer thicknesses were determined. The reservoir stress difference and interlayer thickness were used as two-dimensional coordinate axes to draw the fracture height cross-layer evaluation chart corresponding to the target location.

[0194] Use the fracture height cross-layer evaluation chart to evaluate the cross-layers before fracturing construction in the production block to be built;

[0195] The first fracture propagation analysis result was determined based on the monitoring data from the physical simulation experiment of fracture propagation under different lamination characteristics, and the second fracture propagation analysis result was determined based on the output data from the numerical simulation experiment of multi-cluster fracturing fracture propagation in horizontal wells under different lamination characteristics.

[0196] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. 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. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0197] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.

[0198] The above are merely embodiments of the present invention and are not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principle of the present invention should be included within the scope of the claims of the present invention.

Claims

1. A method for evaluating high-penetration fractures in laminar reservoirs, characterized in that, The method includes: By combining the geomechanical model based on the target location of the study area, the results of the first fracture propagation analysis, and the results of the second fracture propagation analysis, the conditions for fracture-induced bedding opening and cross-layer propagation under different reservoir stress differences and interlayer thicknesses were determined. The reservoir stress difference and interlayer thickness were used as two-dimensional coordinate axes to draw the fracture height cross-layer evaluation chart corresponding to the target location. The aforementioned fracture height and layer-penetration evaluation chart is used to evaluate the layer-penetration before fracturing construction in the production block to be built. Specifically, the first fracture propagation analysis result is determined based on the monitoring data from the physical simulation experiment of fracture propagation under different lamination characteristics, and the second fracture propagation analysis result is determined based on the output data from the numerical simulation experiment of multi-cluster fracturing fracture propagation in horizontal wells under different lamination characteristics.

2. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 1, characterized in that, The process of constructing the geomechanical model includes: Quantitative characterization of the laminar mechanical features at the target location in the study area; A geomechanical model was constructed based on the quantitatively characterized laminar mechanical features.

3. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 2, characterized in that, The laminar mechanical characteristics include interlayer interface mechanical characteristics, rock matrix mechanical characteristics, and geostress characteristics.

4. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 1, characterized in that, The process for determining the results of the first crack propagation analysis includes: Based on the rock samples obtained from the target location in the study area, fracturing specimens with different lamination characteristics were constructed. Physical simulation experiments on crack propagation were carried out on the fracturing specimens to obtain monitoring data of the experimental process. The first crack propagation analysis result is generated based on the monitoring data; Among them, the results of the first crack propagation analysis reflect the laws and main controlling factors of artificial crack competition initiation, penetration, and induced bedding opening under different lamination characteristics.

5. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 1, characterized in that, The process for determining the results of the second crack propagation analysis includes: After characterizing the bedding features of the target location in the study area, a three-dimensional model of crack propagation was constructed. Numerical simulation experiments on the propagation of multiple clusters of hydraulic fracturing fractures in horizontal wells under different lamination characteristics were carried out using the aforementioned three-dimensional model of fracture propagation, and the output data of the experimental process were obtained. A second crack propagation analysis result is generated based on the output data; The second crack propagation analysis results reflect the crack morphology and the law affecting crack height propagation under different laminar flow characteristics.

6. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 5, characterized in that, The process of constructing a three-dimensional model of crack propagation after characterizing the bedding feature parameters of the target location in the study area includes: Based on discrete element mechanics theory, a discrete element mechanics model considering the bedding of the target location in the study area is established; After the injected fracturing fluid passes through the wellbore and perforation, it enters each cluster of fractures. The fluid inflow into each cluster of fractures is controlled by the parameters of the wellbore, perforation, and fractures to implement multi-cluster fracturing. A fluid-structure interaction model for crack propagation is established by coupling rock deformation and fluid flow within cracks. Characterization of bedding feature parameters for strata with bedding development at different scales.

7. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 5 or 6, characterized in that, The bedding characteristic parameters include bedding density, bedding width, bedding spacing, bedding permeability, and bedding mechanical properties.

8. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 7, characterized in that, The characterization of bedding characteristic parameters of strata with bedding development at different scales includes: The total permeability of multiple bedding planes in the target location of the study area is equivalently represented as the permeability of a single bedding plane, thus obtaining the equivalent bedding plane permeability. The bedding density of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding permeability, and the equivalent bedding density is obtained. The equivalent bedding spacing of the strata at the target location in the study area is equivalently characterized based on the equivalent bedding density. The equivalent bedding width is obtained by equivalent characterizing the bedding width of the strata at the target location in the study area based on the equivalent bedding permeability.

9. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 1, characterized in that, The conditions for cross-layer propagation include the minimum displacement during fracturing operations.

10. The method for evaluating high-penetration fractures in laminar reservoirs according to claim 1, characterized in that, In the aforementioned fracture height penetration evaluation chart, based on the differences in reservoir stress difference and interlayer thickness, it is divided into a direct penetration zone, a T-shaped penetration zone, and an inaccessible suppression zone.

11. A device for evaluating high-penetration fractures in laminar reservoirs, characterized in that, The device includes: The fracture height and cross-layer evaluation chart construction module is used to combine and analyze the geomechanical model based on the target location of the study area, the results of the first fracture propagation analysis and the results of the second fracture propagation analysis, to determine the conditions for fracture-induced bedding and cross-layer propagation under different reservoir stress differences and interlayer thicknesses, and to draw the fracture height and cross-layer evaluation chart corresponding to the target location using the reservoir stress difference and interlayer thickness as two-dimensional coordinate axes. The layer-penetration evaluation module is used to evaluate the layer-penetration of the production block before fracturing construction using the fracture height layer-penetration evaluation chart. Specifically, the first fracture propagation analysis result is determined based on the monitoring data from the physical simulation experiment of fracture propagation under different lamination characteristics, and the second fracture propagation analysis result is determined based on the output data from the numerical simulation experiment of multi-cluster fracturing fracture propagation in horizontal wells under different lamination characteristics.

12. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the evaluation method for high-penetration fractures in laminar reservoirs as described in any one of claims 1 to 10.

13. A machine-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the evaluation method for high-penetration fractures in laminar reservoirs as described in any one of claims 1 to 10.

Images