Method and storage medium for determining well layout pattern based on shale gas development

Abstract

The embodiment of the present invention provides a method and storage medium for determining a well layout pattern based on shale gas development, which belongs to the field of shale gas exploration technology. The method includes: obtaining stress state information corresponding to the target layer in the target well layout area; when the stress state information meets the preset stress state condition, determining the fracture stability information of the target layer; when the fracture stability information meets the preset formation stability condition, and the formation trend information of the target layer meets the preset formation trend condition, determining the fault data corresponding to each grid unit in the geological grid model, wherein the geological grid model is generated based on the target well layout area, and the fault data includes the target fault distance and the fault stress angle; based on the fault data, determining the well layout pattern corresponding to the target layer. The well layout design process of the present invention fully combines the stress characteristics and geological structure characteristics of the mountain, thereby effectively improving the degree of adaptation between the formation and the well layout pattern, thereby effectively improving the efficiency of shale gas extraction.

Description

Technical Field

[0001] The present invention relates to the technical field of shale gas exploration, and in particular to a method for determining a well layout pattern based on shale gas development and a storage medium. Background Art

[0002] In the process of realizing the present invention, the inventors of the present application discovered that currently, a single well type is usually used to lay horizontal wells in the mining area along the direction of the minimum horizontal principal stress. However, for areas with complex terrain and structure, a single well type cannot fully adapt to the geological structure, stress characteristics and mining technology requirements, which may affect the mining efficiency. For example, when the reservoir is in a reverse fault stress state, horizontal fractures are easily generated during the fracturing process, thereby affecting the gas well transformation effect and causing casing damage. In addition, for reservoirs with a large formation dip angle, due to the large elevation difference between the heel and toe ends of the horizontal section of the horizontal well, liquid accumulation and sand deposition are easy to occur near the heel end of the up-dip well, and it is also difficult for the down-dip well to return liquid, thereby affecting the production volume of the gas reservoir, resulting in low shale gas production efficiency.

[0003] Other studies have shown that in addition to being affected by the drilling rate of high-quality reservoirs, the single well productivity is also mainly affected by the orientation of the horizontal well trajectory. Rock mechanics experiments and fracture monitoring show that when the horizontal well orientation is perpendicular to the maximum horizontal geostress orientation, the rock is most likely to crack and form a fracture system. However, due to the complex structural characteristics, the superposition of multiple tectonic movements and geostress, and in order to maximize resource utilization, it is difficult to deploy wells strictly in accordance with the orientation of the two being perpendicular to each other. Therefore, the appropriate well type and horizontal well trajectory design are the key to achieving effective utilization of shale gas resources. Summary of the invention

[0004] The purpose of the embodiments of the present invention is to provide a method for determining a well layout pattern based on shale gas development and a storage medium, so as to alleviate the current technical problem of low shale gas extraction efficiency.

[0005] In order to achieve the above-mentioned object, one aspect of an embodiment of the present invention provides a method for determining a well layout pattern based on shale gas development, comprising:

[0006] Obtaining stress state information corresponding to the target layer in the target well layout area;

[0007] When the stress state information meets the preset stress state condition, determining the fracture stability information of the target layer;

[0008] When the fracture stability information satisfies the preset formation stability condition, and the formation trend information of the target layer satisfies the preset formation trend condition, the fault data corresponding to each grid unit in the geological grid model is determined, wherein the geological grid model is generated based on the target well layout area, and the fault data includes the target fault distance and the fault stress angle;

[0009] Based on the fault data, the well layout pattern corresponding to the target layer is determined.

[0010] Optionally, before the step of obtaining the stress state information corresponding to the target layer in the target well arrangement area, the method further includes:

[0011] Obtaining the formation pressure of the target layer in the preset area;

[0012] When the formation pressure meets the pressure saturation condition, the preset area is determined as the target well deployment area.

[0013] Optionally, the step of obtaining stress state information corresponding to the target layer in the target well layout area includes:

[0014] Calculate the stress value of the target layer in the target well layout area to obtain the stress index of the target layer, where the stress value includes the average horizontal maximum principal stress, the average horizontal minimum principal stress and the average vertical stress;

[0015] According to the stress index, the stress state information of the target layer is determined.

[0016] Optionally, when the stress state information satisfies a preset stress state condition, the step of determining the fracture stability information of the target layer includes:

[0017] When the stress state information is of a normal fault type or a strike-slip fault type, determining that the stress state information satisfies a preset stress state condition; and

[0018] A fracture stability index of a target layer is determined, wherein the fracture stability information includes the fracture stability index.

[0019] Optionally, when the fracture stability information satisfies the preset formation stability condition and the formation trend information of the target layer satisfies the preset formation trend condition, before the step of determining the fault data corresponding to each grid unit in the geological grid model, the method further includes:

[0020] Get the formation dip of the target layer.

[0021] Optionally, when the fracture stability information satisfies a preset formation stability condition, and the formation trend information of the target layer satisfies a preset formation trend condition, the step of determining the fault data corresponding to each grid unit in the geological grid model includes:

[0022] When the fracture stability index of the target layer is less than the fracture stability judgment coefficient, it is determined that the fracture stability information meets the preset formation stability condition;

[0023] When the formation dip angle of the target layer is less than or equal to a preset angle, determining that the formation trend information of the target layer satisfies a preset formation trend condition; and

[0024] Based on the fault distribution information, the target fault distance and fault stress angle corresponding to each grid unit in the geological grid model are determined.

[0025] Optionally, the step of determining the target fault distance and fault stress angle corresponding to each grid unit in the geological grid model based on the fault distribution information includes:

[0026] Obtain the shortest distance between the center point of each grid unit in the geological grid model and the distribution position of each fault;

[0027] Comparing the shortest distances to determine a target fault distance, wherein the target fault distance includes the shortest distance with the smallest value; and

[0028] The angle between the minimum horizontal principal stress direction and the fault strike corresponding to each grid unit is determined, wherein the fault stress angle includes the angle between the minimum horizontal principal stress direction and the fault strike.

[0029] Optionally, the step of determining a well layout pattern corresponding to a target layer based on the fault data includes:

[0030] The target fault distance and fault stress angle are calculated to obtain the well placement judgment index of each grid unit;

[0031] Based on the well layout judgment index, the well layout pattern of the target layer corresponding to each grid unit is determined, wherein the well layout pattern includes the well layout direction and the well layout type.

[0032] Optionally, after the step of determining the well layout pattern corresponding to the target layer based on the fault data, the method further includes:

[0033] According to the well layout pattern, a well layout decision diagram is generated.

[0034] On the other hand, the present invention provides a machine-readable storage medium having instructions stored thereon, which, when executed by a processor, enables the processor to implement any of the above-mentioned methods for determining a well layout pattern based on shale gas development in the present application.

[0035] [The present invention has the following beneficial effects] Through the above technical scheme, the well layout pattern matching the target layer is determined according to the stress state information and fracture stability information in the target well layout area, rather than considering a single horizontal well type; in complex structural areas with developed faults, considering the difficulty in implementing the horizontal well network, vertical well or directional well patterns are preferably designed; in addition to considering the well layout pattern matching the target layer, two well layout orientations along the direction of the minimum horizontal principal stress and the direction of the structural contour line are optimized on the basis of the preferred design of the horizontal well type; the well layout design process fully combines the stress characteristics and geological structure characteristics of the mountain, thereby effectively improving the adaptability of the formation and the well layout pattern, thereby effectively improving the shale gas extraction efficiency.

[0036] Other features and advantages of the embodiments of the present invention will be described in detail in the subsequent detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following specific implementations, they are used to explain the embodiments of the present invention, but do not constitute a limitation on the embodiments of the present invention. In the accompanying drawings:

[0038] Figure 1 It is a flow chart of a method for determining a well layout pattern based on shale gas development provided in an embodiment of the present invention.

[0039] Figure 2 Schematic diagram of a target fault distance determination process provided by an embodiment of the present invention.

[0040] Figure 3 It is a schematic diagram of a well layout decision diagram provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0041] The specific implementation of the embodiment of the present invention is described in detail below in conjunction with the accompanying drawings. It should be understood that the specific implementation described here is only used to illustrate and explain the embodiment of the present invention, and is not used to limit the embodiment of the present invention.

[0042] One aspect of an embodiment of the present invention provides a method for determining a well layout pattern based on shale gas development. Figure 1 As shown, Figure 1 : is a flow chart of a method for determining a well layout pattern based on shale gas development provided by an embodiment of the present invention. The specific process may be as follows:

[0043] S101. Obtain stress state information corresponding to the target layer in the target well layout area.

[0044] The target layer is a predetermined layer required to be drilled according to different drilling design tasks, which mostly refers to oil and gas reservoirs; the stress state information is the type of force acting on a unit area inside the medium caused by shale deformation.

[0045] Exploration practice shows that formation pressure will affect drilling efficiency, gas production stability and hydraulic fracturing effect: if the formation pressure is too low, it will increase the difficulty of drilling, thus affecting drilling efficiency. In addition, if wells are laid in areas with too low formation pressure, the bottom hole pressure is too low, which is prone to imbalance, thus affecting gas production rate and stability, and reducing the hydraulic fracturing effect, thus affecting the release and recovery of shale gas. In order to ensure the smooth progress of gas reservoir exploitation, before conducting well layout design, it is necessary to predetermine the target well layout area suitable for well layout and exploitation.

[0046] Specifically, in the embodiment of the present invention, before the above step S101, the following is further included:

[0047] Obtain the formation pressure P0 of the target layer in the preset area;

[0048] When the formation pressure meets the pressure saturation condition, the preset area is determined as the target well deployment area.

[0049] Among them, in the actual exploration process, when the formation pressure P0 of the target layer is greater than or equal to the saturation pressure P1, it is determined that it meets the pressure saturation condition, and the preset area to which the target layer belongs is determined as the target well layout area. Optionally, the formation pressure P0 of the target layer can be obtained through field interpretation data or adjacent well gas test data, and the pressure corresponding to the maximum excess adsorption amount measured in the methane high-pressure adsorption experiment in the target layer sample chamber is used as the saturation pressure P1. For example, the saturation pressure P1 = 10.0 MPa, the formation pressure P0 of the target layer in the first preset area = 8.0 MPa, the formation pressure P0 of the target layer in the second preset area = 10.6 MPa, the formation pressure P0 of the target layer in the third preset area = 11.3 MPa, and the formation pressure P0 of the target layer in the fourth preset area = 12.0 MPa. Since the formation pressure of the target layer in the second preset area, the third preset area and the fourth preset area is greater than the saturation pressure P1, the second preset area, the third preset area and the fourth preset area are determined as the target well layout area.

[0050] Further, in an embodiment of the present invention, the process of obtaining the stress state information corresponding to the target layer in the target well layout area specifically includes: firstly calculating the stress value of the target layer in the target well layout area to obtain the stress index of the target layer, and then determining the stress state information of the target layer according to the stress index. Optionally, the stress value includes the average horizontal maximum principal stress, the average horizontal minimum principal stress and the average vertical stress. In the actual exploration process, the average horizontal maximum principal stress, the average horizontal minimum principal stress and the average vertical stress can be obtained through the one-dimensional rock mechanics interpretation results of a single well. The calculation formula of the stress index Q is shown in formula (1):

[0051]

[0052] Among them, σ H is the average maximum principal stress, σ h is the average minimum principal stress, σ V is the mean vertical stress.

[0053] When the stress index Q≤0, the stress state information of the target layer is a reverse fault type (i.e., one side of the fault plane moves downward relative to the other side, forming a structural morphology in which the lower plate is above the upper plate); when the stress index 1<Q<0, the stress state information of the target layer is a normal fault type (i.e., one side of the fault plane moves upward relative to the other side, forming a structural morphology in which the upper plate is above the lower plate); when the stress index Q≥1, the stress state information of the target layer is a strike-slip fault type (i.e., the two sides of the fault plane move horizontally relative to each other, without obvious vertical displacement).

[0054] S102. When the stress state information satisfies a preset stress state condition, the fracture stability information of the target layer is determined.

[0055] The fracture stability information is used to characterize the fracture stability of the target layer. Specifically, in the embodiment of the present invention, when the stress state information is a normal fault type (i.e., 1<Q<0) or a strike-slip fault type (i.e., Q≥1), it is determined that the stress state information satisfies a preset stress state condition, and the fracture stability index of the target layer is determined, and the fracture stability index is used as the fracture stability information.

[0056] Optionally, the calculation formula of the crack stability index ST is as shown in formula (2):

[0057]

[0058] Where τ is the shear stress acting on the crack, σ n is the effective normal stress acting on the crack, μ is the internal friction coefficient of the crack, C is the internal friction angle of the crack, and the above τ, σ n , μ and C can all be obtained through indoor rock mechanics experiments on rock samples from the target layer in adjacent blocks.

[0059] Furthermore, when the stress state information does not meet the preset stress state condition, that is, when the stress state information of the target layer is of reverse fault type (i.e., Q≤0), the well layout pattern (including well layout direction and well layout type) corresponding to the target layer can be directly determined. Specifically, when the stress state information of the target layer is of reverse fault type, the well layout direction represented by the well layout pattern corresponding to the target layer is determined to be a vertical direction / other non-horizontal direction, and the well layout type is a non-horizontal well type (the specific non-horizontal well type can be selected according to actual geological conditions).

[0060] S103. When the fracture stability information meets the preset formation stability conditions and the formation trend information of the target layer meets the preset formation trend conditions, the fault data corresponding to each grid unit in the geological grid model is determined, wherein the geological grid model is generated based on the target well layout area, and the fault data includes the target fault distance and the fault stress angle.

[0061] Among them, in the actual application process, geological characteristic data such as geological structure, lithology distribution and reservoir properties of the target well layout area can be collected in advance, and the geological characteristic data can be modeled and processed through three-dimensional geological modeling software (for example, Petrel, GeoModeller, etc.) to generate a geological grid model. In addition, each grid unit in the geological grid model can be assigned corresponding attribute parameters (for example, lithology, formation thickness, etc.). The geological grid model is used to describe the geological characteristics of the target well layout area. By establishing and analyzing the geological grid model, exploration personnel can more accurately understand the geological conditions, lithology distribution and reservoir characteristics of the target well layout area, and provide an effective scientific basis for well layout design, fracturing construction and production management in the shale gas extraction process.

[0062] Specifically, the fracture stability judgment coefficient is usually used to evaluate the stability of the fracture, and the formation dip is usually used to evaluate the formation trend. In the actual exploration process, the fracture stability and trend of the formation have a vital impact on the research of well layout design. In an embodiment of the present invention, the formation dip of the target layer is obtained in advance. When the fracture stability index ST of the target layer is less than S0 (fracture stability judgment coefficient), it means that the fracture stability of the target layer is good, that is, it is determined that the fracture stability information meets the preset formation stability condition; when the formation dip of the target layer is θ≤θ0 (preset angle), it means that the formation structure of the target layer is wide or moderate, and it is determined that the formation trend information of the target layer meets the preset formation trend condition. Optionally, the formation dip of the target layer can be obtained by analyzing the seismic interpretation data of the target layer, analyzing the interpretation results of the formation wellbore imaging logging (FMI), or measuring in the well using a formation dip logging instrument.

[0063] Specifically, according to a large number of exploration studies and actual experience, when the fracture stability index gradually increases from 0.5, the fracture stability of the formation gradually deteriorates; conversely, when the fracture stability index gradually decreases from 0.5, the fracture stability of the formation gradually improves. In addition, the inclination angle corresponding to the formation with high and steep structure is usually greater than 20°, so in the embodiment of the present invention, the value of the fracture stability judgment coefficient S0 is set to 0.5, and the preset angle θ0 is set to 20°. It should be noted that in the actual application process, the actual geological conditions and shale characteristics and other factors can also be comprehensively considered to set the values ​​of the fracture stability judgment coefficient S0 and the preset angle θ0, so the values ​​of the fracture stability judgment coefficient S0 and the preset angle θ0 are not specifically limited in the embodiment of the present invention.

[0064] Next, the seismic data interpretation of the target well area is analyzed to obtain the fault distribution information such as the fault distribution position, fault direction and fault distribution in the area, and the target fault distance and fault stress angle corresponding to each grid unit in the geological grid model are determined based on the above fault distribution information. Specifically, first obtain the shortest distance between the center point of each grid unit in the geological grid model and each fault distribution position (i.e., fault line), and compare each shortest distance to determine the target fault distance (i.e., the shortest distance with the smallest value), so as to assign the target fault distance to the center point of the corresponding grid unit, and obtain the planar distribution of the target fault distance corresponding to each unit grid.

[0065] Optionally, a rectangular area is selected in the target well layout area along the due east direction (recorded as the X direction) and the due north direction (recorded as the Y direction), and the grid area is divided along the X direction and the Y direction in the geological grid model. The shortest line segment (representing the shortest distance) is drawn from the center point of each grid unit to all fault lines, and the length of each shortest line segment is recorded as D n (n is the number of fault lines), then the target fault distance D = min (D1, D2, D3, ..., D n ), and assign the target fault distance D to the center point of the corresponding grid unit. Figure 2 As shown, the coordinates of the center point P of the first grid unit are (i, j), the shortest distance between the center point P and the first fault distribution position is D1 = 10m, the shortest distance between the center point P and the second fault distribution position is D2 = 20m, the shortest distance between the center point P and the third fault distribution position is D3 = 15m, the shortest distance between the center point P and the fourth fault distribution position is D4 = 24m, the shortest distance between the center point P and the fifth fault distribution position is D5 = 12m, and the target fault distance D corresponding to the center point P is P =min(D1, D2, D3, D4, D5)=D1=10m, and D P =10m is assigned to the center point P of the first grid unit.

[0066] Furthermore, the angle between the minimum horizontal principal stress direction and the fault strike corresponding to each grid unit is determined, and the angle between the minimum horizontal principal stress direction and the fault strike is used as the fault stress angle, so as to assign the fault stress angle to the corresponding grid unit, and obtain the plane distribution of the fault stress angle corresponding to each unit grid. Optionally, in actual application, the minimum horizontal principal stress direction at each exploration well can be obtained based on the formation borehole imaging (FMI) interpretation results of the target well layout area, and the minimum horizontal principal stress direction corresponding to each grid unit in the geological grid model can be obtained by interpolation processing; the fault strike of each fault line can be obtained based on the seismic interpretation data of the target well layout area, and the fault strike corresponding to each grid unit in the geological grid model can be obtained by interpolation processing.

[0067] Specifically, the calculation formula of the fault stress angle A0 is shown in formula (3):

[0068] A0=|A σh -A f | Formula (3)

[0069] Among them, A σh is the direction of the minimum horizontal principal stress, A f is the fault direction, |A σh -A f | is the absolute value operation.

[0070] Furthermore, when the fracture stability information does not satisfy the preset formation stability condition (i.e., ST≥S0), or when the fracture stability information satisfies the preset formation stability condition (ST<S0) and the formation strike information does not satisfy the preset formation strike condition (θ>θ0), the well layout pattern corresponding to the target layer can be directly determined according to the formation dip. Specifically, when the fracture stability index of the target layer ST≥S0 (i.e., the fracture stability of the target layer is poor), if the formation dip θ≤5° (indicating that the formation structure is wide and gentle), then the well layout direction represented by the well layout pattern corresponding to the target layer is determined to be a vertical direction, and the well layout type is a non-horizontal well type (the specific non-horizontal well type can be selected according to the actual geological conditions); if the formation dip satisfies 5°<θ≤θ0 (indicating that the formation structure is moderate), then the well layout direction represented by the well layout pattern corresponding to the target layer is determined to be a directional well layout, and the well layout type is a non-horizontal well type. The type is a non-horizontal well type (the specific non-horizontal well type can be selected according to the actual geological conditions); if the formation dip angle θ>θ0 (indicating that the formation structure is high and steep), the well layout direction represented by the well layout mode corresponding to the target layer is determined to be directional, and the well layout type is a high-angle well type; when the fracture stability index ST of the target layer is less than S0 (that is, the fracture stability of the target layer is better), if the formation dip angle θ>θ0, the well layout direction represented by the well layout mode corresponding to the target layer is determined to be the structural contour line direction, and the well layout type is a horizontal well type.

[0071] S104. Determine the well layout pattern corresponding to the target layer based on the fault data.

[0072] Among them, the well layout pattern is the decision-making result of the well layout design work, and the well layout pattern includes the well layout direction and well layout type.

[0073] Specifically, the target fault distance and fault stress angle are first calculated to obtain the well layout judgment index of each grid unit, and the well layout pattern of the target layer corresponding to each grid unit is determined based on the well layout judgment index. The calculation formula of the well layout judgment index is shown in formula (4):

[0074] W=δ(A0-θ1)+2×δ(D-D0) Formula (4)

[0075] Among them, δ is defined as: θ1 is the critical angle; D0 is the critical fault spacing. Optionally, in the embodiment of the present invention, the value of the critical angle θ1 is set to 30°, and the value of the critical fault spacing D0 is set to 750m.

[0076] Furthermore, when W = -3 or 1, if the target fault distance D corresponding to the center point of the grid unit is less than d min (Minimum avoidance distance for engineering) is difficult to achieve construction at the engineering level, so the well is not arranged in the target layer corresponding to the grid unit. If the target fault distance D corresponding to the center point of the grid unit is greater than or equal to d min , then a horizontal well is arranged in the target layer corresponding to the grid unit along the direction parallel to the fault (that is, the well arrangement direction is parallel to the fault direction, and the well arrangement type is a horizontal well type). Optionally, d min The value of is 200m; when W=-1, a vertical well or a directional well is arranged in the target layer corresponding to the grid unit (that is, the well arrangement direction is the vertical direction / other non-horizontal direction, and the well arrangement type is a non-horizontal well type, and the specific non-horizontal well type can be selected according to the actual geological conditions); when W=3, a horizontal well is arranged in the target layer corresponding to the grid unit along the direction of the minimum horizontal principal stress (that is, the well arrangement direction is the direction of the minimum horizontal principal stress, and the well arrangement type is a horizontal well type).

[0077] For example, min =200m, if the well placement judgment index W1 of the first grid unit is -3, and the target fault distance D corresponding to the center point of the first grid unit is 100m, then no well is placed in the target layer corresponding to the first grid unit; if the well placement judgment index W2 of the second grid unit is 1, and the target fault distance D corresponding to the center point of the second grid unit is 220m, then a horizontal well is placed in the target layer corresponding to the second grid unit along the direction parallel to the fault; if the well placement judgment index W3 of the third grid unit is 3, then a horizontal well is placed in the target layer corresponding to the third grid unit along the direction of the minimum horizontal principal stress.

[0078] Alternatively, if Figure 3 As shown, a well layout decision diagram can be generated based on the formation pressure, stress state information and well layout pattern corresponding to each target layer obtained in the above process, so as to intuitively display the formation pressure, stress state information and well layout pattern and other information corresponding to each position of the target layer in the target well layout area through the well layout decision diagram, so that the exploration personnel can intuitively understand the best well layout plan corresponding to each position in the target well layout area, and carry out drilling work in a targeted manner, wherein, Figure 3 The dotted line is the formation pressure boundary line, the solid line is the formation stress state boundary line, the dot-dashed line is the fault line, and the solid line with numbers is the structural contour line.

[0079] The method for determining a well layout pattern based on shale gas development provided by the present invention determines a well layout pattern that matches the target layer according to stress state information and fracture stability information in a target well layout area, rather than considering a single horizontal well type; in a complex structural area with developed faults, considering that the implementation of a horizontal well network is difficult, a vertical well or directional well pattern is preferentially designed; in addition to considering a well layout pattern that matches the target layer, two well layout orientations along the direction of the minimum horizontal principal stress and the direction of the structural contour line are optimized and designed on the basis of the preferred design of the horizontal well type; the well layout design process fully combines the stress characteristics and geological structure characteristics of the mountain, thereby effectively improving the degree of adaptation between the formation and the well layout pattern, thereby effectively improving the shale gas extraction efficiency.

[0080] In addition, another aspect of an embodiment of the present invention provides a device for determining a well layout pattern based on shale gas development, which includes a processor and a memory. The steps in the above-mentioned method for determining a well layout pattern based on shale gas development are all stored in the memory as program units, and the processor executes the above-mentioned program units stored in the memory to implement corresponding functions.

[0081] The processor includes a kernel, which retrieves the corresponding program unit from the memory. One or more kernels can be set, and the efficiency of shale gas extraction can be improved by adjusting kernel parameters.

[0082] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0083] An embodiment of the present invention provides a storage medium on which a program is stored. When the program is executed by a processor, the method for determining a well layout pattern based on shale gas development is implemented.

[0084] An embodiment of the present invention provides a processor, which is used to run a program, wherein the program executes the above-mentioned method for determining a well layout pattern based on shale gas development when running.

[0085] The embodiment of the present invention provides a device, the device includes a processor, a memory, and a program stored in the memory and can be run on the processor, and the processor implements the following steps when executing the program: obtaining stress state information corresponding to the target layer in the target well layout area; when the stress state information meets the preset stress state condition, determining the fracture stability information of the target layer; when the fracture stability information meets the preset formation stability condition, and the formation trend information of the target layer meets the preset formation trend condition, determining the fault data corresponding to each grid unit in the geological grid model, wherein the geological grid model is generated based on the target well layout area, and the fault data includes the target fault distance and the fault stress angle; based on the fault data, determining the well layout mode corresponding to the target layer. The device in this article can be a server, a PC, a PAD, a mobile phone, etc.

[0086] The present application also provides a computer program product, which, when executed on a data processing device, is suitable for executing a program that is initialized with the following method steps: obtaining stress state information corresponding to a target layer in a target well layout area; when the stress state information satisfies a preset stress state condition, determining fracture stability information of the target layer; when the fracture stability information satisfies a preset formation stability condition, and the formation trend information of the target layer satisfies a preset formation trend condition, determining fault data corresponding to each grid unit in a geological grid model, wherein the geological grid model is generated based on the target well layout area, and the fault data includes a target fault distance and a fault stress angle; based on the fault data, determining a well layout pattern corresponding to the target layer.

[0087] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may 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 may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0088] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented 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 processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1A device that provides the functions specified in a block or multiple blocks.

[0089] These computer program instructions may also be stored in a computer readable memory capable of directing 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 an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0090] 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, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0091] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0092] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0093] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0094] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0095] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A method for determining a well layout pattern based on shale gas development, characterized in that: include: Obtaining stress state information corresponding to the target layer in the target well layout area; When the stress state information satisfies a preset stress state condition, determining the fracture stability information of the target layer; When the fracture stability information satisfies a preset formation stability condition, and the formation strike information of the target layer satisfies a preset formation strike condition, the fault data corresponding to each grid unit in the geological grid model is determined, wherein the geological grid model is generated based on the target well layout area, and the fault data includes a target fault distance and a fault stress angle; wherein the target fault distance is the minimum value of the shortest distances between the center point of each grid unit in the geological grid model and the distribution position of each fault; and the fault stress angle is the angle between the minimum horizontal principal stress direction corresponding to each grid unit in the geological grid model and the fault strike; Based on the fault data, determining a well layout pattern corresponding to the target layer; Wherein, when the stress state information satisfies a preset stress state condition, the step of determining the fracture stability information of the target layer comprises: When the stress state information is of a normal fault type or a strike-slip fault type, determining that the stress state information satisfies a preset stress state condition; and Determining a fracture stability index of the target layer, wherein the fracture stability information includes the fracture stability index; The crack stability index is determined according to the following formula: Where τ is the crack shear stress, σ n is the effective normal stress of the crack, μ is the friction coefficient of the crack, and C is the friction angle of the crack; The step of determining the well layout pattern corresponding to the target layer based on the fault data includes: Calculating the target fault distance and the fault stress angle to obtain a well placement judgment index for each grid unit; Based on the well layout judgment index, determine the well layout mode of the target layer corresponding to each grid unit, wherein the well layout mode includes a well layout direction and a well layout type; The well layout judgment index is determined according to the following formula: W=δ(A0-θ1)+2×δ(D-D0) Among them, δ is defined as: θ1 is the critical angle, D0 is the critical fault spacing, A0 is the fault stress angle, D is the target fault distance, and W is the well layout judgment index.

2. The method for determining a well layout pattern based on shale gas development according to claim 1, characterized in that: Before the step of obtaining the stress state information corresponding to the target layer in the target well arrangement area, the method further includes: Obtaining the formation pressure of the target layer in the preset area; When the formation pressure meets the pressure saturation condition, the preset area is determined as the target well deployment area.

3. The method for determining a well layout pattern based on shale gas development according to claim 2, characterized in that: The step of obtaining stress state information corresponding to the target layer in the target well layout area includes: Calculating the stress value of the target layer in the target well layout area to obtain the stress index of the target layer, wherein the stress value includes the average horizontal maximum principal stress, the average horizontal minimum principal stress and the average vertical stress; The stress state information of the target layer is determined according to the stress index.

4. The method for determining a well layout pattern based on shale gas development according to claim 1, characterized in that: Before the step of determining the fault data corresponding to each grid unit in the geological grid model when the fracture stability information satisfies the preset formation stability condition and the formation trend information of the target layer satisfies the preset formation trend condition, the method further includes: The formation dip angle of the target layer is obtained.

5. The method for determining a well layout pattern based on shale gas development according to claim 4, characterized in that: When the fracture stability information satisfies the preset formation stability condition, and the formation trend information of the target layer satisfies the preset formation trend condition, the step of determining the fault data corresponding to each grid unit in the geological grid model comprises: When the fracture stability index of the target layer is less than the fracture stability judgment coefficient, it is determined that the fracture stability information meets the preset formation stability condition; When the formation dip angle of the target layer is less than or equal to a preset angle, determining that the formation trend information of the target layer satisfies a preset formation trend condition; and Based on the fault distribution information, the target fault distance and the fault stress angle corresponding to each grid unit in the geological grid model are determined.

6. The method for determining a well layout pattern based on shale gas development according to claim 5, characterized in that: The step of determining the target fault distance and the fault stress angle corresponding to each grid unit in the geological grid model based on the fault distribution information includes: Obtaining the shortest distance between the center point of each grid unit in the geological grid model and the distribution position of each fault; Comparing the shortest distances to determine a target fault distance, wherein the target fault distance includes the shortest distance with the smallest value; and The angle between the minimum horizontal principal stress direction corresponding to each grid unit and the fault strike is determined, wherein the fault stress angle includes the angle between the minimum horizontal principal stress direction and the fault strike.

7. The method for determining a well layout pattern based on shale gas development according to claim 1, characterized in that: After the step of determining the well layout pattern corresponding to the target layer based on the fault data, the method further includes: A well layout decision diagram is generated based on the well layout pattern.

8. A machine-readable storage medium, characterized in that: The machine-readable storage medium stores instructions, and when the instructions are executed by the processor, the processor implements the method for determining a well layout pattern based on shale gas development according to any one of claims 1 to 7.

Images