Experimental Study on Formation Slip under Injection- Production Interregional Pressure Difference Based on the Abnormal Similarity Theory

Laboratory of Enhanced Oil Recovery of Education Ministry, Northeast Petroleum University, Daqing, 163318 Heilongjiang, China School of Electrical Engineering and Information, Northeast Petroleum University, Daqing, 163318 Heilongjiang, China Daqing Oilfield Production Engineering Research Institute, Daqing, 163453 Heilongjiang, China Heilongjiang Provincial Key Laboratory of Oil and Gas Reservoir Stimulation, Daqing, 163453 Heilongjiang, China


Introduction
When an oil field is developed, the average formation pressure difference between adjacent blocks gradually occurs, which leads to a formation pressure difference between adjacent areas. The pressure difference between adjacent areas can make the formation produce nonuniform longitudinal deformation and exert a force on the upper formation. For formations with weak upper horizontal structural surfaces, the pressure difference between the adjacent areas can lead to the formation of cracks and slips at the weak horizontal structural surfaces, leading to shear-type casing damage. Taking the Daqing Oilfield in China as an example, at the bottom of the Nen'er member in the upper oil layers, there is a mud-shale standard layer, which can be used as a standard for identifying other layers because of its obvious logging characteristics. This mud-shale standard layer has a small dip angle and a stable distribution throughout the entire region. There are a large number of low-strength fossils in the standard layer, which constitute the horizontal weak structural surface, and have been damaged in a large area, forming a horizontal crack surface. For example, after long-term high-pressure water injection and polymer injection development in the Daqing Oilfield, the Sazhong Development Zone has formed nine casing-damage concentration areas with a total area of 38.34 km 2 . The casing damage rate in the concentration area is as high as 80%, and most are shear-type casing damage [1][2][3][4]. The casing damage in these areas has the characteristics of layer concentration, region concentration, and time concentration.
Domestic and foreign scholars have studied the mechanism of casing damage for many years, and both have confirmed that the change in reservoir injection and production pressure can cause the crack and slip of the horizontal weak structural surface [5][6][7][8]. At present, the typical research method is to use a finite element model and simulations to determine the influence of reservoir deformation on the degree of cracks, slips, and casing damage [9][10][11]. However, with the long-term overpressure development of oil fields, casing damage is increasingly frequent, and even regional casing damage is formed in the standard mud-shale formations of the Daqing Oilfield [12][13][14]. Because the mechanism of slip casing damage on the mud-shale crack surface is still not clear, and the mechanical model of formation slip and experimental methods have not been formulated, only field experience or conservative methods can be generally used for prevention measures at present. Additionally, a large area of injection reduction and closure of more than 100 injection wells in the casing damage area seriously affects the production of oil fields. The means of prevention and control are not informed, and the pertinence of the problem is poor, resulting in a lack of active monitoring of casing damage. Understanding the mechanism of casing damage is the key point to preventing, controlling, and managing casing damage on weak horizontal structure surfaces. Therefore, establishing a laboratory experimental device that simulates the sliding casing damage on a horizontal crack surface and reveals the mechanism of regional casing damage is critical.
It is important to verify the mechanism of standard layer shear-type casing damage through laboratory experiments [15,16]. However, owing to the limitations of experimental conditions and measurement methods, there are only direct dynamic shear-type casing damage experiments and underground mining ground deformation experiments. Sun [17] used a high-strength expansion agent to simulate formation stress; they inserted the simulation casing into the middle of the device, then inserted the expansion agent into one side of the experimental device, and, through the expansion agent, simulated the formation near the wellbore to shear the casing directly. Yin [18] conducted a direct shear test on a J55 casing with a TS-1 casing damage testing machine independently developed for in situ stress by the Laboratory of China Petroleum Exploration and Development Research Institute. The above two experiments simulated the formation of shear-type casing damage, but the actual casing damage process was greatly simplified in the experiments. In the experiments, the source of shear force was a direct thrust from two sides, but it could not explain the source of shear stress in the actual formation. In addition, there are no other reports on laboratory experimental devices for shear-type casing damage in the petroleum development field.
In an indoor physical simulation experiment of formation deformation in the mining field, Gong et al. [19] established a simulation experiment device for mining subsidence under complex geological conditions by using a similar material model with a scale of 1 : 200, model frame length of 2.1 m, and height of 1 m. The model was divided into 10 layers. It was laid in layers, and the layers were separated using an appropriate amount of mica powder. Zou and Chai [20] established a similar two-dimensional model test-bed using similar materials for simulated formation and measured the displacement state of the simulated formation after mining using laser reflection. However, the above experimental equipment was mainly aimed at underground air excavation to surface deformation and cannot determine the role of regional formation pore pressure difference [21]. The slip of the horizontal weak structural surface was caused by the relative slip between layers inside the formation, which was quite different from the measurement purpose of this experiment. In addition, there are no reports on laboratory test equipment and methods for shear-type casing damage caused by formation slip.
In this study, we propose a set of experimental devices that can simulate and measure formation slip. Based on the abnormal similarity theory with different horizontal and vertical scale similarity ratios, we carried out a formation slip simulation experiment under injection-production interregional pressure difference. The corresponding measurement method and experimental steps of formation slip were designed according to the slip characteristics of the mud-shale crack surface. Based on the abnormal similarity theory, the corresponding relationship between the experimental model and the prototype formation slip was established, and the formation slip at different positions under different conditions was measured and compared with the theoretical calculation results, and the factors influencing the formation slip were analyzed.

Materials and Methods
2.1. Experiment Method Based on the Similarity Model. The formation swelling and deformation caused by pore pressure changes are from the underground target layer to several hundred meters above the ground and several thousand meters in the horizontal direction, which reflects two prominent problems: a big volume is involved [22][23][24][25] and the great difference between horizontal and vertical proportions. The sliding layer is not a water flooding reservoir, but the nonproduction horizontal layer is approximately 70 m above it. The geological conditions are shown in Figure 1.
To solve the above problems, the simulated formation material was used instead of the actual formation, and a similar method was adopted to conduct the simulation experiment. Model experiments with different aspect ratios were designed according to the characteristics of the formation. The dimensional analysis method (π theorem) is applied to derive the similarity relationship between each physical quantity and determine the similarity constant of each physical quantity. In this experiment, the horizontal and vertical dimensions of the geometric parameters were analyzed with different scaling coefficients, and the length dimensions representing the horizontal and vertical directions were taken as different dimensions, ½L r and ½L z , respectively.

Geofluids
According to the definition of each physical quantity, the dimensions of each physical quantity based on the above assumptions were obtained, as shown in Table 1. D is the radial size (m) H is the vertical dimension (m), u is the radial displacement (m), w is the vertical displacement (m), σ r is the radial stress (MPa), σ z is the vertical stress (MPa), ε r is the radial strain (nondimensional), ε z is the vertical strain (nondimensional), E r is the radial Young's modulus (MPa), E z is the vertical Young's modulus (MPa), μ zr is the radial Poisson's ratio (nondimensional), and μ rz is the nondimensional vertical Poisson's ratio.
In this physical simulation experiment, the dimensions ½T and ½M corresponding to the basic dimension are determined in proportion, so only the change in dimension ½M needs to be considered. Therefore, D, H, and E r are taken as the basic dimensions to obtain the dimensionless π terms of other physical quantities, as shown in Table 2.
The horizontal and vertical similarity ratios of the physical model and formation prototype were different. According to the distortion compensation model, the respective π terms are assumed to form a product relation. Then, the relation between the independent π terms is as follows: where π ai is the π term with dimensional parameters. When the parameters of some similar materials fail to meet the similarity ratio requirements in model design and material selection, they distort the model system. π u1 and π u2 are dimensionless parameter π terms, e i is distorted idempotent, and A i is an infinite constant. The prediction coefficient δ i generated by model distortion is expressed as follows: where δ i is the coefficient of prediction, β i is the distortion coefficient of each parameter π term, β i = π im /π ip , and π im and π ip are dimensionless π terms for each parameter. The subscript "p" represents the parameters of the formation prototype, and the subscript "m" represents the parameters of the experimental model. Thus, the similarity ratio relation can be obtained as follows: where C f is the similarity ratio of the parameter f = ðD , H, u,⋯,μ rz Þ and δ i is the prediction coefficient generated by the model distortion.
Both experimental models and stratigraphic prototypes, simulated stratigraphic materials, and stratigraphic rocks can be considered as isotropic materials. Therefore, the following relations are obtained: Then, the similarity ratio of the dimensionless parameters is obtained as follows: In the design of a structure dissimilar model, factors such as type, material, manufacturing conditions, loading capacity, arrangement of measuring points, and equipment conditions should be considered comprehensively in the geometric size of the model. An experimental model that is too small will lead to high requirements for experimental processing and instrument measurement accuracy, and an excessively large experimental model will inconvenience the experiment. Therefore, it is necessary to consider these factors in selecting the experimental size and materials. For this experiment, given the lateral and vertical dimensional similarity ratio and mechanical parameters of the simulated stratum, the similarity ratio of the other parameters and the corresponding prediction coefficients can be obtained.
Considering the elastic modulus, strength, and actual processing capacity of each material comprehensively, soft polyvinyl chloride was selected as the simulated formation material of the similar model. Considering that the actual water injection pressure of the experimental model was generally less than 1.5 MPa, the simulated stress of the formation material is in the linear elastic deformation range. At the same time, under the action of interregional pore pressure, the formation prototype was in the linear elastic deformation range regardless of the oil layer, overburden, or even the surface position and did not involve nonlinear yield deformation and failure. Therefore, for the experiments, this polymer material could be used as a simulated stratum instead of rocks.
Thus, the experimental model and formation prototype parameters were obtained, as shown in Table 3.
According to the similarity theory, the similarity ratio and prediction coefficient of each parameter between the actual stratigraphic prototype and the experimental model of the upper and lower layers are shown in Tables 4 and 5, respectively.
The values of the corresponding parameters of the actual formation can be obtained by multiplying the parameters of the simulated formation in the experimental model with the corresponding parameter similarity ratio.

Equipment.
The crack and slip of the horizontal weak structural surface at the bottom of the Nen'er member is caused by the vertical deformation of the bottom oil layer, which is characterized by no stress at the top and deformation at the bottom under the action of formation pressure. However, deformation through the formation pore pressure requires the accumulation of a large thickness, which cannot be achieved by laboratory experiments to produce the bottom deformation of the formation slip. In a previous finite element study, bottom water injection pressurization pressure was similar to that of interregional formation pore pressure. For this reason, we used bottom water injection pressurization to simulate the interregional formation pressure difference. However, this type of experimental equipment is unprecedented and requires additional design and manufacturing. To experimentally achieve formation slip deformation with different radial sizes, a flange adapter was used, which could correspond to the experimental materials with diameters of 250 mm, 300 mm, and 400 mm. In this experiment, a two-dimensional axisymmetric model was used. To ensure that the deformation of the equipment is very small when simulating formation deformation, steel equipment with a height of 1000 mm and a wall thickness of 40 mm that can withstand a pressure of 5 MPa without leakage was used. The maximum bearing capacity of the designed physical simulation device of a weak structure surface slip under the action of the interlayer pressure difference was 5 MPa, and the diagram of its structure is shown in Figure 2.
In the experiment, a combination of a manual pump and a constant-flux pump was used to inject liquid into the device to achieve pressurization. The maximum allowable pressure of the constant-flux pump was 42 MPa, the minimum flow rate was 0.001 mL/s, and the maximum flow rate was 1.000 mL/s. First, the manual pump was used for rapid pressurization. Then, when all the liquid in the pump was pumped into the equipment to reach the experimental pressure, the constant-flux pump was used for slow and accurate pressurization until the pressure required by the experiment was reached. The manual liquid booster pump, advection pump, and formation deformation test device were connected through six channels, and a pressure gauge was installed on the six channels to measure the internal liquid pressure ( Figure 2).

Method for Measuring Formation Slip.
By connecting the water injection pressurization system with the main experimental equipment, simulated formation expansion can be achieved. The vertical displacement of the formation was accurately measured with a displacement measurement device with an accuracy of 0.001 mm and a range of 15.000 mm. For continuous and stable measurements of the formation with the displacement measurement device, a custom steel plate was made to fix it. The middle and both ends of the customized steel plate were equipped with a drill hole. The drill hole in the middle was used to install the 4 Geofluids displacement measurement device. The two ends could be fitted to the flange of the formation-deformation measuring device.
The soft adhesive of epoxy resin selected in the formation slip experiment, with an operation time of 4 h and solidification time of 12 h, was used to simulate the casing deformation and slip under the formation-deformation condition. The epoxy resin adhesive was formed by gradual thickening and condensation after mixing liquids A and B. The hardness and softness of the epoxy resin could also be adjusted according to the different formulations and solidification times. In this experiment, a soft adhesive with a solidification time of 12 h was selected. Both liquids A and B were of this type of epoxy resin adhesive and were low-viscosity solutions. After mixing, the viscosity gradually increases; good fluidity was maintained in the first 30 min, and the initial setting time was 4 h. From 4 h to 10 h, it became a solid, and the mobility almost disappeared, but the plasticity was strong; that is, after the deformation under the action of an external force, the shape could not be recovered. After 10 h, the plastic characteristics gradually trans-formed into elastic characteristics. After 12 h, the epoxy resin adhesive was set with strong elasticity and a certain memory effect; that is, after the external force was applied to deform it, the deformation could be recovered by removing the external force. The AB epoxy resin adhesive maintained a long fluidity time. In the slip experiment of a weak horizontal structural surface, the amount of expansion of the simulated material should be kept constant; thus, a material with moderate solidification time was needed to simulate the deformation condition of the casing in maintaining formation uplift. However, the preparation process took a long time, so the epoxy resin needed to maintain a certain period of fluidity before solidification. The solidified epoxy resin adhesive maintained its deformed state during solidification and returned to its original state when the force was removed. The epoxy resin adhesive is a fluid before curing and can deform with the internal deformation of the casing. Soft epoxy adhesives have good scalability after curing and are not easily damaged by large deformations. This type of epoxy resin deforms during disassembly, but the deformation is within the elastic range. After the epoxy resin adhesive was Vertical Young's modulus Vertical Poisson's ratio Vertical Poisson's ratio π 7 = μ rz /u 2 w −2   5 Geofluids removed, it returned to the shape of the formation slip and did not affect the slip measurement results.
The epoxy resin was solidified in the hole, and after removing all external forces, the epoxy resin deformed to the solidified shape, that is, the hole position slippage state due to the formation expansion. This method can be used to easily observe the variation characteristics of slip at different planes for formation expansion induced by hydraulic fracturing. To avoid the mutual interference of the stress concentration near the hole, the spiral layout mode was selected after observing the hole layout. The holes of the simulated material were all distributed on the helix. This ensured that the spacing between the holes was the same and the straight distance between the holes was not too short. In the experiment, the distance between the actual holes was extended to an average of 34.6 mm by spiral distribution drilling. The shortest distance between two holes was 25.4 mm. The stress between holes with a diameter of 2 mm was very small.
If the method of injecting epoxy resin directly into the simulated formation was adopted, the epoxy resin glue would be completely cemented in the inner wall of the observed hole, and the epoxy resin glue would not be removable or would be seriously damaged after removal. Therefore, a method of inserting a polytetrafluoroethylene (PTFE) pipe into the hole to simulate the casing was designed. PTFE has a strong resistance to organic solvents and is not easy to glue. The epoxy resin adhesive did not stick to the PTFE after solidification and could be easily removed. In addition, in order to prevent the leakage of epoxy resin glue during injection and the experiment, a length of 1-2 mm at the bottom of the simulated casing of the PTFE pipe was sealed in advance. Because the epoxy resin adhesive was transparent, dye was added to the epoxy resin adhesive to change its color in order to facilitate observation in the experiment; the dye was very low in content, with only 0.1% concentration dying the epoxy resin. After the formation slip experiment was completed, the morphology of the epoxy resin glue was the internal morphology of the observation hole. Owing to the small amount of deformation in the observation of the inner shape of the hole, there would be a large error in the direct measurement using the measuring tool; thus, the deformation could not be measured directly. In the actual measurement process, the deformation was obtained by enlarging a picture of the epoxy resin for accurate measurement and then reducing the corresponding multiple, as shown in Figure 3.
Through the magnification measurement method, the transverse shear deformation at the position of the crackslip surface in the observation hole was obtained. Combined with the relationship between the formation slip and transverse shear deformation of the observation hole, the interlayer slip at the observation hole position was calculated, as shown in the following equation:   6 Geofluids is the reaming rate of the bit (for PVC, an easy-to-drill material, the reaming rate is generally 5%), and D c is the simulated casing outer diameter (mm).

Experimental Procedure
(1) The flange conversion joint corresponding to the experimental size of the experimental material was installed on the experimental equipment, and the water tank was filled with water in advance (2) The simulated material was fixed to an experimental device. The simulated stratum was placed in the experimental device according to the designed sequence, and the bottom of the simulated stratum was sealed to prevent liquid leakage. The top bolt was tightened to fix the simulated formation boundary. Subsequently, the sealing and connectivity of the experimental device were examined (3) The boreholes were marked on the simulated formation device according to the designed position, and the simulated boreholes were drilled with an electric drill, with the drilling depth guaranteed to meet the experimental requirements. This is shown in Figure 4 (4) Timing immediately started at the configuration of the epoxy resin adhesive. The prepared simulated casing was inserted into the simulated well, and the configured epoxy resin adhesive was injected into the simulated casing with a syringe installed with a long glue injection needle. The casing ensures that the epoxy resin adhesive will not flow into the microtensioned layer when the formation slips. The thickness of the casing was 0.1 mm, which hardly changed when the formation was deformed, and it did not affect the experimental results. The viscosity of the epoxy resin adhesive increased sharply. After 15 min, although the epoxy resin adhesive is still a fluid, it cannot be injected into the casing through a syringe, resulting in the failure of the experiment. Therefore, all epoxy resins should be injected into the simulated casing within 15 min (5) The displacement measuring device was installed on the formation simulation device to detect the swelling displacement of the simulated material (6) A manual pump was used to inject water into the pressurizing device and to apply pressure to the simulated material for recording the longitudinal deformation and water injection pressure data of the formation. The changes in value of the pressure gauge and displacement measuring device were observed. The water injection pressure in the experiment is different from the pore pressure in the actual formation; therefore, it cannot be directly used as the pore pressure of the formation. It must be equivalent to the pore pressure of the reservoir required by the simulation of the uplift of the formation center in the experiment and the corresponding uplift of the actual formation. During this period, it is necessary to find the corresponding relationship between the pore pressure variation and the experimental water injection pressure by numerical simulation. When the displacement measurement device reached the experimental design of longitudinal deformation, the pressure was stopped to keep the pressure unchanged. The above steps were carried out within 1 h after the epoxy resin adhesive was mixed, that is, before the initial setting of the epoxy resin adhesive (7) A constant-flux pump was used to keep the pressure value of the bottom booster unchanged. Then, the epoxy resin adhesive took 12 h to set (8) After 12 h, the drainage valve was opened to relieve pressure. Then, the pressure of the experimental equipment was removed, and the simulated casing in the simulated formation was removed. After the simulated casing was cut open, the epoxy resin adhesive reflecting the formation slip shape of each simulated well was removed and pasted on the coordinate paper with formation information to complete the slip measurement

Experimental Result
The experimental device we developed was used to carry out several groups of formation slip experiments, and the formation slip distribution was measured under different pressures, simulated formation diameters, reservoir-buried thicknesses, and crack surface depths, and the internal morphology of the simulated well after deformation was obtained, as shown in Figure 5.
The mechanical parameters of the experimental model and the formation prototype were based on the data in Table 3. We conducted a total of 14 groups of experiments, nine of which were effective. Nine valid sets of experimental data are listed in Table 6.
The slip distributions on the crack surface between the top and bottom formations under different water injection pressures (serial numbers 1, 2, and 3) were compared for a simulated formation diameter of 250 mm, as shown in Figure 6.
The slip of the crack surface is a relative displacement, that is, the difference between the radial displacements of the upper and lower layers of the horizontal crack surface. In the simulation experiment, the displacement of the upper and lower parts of the crack surface could not be measured, and the slip could only be determined by observing the solidification shape of the epoxy resin adhesive in the borehole. The experimental results showed that after the deformation of the simulated formation bottom, the simulated formation above was deformed by pressure from the bottom. Although the layer in the middle of the simulated formation could maintain the longitudinal continuity of the stratum, it had a relative displacement in the horizontal direction. Although the surrounding area of the simulated formation was completely fixed and the center of the formation did not slip under symmetric action, there was also a slip in the middle. It can be seen from Figure 6 that under different vertical displacements of the formation center, the variation trend of the interlayer slip was the same, which first increased and then decreased with an increase in distance from the center. By comparing the slippage of different vertical displacements at the center of the formation, it was found that the slippage on the horizontal crack surface of the formation had a linear relationship with the longitudinal displacements at the center.
The actual slip of the formation can be obtained by multiplying the slip value by the transverse proportionality coefficient. The interregional formation pressure difference did not correspond to the water injection pressure at the bottom of the test equipment. The pressure difference between the adjacent regions of the actual formation can be calculated according to its effect on the formation. First, the corresponding relationship between the formation pressure difference and the vertical deformation of the surface center was calculated based on the finite element simulation. Then, according to the simulated uplift height multiplied by the longitudinal length ratio, the actual maximum vertical deformation height of the formation can be obtained to determine the regional formation pressure difference of the actual formation. Through calculations, it was determined that in the actual formation corresponding to the experimental conditions, the slip amounts at different positions corresponding to the interregional pressure difference (serial numbers 1, 2, and 3) were as in Figure 7.
Through scale conversion, the size of the experimental model was converted to the actual size of the formation, as shown in Figure 7. The origin of the r coordinate was the center of the abnormally high-pressure area, and the position r = 1500 m was the center of the low-pressure area. The slip of the crack surface occurred between the highpressure and low-pressure zones, and the slip first increased and then decreased. In actual oil field production, formation slip will shear the wellbore on the slip surface, resulting in shear-casing damage. In the Daqing Oilfield, the crack surface was uncemented to prevent shear-casing damage, leaving room for formation slip. When the formation slip exceeded 60 mm, the formation would cause the shear deformation of the casing and affect the production of the oil and water wells. Based on this amount of slip, the maximum interregional formation pressure difference was 0.69 MPa for no casing damage occurred on this field.
The range of the interregional formation pressure difference also affects the slip on the horizontal crack surface of the horizontal weak structure layer. Therefore, experiments with different diameters were performed to simulate   9 Geofluids formation slip. The diameters of the experiment were 250 mm, 300 mm, and 400 mm (serial numbers 1, 4, and 5). The interlayer slip was measured and analyzed where all the longitudinal displacements of the formation center were 4.20 mm, and it was converted into the actual formation conditions under the formation pressure difference in the same area. The results are presented in Figure 8.
After transformation, the formation pressure difference between regions was 1.59 MPa. With an increase in the range of interregional formation pressure, the maximum value of slippage on the mud-shale crack surface decreased slightly and was extrapolated, and the influence range of the interzonal slippage increased. The top depth of the oil layer was different in different blocks, so the variation degree of slippage on the mud-shale crack surface would also change accordingly.
To study the effect of the reservoir depth on the interlayer slip, different simulated formation thickness experiments were carried out to analyze the formation slip under different standard depth conditions. The simulated material in this group was 400 mm in diameter, and the total thicknesses of the simulated formations were 24, 20, and 16 mm. The corresponding total thicknesses of the actual formations were 960 m, 800 m, and 640 m (serial numbers 7, 5, and 6). The relative position of the mud-shale standard layer remained unchanged, and the thickness ratio of the top and bottom layers was 3 : 1. Through the transformation of the measurement results to the actual formation, the slip distribution on the horizontal crack surface of the standard layer was calculated and obtained for an interregional formation pressure difference of 1.59 MPa and different top depths of the oil layers, as shown in Figure 9.
In terms of numerical value and trend, with the increase in reservoir depth, the maximum slip on the mud-shale crack surface increased; however, the increase was small. Within the range of the reservoir top depth in the Daqing Oilfield, the slip varied by 25% at most. The maximum value of slip and the distribution law of slip showed no obvious changes in the three groups of experiments.
To determine the influence of depth on the slippage of the mud-shale crack surface, 400 mm was selected to simulate the formation and ensure that the total thickness of the simulated formation was uniformly 20 mm. Three combined conditions of 13 mm: 7 mm, 15 mm: 5 mm, and 17 mm : 3 mm thickness of the upper and lower layers were studied, and the corresponding horizontal crack surface depths were 520 m, 600 m, and 680 m, respectively (serial numbers 8, 5, and 9). Through the transformation of the measured results to the actual formation, the slip distribution on the horizontal crack surface of the standard layer with different top depths of oil layers was calculated and obtained when the interregional formation pressure difference was 1.59 MPa, as shown in Figure 10.
The calculation results showed that the interlayer slip increased with the upward movement of the relative position of the standard layer. When the mud-shale standard layer was located in the middle of the oil layer and the surface, the maximum slip was caused by the interregional formation pressure difference. This law was very similar to the distribution law of shear stress; that is, under the action of interregional pressure difference, the formation shear stress accumulated from the bottom, reached a maximum in the middle of the formation, and became zero on the surface. In addition, from the calculation results, for different interregional formation pressure differences, action ranges, formation thicknesses, and crack surface depths, the comparison between the experimentally measured data and the theoretically calculated results was in good agreement. In the simulation experiment, a mud-shale standard interlayer slip was observed, the amount of interlayer slip under different conditions was measured, and the accuracy of the theoretical model of formation slip under the effect of interregional formation pressure difference was verified.
The experimental analysis showed that the interregional formation pressure difference between the bottom oil layers was the driving force for the slip of the mudshale standard layer in the formation, and the formation pore pressure did not need to act directly on the