Casing reinforcement parameter determination method, device, computer equipment and storage medium

Abstract

The present application is about a method, device, computer equipment and storage medium for determining casing reinforcement parameters, and specifically relates to the field of energy extraction. The method comprises: obtaining a three-dimensional geological model and a three-dimensional geostress field model of a target area; obtaining the deformation risk position of a target casing in the target area according to the three-dimensional geological model and the three-dimensional geostress field model; obtaining the maximum stress borne by the deformation risk position according to the historical data of the target area; and determining the casing reinforcement parameters corresponding to the deformation risk position according to the maximum stress borne by the deformation risk position. The above scheme determines the position where the casing may deform by establishing a three-dimensional geological model and a three-dimensional stress field model of the target area, and determines how to reinforce the casing at this position according to the maximum stress that the deformation risk position may bear at this position, thereby reducing the risk of deformation of the casing due to excessive stress and improving the extraction efficiency of shale oil and gas.

Description

Technical Field

[0001] The present application relates to the field of energy extraction, and in particular to a casing reinforcement method, device, computer equipment and storage medium. Background Art

[0002] With the development of oil extraction technology, the proportion of shale oil extraction in oil extraction is gradually increasing. Shale oil refers to the oil resources contained in shale formations dominated by shale. It includes oil in the pores and cracks of mud shale, as well as oil resources in the adjacent layers and interlayers of dense carbonate or clastic rocks in mud shale formations.

[0003] In the related technologies, the effective development methods are usually horizontal wells and staged fracturing technology. The staged fracturing technology of horizontal wells is that when the fracturing fluid enters the reservoir through the perforated holes at a high speed, the perforated friction resistance will be generated and increase with the increase of pumping displacement, driving the bottom hole pressure to rise. Once the bottom hole pressure exceeds the fracture pressure of multiple fracturing layers, cracks will be pressed in each layer, requiring the fracture pressure of each layer to be basically close, which can be adjusted by the perforated friction resistance.

[0004] However, in the related technology, the problem of casing deformation occurs during the fracturing process, which destroys the integrity of the well, shortens the life cycle of the well, and leads to low shale oil and gas extraction efficiency. Summary of the invention

[0005] The embodiments of the present application provide a method, device, computer equipment and storage medium for determining casing reinforcement parameters, which can reduce the risk of casing deformation due to excessive stress and improve the extraction efficiency of shale oil and gas. The technical solution is as follows:

[0006] In one aspect, a method for determining casing reinforcement parameters is provided, the method comprising:

[0007] Acquire a three-dimensional geological model and a three-dimensional geostress field model of the target area; the three-dimensional geological model is used to indicate the distribution of faults or cracks corresponding to the target area; the three-dimensional geostress field model is used to indicate the stress distribution of the target area;

[0008] According to the three-dimensional geological model and the three-dimensional geostress field model, a deformation risk position of the target casing in the target area is obtained; the deformation risk position is a position among various positions of the target casing where the probability of deformation is greater than a deformation threshold;

[0009] According to the historical data of the target area, the maximum stress borne by the deformation risk position is obtained; the historical data includes the deformation amount of the casing deformation occurring in the target area;

[0010] According to the maximum stress borne by the deformation risk position, a casing reinforcement parameter corresponding to the deformation risk position is determined, and the casing reinforcement parameter is used to indicate the casing reinforcement strength of the deformation risk position.

[0011] In another aspect, a casing reinforcement parameter determination device is provided, the device comprising:

[0012] A model acquisition module, for a three-dimensional geological model and a three-dimensional geostress field model of a target area; the three-dimensional geological model is used to indicate the distribution of faults or cracks corresponding to the target area; the three-dimensional geostress field model is used to indicate the stress distribution of the target area;

[0013] A deformation position acquisition module is used to acquire a deformation risk position of a target casing in the target area according to the three-dimensional geological model and the three-dimensional geostress field model; the deformation risk position is a position of the target casing where the probability of deformation is greater than a deformation threshold among various positions of the target casing;

[0014] A maximum stress acquisition module, used to acquire the maximum stress borne by the deformation risk position according to the historical data of the target area; the historical data includes the deformation amount of the casing deformation occurring in the target area;

[0015] The casing reinforcement parameter determination module is used to determine the casing reinforcement parameter corresponding to the deformation risk position according to the maximum stress borne by the deformation risk position, and the casing reinforcement parameter is used to indicate the casing reinforcement strength of the deformation risk position.

[0016] In a possible implementation, the deformation position acquisition module includes:

[0017] A spatial parameter acquisition unit, used to determine the spatial parameters of each fault or crack in the target area according to the three-dimensional geological model; the spatial parameters include at least one of the position information and size information of each fault or crack;

[0018] A potential activation layer acquisition unit is used to acquire a potential activated fault from each fault or crack in the target area according to the three-dimensional geostress field model and the spatial parameters of each fault or crack; the potential activated fault is a fault or crack with a non-zero slip probability;

[0019] a slip probability determination unit, configured to determine the slip probability of the potential activated fault according to the pore pressure increment of the potential activated fault; the pore pressure increment is the increase in pore pressure at the potential activated fault when the hydraulic fracturing operation is performed;

[0020] The risk position acquisition unit is used to acquire the deformation risk position of the target casing in the target area according to the slip probability of the potentially activated fault and the deformation threshold.

[0021] In a possible implementation, the risk location acquisition unit is used to:

[0022] Acquire the potential activated faults that have intersections with the target casing in the target area among the potential activated faults as candidate deformation faults;

[0023] The intersection of the candidate deformation faults whose sliding probability is greater than the deformation threshold among the candidate deformation faults and the target casing in the target area is obtained as the deformation risk position of the target casing in the target area.

[0024] In a possible implementation, the historical data also includes microseismic data;

[0025] The maximum stress acquisition module further includes:

[0026] a sliding distance acquisition unit, configured to acquire a sliding distance corresponding to the deformation amount of the sample casing according to the deformation amount of the sample casing when the sample casing in the target area is deformed and microseismic data corresponding to the deformation amount of the sample casing; the sample casing is a casing of the same type as the target casing;

[0027] The maximum stress acquisition unit is used to acquire the maximum stress borne by the deformation risk position of the target casing according to the sliding distance corresponding to the deformation amount of the sample casing.

[0028] In a possible implementation, the maximum stress acquisition unit is further used to:

[0029] A fault stress temperature-pressure coupling model of the target area is constructed according to the geological parameters and structural parameters of the target area; the geological parameters include at least one of the temperature parameters, pressure parameters and rock formation parameters of the target area; the structural parameters include at least one of the casing parameters and the parameters of the external reinforcement material of the casing; the fault stress temperature-pressure coupling model is used to indicate the relationship between the sliding distance of the fault or crack in the target area and the stress on the casing;

[0030] The maximum stress borne by the deformation risk position of the target casing is obtained according to the sliding distance and the fault stress temperature-pressure coupling model of the target area.

[0031] In a possible implementation manner, the device further includes:

[0032] A reinforcement threshold acquisition module, used to acquire a reinforcement stress threshold corresponding to the casing reinforcement parameter; the reinforcement stress threshold is used to indicate the maximum stress borne by the target casing after reinforcement according to the casing reinforcement parameter;

[0033] The reinforcement parameter updating module is used to update the casing reinforcement parameter according to the maximum stress borne by the deformation risk position when the maximum stress borne by the deformation risk position is greater than the reinforcement stress threshold.

[0034] In a possible implementation, the casing reinforcement parameter includes at least one of the performance, weaving method, and dosage thickness of the reinforced casing material.

[0035] On the other hand, a computer device is provided, which includes a processor and a memory, wherein the memory stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by the processor to implement the above-mentioned casing reinforcement parameter determination method.

[0036] On the other hand, a computer-readable storage medium is provided, wherein at least one instruction, at least one program, a code set or an instruction set is stored in the storage medium, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by a processor to implement the above-mentioned casing reinforcement parameter determination method.

[0037] In another aspect, a computer program product or a computer program is provided, the computer program product or the computer program comprising computer instructions, the computer instructions being stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, so that the computer device executes the above-mentioned casing reinforcement parameter determination method.

[0038] The technical solution provided by this application may have the following beneficial effects:

[0039] A three-dimensional geological model and a three-dimensional geostress model of the target area are constructed according to the parameter information of the target area, and the sliding probability and sliding distance of the cracks or faults in the target area are determined according to the three-dimensional geological model and the three-dimensional geostress model. The deformation risk position of the casing in the target area is determined according to the sliding probability, and the maximum stress corresponding to the deformation risk position is determined according to the sliding distance, and the parameters for reinforcing the casing at the deformation risk position are determined according to the magnitude of the maximum stress. Through the above scheme, the position where the casing may deform is determined, and at this position, according to the maximum stress that the deformation risk position may withstand, it is determined how to reinforce the casing at this position, thereby reducing the risk of deformation of the casing due to excessive stress and improving the efficiency of shale oil and gas extraction.

[0040] It should be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0042] Figure 1 The system structure diagram of a casing reinforcement parameter determination system is shown according to an exemplary embodiment.

[0043] Figure 2 It is a flow chart of a method for determining casing reinforcement parameters according to an exemplary embodiment.

[0044] Figure 3 The present invention is a flow chart of a method for determining casing reinforcement parameters according to an exemplary embodiment.

[0045] Figure 4 Shows Figure 3 The illustrated embodiment relates to a flow chart of a method for determining casing reinforcement parameters.

[0046] Figure 5 Shows Figure 3 The illustrated embodiment relates to a flow chart of a method for determining a set change risk.

[0047] Figure 6 Shows Figure 3 The illustrated embodiment relates to a flow chart of a method for calculating casing deformation force.

[0048] Figure 7 Shows Figure 3 The illustrated embodiment is a flow chart of a method for determining reinforcement parameters.

[0049] Figure 8 Shows Figure 3 A schematic diagram of casing reinforcement according to an exemplary embodiment is shown.

[0050] Fig. 9 It is a flow chart of a method for determining casing reinforcement parameters according to an exemplary embodiment.

[0051] Fig.10 It is a structural block diagram of a casing reinforcement parameter determination device according to an exemplary embodiment.

[0052] Fig.11 The diagram is a schematic structural diagram of a computer device according to an exemplary embodiment. DETAILED DESCRIPTION

[0053] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application. Instead, they are merely examples of devices and methods consistent with some aspects of the present application as detailed in the appended claims.

[0054] It should be understood that the "several" mentioned in this article refers to one or more, and "multiple" refers to two or more. "And / or" describes the association relationship of the associated objects, indicating that there can be three relationships. For example, A and / or B can mean: A exists alone, A and B exist at the same time, and B exists alone. The character " / " generally indicates that the associated objects before and after are in an "or" relationship.

[0055] Figure 1 1 is a system configuration diagram of a casing reinforcement parameter determination system according to an exemplary embodiment. The system includes: a feature information acquisition device 110 , a reinforcement parameter determination device 120 and a sensor device 130 .

[0056] The characteristic information collection device 110 can be used to collect seismic, well logging, and core data obtained by the sensor device, or the characteristic information collection device 110 can also be used to collect formation, fault, casing, and cement ring data obtained by the sensor device 130. For example, the characteristic information collection device 110 can be a terminal device with a function of collecting characteristic information, for example, the terminal can be a mobile phone, a tablet computer, smart glasses, a smart watch, a laptop portable computer, a desktop computer, etc. The number of characteristic information collection devices 110 is not limited.

[0057] The reinforcement parameter determination device 120 may be a terminal device that has a function of processing characteristic information to determine casing reinforcement parameters, for example, the terminal may be a mobile phone, a tablet computer, smart glasses, a smart watch, a laptop computer, a desktop computer, etc. The number of reinforcement parameter determination devices 120 is not limited.

[0058] Among them, the feature information collection device 110 and the reinforcement parameter determination device 120 are connected through a communication network and can transmit information.

[0059] Alternatively, information can be transmitted between the feature information collection device 110 and the reinforcement parameter determination device 120 through the Internet. For example, the feature information collection device 110 and the reinforcement parameter determination device 120 can respectively install applications with information transmission functions, and use the server corresponding to the application in the Internet as a relay to transmit feature information.

[0060] The characteristic information collection device 110 is connected to the reinforcement parameter determination device 120 via a communication network. Optionally, the communication network is a wired network or a wireless network.

[0061] In the embodiment of the present application, the characteristic information acquisition device 110 can send the characteristic information data collected 130 from the sensor device to the reinforcement parameter determination device 120, and the reinforcement parameter determination device 120 further processes the characteristic information to determine the reinforcement parameters of the casing.

[0062] Optionally, the above-mentioned wireless network or wired network uses standard communication technology and / or protocol. The network is usually the Internet, but it can also be any network, including but not limited to any combination of local area network (LAN), metropolitan area network (MAN), wide area network (WAN), mobile, wired or wireless network, private network or virtual private network. In some embodiments, the data exchanged through the network is represented by technology and / or format including Hyper Text Mark-up Language (HTML), Extensible Markup Language (XML), etc. In addition, conventional encryption technologies such as Secure Socket Layer (SSL), Transport Layer Security (TLS), Virtual Private Network (VPN), Internet Protocol Security (IPsec) can also be used to encrypt all or some links. In other embodiments, customized and / or dedicated data communication technology can also be used to replace or supplement the above-mentioned data communication technology.

[0063] Please refer to Figure 2 , which is a flow chart of a method for determining casing reinforcement parameters according to an exemplary embodiment. The method can be executed by a computer device, which can be Figure 1 The reinforcement parameter determination device 120 in the illustrated embodiment. Figure 2 As shown, the process of the method for determining casing reinforcement parameters may include the following steps:

[0064] Step 201, obtain a three-dimensional geological model and a three-dimensional geostress field model of a target area; the three-dimensional geological model is used to indicate the distribution of faults or cracks corresponding to the target area; the three-dimensional geostress field model is used to indicate the stress distribution of the target area.

[0065] In a possible implementation, the three-dimensional geological model of the target area is a three-dimensional fracture geological model constructed based on ant tracking technology. That is, based on the geological structural characteristics of the target area, the faults identified by ant body technology are taken as the main body, and the faults interpreted by microseismic signals are used as supplements to comprehensively observe the fault occurrence and establish a fault model.

[0066] The occurrence is used to indicate the dip information, inclination information and strike information corresponding to each point on the fault or fracture.

[0067] In a possible implementation, the geostress and pore pressure of the target area can be obtained based on conventional logging, imaging logging, and small-scale fracturing test data in the target area to construct a three-dimensional geostress field model of the target area.

[0068] Step 202, based on the three-dimensional geological model and the three-dimensional geostress field model, obtain the deformation risk position of the target casing in the target area; the deformation risk position is a position of the target casing where the probability of deformation is greater than the deformation threshold.

[0069] Step 203: obtaining the maximum stress borne by the deformation risk position according to the historical data of the target area; the historical data includes the deformation amount of the casing deformation occurring in the target area.

[0070] Step 204: determining a casing reinforcement parameter corresponding to the deformation risk position according to the maximum stress borne by the deformation risk position, wherein the casing reinforcement parameter is used to indicate the casing reinforcement strength of the deformation risk position.

[0071] In a possible implementation, corresponding reinforcement materials are installed at deformation risk positions of the target casing according to materials and reinforcement methods corresponding to casing reinforcement parameters to improve the ability to withstand stress.

[0072] The reinforcing material may be a high tensile strength and toughness material, such as carbon fiber, glass fiber, aramid fiber, alumina fiber and silicon carbide fiber, or a combination of fiber composite material and resin or geopolymer mortar.

[0073] In summary, in the scheme shown in the embodiment of the present application, a three-dimensional geological model and a three-dimensional geostress model of the target area are constructed according to the parameter information of the target area, and the sliding probability and sliding distance of the cracks or faults in the target area are determined according to the three-dimensional geological model and the three-dimensional geostress model, the deformation risk position of the casing in the target area is determined according to the sliding probability, the maximum stress corresponding to the deformation risk position is determined according to the sliding distance, and the parameters for reinforcing the casing at the deformation risk position are determined according to the magnitude of the maximum stress. Through the above scheme, the position where the casing may deform is determined, and at this position, according to the maximum stress that the deformation risk position may withstand, it is determined how to reinforce the casing at this position, thereby reducing the risk of deformation of the casing due to excessive stress and improving the efficiency of shale oil and gas extraction.

[0074] Please refer to Figure 3 , which is a method flow chart of a method for determining casing reinforcement parameters according to an exemplary embodiment. The method can be executed by a computer device, which can be Figure 1 The reinforcement parameter determination device 120 in the illustrated embodiment. Figure 3 As shown, the method for determining casing reinforcement parameters may include the following steps:

[0075] Step 301: Obtain a three-dimensional geological model and a three-dimensional geostress field model of a target area.

[0076] In a possible implementation, the three-dimensional geological model is a three-dimensional quantitative stochastic model constructed according to geological parameters of the target area.

[0077] The geological parameters include at least one of the seismic parameters, well logging and core parameters, and conceptual model parameters of the target area. The seismic parameters include at least one of seismic horizons, faults, seismic phases, rock types, and rock properties; the well logging and core parameters include well-connected profiles, lithology, lithofacies, and rock properties; permeability; at least one of the oil-gas-water interface, histogram, scatter plot, and vertical semi-variogram; the conceptual model parameters include at least one of the sedimentary phase model, sedimentary body superposition relationship, mudstone distribution characteristics, percentage and histogram corresponding to the size of sediments, and spatial continuity lateral semi-variogram.

[0078] The three-dimensional geological model can be a three-dimensional quantitative random model generated based on the above geological parameters, that is, based on geological, well logging, geophysical data and various interpretation results or conceptual models.

[0079] In a possible implementation, the three-dimensional geostress field model is used to indicate the initial stress distribution of the target area.

[0080] Among them, the initial stress values ​​at various positions in the three-dimensional geostress field model can be calculated according to the displacement back analysis method, that is, a mechanical model is constructed according to the geometric structure of the rock mass, the mechanical parameters of the rock mass, the structural size parameters of the rock mass, etc., and the various parameters in the mechanical model are determined according to the measured values ​​of the rock and soil displacement caused by the microseismicity caused by the hydraulic fracturing operation. Finally, the mechanical model is analyzed and calculated to obtain the three-dimensional geostress field model of the target area.

[0081] Alternatively, the initial stress values ​​at each position in the three-dimensional geostress field model can be obtained by a stress function method, that is, a uniform stress field is assumed in the target area, and the three-dimensional geostress field model is obtained according to the stress function by making the calculated stress value match the measured value at a given point.

[0082] In a possible implementation, the three-dimensional geological model and the three-dimensional geostress field model have the same three-dimensional coordinate system, so each point in the target area has the same three-dimensional coordinates in the three-dimensional geological model and the three-dimensional geostress field model.

[0083] Step 302: Determine the spatial parameters of each fault or fracture in the target area according to the three-dimensional geological model.

[0084] The spatial parameter includes at least one of the position information and size information of each fault or crack.

[0085] In a possible implementation, the three-dimensional coordinates corresponding to each point on each fault or fracture in the target area and the occurrence information corresponding to each point are determined based on the three-dimensional geological model.

[0086] In a possible implementation, the spatial parameters of the fault or crack may be identified by using ant body technology.

[0087] In a possible implementation, the spatial parameters of the fault or crack may be acquired based on microseismic data.

[0088] For example, the spatial parameters of faults or cracks can be obtained through microseismic data based on coherence volume technology. When geological factors such as stratum lithology and characteristics change laterally, seismic waves will inevitably change, which will further cause changes in the value of seismic waves. Therefore, the development of cracks and faults can be intuitively displayed through the slice diagram of seismic coherence volume. Usually, the longer linear or large curvature radius curves on the slice diagram are fault displays, and the shorter ones are crack displays.

[0089] Step 303: Acquire potential activated faults in each fault or fracture in the target area according to the three-dimensional geostress field model and the spatial parameters of each fault or fracture.

[0090] The potentially activated fault is a fault or crack whose slip probability is not zero.

[0091] In a possible implementation, initial stress parameters corresponding to each point on each fault or crack are obtained according to the three-dimensional geostress field.

[0092] The initial stress parameters may include at least one of horizontal maximum principal stress, horizontal minimum principal stress, vertical stress, pore pressure and horizontal maximum principal stress orientation.

[0093] In a possible implementation, a potential activated fault is obtained from each fault or fracture in the target area according to the vector sum of the initial stress corresponding to each point on the each fault or fracture.

[0094] When the vector sum of the initial stresses corresponding to each point on a fault or crack is greater than a threshold, it means that the stress received by the fault or crack is close to the stress threshold at which the fault or crack slips. At this time, hydraulic fracturing may increase the pore pressure, or micro-earthquakes caused by other geological activities may cause the stress received by the fault or crack to be greater than the stress threshold at which the fault or crack slips. Therefore, when the vector sum of the initial stresses corresponding to each point on the fault or crack is greater than the stress threshold corresponding to the fault or crack, the fault or crack is obtained as a potentially activated fault.

[0095] Among them, the stress threshold corresponding to the fault or crack (or called critical stress fault) can be determined according to the spatial parameters corresponding to the fault or crack. For example, the smaller the size of the fault or crack, the greater the stress required to destroy the original geological structure to enable the fault to slide, and therefore the greater the stress threshold corresponding to the fault or crack.

[0096] In a possible implementation, the three-dimensional coordinates, initial stress parameters and occurrence information of each point on each fault or fracture are obtained according to the spatial parameters of the fault or fracture, and the potential activity index of each point on the fault or fracture is determined. The potential activity index is used to indicate the potential mechanical activity of each point on the fault or fracture.

[0097] In a possible implementation, based on whether the potential activity index of each point on the fault or fracture satisfies a specified condition, the fault or fracture is obtained as a potential activated layer.

[0098] For example, when the sum of the potential activity indexes of all points on the fault or crack is greater than the activity threshold, or the number of points on the fault or crack whose potential activity index is greater than the activity point threshold is greater than the proportion threshold, it means that the overall potential activity of the fault or crack is high, and the fault or crack may be active (i.e., fault slip), so the fault or crack is obtained as a potential activated layer.

[0099] Step 304: Determine the slip probability of the potential activated fault according to the pore pressure increment of the potential activated fault.

[0100] The pore pressure increment is the increase in pore pressure at the potentially activated fault when the hydraulic fracturing operation is performed.

[0101] In a possible implementation, the pore pressure increment is generated after hydraulic fracturing is performed on the target area.

[0102] Among them, pore pressure is the pressure of the fluid in the rock pores in the target area. The pore pressure at a certain location in the target area is related to the depth of the location. However, when hydraulic fracturing operations are performed on the target area, the impermeable layer in the target area blocks the free flow of pore fluid in the aquifer in the target area. The pressure in the pores of the rock will increase, that is, a certain pore pressure increment will be generated.

[0103] In a possible implementation, a sliding probability function is constructed according to the pore pressure increment and the sliding probability of the potential activated layer; the sliding probability function is used to indicate the relationship between the sliding probability of the activated layer and the pore pressure increment.

[0104] In one possible implementation, the slip probability model is constructed according to the fault data of the target area and the three-dimensional geostress field model of the target area, based on the fault slip hypothesis and the quantitative risk analysis (QRA) method. The slip probability model is a function of the pore pressure increment and the slip probability of the potential activation layer.

[0105] First, the three-dimensional geostress field model is obtained. According to the uncertainty of geological parameters in the three-dimensional geostress field model, the uniformly distributed geological parameters are randomly sampled according to the uncertainty to generate a geomechanical probability model. Then, according to the distribution of the disturbance value of the pore pressure (critical pore pressure) required for fault sliding in the probability model, the function of the sliding probability increasing with the pore pressure is obtained, that is, the above-mentioned sliding probability model.

[0106] In a possible implementation, the sliding probability of the potential activated layer in the target area is obtained according to the sliding probability function and the pore pressure increment corresponding to hydraulic fracturing in the target area.

[0107] The pore pressure increment may be determined based on the pore pressure coefficient of the target area and the stress generated by the hydraulic fracturing operation on the target area. The pore pressure coefficient is used to indicate the ratio of the pore pressure increment of the target area to the deviatoric stress on the target area.

[0108] Step 305: acquiring a deformation risk position of a target casing in the target area according to the slip probability of the potentially activated fault and the deformation threshold.

[0109] The deformation threshold is used to indicate the maximum value of the sliding probability of the potential activated fault allowed in the target area, that is, when the sliding probability of the potential activated fault is greater than the deformation threshold, the potential activated fault can be regarded as an unsafe fault or crack that may slide, and thus the position in the target casing in the target area is determined as a deformation risk position with deformation risk.

[0110] In a possible implementation, among the potential activated faults, the potential activated faults that have an intersection with the target casing in the target area are obtained as candidate deformation faults; and among the candidate deformation faults, the intersection of the candidate deformation faults having a sliding probability greater than a deformation threshold and the target casing in the target area is obtained as the deformation risk position of the target casing in the target area.

[0111] In a possible implementation, the intersection of the target casing in the target area and the potentially activated fault is obtained as a candidate deformation point; and among the candidate deformation faults corresponding to each candidate deformation point, the candidate deformation fault whose sliding probability is greater than the deformation threshold is obtained as the target deformation fault.

[0112] That is, among the potentially activated faults, what needs to be considered are the cracks and faults intersecting with the target casing in the target area. When the cracks or faults intersecting with the target casing in the target area slip, there is a high probability that the stress will be directly transmitted to the position where the target casing intersects with the cracks or faults, thereby causing casing deformation. The cracks or faults that do not intersect with the target casing slip, usually with a small sliding distance, and their stress cannot be transmitted to the target casing in the target area, or the stress transmitted to the target casing in the target area is small and cannot cause casing deformation.

[0113] Therefore, in the potentially activated fault, the intersection position corresponding to the target deformation fault that intersects with the target casing and has a sliding probability greater than the deformation threshold is obtained as the deformation risk position.

[0114] In a possible implementation, the deformation threshold may include a first deformation threshold and a second deformation threshold, wherein the first deformation threshold may be greater than the second deformation threshold.

[0115] The intersection of the candidate deformation fault whose sliding probability is greater than the second deformation threshold among the candidate deformation faults and the target casing in the target area is obtained as the deformation medium risk position of the target casing in the target area. The intersection of the candidate deformation fault whose sliding probability is greater than the first deformation threshold among the candidate deformation faults and the target casing in the target area is obtained as the deformation high risk position of the target casing in the target area.

[0116] In a possible implementation, a potential activated fault having a point whose distance from a target casing in a target area is less than a sliding threshold and whose sliding probability is greater than a deformation threshold is obtained as a target deformation fault. The sliding threshold is used to indicate the maximum value of the allowable sliding distance of the potential activated fault.

[0117] When there is a point in the potentially activated fault whose distance to the target casing in the target area is less than the sliding threshold, that is, the target casing is close to the potentially activated fault, although the target casing does not intersect with the potentially activated fault, the gravitational force generated when the potentially activated fault slides can still be transmitted to the target casing, so it may still cause deformation of the target pipeline.

[0118] In a possible implementation, the sliding threshold may be determined based on geological parameters of the target area and size information of the target deformation fault.

[0119] In a possible implementation, the deformation risk position of the target pipeline is determined according to the point in the target deformation fault that is the shortest distance from the target casing and the size information of the target deformation fault.

[0120] For example, the point with the smallest distance to the target casing is taken as the center of the circle, the radius is determined according to the size information of the target deformation fault, the deformation risk circle is obtained, and all target pipelines within the deformation risk circle are determined as the deformation risk positions of the target pipelines.

[0121] Step 306, according to the deformation amount of the sample casing in the target area when the sample casing is deformed and the microseismic data corresponding to the deformation amount of the sample casing, obtain the sliding distance corresponding to the deformation amount of the sample casing; the sample casing is a casing of the same type as the target casing.

[0122] Among them, historical data also includes microseismic data.

[0123] In a possible implementation, the sample casing may be a casing in the target area whose distance from the target casing is less than a geological threshold.

[0124] That is, the sliding conditions of the potential activated faults in the target area can be obtained through the deformation conditions of the sample casings that are closer to the target casing in the area corresponding to the target casing, and the microseismic data corresponding to the deformation conditions. Since the distance between the target casing and the sample casing is less than the geological threshold, and the geological properties in a certain area are generally the same, the sliding probability and sliding distance of the activated faults corresponding to the target casing in the target area can be simulated through the historical data of the sample casing, and the corresponding relationship between the sliding distance of the faults that slip in each potential activated fault and the microseismic data can be determined.

[0125] Step 307: Obtain the maximum stress borne by the deformation risk position of the target casing according to the sliding distance corresponding to the deformation amount of the sample casing.

[0126] In a possible implementation, the maximum stress among the stresses corresponding to the sample casing deformation is obtained according to the sliding distance corresponding to the sample casing deformation, and the maximum stress corresponding to the sample casing deformation is obtained as the maximum stress borne by the deformation risk position of the target casing.

[0127] In a possible implementation, a fault stress temperature-pressure coupling model of the target area is constructed according to the geological parameters and structural parameters of the target area; the geological parameters include at least one of the temperature parameters, pressure parameters and rock formation parameters of the target area; the structural parameters include at least one of the casing parameters and the parameters of the external reinforcement material of the casing; the fault stress temperature-pressure coupling model is used to indicate the relationship between the sliding distance of the fault or crack in the target area and the stress on the casing.

[0128] According to the sliding distance corresponding to the deformation amount of the sample casing and the fault stress temperature-pressure coupling model of the target area, the maximum stress borne by the deformation risk position of the target casing is obtained.

[0129] In a possible implementation, the fault stress temperature-pressure coupling model may be a casing-cement sheath-formation combination temperature-pressure coupling finite element model determined according to the temperature parameters, pressure parameters, geological parameters, casing parameters, and cement sheath parameters corresponding to the casing in the target area.

[0130] The temperature-pressure coupling finite element model of the casing-cement sheath-formation combination is established according to the parameter information of the casing, cement sheath and formation, and based on the fluid-solid thermal coupling effect. The model can obtain the stress value of the casing caused by the fault sliding distance displacement according to the fault sliding distance displacement, the inner diameter of the casing, the wall thickness, the cement sheath parameters, and the temperature and pressure of the target area.

[0131] In a possible implementation, the sliding distance corresponding to the sample casing deformation amount is used to obtain the maximum stress borne by the target casing deformation risk position, and the sliding distance corresponding to the sample casing deformation amount and the corresponding maximum stress are used to obtain the corresponding relationship between the sliding distance and the maximum stress; based on the corresponding relationship and the estimated sliding distance of the crack or fault corresponding to the deformation risk position, the maximum stress corresponding to the deformation risk position is obtained.

[0132] In a possible implementation, the spatial parameters of the fracture or fault are substituted into a temperature-pressure coupled finite element model of the casing-cement sheath-formation assembly to obtain an estimated sliding distance of the fracture or fault.

[0133] Step 308: Determine the casing reinforcement parameter corresponding to the deformation risk position according to the maximum stress borne by the deformation risk position.

[0134] In a possible implementation, the casing reinforcement parameter includes at least one of the performance, weaving method, and dosage thickness of the reinforced casing material.

[0135] In a possible implementation, the casing reinforcement parameter is also used to indicate the size of the casing reinforcement area.

[0136] For example, when the casing is reinforced according to the casing reinforcement parameter, the larger the maximum stress corresponding to the casing reinforcement parameter, the more likely the maximum stress is to affect the area around the deformation risk area. Therefore, increasing the reinforcement area corresponding to the casing reinforcement parameter can improve the safety of casing reinforcement.

[0137] In a possible implementation, a reinforcement stress threshold corresponding to the casing reinforcement parameter is obtained; the reinforcement stress threshold is used to indicate the maximum stress borne by the target casing after reinforcement according to the casing reinforcement parameter; when the maximum stress borne by the deformation risk position is greater than the reinforcement stress threshold, the casing reinforcement parameter is updated according to the maximum stress borne by the deformation risk position.

[0138] After the casing is reinforced according to the casing reinforcement parameters, it can be determined whether the reinforced casing can withstand the maximum stress and not deform according to the casing reinforcement parameters and the parameters corresponding to the casing and the above-mentioned fault stress temperature-pressure coupling model. When the reinforced casing cannot withstand the maximum stress and causes deformation, the casing reinforcement parameters can be updated according to the difference between the maximum stress and the maximum stress borne by the reinforced casing when it is deformed, so as to ensure that the casing reinforced according to the casing reinforcement parameters can withstand the maximum stress in the area.

[0139] Please refer to Figure 4 , which shows a flow chart of a method for determining casing reinforcement parameters involved in an embodiment of the present application. Figure 4As shown, S401, using the three-dimensional geological model and three-dimensional geostress field of the target work area, combined with the geomechanical conditions of fracture / fault sliding and the impact of fracturing construction, the fault sliding risk is evaluated to obtain the casing deformation risk degree of the intersection of the fault / crack and the well trajectory.

[0140] Please refer to Figure 5 , which shows a flow chart of a method for determining a change risk involved in an embodiment of the present application. Figure 5 As shown:

[0141] S501, using the geological model of the target work area, identifying reservoir-related faults / fractures and their distribution in the target work area.

[0142] S502, using the three-dimensional geostress field of the target work area and the three-dimensional Moore circle diagram to determine the stress state of the fault and preliminarily screen the potential activated faults.

[0143] S503: Based on the geomechanical parameters, a function of fault slip probability and pore pressure is constructed, and the fault slip risk is evaluated by combining the pore pressure increment caused by hydraulic fracturing in the block.

[0144] S504, based on the high, medium and low risks of fault sliding, infer the risk of casing deformation at the intersection of the fault and the well trajectory.

[0145] like Figure 4 As shown, S402, using the temperature-pressure coupling finite element model of the casing-cement sheath-formation combination of the target work area, according to the fault slip distance of the work area, the stress value acting on the casing under the slip distance is calculated to obtain the maximum force causing the casing deformation.

[0146] Please refer to Figure 6 , which shows a flow chart of a method for calculating casing deformation force involved in an embodiment of the present application. Figure 6 As shown:

[0147] S601: Using the formation, fault, casing and cement sheath data of the target work area, a temperature-pressure coupling finite element model of the casing-cement sheath-formation combination after cementing is established.

[0148] S602 uses the on-site microseismic data and casing deformation point data to invert the fault sliding distance when casing deformation occurs in the well section.

[0149] S603, according to the elastic-plastic theory, the stress value acting on the casing at the corresponding sliding distance is calculated from the interface stress and displacement continuity conditions.

[0150] S604, compare and obtain the maximum force of casing deformation caused by fault sliding in the work area.

[0151] like Figure 4As shown, S403, high tensile strength and toughness material is installed outside the conventional casing to verify the ability of the casing with high tensile strength and toughness material to resist external forces. The casing will be able to resist the sliding force of the fault in the work area and is used in high and medium risk casing change well sections.

[0152] Please refer to Figure 7 , which shows a flow chart of a reinforcement parameter determination method involved in an embodiment of the present application. Figure 7 As shown, it describes the steps of combining casing that can resist the sliding force of faults in the work area according to an exemplary embodiment of the present application and applying it in high- and medium-risk casing change well sections.

[0153] S701, install high tensile strength and toughness material outside the conventional casing to verify whether the casing installed with high tensile strength and toughness material has the ability to resist external forces.

[0154] S702, using experimental and numerical simulation methods, verify whether the externally installed casing made of tensile strength and ductile material can resist the maximum force exerted on the casing by fault sliding in the work area.

[0155] S703, if the casing made of high tensile strength and toughness material does not have the ability to resist external forces, adjust the material properties, weaving method, dosage thickness and other parameters.

[0156] S704: If the casing made of high tensile strength and toughness material has the ability to resist external forces, the reinforcement scheme corresponding to the reinforcement parameters is applied to the high and medium risk casing change well sections.

[0157] Please refer to Figure 8 , which shows a schematic diagram of casing reinforcement involved in an exemplary embodiment of the present application. The externally installed high-tensile strength and toughness material casing is used in high- and medium-risk casing change well sections, and the application length of a single casing change point is 10-40m. The above-mentioned high-tensile strength and toughness material can be carbon fiber, glass fiber, aramid fiber, alumina fiber and silicon carbide fiber, or a combination of fiber composite materials and resin, geopolymer mortar. There is no restriction on the method of obtaining the high-tensile strength and toughness material in this application.

[0158] The high tensile strength toughness material can be fixed to the outside of the conventional casing in various ways such as winding, wrapping, and inlaying by conventional or unconventional weaving methods.

[0159] like Figure 8As shown, in the fracture development zones 804 and 805 in the shale formation 803, the high-risk sliding fault 804 is activated by hydraulic fracturing, and a certain amount of displacement occurs. When the sliding generates a local load such as a "point" or "line" acting on the oil layer casing 802, the high tensile strength toughness material 806 extends and disperses the force on the oil layer casing 802, so that the casing 802 shares the external force in the form of a "surface" or "body" and converts it into a regional load. The surface force area of ​​the casing 802 is increased, the force per unit area is reduced, and the casing 802 is prevented from being damaged or the deformation of the casing 802 is controlled within a range that does not affect production needs.

[0160] In summary, in the scheme shown in the embodiment of the present application, a three-dimensional geological model and a three-dimensional geostress model of the target area are constructed according to the parameter information of the target area, and the sliding probability and sliding distance of the cracks or faults in the target area are determined according to the three-dimensional geological model and the three-dimensional geostress model, the deformation risk position of the casing in the target area is determined according to the sliding probability, the maximum stress corresponding to the deformation risk position is determined according to the sliding distance, and the parameters for reinforcing the casing at the deformation risk position are determined according to the magnitude of the maximum stress. Through the above scheme, the position where the casing may deform is determined, and at this position, according to the maximum stress that the deformation risk position may withstand, it is determined how to reinforce the casing at this position, thereby reducing the risk of deformation of the casing due to excessive stress and improving the efficiency of shale oil and gas extraction.

[0161] Please refer to Fig. 9 , which is a flow chart of a casing reinforcement parameter determination method according to an exemplary embodiment. The casing reinforcement parameter determination method is jointly executed by a feature information acquisition device 900 and a reinforcement parameter determination device 910, wherein the feature information acquisition device 900 can receive corresponding feature information from a sensor device, and the reinforcement parameter determination device 910 can be a terminal or a server. Fig. 9 As shown, the reinforcement parameter determination device 910 sends an information acquisition request to the feature information acquisition device 900, obtains the feature information on the feature information acquisition device, such as geological parameters, etc., and constructs a three-dimensional geological model 911 and a three-dimensional ground stress model 912 corresponding to the target area based on the feature information.

[0162] The reinforcement parameter determination device 910 obtains the spatial parameters of the faults or cracks 913 corresponding to the target area according to the three-dimensional geological model 912 corresponding to the target area, determines the size information and position information of the faults or cracks 913 in the target area, and screens out potential activated faults 914 that may slip among the faults or cracks in the target area according to the spatial information of the faults or cracks and the three-dimensional geostress model 912 of the target area.

[0163] The reinforcement parameter determination device 910 then simulates the hydraulic fracturing operation on the target area, obtains the pore pressure increment 915 corresponding to the target area, and determines the sliding probability 916 of each potential activated fault under the pore pressure increment based on the potential activated fault and the pore pressure increment, and determines the deformation risk position 917 of the target casing at each position where casing deformation may occur based on the sliding probability 916.

[0164] The reinforcement parameter determination device 910 obtains the historical deformation data 918 of the target area from the characteristic information acquisition device 900. The historical deformation data 918 at least includes the sample casing deformation variable and microseismic data. According to the relationship between the sample casing deformation variable and the microseismic data, the sliding distance 919 of the corresponding crack or fault when the sample casing is deformed is determined. According to the sliding distance 919, according to the elastic-plastic theory, the stress value acting on the sample casing at the corresponding sliding distance is calculated by the interface stress and displacement continuity conditions, and the largest of the stress values ​​is used as the maximum stress 920 of the target casing.

[0165] And according to the maximum stress 920, the casing reinforcement parameter 921 when performing the casing reinforcement operation at the deformation risk position 917 is determined, and a simulation test is performed on the reinforced casing to determine whether it can withstand the maximum stress 920. When it cannot withstand the maximum stress, the casing reinforcement parameter 921 is updated based on the maximum stress 920.

[0166] Fig.10 is a structural block diagram of a casing reinforcement parameter determination device according to an exemplary embodiment. The casing reinforcement parameter determination device can be implemented Figure 2 or Figure 3 All or part of the steps in the method provided in the illustrated embodiment. The casing reinforcement parameter determination device may include:

[0167] Model acquisition module 1001, for a three-dimensional geological model and a three-dimensional geostress field model of a target area; the three-dimensional geological model is used to indicate the distribution of faults or cracks corresponding to the target area; the three-dimensional geostress field model is used to indicate the stress distribution of the target area;

[0168] The deformation position acquisition module 1002 is used to acquire the deformation risk position of the target casing in the target area according to the three-dimensional geological model and the three-dimensional geostress field model; the deformation risk position is a position of the target casing where the probability of deformation is greater than a deformation threshold among various positions of the target casing;

[0169] The maximum stress acquisition module 1003 is used to acquire the maximum stress borne by the deformation risk position according to the historical data of the target area; the historical data includes the deformation amount of the casing deformation in the target area;

[0170] The casing reinforcement parameter determination module 1004 is used to determine the casing reinforcement parameter corresponding to the deformation risk position according to the maximum stress borne by the deformation risk position, and the casing reinforcement parameter is used to indicate the casing reinforcement strength of the deformation risk position.

[0171] In a possible implementation, the deformation position acquisition module 1002 includes:

[0172] A spatial parameter acquisition unit, used to determine the spatial parameters of each fault or crack in the target area according to the three-dimensional geological model; the spatial parameters include at least one of the position information and size information of each fault or crack;

[0173] A potential activation layer acquisition unit is used to acquire a potential activated fault from each fault or crack in the target area according to the three-dimensional geostress field model and the spatial parameters of each fault or crack; the potential activated fault is a fault or crack with a non-zero slip probability;

[0174] a slip probability determination unit, configured to determine the slip probability of the potential activated fault according to the pore pressure increment of the potential activated fault; the pore pressure increment is the increase in pore pressure at the potential activated fault when the hydraulic fracturing operation is performed;

[0175] The risk position acquisition unit is used to acquire the deformation risk position of the target casing in the target area according to the slip probability of the potentially activated fault and the deformation threshold.

[0176] In a possible implementation, the risk location acquisition unit is used to:

[0177] Acquire the potential activated faults that have intersections with the target casing in the target area among the potential activated faults as candidate deformation faults;

[0178] The intersection of the candidate deformation faults whose sliding probability is greater than the deformation threshold among the candidate deformation faults and the target casing in the target area is obtained as the deformation risk position of the target casing in the target area.

[0179] In a possible implementation, the historical data also includes microseismic data;

[0180] The maximum stress acquisition module 1003 further includes:

[0181] a sliding distance acquisition unit, configured to acquire a sliding distance corresponding to the deformation amount of the sample casing according to the deformation amount of the sample casing when the sample casing in the target area is deformed and microseismic data corresponding to the deformation amount of the sample casing; the sample casing is a casing of the same type as the target casing;

[0182] The maximum stress acquisition unit is used to acquire the maximum stress borne by the deformation risk position of the target casing according to the sliding distance corresponding to the deformation amount of the sample casing.

[0183] In a possible implementation, the maximum stress acquisition unit is further used to:

[0184] A fault stress temperature-pressure coupling model of the target area is constructed according to the geological parameters and structural parameters of the target area; the geological parameters include at least one of the temperature parameters, pressure parameters and rock formation parameters of the target area; the structural parameters include at least one of the casing parameters and the parameters of the external reinforcement material of the casing; the fault stress temperature-pressure coupling model is used to indicate the relationship between the sliding distance of the fault or crack in the target area and the stress on the casing;

[0185] The maximum stress borne by the deformation risk position of the target casing is obtained according to the sliding distance and the fault stress temperature-pressure coupling model of the target area.

[0186] In a possible implementation manner, the device further includes:

[0187] A reinforcement threshold acquisition module, used to acquire a reinforcement stress threshold corresponding to the casing reinforcement parameter; the reinforcement stress threshold is used to indicate the maximum stress borne by the target casing after reinforcement according to the casing reinforcement parameter;

[0188] The reinforcement parameter updating module is used to update the casing reinforcement parameter according to the maximum stress borne by the deformation risk position when the maximum stress borne by the deformation risk position is greater than the reinforcement stress threshold.

[0189] In a possible implementation, the casing reinforcement parameter includes at least one of the performance, weaving method, and dosage thickness of the reinforced casing material.

[0190] In summary, in the scheme shown in the embodiment of the present application, a three-dimensional geological model and a three-dimensional geostress model of the target area are constructed according to the parameter information of the target area, and the sliding probability and sliding distance of the cracks or faults in the target area are determined according to the three-dimensional geological model and the three-dimensional geostress model, the deformation risk position of the casing in the target area is determined according to the sliding probability, the maximum stress corresponding to the deformation risk position is determined according to the sliding distance, and the parameters for reinforcing the casing at the deformation risk position are determined according to the magnitude of the maximum stress. Through the above scheme, the position where the casing may deform is determined, and at this position, according to the maximum stress that the deformation risk position may withstand, it is determined how to reinforce the casing at this position, thereby reducing the risk of deformation of the casing due to excessive stress and improving the efficiency of shale oil and gas extraction.

[0191] Fig.11 1 is a schematic diagram of the structure of a computer device according to an exemplary embodiment. The computer device can be implemented as a model search device and / or an image segmentation device in each of the above-mentioned method embodiments. The computer device 1100 includes a central processing unit (CPU) 1101, a system memory 1104 including a random access memory (RAM) 1102 and a read-only memory (ROM) 1103, and a system bus 1105 connecting the system memory 1104 and the central processing unit 1101. The computer device 1100 also includes a basic input / output system 1106 that helps transmit information between various devices in the computer, and a large-capacity storage device 1107 for storing an operating system 1113, an application program 1114 and other program modules 1115.

[0192] The mass storage device 1107 is connected to the central processing unit 1101 through a mass storage controller (not shown) connected to the system bus 1105. The mass storage device 1107 and its associated computer readable medium provide non-volatile storage for the computer device 1100. That is, the mass storage device 1107 may include a computer readable medium (not shown) such as a hard disk or a Compact Disc Read-Only Memory (CD-ROM) drive.

[0193] Without loss of generality, the computer-readable medium may include computer storage media and communication media. Computer storage media include volatile and non-volatile, removable and non-removable media implemented by any method or technology for storing information such as computer-readable instructions, data structures, program modules or other data. Computer storage media include RAM, ROM, flash memory or other solid-state storage technologies, CD-ROM, or other optical storage, cassettes, magnetic tapes, disk storage or other magnetic storage devices. Of course, those skilled in the art will appreciate that the computer storage media is not limited to the above. The above-mentioned system memory 1104 and mass storage device 1107 can be collectively referred to as memory.

[0194] The computer device 1100 can be connected to the Internet or other network devices through a network interface unit 1111 connected to the system bus 1105 .

[0195] The memory also includes one or more programs, which are stored in the memory. The CPU 1101 implements the Figure 3 , Figure 4 or Fig. 9 All or part of the steps of the method shown.

[0196] In an exemplary embodiment, a non-temporary computer-readable storage medium including instructions is also provided, such as a memory including a computer program (instructions), and the above program (instructions) can be executed by a processor of a computer device to complete the method shown in each embodiment of the present application, the method performed by a server or a user terminal. For example, the non-temporary computer-readable storage medium can be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, and an optical data storage device, etc.

[0197] In an exemplary embodiment, a computer program product or a computer program is also provided, the computer program product or the computer program includes computer instructions, the computer instructions are stored in a computer-readable storage medium. The processor of the computer device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, so that the computer device executes the methods shown in the above embodiments.

[0198] Those skilled in the art will readily appreciate other embodiments of the present application after considering the specification and practicing the invention disclosed herein. The present application is intended to cover any modification, use or adaptation of the present application, which follows the general principles of the present application and includes common knowledge or customary techniques in the art that are not disclosed in the present application. The specification and examples are intended to be exemplary only, and the true scope and spirit of the present application are indicated by the claims.

[0199] It should be understood that the present application is not limited to the precise structures that have been described above and shown in the drawings, and that various modifications and changes may be made without departing from the scope thereof. The scope of the present application is limited only by the appended claims.

Claims

1. A method for determining casing reinforcement parameters, characterized in that: The method comprises: Acquire a three-dimensional geological model and a three-dimensional geostress field model of the target area; the three-dimensional geological model is used to indicate the distribution of faults or cracks corresponding to the target area; the three-dimensional geostress field model is used to indicate the stress distribution of the target area; According to the three-dimensional geological model and the three-dimensional geostress field model, a deformation risk position of the target casing in the target area is obtained; the deformation risk position is a position among various positions of the target casing where the probability of deformation is greater than a deformation threshold; According to the historical data of the target area, the maximum stress borne by the deformation risk position is obtained; the historical data includes the deformation amount of the casing deformation occurring in the target area; Determine, according to the maximum stress borne by the deformation risk position, a casing reinforcement parameter corresponding to the deformation risk position, wherein the casing reinforcement parameter is used to indicate the casing reinforcement strength at the deformation risk position; Wherein, obtaining the deformation risk position of the target casing in the target area according to the three-dimensional geological model and the three-dimensional geostress field model includes: Determine the spatial parameters of each fault or fracture in the target area according to the three-dimensional geological model; the spatial parameters include at least one of the position information and size information of each fault or fracture; According to the three-dimensional geostress field model and the spatial parameters of each fault or fracture, a potential activated fault is obtained from each fault or fracture in the target area; the potential activated fault is a fault or fracture with a non-zero slip probability; Determining the slip probability of the potential activated fault according to the pore pressure increment of the potential activated fault; the pore pressure increment is the increase in pore pressure at the potential activated fault when the hydraulic fracturing operation is performed; Acquire the potential activated faults that have intersections with the target casing in the target area among the potential activated faults as candidate deformation faults; The intersection of the candidate deformation faults whose sliding probability is greater than the deformation threshold among the candidate deformation faults and the target casing in the target area is obtained as the deformation risk position of the target casing in the target area.

2. The method according to claim 1, characterized in that The historical data also includes microseismic data; The maximum stress borne by the deformation risk position of the target casing is obtained according to the historical data of the target area; the historical data includes the deformation amount of the casing deformation in the target area, including: According to the deformation amount of the sample casing in the target area when the sample casing is deformed, and the microseismic data corresponding to the deformation amount of the sample casing, the sliding distance corresponding to the deformation amount of the sample casing is obtained; the sample casing is a casing of the same type as the target casing; According to the sliding distance corresponding to the deformation amount of the sample casing, the maximum stress borne by the deformation risk position of the target casing is obtained.

3. The method according to claim 2, characterized in that The step of obtaining the maximum stress borne by the deformation risk position of the target casing according to the sliding distance includes: A fault stress temperature-pressure coupling model of the target area is constructed according to the geological parameters and structural parameters of the target area; the geological parameters include at least one of the temperature parameters, pressure parameters and rock formation parameters of the target area; the structural parameters include at least one of the casing parameters and the parameters of the external reinforcement material of the casing; the fault stress temperature-pressure coupling model is used to indicate the relationship between the sliding distance of the fault or crack in the target area and the stress on the casing; The maximum stress borne by the deformation risk position of the target casing is obtained according to the sliding distance and the fault stress temperature-pressure coupling model of the target area.

4. The method according to claim 1, characterized in that: The method further comprises: Obtaining a reinforcement stress threshold corresponding to the casing reinforcement parameter; the reinforcement stress threshold is used to indicate the maximum stress borne by the target casing after being reinforced according to the casing reinforcement parameter; When the maximum stress borne by the deformation risk position is greater than the reinforcement stress threshold, the casing reinforcement parameter is updated according to the maximum stress borne by the deformation risk position.

5. The method according to claim 1, characterized in that The casing reinforcement parameters include at least one of the performance, weaving method, dosage and thickness of the reinforced casing material.

6. A casing reinforcement parameter determination device, characterized in that: The device comprises: A model acquisition module, for a three-dimensional geological model and a three-dimensional geostress field model of a target area; the three-dimensional geological model is used to indicate the distribution of faults or cracks corresponding to the target area; the three-dimensional geostress field model is used to indicate the stress distribution of the target area; A deformation position acquisition module is used to acquire a deformation risk position of a target casing in the target area according to the three-dimensional geological model and the three-dimensional geostress field model; the deformation risk position is a position of the target casing where the probability of deformation is greater than a deformation threshold among various positions of the target casing; A maximum stress acquisition module, used to acquire the maximum stress borne by the deformation risk position according to the historical data of the target area; the historical data includes the deformation amount of the casing deformation occurring in the target area; A casing reinforcement parameter determination module, used to determine the casing reinforcement parameter corresponding to the deformation risk position according to the maximum stress borne by the deformation risk position, wherein the casing reinforcement parameter is used to indicate the casing reinforcement strength of the deformation risk position; Wherein, the deformation position acquisition module includes: A spatial parameter acquisition unit, used to determine the spatial parameters of each fault or crack in the target area according to the three-dimensional geological model; the spatial parameters include at least one of the position information and size information of each fault or crack; A potential activation layer acquisition unit is used to acquire a potential activated fault from each fault or crack in the target area according to the three-dimensional geostress field model and the spatial parameters of each fault or crack; the potential activated fault is a fault or crack with a non-zero slip probability; a slip probability determination unit, configured to determine the slip probability of the potential activated fault according to the pore pressure increment of the potential activated fault; the pore pressure increment is the increase in pore pressure at the potential activated fault when the hydraulic fracturing operation is performed; The risk position acquisition unit is used to acquire the potential activated faults among the potential activated faults that have an intersection with the target casing in the target area as candidate deformation faults; and acquire the intersection of the candidate deformation faults whose sliding probability is greater than the deformation threshold among the candidate deformation faults and the target casing in the target area as the deformation risk position of the target casing in the target area.

7. A computer device, characterized in that: The computer device includes a processor and a memory, wherein the memory stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by the processor to implement the casing reinforcement parameter determination method according to any one of claims 1 to 5.

8. A computer-readable storage medium, characterized in that: The storage medium stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, the at least one program, the code set or the instruction set is loaded and executed by the processor to implement the casing reinforcement parameter determination method according to any one of claims 1 to 5.

Images