Reservoir fracture effectiveness grade evaluation method

Abstract

The invention discloses a reservoir fracture effectiveness grade evaluation method. The method comprises the following specific steps: preparing work area construction data and calibration well data; drawing a construction plane contour map according to the construction data of the work area, and interpolating to obtain a construction data volume; substituting the constructed data volume into a curvature algorithm to obtain a curvature data volume, and extracting a corresponding curvature value from the curvature data volume according to the calibration well data; performing correlation analysis on the calibrated well data and the extracted curvature value to obtain a fitting formula; and substituting the yield standard into the fitting formula to obtain a curvature range, and establishing a reservoir fracture effectiveness grade evaluation standard for evaluation. According to the reservoir fracture effectiveness grade evaluation method, fracture effectiveness evaluation parameter values of different types of gas wells such as a low-yield well, an industrial airflow well and a high-yield well are determined through establishment of fracture effectiveness evaluation parameters, the reservoir fracture effectiveness grade is evaluated on the surface, exploration deployment is guided, and the gas testing yield is improved.

Description

Technical Field

[0001] The invention belongs to the technical field of oil and gas exploration and development, and relates to a reservoir fracture effectiveness grade evaluation method. Background Art

[0002] For fractured carbonate gas reservoirs, the effectiveness of reservoir fractures is closely related to the level of natural gas production. At the same time, fracture effectiveness is closely related to oil and gas migration, accumulation and reserves, and is one of the main evaluation parameters of reservoir "sweet spots". Therefore, reservoir fracture effectiveness level evaluation is one of the important tasks in oil and gas exploration and development. At present, the commonly used fracture effectiveness level evaluation methods mainly include: geostatistical method, well logging evaluation method, seismic method, etc. The geostatistical method directly evaluates the effectiveness of fractures through field profiles and cores. This technology is simple, direct and efficient, but the field profile is easily affected by external forces such as weathering, and the acquisition of cores is affected by drilling costs and has a low coverage. The well logging evaluation method can also only characterize the effectiveness of formation fractures around the well wall, and the interpretation results are somewhat subjective. For fractured gas reservoirs, the better the fracture effectiveness, the higher the natural gas production. Although seismic prediction has advantages in identifying fracture distribution on the surface, it is currently difficult to evaluate the effectiveness of fractures. Summary of the invention

[0003] The purpose of the present invention is to provide a reservoir fracture effectiveness grade evaluation method, which solves the problems of low coverage and difficult evaluation in the fracture effectiveness evaluation in the prior art.

[0004] The technical solution adopted by the present invention is a reservoir fracture effectiveness grade evaluation method, which is specifically implemented according to the following steps:

[0005] Step 1: Prepare the structural data of the work area and the calibration well data;

[0006] Step 2: Draw a structural plane contour map based on the work area structural data in step 1, and interpolate to obtain a structural data volume;

[0007] Step 3, substituting the structural data volume obtained in step 2 into the curvature algorithm to obtain a curvature data volume, and extracting the corresponding curvature value from the curvature data volume in step 2 according to the calibration well data in step 1;

[0008] Step 4: Perform correlation analysis on the calibration well data in step 1 and the curvature value extracted in step 3 to obtain a fitting formula;

[0009] Step 5: Substitute the production standard into the fitting formula of step 4 to obtain the curvature range, and establish the reservoir fracture effectiveness grade evaluation standard according to the curvature range for evaluation.

[0010] The present invention is also characterized in that:

[0011] In step 1, the structural data of the work area are obtained from drilling analysis data or seismic interpretation data, including horizontal and vertical coordinates, structural values, and structural data of the work area.

[0012] In step 1, the calibration well data includes the horizontal and vertical coordinates of the well location, the test gas production P gas The test gas production includes low-yield, industrial gas flow and high-yield test gas production.

[0013] Step 2 is implemented as follows:

[0014] Step 2.1, draw a structural plane contour map based on the work area structural data in step 1;

[0015] Step 2.2: Use geological modeling software to interpolate the plane contour map of step 2.1 to obtain the structural data volume of the work area.

[0016] Step 2.1 is specifically as follows: using industrial drawing software, draw a structural plane contour map based on the work area structural data of step 1, check whether there are singular values ​​in the structural plane contour map, and if so, correct the layered data until there are no singular values.

[0017] Step 3 is implemented as follows:

[0018] Step 3.1: Substitute the structural data obtained in step 2 into the geological modeling software to obtain the maximum principal curvature K of the structure. max , minimum principal curvature value K min And the directional coefficient d, if K max , K min If both are greater than 0, then d = 1; otherwise d = -1;

[0019] Step 3.2: According to E f =d*K min *K max Calculate the curvature data volume E f ;

[0020] Step 3.3: According to the ordinate and abscissa coordinates (x, y) of the well calibrated in step 1, the curvature data volume E in step 3.2 is obtained. f Extract the corresponding curvature value E f (x, y).

[0021] Step 4 is as follows: the test gas production P of the well calibrated in step 1 is gas The curvature value E extracted in step 3 f (x, y) is subjected to correlation analysis to obtain the functional relationship between the two, that is, the fitting formula.

[0022] Step 5 is as follows: Substitute the production standards of low-yield wells, industrial gas flow wells, and high-yield wells of different oil companies into the fitting formula of step 4 to obtain the E corresponding to different production standards. f Value, E f The range of values ​​is the curvature range. According to the curvature range, a reservoir fracture effectiveness grade evaluation standard is established. Then, according to the reservoir fracture effectiveness grade evaluation standard, a favorable area range is established for evaluation.

[0023] The beneficial effects of the present invention are:

[0024] By establishing fracture effectiveness evaluation parameters, the fracture effectiveness evaluation parameter values ​​of different types of gas wells such as low-yield wells, industrial gas flow wells, and high-yield wells are determined, and the reservoir fracture effectiveness level is evaluated on the surface, which can guide exploration deployment and increase test gas production. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1 is a flow chart of a reservoir fracture effectiveness grade evaluation method of the present invention;

[0026] Figure 2 It is a structural plane contour map of the reservoir fracture effectiveness grade evaluation method of the present invention;

[0027] Figure 3 is the test gas output P in the method of the present invention gas and the evaluation parameter E f Function relationship fitting graph;

[0028] Figure 4 It is a contour map of the reservoir fracture effectiveness grade evaluation plane map in the method of the present invention. DETAILED DESCRIPTION

[0029] The present invention is described in detail below with reference to the accompanying drawings and specific embodiments.

[0030] The reservoir fracture effectiveness grade evaluation method of the present invention has the following process: Figure 1 As shown, the specific implementation steps are as follows:

[0031] Step 1: Prepare the structural data of the work area and the calibration well data;

[0032] The structural data of the work area mainly include horizontal and vertical coordinates and structural values. The structural data of the work area is obtained from drilling analysis data or seismic interpretation data;

[0033] The calibration well data in the work area mainly include the horizontal and vertical coordinates of the well location, the test gas production P gas , the test gas production includes low-yield, industrial gas flow, and high-yield test gas production;

[0034] Step 2: Draw a structural plane contour map based on the work area structural data in step 1, and interpolate to obtain a structural data volume;

[0035] Step 2.1, draw the contour map of the construction plane;

[0036] Use industrial drawing software, such as Geomap3.6, to draw a construction plane contour map based on the construction data of the work area in step 1, and check whether there are singular values ​​in the construction plane contour map, such as Figure 2 As shown, if there is, correct the hierarchical data, and individual singular values ​​can be directly deleted until there are no singular values;

[0037] Step 2.2, using geological modeling software, such as Petrol software, interpolate the plane contour map of step 2.1 to obtain the structural data volume of the work area;

[0038] Step 3: Substitute the structural data volume obtained in step 2 into the curvature algorithm to obtain the curvature data volume, and extract the corresponding curvature value E from the curvature data volume in step 2 according to the horizontal and vertical coordinates (x, y) of the calibrated well in step 1. f (x, y);

[0039] Step 3.1: Substitute the structural data obtained in step 2 into geological modeling software, such as Petrol software or Geoeast, to obtain the maximum principal curvature K of the structure. max , minimum principal curvature value K min And the directional coefficient d, if K max , K min are all greater than 0, d=1; in other cases, d=-1.

[0040] Step 3.2: According to E f =d*K min *K max Calculate the curvature data volume E f .

[0041] Step 3.3: According to the ordinate and abscissa coordinates (x, y) of the well calibrated in step 1, the curvature data volume E in step 3.2 is obtained. f Extract the corresponding curvature value E f (x, y);

[0042] Step 4: Set the test gas production of the calibration well in step 1 to P gas The curvature value E extracted in step 3 f (x, y) performs correlation analysis to obtain the fitting formula;

[0043] Specifically: the gas production P of the calibration well in step 1 gas and the extracted curvature value E f (x, y) for correlation analysis, such as Figure 3 As shown, the functional relationship between the two is obtained, that is, the fitting formula is established.

[0044] Step 5: Substitute the production standard into the fitting formula of step 4 to obtain the curvature range, and establish the reservoir fracture effectiveness grade evaluation standard according to the curvature range.

[0045] Step 5 is as follows:

[0046] Substitute the production standards of low-yield wells, industrial gas flow wells, and high-yield wells of different oil companies into the fitting formula in step 4 to obtain the E corresponding to different production standards. f Value, E f The range of values ​​is the curvature range. According to the curvature range, a reservoir fracture effectiveness grade evaluation standard is established. Then, according to the reservoir fracture effectiveness grade evaluation standard, a favorable area range is established for evaluation.

[0047] The reservoir fracture effectiveness grade evaluation method of the present invention is applicable to geological oil and gas field exploration, and the scope of the favorable area corresponding to the reservoir fracture effectiveness grade evaluation standard is as follows: Figure 4 As shown, exploration deployment is carried out in combination with the geological background of reservoir formation.

[0048] The technical principle of the reservoir fracture effectiveness grade evaluation method of the present invention is:

[0049] For fractured gas reservoirs, in general, the greater the deformation (reservoir), that is, the greater the curvature, the greater the fracture opening, the better the fracture effectiveness, and the curvature value E f Describe the effectiveness of fractures. Under similar reservoir geological conditions such as source rock, formation thickness or reservoir thickness, and capping conditions, the better the fracture effectiveness, the higher the single well production, and the curvature value E is established. f and test gas output P gas The functional relationship between them, combined with the requirements of low-yield wells, industrial gas flow wells, and high-yield wells for test gas production, can evaluate the fracture effectiveness level.

[0050] The reservoir fracture effectiveness grade evaluation method of the present invention can be applied to geological oil and gas field exploration, achieving the purpose of finding favorable gas-bearing geological targets.

[0051] Example 1

[0052] The reservoir fracture effectiveness rating evaluation method provided in this embodiment has the following steps:

[0053] 1. Prepare structural data of the work area and calibration well data;

[0054] 1) Arrange the structural data of the work area. The structural data is obtained by subtracting the core elevation data from the top surface depth in Table 1 below the drilling geological stratification data.

[0055] Table 1 Drilling geological stratification data table

[0056] well name Horizontal axis Vertical axis Top Depth Heart Repair Poster S210 19271020.9 4278873.2 2871.3 986.2 Q21 19271048.6 4259132.1 2394.6 482.3 Q9 19271051.1 4251412.3 2554.6 697.1 S258 19271066.5 4255321.6 2763.7 897.0 M120 19271071.1 4239256.1 2596.1 785.5 S261 19271078.2 4213181.1 2716.0 865.0 S265 19271082.3 4198005.6 2719.4 882.8 S278 19271089.1 4261077.3 2653.5 823.3 M290 19271111.6 4209763.2 2690.2 880.1 M263 19271125.6 4160455.8 2675.3 886.0 L6 19271090.1 4202462.9 2644.5 854.0 L15 19271091.5 4195941.3 2621.4 786.0

[0057] 2) Arrange the calibration well data of the work area, including the horizontal and vertical coordinates of the well location and the test gas production, as shown in Table 2.

[0058] Table 2 Data table of calibration wells in the work area

[0059]

[0060]

[0061] 2. Draw the structural plane contour map through the structural data, and interpolate to obtain the structural data volume;

[0062] 1) Draw the contour map of the construction plane. Use Geomap3.6 drawing software to draw the contour map of the construction plane, such as Figure 2 Check whether there are singular values. If there are, correct the stratified data and delete the individual singular values ​​directly.

[0063] 2) Obtain the work area structure data volume through interpolation of Petrol software;

[0064] 3. Substitute the interpolated structural data volume into the curvature algorithm to obtain the curvature data volume, and extract the corresponding curvature value E from the curvature data volume according to the horizontal and vertical coordinates (x, y) of the calibration well. f (x, y);

[0065] 1) Import the constructed data volume obtained in the previous step into the Petrol software to obtain the maximum principal curvature K of the construction max , minimum principal curvature value K min and direction factor d(K max , K min All are greater than 0: d = 1; other cases: d = -1), as shown in Table 3 below.

[0066] Table 3 Curvature parameter K of the working area max , K min And detailed table

[0067]

[0068]

[0069] 2) Let E f =d*K min *K max , E f That is, the curvature data volume is calculated, as shown in Table 4.

[0070] Table 4 Work Area E f Parameters table

[0071] well name Horizontal axis Vertical axis <![CDATA[E f ]]> S376 18743843.0 4278874.3 3.1E-14 S169 18730125.9 4259152.7 2.8E-14 S2 19271071.9 4251492.3 4.1E-14 S170 18747703.9 4255305.6 -3.2E-14 S172 18743253.9 4239253.7 -5.1E-14 S173 18745146.0 4213188.3 3.1E-14 S174 18733262.3 4198002.5 6.2E-14 S18 18760386.1 4261076.3 5.3E-14 S190 18733270.3 4209767.6 -2E-14 S191 18744021.2 4160455.3 5E-15 S194 18755213.7 4202464.9 1.8E-14 S196 18758737.9 4195949.3 4E-14

[0072] 3) Extract the corresponding curvature value E according to the vertical and horizontal coordinates (x, y) of the calibration well f (x, y) values, as shown in Table 5 below;

[0073] Table 5 E corresponding to the production of calibration wells f Value table

[0074]

[0075]

[0076] 4. Compare the test gas production of the calibration well (low production, industrial gas flow, high production) with the extracted curvature value E f (x, y) to perform correlation analysis and determine the functional relationship between the two, that is, to establish a fitting formula;

[0077] The E of the calibrated wells in Table 5 f Its test gas production P gas Do correlation analysis, such as Figure 3 As shown, the functional relationship between the two is determined and the fitting formula is established as follows:

[0078] E f =1.1949*P gas -0.7987 (1);

[0079] This functional relationship continues to be used in the working areas of Examples 2, 3, and 4.

[0080] 5. Establish evaluation criteria.

[0081] 1) The production standards of low-yield wells, industrial gas flow wells, and high-yield wells of different oil companies are used (this example uses the gas well classification standard of PetroChina Changqing Oilfield Branch: low-yield wells < 2.0*10 4 m 3 / d, 4.0*10 4 m 3 / d>Industrial air flow well ≥2.0*10 4 m 3 / d, high-yield wells ≥4.0*10 4 m 3 / d) is substituted into the above fitting formula (1), the corresponding curvature range is determined, and the evaluation standard is established, as shown in Table 6. This evaluation standard is continued to be used in the working areas of Examples 2, 3, and 4.

[0082] Table 6 Evaluation table of reservoir fracture effectiveness level in the work area

[0083] Industry Standards Low production wells Industrial air flow well High production wells <![CDATA[Output (10*10 4 m 3 / d)]]> <![CDATA[P gas <2.0]]> <![CDATA[4.0>P gas ≥2.0]]> <![CDATA[P gas ≥4.0]]> <![CDATA[E f *10 14 ]]> <![CDATA[E f <1.59]]> <![CDATA[3.98>AND f ≥1.59]]> <![CDATA[E f ≥3.98 <!-- 6 -->]]> Rating Category III Category II Category I

[0084] 2) Establish the scope of favorable area and establish evaluation criteria according to the evaluation criteria.

[0085] Based on the above work, E f Areas <1.59 are defined as Class III, and areas 3.98>E f The area with a value of ≥1.59 is defined as Class II. f Areas with a value of ≥3.98 are defined as Class I. Among them, Class I is the priority area for implementation in the case of other types of geological conditions.

[0086] Draw E f The plane map and contour map are used to deploy exploration well locations in combination with the geological background of reservoir formation.

[0087] According to the E corresponding to Class I, Class II, and Class III f The value is used to draw the contour map of the reservoir fracture effectiveness level evaluation plane map at intervals, such as Figure 4 As shown, that is: E f The area delineated by ≥3.98 is the Class I favorable area, where the crack effectiveness is the greatest; 3.98>E f The area delineated by Ef ≥1.59 is a Class II favorable area, where the fracture effectiveness is second; the area delineated by Ef <1.59 is a Class III favorable area, where the fracture effectiveness is poor.

[0088] According to the results of the effectiveness grade evaluation, YT3H1 was deployed and implemented in the Class I area. The horizontal section of the well is 759 meters long, 664 meters of fractured reservoirs are drilled, the reservoir drilling rate is 87.5%, and the unimpeded flow rate obtained in the gas test is 629,700 cubic meters per day. The present invention applies the fracture effectiveness grade evaluation to geological oil and gas field exploration, achieving the purpose of finding favorable gas-bearing geological targets.

[0089] Example 2

[0090] The reservoir fracture effectiveness rating evaluation method provided in this embodiment has the following steps:

[0091] 1. Prepare structural data of the work area and calibration well data;

[0092] 1) Arrange the structural data of the work area. The structural data is obtained by subtracting the core elevation data from the top surface depth in Table 7 below through the drilling geological stratification data.

[0093] Table 7 Drilling geological stratification data table

[0094]

[0095]

[0096] 2. Draw the structural plane contour map through the structural data, and interpolate to obtain the structural data volume;

[0097] 1) Draw the contour map of the construction plane. Use Geomap3.6 to draw the contour map of the construction plane and check whether there are singular values. If there are, correct the layered data. Individual singular values ​​can be directly deleted until there are no singular values.

[0098] 2) Obtain the work area structure data volume through interpolation of Petrol software;

[0099] 3. Substitute the interpolated constructed data volume into the curvature algorithm to obtain the curvature data volume;

[0100] 1) Import the constructed data volume obtained in the previous step into Geoeast software to obtain the maximum principal curvature K of the construction max , minimum principal curvature value K min and direction factor d(K max , K min All are greater than 0: d = 1; other cases: d = -1), as shown in Table 8.

[0101] Table 8 Curvature parameter K of working area max , K min and d table

[0102] well name Horizontal axis Vertical axis <![CDATA[K max ]]> <![CDATA[K min ]]> d CH3 19434012.2 4117162.7 -6.70E-07 5.00E-08 -1 CH5 18684473.2 4323157.7 5.20E-06 6.00E-09 1 Q5 19248439.3 3956425.1 2.80E-07 3.00E-07 1 Q9 18754235.3 3987366.2 3.50E-06 3.10E-08 1 B12 18722625.5 4040379.3 2.30E-07 3.00E-07 1 BX6 18713010.3 4020330.0 1.30E-06 2.00E-10 1 BX7 18707672.8 3965259.1 2.00E-07 3.00E-08 1 BX8 18707662.1 3965266.2 2.00E-06 1.00E-08 1 B90 18687223.2 3971801.5 1.50E-06 2.50E-08 1 CH10 18687003.8 4133420.1 -1.00E-06 9.00E-09 -1 CH20 19336365.3 4176716.2 5.00E-06 4.00E-09 1 CH25 19434029.6 4117150.5 5.20E-06 -1.50E-08 -1

[0103] 2) Let E f =d*K min *K max , E f That is, the curvature data volume is calculated, as shown in Table 9.

[0104] Table 9 Work Area E f Parameters table

[0105] well name Horizontal axis Vertical axis <![CDATA[E f *10^14]]> CH3 19434012.2 4117162.7 -3.35 CH5 18684473.2 4323157.7 3.12 Q5 19248439.3 3956425.1 8.40 Q9 18754235.3 3987366.2 10.85 B12 18722625.5 4040379.3 6.90 BX6 18713010.3 4020330.0 0.03 BX7 18707672.8 3965259.1 0.60 BX8 18707662.1 3965266.2 2.00 B90 18687223.2 3971801.5 3.75 CH10 18687003.8 4133420.1 -0.90 CH20 19336365.3 4176716.2 2.01 CH25 19434029.6 4117150.5 -7.80

[0106] 4. Compare the test gas production of the calibration well (low production, industrial gas flow, high production) with the extracted curvature value E f (x, y) do correlation analysis to determine the functional relationship between the two (establish a fitting formula);

[0107] The fitting formula is shown in formula (1) in Example 1.

[0108] 5. Establish evaluation criteria.

[0109] 1) The evaluation criteria are shown in Table 6;

[0110] 2) Establish the scope of favorable area and establish evaluation criteria according to the evaluation criteria.

[0111] In the standard of Table 6, E f Areas <1.59 are defined as Class III, and areas 3.98>E f The area with a value of ≥1.59 is defined as Class II.f Areas with a value of ≥3.98 are defined as Class I. Among them, Class I is the priority area for implementation in the case of other types of geological conditions.

[0112] According to Table 9, the parameter value E is evaluated f The value is determined and exploration well locations are deployed in combination with the geological background of reservoir formation.

[0113] According to the E corresponding to Class I, Class II, and Class III f According to the value standard, the area near the Q9 well in the area of ​​this embodiment is a Class I favorable area with the greatest fracture effectiveness. In addition, the hydrocarbon source in this area is sufficient and the capping conditions are good. Therefore, the QT1H well was deployed and implemented in the Q9 adjacent area. The horizontal section of the well is 1,500 meters long, 1,301 meters of fractured reservoirs are drilled, the reservoir drilling rate is 86.7%, and the unobstructed flow rate obtained by gas testing is 513,000 cubic meters per day. The present invention applies the fracture effectiveness grade evaluation to geological oil and gas field exploration, achieving the purpose of finding favorable gas-bearing geological targets.

[0114] Example 3

[0115] The reservoir fracture effectiveness rating evaluation method provided in this embodiment has the following steps:

[0116] 1. Prepare structural data of the work area and calibration well data;

[0117] 1) Arrange the structural data of the work area. The structural data is obtained by subtracting the core elevation data from the top surface depth in the drilling geological layering data table 10.

[0118] Table 10 Drilling geological stratification data table (partial)

[0119] well name Horizontal axis Vertical axis Top Depth Heart Repair Poster E15 19385468.6 3955957.7 3770.0 1320.9 E21 19309657.8 4182431.7 4054.0 1336.2 E27 19388555.8 4165208.3 3779.6 1412.9 E38 19363690.1 4169604.7 3727.0 1394.7 E49 19471726.1 4159364.5 3852.4 1462.2 E52 19384340.8 4185230.6 2962.5 1262.8 E60 19379006.4 4196221.9 3680.8 1487.7 E61 19382837.1 4202061.4 3798.8 1276.5 E72 19403360.0 4198185.0 3830.0 1219.3 E83 19430442.5 4382930.0 3939.4 1203.2 E24 19354159.9 4166817.4 3870.0 1394.3 E25 18691160.9 4143305.4 3690.8 1274.9

[0120] 2. Draw the structural plane contour map through the structural data, and interpolate to obtain the structural data volume;

[0121] 1) Draw the contour map of the construction plane. Use Geomap3.6 to draw the contour map of the construction plane and check whether there are singular values. If there are, correct the layered data and the individual singular values ​​can be deleted directly.

[0122] 2) Obtain the work area structure data volume through interpolation of Petrol software;

[0123] 3. Substitute the interpolated constructed data volume into the curvature algorithm to obtain the curvature data volume;

[0124] 1) Import the constructed data volume obtained in the previous step into the Petrol software to obtain the maximum principal curvature K of the construction max , minimum principal curvature value K min and direction factor d(K max , K minAll are greater than 0: d = 1; other cases: d = -1), as shown in Table 11.

[0125] Table 11 Curvature parameter K of the working area max , K min and d table

[0126] well name Horizontal axis Vertical axis <![CDATA[K max ]]> <![CDATA[K min ]]> d E18 19314200.7 4241710.4 -3.20E-07 -1.10E-07 -1 E25 19359710.9 4287628.8 2.50E-06 6.90E-09 1 E32 19367597.8 4318364.0 3.90E-07 1.00E-07 1 E46 19254535.9 4158181.3 3.30E-05 -2.00E-09 -1 E50 19264923.1 4146140.0 3.90E-07 -1.60E-07 -1 E53 19305068.9 4198135.9 1.80E-06 2.50E-09 1 E65 19285136.5 4003379.4 5.00E-07 2.10E-07 1 E68 19248612.8 4049723.6 3.20E-06 3.00E-08 1 E71 18753610.5 4110034.7 1.50E-06 2.50E-08 1 E90 19271009.4 4123235.7 6.00E-07 1.90E-08 -1 E102 18688768.0 4235447.3 3.90E-06 5.00E-09 1 E110 18698761.5 4252606.5 3.00E-05 1.50E-09 1

[0127] 2) Let E f =d*K min *K max , E f That is, the curvature data volume is calculated, as shown in Table 12.

[0128] Table 12 Work Area E f Parameters table

[0129] well name Horizontal axis Vertical axis <![CDATA[E f *10^14]]> CH3 19434012.2 4117162.7 -3.52 CH5 18684473.2 4323157.7 1.73 Q5 19248439.3 3956425.1 3.90 Q9 18754235.3 3987366.2 -6.60 B12 18722625.5 4040379.3 -6.24 BX6 18713010.3 4020330.0 0.45 BX7 18707672.8 3965259.1 10.50 BX8 18707662.1 3965266.2 9.60 B90 18687223.2 3971801.5 3.75 CH10 18687003.8 4133420.1 -1.14 CH20 19336365.3 4176716.2 1.95 CH25 19434029.6 4117150.5 4.50

[0130] 4. Compare the test gas production of the calibration well (low production, industrial gas flow, high production) with the extracted curvature value E f (x, y) do correlation analysis to determine the functional relationship between the two (establish a fitting formula);

[0131] The fitting formula is shown in formula (1).

[0132] 5. Establish evaluation criteria.

[0133] 1) The evaluation criteria are shown in Table 6;

[0134] 2) Establish the scope of favorable area according to the evaluation criteria and conduct evaluation.

[0135] In the standard of Table 6, E f Areas <1.59 are defined as Class III, and areas 3.98>E f The area with a value of ≥1.59 is defined as Class II. f Areas with a value of ≥3.98 are defined as Class I. Among them, Class I is the priority area for implementation in the case of other types of geological conditions.

[0136] According to Table 12, the parameter value E is evaluated f The value is determined and exploration well locations are deployed in combination with the geological background of reservoir formation.

[0137] According to the E corresponding to Class I, Class II, and Class III fAccording to the value standard, the area near the BX7 well in the area of ​​this embodiment is a Class I favorable area with the greatest fracture effectiveness. In addition, the hydrocarbon source in this area is sufficient and the capping conditions are good. Therefore, the BT1H well was deployed and implemented in the BX7 area. The horizontal section of the well is 1,500 meters long, 1,289 meters of fractured reservoirs are drilled, the reservoir drilling rate is 85.9%, and the unimpeded flow rate obtained by gas testing is 423,000 cubic meters per day. The present invention applies the fracture effectiveness grade evaluation to geological oil and gas field exploration, achieving the purpose of finding favorable gas-bearing geological targets.

[0138] Example 4

[0139] The reservoir fracture effectiveness rating evaluation method provided in this embodiment has the following steps:

[0140] 1. Prepare structural data of the work area and calibration well data;

[0141] 1) Arrange the structural data of the work area. The structural data is obtained by subtracting the core elevation data from the top surface depth in the drilling geological layering data table 13.

[0142] Table 13 Drilling geological stratification data table

[0143]

[0144]

[0145] 2. Draw the structural plane contour map through the structural data, and interpolate to obtain the structural data volume;

[0146] 1) Draw the contour map of the construction plane. Use industrial drawing software to draw the contour map of the construction plane and check whether there are singular values. If there are, correct the individual singular values ​​of the layered data and delete them directly.

[0147] 2) Obtain the work area structure data volume through interpolation of Petrol software;

[0148] 3. Substitute the interpolated constructed data volume into the curvature algorithm to obtain the curvature data volume;

[0149] 1) Import the constructed data volume obtained in the previous step into the Petrol software to obtain the maximum principal curvature K of the construction max , minimum principal curvature value K min and direction factor d(K max , K min All are greater than 0: d = 1; other cases: d = -1), as shown in Table 11.

[0150] Table 14 Curvature parameter K of working area max , K min and d table

[0151] well name Horizontal axis Vertical axis <![CDATA[K max ]]> <![CDATA[K min ]]> d QT5 18694202.1 4072344.5 5.50E-07 3.00E-08 1 J85 18734785.2 4018606.8 3.20E-06 5.90E-09 1 J87 18734794.7 4018606.8 4.30E-07 2.00E-07 1 J84 18753019.2 3964947.4 -1.60E-05 -1.00E-09 -1 J83 18753629.7 4048098.5 5.10E-07 -3.70E-07 -1 J81 18768564.3 3958577.6 5.30E-07 1.80E-09 1 J89 187685672.1 3958510.6 2.00E-07 3.20E-07 1 J78 18715185.1 3955961.6 2.80E-06 3.00E-08 1 J75 18719016.0 3989426.4 3.50E-06 5.10E-08 1 J72 19236268.9 4027214.1 -1.00E-06 -3.00E-08 -1 J69 18691718.0 4010474.9 -4.50E-06 4.00E-09 -1 J58 19233624.7 3975657.3 2.60E-05 -2.00E-09 -1

[0152] 2) Let E f =d*K min *K max , E f That is, the curvature data volume is calculated, as shown in Table 15.

[0153] Table 15 Work Area E f Parameters table

[0154] well name Horizontal axis Vertical axis <![CDATA[E f *10^14]]> QT5 18694202.1 4072344.5 1.65 J85 18734785.2 4018606.8 1.89 J87 18734794.7 4018606.8 8.60 J84 18753019.2 3964947.4 -1.60 J83 18753629.7 4048098.5 -7.77 J81 18768564.3 3958577.6 0.10 J89 187685672.1 3958510.6 6.40 J78 18715185.1 3955961.6 8.40 J75 18719016.0 3989426.4 2.80 J72 19236268.9 4027214.1 -3.00 J69 18691718.0 4010474.9 -1.80 J58 19233624.7 3975657.3 -5.20

[0155] 4. Compare the test gas production of the calibration well (low production, industrial gas flow, high production) with the extracted curvature value E f (x, y) do correlation analysis to determine the functional relationship between the two (establish a fitting formula);

[0156] The fitting formula is shown in formula (1).

[0157] 5. Establish evaluation criteria.

[0158] 1) The evaluation criteria are shown in Table 6;

[0159] 2) Establish the scope of favorable area for evaluation according to the evaluation criteria.

[0160] In the standard of Table 6, E f Areas <1.59 are defined as Class III, and areas 3.98>E f The area with a value of ≥1.59 is defined as Class II. f Areas with a value of ≥3.98 are defined as Class I. Among them, Class I is the priority area for implementation in the case of other types of geological conditions.

[0161] Evaluate the parameter value E according to Table 15 f The value is determined and exploration well locations are deployed in combination with the geological background of reservoir formation.

[0162] According to the E corresponding to Class I, Class II, and Class III f According to the value standard, the area near the J78 well in the area of ​​this embodiment is a Class I favorable area with the greatest fracture effectiveness. In addition, the hydrocarbon source in this area is sufficient and the sealing conditions are good. Therefore, the QT1H well was deployed and implemented in the J78 area. The horizontal section of the well is 1,500 meters long, and the fractured reservoir is drilled for 1,159 meters, with a reservoir drilling rate of 77.3%. The unobstructed flow rate obtained by gas testing is 112,000 cubic meters per day.

Claims

1. A reservoir fracture effectiveness rating evaluation method, characterized in that: Follow the steps below to implement it: Step 1: Prepare the structural data of the work area and the calibration well data; Step 2: Draw a structural plane contour map based on the work area structural data in step 1, and interpolate to obtain a structural data volume; Step 3, substituting the structural data volume obtained in step 2 into the curvature algorithm to obtain a curvature data volume, and extracting the corresponding curvature value from the curvature data volume in step 2 according to the calibration well data in step 1; Step 4: Perform correlation analysis on the calibration well data in step 1 and the curvature value extracted in step 3 to obtain a fitting formula; Step 5: Substitute the production standard into the fitting formula of step 4 to obtain the curvature range, and establish the reservoir fracture effectiveness grade evaluation standard according to the curvature range for evaluation.

2. The reservoir fracture effectiveness grade evaluation method according to claim 1, characterized in that: In step 1, the structural data of the work area are obtained from drilling analysis data or seismic interpretation data, including horizontal and vertical coordinates, structural values, and structural data of the work area.

3. The reservoir fracture effectiveness grade evaluation method according to claim 2, characterized in that: In step 1, the calibration well data includes the horizontal and vertical coordinates of the well location, the test gas production P gas The test gas production includes low-yield, industrial gas flow and high-yield test gas production.

4. The reservoir fracture effectiveness grade evaluation method according to claim 3, characterized in that: The step 2 is specifically implemented according to the following steps: Step 2.1, draw a structural plane contour map based on the work area structural data in step 1; Step 2.2: Use geological modeling software to interpolate the plane contour map of step 2.1 to obtain the structural data volume of the work area.

5. The reservoir fracture effectiveness grade evaluation method according to claim 4, characterized in that: The step 2.1 is specifically as follows: using industrial drawing software, drawing a structural plane contour map according to the work area structural data of step 1, checking whether there are singular values ​​in the structural plane contour map, and if so, correcting the layered data until there are no singular values.

6. The reservoir fracture effectiveness grade evaluation method according to claim 4, characterized in that: The step 3 is specifically implemented according to the following steps: Step 3.1: Substitute the structural data obtained in step 2 into the geological modeling software to obtain the maximum principal curvature K of the structure. max , minimum principal curvature value K min And the directional coefficient d, if K max , K min If both are greater than 0, then d = 1; otherwise d = -1; Step 3.2: According to E f =d*K min *K max Calculate the curvature data volume E f ; Step 3.3: According to the ordinate and abscissa coordinates (x, y) of the well calibrated in step 1, the curvature data volume E in step 3.2 is obtained. f Extract the corresponding curvature value E f (x, y).

7. The reservoir fracture effectiveness grade evaluation method according to claim 6, characterized in that: The step 4 is specifically as follows: the test gas production P of the well calibrated in step 1 is gas The curvature value E extracted in step 3 f (x, y) is subjected to correlation analysis to obtain the functional relationship between the two, that is, the fitting formula.

8. The reservoir fracture effectiveness grade evaluation method according to claim 7, characterized in that: The step 5 is specifically as follows: Substitute the production standards of low-yield wells, industrial gas flow wells, and high-yield wells of different oil companies into the fitting formula of step 4 to obtain E corresponding to different production standards f Value, E f The range of values ​​is the curvature range. According to the curvature range, a reservoir fracture effectiveness grade evaluation standard is established. Then, according to the reservoir fracture effectiveness grade evaluation standard, a favorable area range is established for evaluation.

Images