Disclosure of Invention
Aiming at the defects existing in the prior art, the invention aims to provide a method for identifying a horizontal well fracturing engineering dessert based on MSE, which solves the problem that the horizontal well fracturing engineering dessert can be accurately and rapidly identified under the condition of lacking complete logging data, and the geological engineering integration is promoted.
In order to achieve the above purpose, the technical solution adopted by the invention is as follows:
A method of identifying a hydraulic fracturing engineering dessert based on MSE, comprising the steps of:
S1, acquiring raw data required for calculating mechanical specific energy from well drilling history data, acquiring gamma well logging data from well logging data, and acquiring drilling time data from the well drilling history data.
S2, selecting a Fred E.Dupriest model as a basic formula for calculating the mechanical specific energy, substituting the original data into the basic formula to calculate the mechanical specific energy data of stratum rocks with different depths along the well depth direction.
S3, carrying out noise reduction treatment on the mechanical specific energy data, the gamma logging data and the drilling data by adopting a moving average method to respectively obtain MSE values, GR values or drilling times of different depths along the well depth direction.
S4, dividing the MSE value by the GR value at the same well depth to obtain MSE/GR values, drawing stratum rock MSE/GR curves at different well depths, and thus fusing mechanical specific energy data with gamma logging data to construct an MSE+GR model, and primarily identifying the region with the MSE/GR value not less than 15 in the stratum rock MSE/GR curves as an engineering dessert.
S5, dividing the MSE/GR value by the drilling time of the same well depth to obtain an MSE/(GR time of drilling), drawing fusion graphs of MSE data, gamma logging data and drilling time of stratum rocks at different well depths, further fusing an MSE+GR model and the drilling time data, and preliminarily determining the identification result of the engineering dessert in the region with the MSE/(GR time of drilling) value of more than 5 in the fusion graphs.
S6, preliminarily verifying the identification result of the engineering dessert through the fracturing curve, extracting fracturing curve data from fracturing data, counting the cracking pressure and the pressure change condition after stratum fracture, judging whether the fracturing data correspond to the engineering dessert, judging whether the cracking pressure is not less than 50MPa and the pressure is not less than 7.5MPa after the cracking, and conforming to the characteristics of the engineering dessert.
If the fracturing data does not correspond to the engineering dessert, the analysis and calculation results of the rock cuttings mineral composition and the brittleness index are adopted to verify the identification result of the engineering dessert, the well section with the rock cuttings brittleness index larger than 0.7 can be determined as the engineering dessert, and the identified dessert area is verified to have larger rock hardness and brittleness, so that a complex fracture network is easy to form.
Further, the raw data required for calculating the mechanical specific energy in S1 includes drilling depth, weight on bit, bit size, rotation speed, torque and mechanical rotation speed, and the fracturing curve of each fracturing segment along the well depth direction is obtained from the drilling well history data.
Further, the calculation formula of the Fred E.Dupriest model is adopted in S2,
;
In the formula,Is the mechanical specific energy; Is weight on bit; is the area of the drill bit; is the rotation speed; is torque; is the rate of penetration.
Further, the MSE value of a certain depth obtained after the noise reduction processing in S3 refers to an average value of MSE values at 5 adjacent depths above and below the depth, respectively.
The GR value at a certain depth obtained after the noise reduction processing in S3 means an average value of GR values at 5 adjacent depths above and below the depth, respectively.
The drill time at a certain depth obtained after the noise reduction processing in S3 is an average value of drill times at 5 adjacent depths above and below the depth, respectively.
Further, the step of preliminarily verifying the engineering dessert by the fracturing curve in S6 comprises the following steps:
and S61, analyzing each fracturing segment by combining the fracturing curve extracted from the fracturing data, analyzing the region with MSE/(GR) value more than 5 as engineering dessert, and identifying the region with MSE/GR value not less than 15 in the fracturing segment range as the engineering dessert fracturing segment.
S62, counting the stratum fracture pressure and the pressure drop condition after stratum fracture in the fracture curve, further screening out a fracture section with the fracture starting pressure not less than 50MPa and the pressure drop after fracture not less than 7.5MPa from the fracture sections with the discharge capacity of 4.5-5.5m 3/min, conforming to the characteristics of engineering dessert, and if the fracture section with the pressure drop after fracture is less than 7.5MPa, not serving as the engineering dessert.
Further, sampling a rock fragment sample at each well section of the well, carrying out mineral composition analysis on the rock fragment sample to obtain the content of quartz minerals, the content of carbonate minerals and the content of clay minerals of the rock fragment sample at each well section, calculating the rock brittleness index of each well section, wherein the calculation formula of the rock brittleness index is as follows,
;
In the formula,Is the content of quartz minerals; Is the content of carbonate mineral; is the clay mineral content, and BI 3 is the brittleness index.
Further, the calculated results of the rock brittleness index of each well section are marked in the fusion graph, and the well section with the rock brittleness index greater than 0.7 corresponds to the area with the MSE/(GR) value greater than 5, and the engineering dessert is determined.
Compared with the prior art, the invention has the beneficial technical effects that:
1. the invention can quickly and accurately identify engineering desserts by gamma logging data in logging data without depending on other data in the logging data, improves engineering dessert identification efficiency and promotes geological engineering integration.
2. Under the condition of incomplete logging data, the invention provides a method for identifying the engineering dessert of the horizontal well through an MSE+GR model, gamma logging data can be obtained after well drilling is completed, a fusion diagram is obtained by combining the gamma logging data with the well history data, the engineering dessert is preliminarily determined, and then a well section conforming to the engineering dessert is rapidly judged by combining a pressure curve.
3. The invention provides a method for verifying the recognition result of the engineering dessert through pressure curve and rock brittleness analysis for the first time so as to verify the accuracy of the horizontal well engineering dessert recognition.
Detailed Description
The invention is described in detail below with reference to the attached drawing figures:
example 1, in connection with fig. 9, a method of identifying a hydraulic fracturing engineering dessert based on MSE, comprising the steps of:
S1, acquiring raw data required for calculating mechanical specific energy from well drilling history data, acquiring gamma well logging data from well logging data, and acquiring drilling time data from the well drilling history data.
S2, selecting a Fred E.Dupriest model as a basic formula for calculating the mechanical specific energy, substituting the original data into the basic formula to calculate the mechanical specific energy data of stratum rocks with different depths along the well depth direction.
The raw data required for calculating the mechanical specific energy comprises drilling depth, weight on bit, bit size, rotating speed, torque and mechanical rotating speed, and the fracturing curve of each fracturing segment along the well depth direction is obtained from the drilling well history data.
The Fred e.dupriest model is used as the mechanical specific energy calculation model, the mechanical specific energy calculation formula,
;
In the formula (i),Is the mechanical specific energy; Is weight on bit; is the area of the drill bit; is the rotation speed; is torque; is the rate of penetration.
The calculated MSE and GR data are plotted in the same graph (shown in FIG. 1) to obtain corresponding graphs of MSE and GR for different well depths, and the well depths extracted from the well history data are plotted with the time-of-drilling data to obtain a time-of-drilling graph (shown in FIG. 2).
S3, carrying out noise reduction treatment on the mechanical specific energy data, the gamma logging data and the drilling data by adopting a moving average method to respectively obtain MSE values, GR values or drilling times of different depths along the well depth direction. The correction mechanical specific energy calculation model is used for carrying out noise reduction treatment on MSE data, so that the influence caused by abnormal fluctuation is eliminated, the change rules of MSE, GR and drilling are conveniently observed, and the influence of the abnormal fluctuation is weakened by adopting a moving average method.
The MSE value at a certain depth obtained after the noise reduction processing refers to an average value of MSE values at 5 adjacent depths above and below the depth, respectively. The GR value at a certain depth after the noise reduction processing is an average value of GR values at 5 adjacent depths above and below the depth, and the MSE and GR data after the noise reduction processing are plotted in the same graph (shown in fig. 3). The drilling time of a certain depth obtained after the noise reduction treatment is an average value of drilling times of 5 adjacent depths above and below the depth, and drilling time data after the noise reduction treatment is plotted to obtain a drilling time curve correction chart (shown in fig. 4).
Comparing MSE, GR after moving average denoising and well depth with the original data graph, it can be seen that the data fluctuation is obviously frequent before denoising, engineering desserts are not easy to observe and identify, the data after denoising obviously becomes easy to observe, and errors caused by other influences are removed, so that the next processing and analysis are convenient.
S4, dividing the MSE value by the GR value at the same well depth to obtain MSE/GR values, drawing stratum rock MSE/GR curves (shown in figure 5) at different well depths, and thus fusing mechanical specific energy data with gamma logging data to construct an MSE+GR model, primarily identifying an area with the MSE/GR value not smaller than 15 in the stratum rock MSE/GR curves as an engineering dessert, wherein the higher the mechanical specific energy is, the lower the gamma logging curve value is, the higher the stratum hardness brittleness is, the lower the clay content is, and the crack opening degree is easier to keep and a complex fracture network is formed.
S5, dividing the MSE/GR value by the drilling time of the same well depth to obtain MSE/(GR time), drawing fusion graphs of MSE data, gamma logging data and drilling time of stratum rocks at different well depths, further fusing MSE+GR models with the drilling time data, preliminarily determining the identification result of the engineering dessert in the fusion graph in the region with the MSE/(GR time) value larger than 5, further fusing the MSE+GR models with the drilling time data, further reducing the position of an MSE/GR curve, enabling the horizontal axis to be approximately used as a datum line for identifying the engineering dessert, and eliminating the influence of small peaks, so that the engineering dessert is easier to identify, as shown in FIG. 6.
S6, preliminarily verifying the identification result of the engineering dessert through a fracturing curve, extracting fracturing curve data from fracturing data, counting the cracking pressure and the pressure change condition after stratum fracture, judging whether the fracturing data correspond to the engineering dessert, judging whether the cracking pressure is not less than 50MPa, and judging that the pressure after cracking is reduced by not less than 7.5MPa, wherein the characteristics of the engineering dessert are met.
The method for preliminarily verifying engineering desserts by using the fracturing curve comprises the following steps:
s11, analyzing each fracturing segment by combining the fracturing curve extracted from the fracturing data, analyzing the region with MSE/(GR) value not smaller than 15 as engineering dessert, and identifying the well segment with MSE/GR value not smaller than 15 in the fracturing segment range as the engineering dessert fracturing segment.
S12, counting the stratum fracture pressure and the pressure drop condition after stratum fracture in the fracture curve, and further screening out the fracture sections with the fracture starting pressure not less than 50MPa and the pressure drop after fracture not less than 7.5MPa from the fracture sections with the displacement of 4.5-5.5m 3/min, wherein the fracture sections with the pressure drop after fracture not less than 7.5MPa accord with the characteristics of engineering dessert, and the fracture sections with the pressure drop after fracture not less than 7.5MPa are not used as engineering dessert.
If the fracturing data does not correspond to the engineering dessert, the analysis and calculation results of the rock cuttings mineral composition and the brittleness index are adopted to verify the identification result of the engineering dessert, the section with the rock cuttings brittleness index exceeding 0.7 is determined as the engineering dessert, and the identified dessert area is verified to have higher rock hardness and brittleness, so that a complex seam network is easy to form.
Sampling a rock fragment sample at each well section of the well, carrying out mineral composition analysis on the rock fragment sample to obtain the content of quartz minerals, the content of carbonate minerals and the content of clay minerals of the rock fragment sample at each well section, analyzing, calculating the rock brittleness index of each well section, calculating the rock brittleness index according to the following calculation formula,
;
In the formula,Is the content of quartz minerals; Is the content of carbonate mineral; the clay mineral content, and BI 3 is brittleness index and dimensionless.
And marking the calculated results of the rock brittleness indexes of the well sections in the fusion graph, wherein the well sections with the rock brittleness indexes exceeding 0.7 correspond to the sections with MSE/(GR) drilling time values greater than 5, and determining the sections as engineering desserts.
In embodiment 2, the conventional engineering dessert identification method often has the defects of complicated parameter selection, large error of identification result and the like, so that the fracturing cost is increased. The invention accurately and rapidly identifies the dessert of the tight sandstone horizontal well engineering under the condition of incomplete logging data based on the mechanical specific energy theory, and promotes the integration of geological engineering.
Taking the G31-P7 well as an example, data required for calculating mechanical specific energy is extracted from site data, the mechanical specific energy is calculated, and the mechanical specific energy is fused with gamma logging data and drilling data to obtain an MSE/(GR time-of-drilling) fusion map of the G31-P7 well. The G31-P7 well is subjected to 19-section fracturing, the engineering dessert is analyzed by taking the MSE/(GR) value which is obviously higher (more than 5) as the engineering dessert, and the 4 th, 6 th, 8 th, 10 th, 12 th, 13 th, 17 th, 18 th and 19 th fracturing sections are provided with areas with MSE/GR values not less than 15, and the areas are initially identified as the engineering dessert fracturing sections. The fracture pressure and pressure drop after fracture of the formation in the fracture curve were counted as shown in table 1.
TABLE 1 formation fracturing pressure at different fracturing stages and pressure drop after fracture
Further, as can be seen from table 1, at similar displacement, the cracking pressure of the fracturing sections of 2, 3, 4, 5, 7, 12, 16, 18 and the like is not less than 50MPa, the pressure drop after cracking is not less than 7.5MPa, and the characteristics of engineering dessert are met.
After the MSE/(GR) fusion map is combined, engineering dessert identification analysis can be carried out, so that MSE/(GR) values of 2, 3, 4, 5, 7, 12, 16, 18 and 19 well sections can be obtained preliminarily, and by analyzing the fracturing curve, the fracturing sections with the discharge capacity of 4.5-5.5 m 3/min can be screened out, wherein the MSE/(GR) values of the 2, 3, 4, 5, 7, 12, 16 and 18 well sections are higher. Is approximately consistent with the depth of the stratum of the earlier predicted engineering dessert. And the cracking pressure of the fracturing section is not less than 50MPa, the pressure drop after cracking is not less than 7.5MPa, and the characteristics of engineering desserts are met. The pressure of the 19 th section after cracking is reduced to 7.5MPa, which is probably because the lower drilling time in the fracturing process leads to larger numerical value, and no obvious pressure drop exists after the rock is cracked for a long time, so the 19 th well section is not used as an engineering dessert.
The fracturing data basically corresponds to the engineering dessert identification result, and the effectiveness of the proposed engineering dessert identification method is preliminarily verified. However, the fracturing data does not completely correspond to the engineering dessert identification result, because the fracturing data is influenced by the formation factors and the engineering factors, so that the engineering dessert identification method is not comprehensive only by the fracturing data, and the engineering dessert identification method needs to be further verified through the rock debris analysis result.
The method for verifying the engineering dessert identification result by using the rock brittleness evaluation method is to avoid limitation of fracturing curve data, and uses rock debris mineral composition and brittleness index analysis calculation results to verify the engineering dessert identification result.
For the G31-P7 well, a cuttings sampling wellbore section was designed, as shown in Table 2 below. Wherein, there is a well section with higher MSE and a well section with lower MSE, so that the sampling is representative.
Table 2 cuttings sampling wellbore section design
Further, 4 cuttings samples were successfully obtained, mineral composition analysis was performed on these samples, and the rock brittleness index (dimensionless) of four cuttings samples was obtained according to the calculation formula of the rock brittleness index, as shown in table 3 below, the four cuttings samples were in the order of 5366, 4882, 4888 and 5534 from the top to the bottom, wherein the first three cuttings brittleness indexes were all greater than 0.7, reflecting that the area was relatively hard and brittle, and could be identified as an engineering dessert.
Table 3 results of successful sampling wellbore section test of cuttings
The results of the cuttings brittleness index analysis are noted in a fusion plot of MSE and gamma and drill time data, as shown in fig. 8 below. It can be seen that the region with a higher chip analysis brittleness index corresponds to the region with a higher MSE/(GR) during drilling, thereby verifying that the dessert region identified by the identifying method Cheng Tiandian has higher rock hardness and is easy to form a complex stitch net.
The parts not described in the invention can be realized by adopting or referring to the prior art.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the description of the present invention, it should be understood that the directions or positional relationships indicated by the terms "upper", "lower", "front", "rear", "left", "right", etc., are based on the directions or positional relationships shown in the drawings, are merely for convenience of describing the present invention and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present invention.
It should be understood that the above description is not intended to limit the invention to the particular embodiments disclosed, but to limit the invention to the particular embodiments disclosed, and that the invention is not limited to the particular embodiments disclosed, but is intended to cover modifications, adaptations, additions and alternatives falling within the spirit and scope of the invention.