Model well-based electric imaging logging fracture width calculation method

Abstract

The present invention relates to the technical field of fractured reservoir research. Disclosed is a model well-based electric imaging logging fracture width calculation method, comprising: by means of a fracture width gauge, obtaining the fracture width of a model well; preparing in batches a slurry in the model well into slurry having different resistivities, so as to obtain slurry resistivity data and obtain the electric imaging logging fracture widths of the model well corresponding to the different slurry of different batches; according to the different slurry resistivity data of the different batches and the corresponding electric imaging logging fracture widths, obtaining stratum contrast ratios of different batches; according to the fracture widths obtained by the fracture width gauge, the electric imaging logging fracture widths and the stratum contrast ratios, obtaining a fracture width correction coefficient of any well; and, according to the fracture width correction coefficient of any well, obtaining a corrected electric imaging logging fracture width of an actual well. The present invention can correct and accurately calculate the electric imaging logging fracture widths of actual wells, thus providing reliable guidance for evaluating the effectiveness and productivity of fractured reservoirs.

Description

A method for calculating fracture width using electrical imaging logging based on a model well Technical Field

[0001] The invention relates to the technical field of fracture reservoir research, and in particular to a method for calculating fracture width based on electrical imaging logging of a model well. Background Art

[0002] Fracture width refers to the average opening of a fracture near the wellbore wall. To date, there is no reliable method for directly measuring fracture width in rock under formation conditions in the wellbore. Using electrical imaging logging to identify fractures and calculate their width is a common fracture evaluation method. Currently, the most widely used electrical imaging logging instrument worldwide is Schlumberger's FMI. Fracture width is calculated by integrating the area of ​​abnormal current caused by the presence of fractures in a high-resistance background. The fracture width calculation formula is:

[0003] W=cAR m b R xo 1-b (1)

[0004] Where W is the crack width, mm; b and c are dimensionless instrument parameters; Rm is the mud resistivity, Ω·m; Rxo is the flushing zone resistivity, Ω·m; A is the abnormal current area caused by the crack, uA.mm / V.

[0005] Different oil and gas service companies use electrical imaging logging instruments with different instrument coefficients, b and c. Schlumberger, through numerical simulation, determined the instrument coefficients for its FMI electrical imaging logging instrument: b = 0.863, c = 0.48. Rm is the mud resistivity, which can be directly measured or obtained by mud sampling and analysis. Schlumberger maintains commercial confidentiality regarding the determination of the A and Rxo parameters in Equation (1) and has not disclosed them to the public. Currently, Schlumberger has developed the Techlog logging software platform, which can process data from its own FMI electrical imaging logging instrument according to Equation (1) and determine fracture widths. However, there is no accurate data to verify the accuracy of Equation (1) and the reliability of its calculation results.

[0006] Fractured reservoirs develop dual porosity spaces, namely pores and fractures. Fractures are important storage spaces and seepage channels in buried-hill fractured reservoirs, and play a key role in fluid migration and production. Electrical imaging logging is widely used both domestically and internationally for qualitative evaluation of fractures and calculation of fracture parameters (such as fracture width). The wider a fracture is, the stronger its seepage capacity is. Accurate calculation of fracture width is of great significance for evaluating the effectiveness of fractured reservoirs and understanding their productivity.

[0007] Therefore, a method that can accurately calculate the width of electrical imaging logging fractures is urgently needed to provide reliable guidance for the effectiveness and productivity evaluation of fractured reservoirs. SUMMARY OF THE INVENTION

[0008] The present invention provides a method for calculating the width of electrical imaging logging fractures based on a model well. The method uses a fracture width meter and an electrical imaging logging instrument to measure the fracture width of the model well, compares the variation relationship between the fracture width obtained by the fracture width meter and the fracture width calculated by electrical imaging logging, and combines formation contrast to obtain a fracture width correction coefficient for any well. The method can correct and accurately calculate the electrical imaging logging fracture width of an actual well, providing reliable guidance for the effectiveness and productivity evaluation of fractured reservoirs. Technical issues

[0009] The present invention provides a method for calculating the width of electrical imaging logging fractures based on a model well, which is used to solve the defect that the accuracy of the width of electrical imaging logging fractures obtained by the prior art needs to be improved. Technical Solutions

[0010] The present invention provides a method for calculating fracture width based on electrical imaging logging of a model well, comprising:

[0011] The fracture width of the model well is obtained by using a fracture width meter;

[0012] The mud in the model well is prepared into muds with different resistivities in batches, mud resistivity data is obtained, and the electrical imaging logging fracture width of the model well corresponding to the muds of different batches is obtained;

[0013] According to the different mud resistivity data and the corresponding electrical imaging logging fracture width, the formation contrast of different times is obtained;

[0014] According to the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a correction coefficient of the fracture width of the electrical imaging logging of the model well is obtained;

[0015] The fracture width correction coefficient of any well is obtained based on the electrical imaging logging fracture width correction coefficient and formation contrast of the model well;

[0016] According to the fracture width correction coefficient of any well, the corrected electrical imaging logging fracture width of the actual well is obtained.

[0017] In one embodiment, the fracture width of the model well obtained by the fracture width meter is a physical fracture width measurement value of the model well, and the fracture width of the model well can be directly obtained from this data. The fracture width of the model well obtained by electrical imaging logging can be obtained by processing the fracture width of the model well obtained by electrical imaging logging using various techniques known in the art, such as FMI electrical imaging technology, computer vision-based automatic measurement technology, such as a fracture width detection algorithm based on OpenCV, etc.

[0018] According to a method for calculating the fracture width of an electrical imaging logging model well provided by the present invention, the fracture width correction coefficient of the electrical imaging logging model well is obtained based on the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging model well, including:

[0019] According to the variation law between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a fitting relationship is obtained;

[0020] According to the fitting relationship, the electrical imaging logging fracture width correction coefficient of the model well is obtained.

[0021] According to a method for calculating fracture width based on electrical imaging logging of a model well provided by the present invention, a fitting relationship is obtained based on the variation pattern between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, including:

[0022] According to the variation law between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, the first fitting relationship is obtained;

[0023] Performing a first optimization process on the first fitting relationship to obtain a second fitting relationship;

[0024] Among them, the first fitting relationship is:

[0025] Log10(FVA_C)=s i Log10(FVA_Ri)+n i0 ,

[0026] In the first fitting equation, FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the model well at the time i of electrical imaging logging, s i Indicates the slope of the first fitting relationship for the i-th time, n i0 represents the intercept of the first fitting relationship of the i-th time;

[0027] And, wherein the second fitting relationship is:

[0028] Log10(FVA_C)=1·Log10(FVA_Ri)+ni ,

[0029] In the second fitting equation, FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the model well at the time i of electrical imaging logging, s i It represents the slope of the first fitting relationship of the i-th time, and its value is 1, n i Represents the intercept of the second fitting relationship of the i-th order.

[0030] According to a method for calculating the fracture width of electrical imaging logging based on a model well provided by the present invention, the electrical imaging logging fracture width correction coefficient of the model well is obtained according to the fitting relationship, specifically:

[0031] Performing a second optimization process on the second fitting relationship to obtain a relationship expression between the fracture width obtained by the fracture width meter and the fracture width obtained by the electrical imaging logging, thereby obtaining a correction coefficient for the fracture width of the electrical imaging logging of the model well;

[0032] The relational expression is:

[0033] FVA_C=10 ni ·FVA_Ri,

[0034] In the relational expression, FVA_C represents the fracture width of the model well obtained by the fracture width gauge, and FVA_R i represents the width of the electrical imaging logging fracture of the model well i, 10 ni represents the correction coefficient of the fracture width of the i-th electrical imaging logging in the model well, d i =10 ni , i=1,2,...,n.

[0035] According to a method for calculating fracture width based on electrical imaging logging of a model well provided by the present invention, the fracture width correction coefficient of any well is expressed as follows:

[0036] d=a+f·(Rm / Rxo),

[0037] In the expression of the fracture width correction coefficient of any well, d represents the fracture width correction coefficient of the model well, Rm represents the mud resistivity of the model well, Rxo represents the resistivity of the flushing zone of the model well, Rm / Rxo represents the formation contrast of the model well, a and f represent the parameters of the fracture width correction coefficient of any well. The expression of the fracture width correction coefficient of any well is used to fit the data points of the model well to obtain a and f, where d is obtained by the fracture width correction coefficient d of n electrical imaging logging. i Rm and Rxo are composed of Rm during n times of electrical imaging logging. i 、Rxo iComposition, i=1,2,...,n.

[0038] According to a method for calculating electrical imaging logging fracture width based on a model well provided by the present invention, the method obtains the corrected electrical imaging logging fracture width of an actual well according to the fracture width correction coefficient of an arbitrary well, comprising:

[0039] According to the fracture width correction coefficient of any well, the fracture width correction coefficient of the actual well is obtained;

[0040] According to the actual well fracture width correction coefficient, the actual well fracture width correction electrical imaging logging is obtained;

[0041] The expression of the fracture width correction coefficient of the actual well is:

[0042] d'=a+f·(Rm' / Rxo'),

[0043] In the expression of the fracture width correction coefficient of the actual well, d' represents the fracture width correction coefficient of the actual well, a and f represent the parameters of the fracture width correction coefficient of any well, Rm' represents the mud resistivity of the actual well, and Rxo' represents the resistivity of the flushing zone of the actual well.

[0044] According to a method for calculating electrical imaging logging fracture width based on a model well provided by the present invention, obtaining the corrected electrical imaging logging fracture width of the actual well according to the fracture width correction coefficient of the actual well includes:

[0045] Obtain the mud resistivity at the actual well fracture;

[0046] According to the actual well electrical imaging logging fracture width, the resistivity of the flushing zone at the actual well electrical imaging logging fracture is obtained;

[0047] According to the mud resistivity of the actual well, the resistivity of the flushing zone, the electrical imaging logging fracture width and the fracture width correction coefficient, the corrected electrical imaging logging fracture width of the actual well is obtained;

[0048] The expression of the corrected electrical imaging logging fracture width of the actual well is:

[0049] W = d'·W0,

[0050] In the expression of the corrected electrical imaging logging fracture width of the actual well, W represents the corrected electrical imaging logging fracture width of the actual well, W0 represents the electrical imaging logging fracture width of the actual well, and d' represents the fracture width correction coefficient of the actual well.

[0051] The present invention also provides a model well-based electrical imaging logging fracture width calculation system, comprising:

[0052] The first data acquisition module is used to obtain the fracture width of the model well by using a fracture width meter;

[0053] The second data acquisition module is used to: prepare the mud in the model well into muds with different resistivities in batches, obtain mud resistivity data and obtain the electrical imaging logging fracture width of the model well corresponding to the muds of different batches;

[0054] The first data processing module is used to obtain different formation contrasts based on different mud resistivity data and corresponding electrical imaging logging fracture widths;

[0055] The second data processing module is used to obtain a correction coefficient for the electrical imaging logging fracture width of the model well based on the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width of the model well;

[0056] The third data processing module is used to obtain the fracture width correction coefficient of any well based on the electrical imaging logging fracture width correction coefficient and formation contrast of the model well;

[0057] The correction module is used to obtain the corrected electrical imaging logging fracture width of an actual well according to the fracture width correction coefficient of any well.

[0058] The present invention also provides an electronic device comprising a processor and a memory storing a computer program, wherein the processor implements any of the above-mentioned methods for calculating fracture width based on electrical imaging logging of a model well when executing the computer program.

[0059] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the computer program implements any of the above-mentioned methods for calculating fracture width based on electrical imaging logging of a model well.

[0060] The present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute any of the above-mentioned methods for calculating fracture width based on electrical imaging logging of model wells. Beneficial effects

[0061] The present invention provides a method for calculating the width of electrical imaging logging fractures based on a model well. The method uses a fracture width meter and an electrical imaging logging instrument to measure the fracture width of the model well, compares the variation relationship between the fracture width obtained by the fracture width meter and the fracture width calculated by electrical imaging logging, and combines formation contrast to obtain a fracture width correction coefficient for any well. The method can correct and accurately calculate the electrical imaging logging fracture width of an actual well, providing reliable guidance for the effectiveness and productivity evaluation of fractured reservoirs. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] In order to more clearly illustrate the technical solutions in the present invention or the prior art, a brief introduction will be given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0063] FIG1 is a flow chart of a method for calculating fracture width based on electrical imaging logging of a model well provided by the present invention.

[0064] Figure 2(a) shows the rock physics model (1 m × 1 m) and the fracture, and Figure 2(b) shows the fracture width measured by the fracture width meter.

[0065] FIG3 shows an electrical imaging logging image obtained by processing electrical imaging logging data through software.

[0066] FIG4 shows the fracture trace width of the electrical imaging image of the model well and the fracture width calculation results.

[0067] FIG5 shows the crack width measurement results of the crack width gauge before and after two times.

[0068] FIG6 shows a cross-plot of the fracture width obtained by electrical imaging logging and the fracture width obtained by the fracture width meter.

[0069] FIG7 is a schematic diagram showing the background resistivity of the tight section of the shallow lateral resistivity in the model well when there is no fracture influence.

[0070] FIG8 shows the electrical imaging logging fracture parameters and the corrected electrical imaging logging fracture width of Well B-2Sa.

[0071] FIG9 is a schematic structural diagram of a model well-based electrical imaging logging fracture width calculation system provided by the present invention.

[0072] FIG10 is a schematic structural diagram of an electronic device provided by the present invention. Modes for Carrying Out the Invention

[0073] In order to make the purpose, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below in conjunction with the drawings in the present invention. Obviously, the embodiments described are part of the embodiments of the present invention, not all of the embodiments, and they should not be understood as limitations on the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention. In the description of the present invention, it should be understood that the terms used are only for descriptive purposes and cannot be understood as indicating or implying relative importance.

[0074] The following describes the method for calculating fracture width based on electrical imaging logging of a model well provided by the present invention in conjunction with FIG. 1 to FIG. 10 .

[0075] FIG1 is a flow chart of a method for calculating fracture width using electrical imaging logging based on a model well provided by the present invention. Referring to FIG1 , a method for calculating fracture width using electrical imaging logging based on a model well provided by the present invention may include:

[0076] Step S110, obtaining the fracture width of the model well by using a fracture width meter;

[0077] Step S120: preparing the mud in the model well into muds with different resistivities in batches, obtaining mud resistivity data and obtaining electrical imaging logging fracture widths of the model well corresponding to the muds of different batches;

[0078] Step S130, obtaining different formation contrasts based on different mud resistivity data and corresponding electrical imaging logging fracture widths;

[0079] Step S140, obtaining a correction coefficient for the electrical imaging logging fracture width of the model well based on the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width of the model well;

[0080] Step S150, obtaining a fracture width correction coefficient of any well based on the electrical imaging logging fracture width correction coefficient and formation contrast of the model well;

[0081] Step S160: Obtain the corrected electrical imaging logging fracture width of the actual well according to the fracture width correction coefficient of the arbitrary well.

[0082] It should be noted that, in step S110, the fracture width of the rock physics model well can be measured using a fracture width gauge. The measurement accuracy of the fracture width gauge is 1 μm. Multiple data points can be measured for one fracture and the average of the multiple fracture width values ​​can be calculated as the average width of the fracture.

[0083] It should be noted that step S120 can be performed by adding a certain amount of industrial salt to the mud in the model well five times to prepare five different muds with different resistivities. Electrical imaging logging data and mud resistivity Rm data for the model well after each addition of industrial salt are recorded. The electrical imaging logging data can be acquired using existing electrical imaging logging equipment, such as FMI, ERMI, MCI, STAR, STAR-II, and EMI. Software is then used to process the electrical imaging logging data to generate electrical imaging logging images. The visual fracture widths of the five electrical imaging logging images corresponding to the five different mud resistivities are then compared to determine the logging response of formation contrast Rm / Rxo (i.e., the ratio of mud resistivity Rm to flush zone resistivity Rxo). Furthermore, the electrical imaging logging data can be processed using existing Techlog software to calculate the electrical imaging logging fracture width.

[0084] It should be noted that the fracture width of the model well can be calculated by five times of electrical imaging logging. The fracture uniformity can be judged according to the color brightness of the fracture trajectory of the model well image and the difference in the fracture width measurement results of the two previous fracture width gauges. The data points with reliable quality obtained by the fracture width gauge measurement corresponding to relatively uniform fractures are preferably used for subsequent data modeling to improve the accuracy of the output results of the present invention.

[0085] In one embodiment, step S150 may include:

[0086] Step S1501: Obtain a fitting relationship based on the variation pattern between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well;

[0087] Step S1502: Obtain the electrical imaging logging fracture width correction coefficient of the model well according to the fitting relationship.

[0088] In one embodiment, step S1501 may include:

[0089] According to the variation law between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, the first fitting relationship is obtained;

[0090] Performing a first optimization process on the first fitting relationship to obtain a second fitting relationship;

[0091] Among them, the first fitting relationship is:

[0092] Log10(FVA_C)=s i Log10(FVA_Ri)+n i0 (2)

[0093] In formula (2), FVA_C represents the fracture width of the model well obtained by the fracture width meter, mm, and FVA_R i Indicates the width of the electrical imaging logging fracture of the model well i (i=1,2,...,n), mm, s i It represents the slope of the first fitting relationship of the i-th (i=1,2,...,n) order, dimensionless, n i0 The intercept of the first fitting equation for the i-th order (i=1,2,...,n) is dimensionless. Log10(FVA_C) represents the base-10 logarithm of FVA_C, and the others are similar.

[0094] It should be noted that step S1501 can obtain the fitting relationship (2) by comparing the change pattern between the fracture width of the model well 5 times of electrical imaging logging and the reliable fracture width measured by the fracture width meter, and then through mathematical operation, let s i =1, and the second fitting relationship (3) is obtained.

[0095] Among them, the second fitting relationship is:

[0096] Log10(FVA_C)=1·Log10(FVA_Ri)+n i (3)

[0097] In formula (3), FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the model well at the time i of electrical imaging logging, s i It represents the slope of the first fitting relationship of the i-th (i=1,2,...,n) time, and its value is 1, n i It represents the intercept of the second fitting relationship of the i-th (i=1,2,...,n) order, dimensionless.

[0098] Furthermore, in step S1502, the second fitting relationship can be optimized to obtain the relationship expression (4) between the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width, thereby obtaining the electrical imaging logging fracture width correction coefficient of the model well;

[0099] The relational expression is:

[0100] FVA_C=10 ni ·FVA_Ri (4)

[0101] In the relational expression, FVA_C represents the fracture width of the model well obtained by the fracture width gauge, and FVA_R i represents the width of the electrical imaging logging fracture of the model well i, 10 nirepresents the correction coefficient of the fracture width of the i-th electrical imaging logging in the model well, d i =10 ni , i=1,2,...,n.

[0102] According to the relational expression (4), the electrical imaging logging fracture width correction coefficient d can be obtained when the i-th electrical imaging logging fracture width of the model well is corrected to the accurate size. i The value of d i =10 ni , i=1,2,...,n.

[0103] Furthermore, based on the relationship between the fracture width correction coefficient d of the electrical imaging logging of the model well and the formation contrast (Rm / Rxo), the calculation model of the fracture width correction coefficient of any well can be obtained. The expression (5) of the fracture width correction coefficient of any well is:

[0104] d=a+f·(Rm / Rxo) (5)

[0105] In formula (5), d represents the fracture width correction coefficient of the model well, which is dimensionless and is taken as d=10 ni , Rm represents the mud resistivity of the model well, Ω.m, Rxo represents the flushing zone resistivity of the model well, Ω.m, Rm / Rxo represents the formation contrast of the model well, dimensionless, a and f represent the parameters of the fracture width correction coefficient for any well. According to Equation (5), the values ​​of a and f can be obtained. Among them, the mud resistivity Rm can be directly obtained by measurement, and the flushing zone resistivity Rxo is the background resistivity when there is no fracture influence. It can be directly obtained by reading the value of the shallow lateral resistivity RS in the tight section without fractures in the model well.

[0106] In formula (5), d and Rm / Rxo are obtained based on the model well, and d is the correction coefficient d of the fracture width of n electrical imaging logging. i Rm and Rxo are composed of Rm during n times of electrical imaging logging. i 、Rxo i The composition, i = 1, 2, ..., n, Rm and Rxo are obtained by reading the values, and then the coefficients a and f are obtained by fitting according to formula (5), which has the property of being applicable to any well. In other words, expression (5) in the model well is applicable to any well.

[0107] In one embodiment, step S160 may include:

[0108] Step S1601: Obtain the fracture width correction coefficient of the actual well based on the fracture width correction coefficient of any well;

[0109] Step S1602: Obtain the corrected electrical imaging logging fracture width of the actual well according to the fracture width correction coefficient of the actual well.

[0110] Specifically, step S1601 can substitute the parameters a and f obtained from equation (5) into the calculation model of the fracture width correction coefficient of the actual well, and its expression (6) is:

[0111] d'=a+f·(Rm' / Rxo') (6)

[0112] In formula (6), d' represents the fracture width correction coefficient of the actual well, which is dimensionless, a and f represent the parameters of the fracture width correction coefficient of any well, Rm' represents the mud resistivity of the actual well, Ω.m, and Rxo' represents the resistivity of the flushing zone of the actual well, Ω.m.

[0113] In one embodiment, step S1602 may include:

[0114] Obtain the mud resistivity at the actual well fracture;

[0115] According to the actual well electrical imaging logging fracture width, the resistivity of the flushing zone at the actual well electrical imaging logging fracture is obtained;

[0116] According to the mud resistivity, flushing zone resistivity, electrical imaging logging fracture width and fracture width correction coefficient of the actual well, the corrected electrical imaging logging fracture width of the actual well is obtained.

[0117] It should be noted that the mud resistivity can be directly measured by existing measuring instruments, while the resistivity of the flushing zone can be obtained by calculating W1 and W2 of each fracture based on the electrical imaging logging fracture width of the actual well using formula (1).

[0118] Specifically, let Rm=1,

[0119] When b=0 and c=1, the calculated crack width W1 is:

[0120] W1=A×Rxo (7)

[0121] When b=1 and c=1, the calculated crack width W2 is:

[0122] W2=A (8)

[0123] Then, the resistivity Rxo of the flushing zone at each fracture position in the electrical imaging logging of the actual well is:

[0124] Rxo=W1 / W2 (9)

[0125] Then, the corrected electrical imaging logging fracture width of the actual well is obtained by combining the mud resistivity, flushing zone resistivity, electrical imaging logging fracture width, and fracture width correction coefficient of the actual well. The expression (10) for the corrected electrical imaging logging fracture width of the actual well is:

[0126] W = d'·W0 (10)

[0127] In formula (10), W represents the corrected electrical imaging logging fracture width of the actual well obtained by the present invention, W0 represents the electrical imaging logging fracture width obtained by formula (1), and d' represents the fracture width correction coefficient of the actual well, which is dimensionless.

[0128] The method for calculating fracture width based on electrical imaging logging of a model well provided by the present invention will be further described below through specific embodiments.

[0129] Example 1: Take the B-2Sa well in a granite buried-hill fractured gas reservoir as an example.

[0130] 1) Fracture widths in a granite rock physics model well were measured using a fracture width gauge with a measurement accuracy of 1 μm. A two-dimensional plan view of the rock physics model well (1 m × 1 m × 1 m) is shown in Figure 2(a), and a schematic diagram of the fracture widths measured by the fracture width gauge is shown in Figure 2(b). The fracture widths at positions 1 and 2 in Figure 2(a) are measured to be 0.2 mm and 0.23 mm, respectively, in Figure 2(b). Multiple data points can be measured for a single fracture, and the average of these fracture widths is calculated as the average width of the fracture. The average widths of individual fractures are shown in Table 1. Table 1 Notes: Blank spaces in Table 1 indicate unmeasured data. The last column, "Fracture Width Quality," has a value of 1 for good quality and 0 for poor quality.

[0131]

[0132] Figure 2(a) shows a 360-degree rock physics model and fractures around the granite wellbore wall. Positions 1 and 2 in the figure indicate fracture width measurement points. Figure 2(b) shows the fracture width measurement results at positions 1 and 2 using a fracture width gauge. The fracture widths at positions 1 and 2 in Figure 2(a) are 0.2 mm and 0.23 mm, respectively, as measured in Figure 2(b).

[0133] 2) Five types of mud with different resistivities were prepared by adding a certain amount of industrial salt to the mud in the model well five times. The electrical imaging logging data and mud resistivity Rm data of the model well were recorded after each addition of industrial salt. The raw data were processed by software to obtain electrical imaging logging images, as shown in Figure 3. Based on this, the electrical imaging logging fracture width of the model well can be obtained.

[0134] In Figure 3, the first track represents the formation depth of the rock physics model; the second track represents the rock block number; the third track represents the first electrical imaging logging image, with the corresponding mud resistivity Rm = 1.167 Ω·m; the fourth track represents the second electrical imaging logging image, with the corresponding Rm = 0.932 Ω·m; the fifth track represents the third electrical imaging logging image, with the corresponding Rm = 0.570 Ω·m; the sixth track represents the fourth electrical imaging logging image, with the corresponding Rm = 0.482 Ω·m; and the seventh track represents the fifth electrical imaging logging image, with the corresponding Rm = 0.282 Ω·m.

[0135] In the rock physics model well in Figure 2 (a), each rock block (1m×1m×1m) is dense granite. After drilling the wellbore, cracks were created by artificial fracture creation. Therefore, the background resistivity Rxo of each rock block without the influence of cracks can be considered to be the same. Comparing the visual width of the cracks in the five electrical imaging logging images corresponding to five different mud resistivities Rm in Figure (3), it is found that the smaller the formation contrast (the smaller the difference between Rm and Rxo), the larger the Rm / Rxo value, and the wider the cracks appear on the electrical imaging logging image.

[0136] 3) Use software to calculate the fracture widths of the model well using five electrical imaging logs according to equation (1). The fracture uniformity is determined by the color brightness of the fracture trace widths in the model well image (as shown in Figure 4) and the difference between the fracture width measurements taken by the fracture width gauge (as shown in Figure 5).

[0137] In Figure 4, the first track represents the formation depth of the rock physics model; the second track represents the electrical imaging log image; the third track represents the fracture trace width map on the electrical imaging log image; and the fourth track represents the fracture width results. The calculated fracture width represents the electrical imaging log fracture width calculated using software, while the model-measured fracture width represents the fracture width obtained using a fracture width gauge. Based on the changes in brightness along the same fracture trace in the second track and the changes in fracture width along the same fracture trace in the third track in Figure 4, combined with the difference between the calculated and model-measured fracture widths in track 4 (whether the data points overlap), it can be concluded that the fractures in frame 2 are more uniform than those in frame 1, and the fracture width measurements are more reliable.

[0138] In Figure 5, the first track represents the stratigraphic depth of the rock physics model; the second track represents the rock block number; the third track represents the electrical imaging log image; the fourth track, with diamonds, represents the first fracture width measurement of the rock physics model, while the rods represent the fracture width calculated using the electrical imaging log; the fifth track, with circles, represents the second fracture width measurement, while the rods represent the fracture width calculated using the electrical imaging log. The uniformity of the fracture color in the third track's electrical imaging image, combined with the fracture width measurements from the fourth and fifth tracks and the fracture width calculated using the logging, indicates that the fractures between blocks 1 and 2 are more uniform, making the fracture width measurements more reliable.

[0139] According to the above method, the fracture widths FVA_Ri (i = 1, 2, ..., 5) calculated from five electrical imaging logs corresponding to relatively uniform fractures and nine reliable data points measured by the fracture width meter are selected for subsequent modeling, as shown in Table 2.

[0140]

[0141] Table 2 Description: TDEP represents the fracture depth; FVA_Ri represents the fracture width calculated by the i-th electrical imaging logging, i = 1, 2, ..., 5; FVA_C represents the fracture width measured by the fracture width meter.

[0142] 4) Comparing the fracture widths calculated from the five electrical imaging logs in the model well in step 3) with the reliable fracture widths measured using a fracture width gauge, we found a good linear relationship between the two on a logarithmic scale. Through mathematical calculations, we set the slope of each linear equation to 1, as shown in Figure 6.

[0143] In Figure 6, the horizontal axis represents the fracture width calculated from five electrical imaging logs, and the vertical axis represents the fracture width measured by the fracture width gauge. The intersection of the fracture width calculated from the i-th electrical imaging log and the width measured by the fracture width gauge is represented by a circle, a cross, a square, and an inverted triangle, respectively, where i = 1, 2, ..., 5. As can be seen from Figure 6, the fracture width calculated from the electrical imaging log generally increases as formation contrast increases.

[0144] Based on the fracture width calculated by the single electrical imaging logging in Figure 6 and the fracture width measured by the fracture width meter, a fitting relationship was established and the fracture width correction coefficient d was calculated, as shown in Table 3.

[0145]

[0146] Table 3 Description: Electrical imaging logging pass column, i-th time represents the i-th electrical imaging logging, i=1, 2, ..., 5; original fitting equation (first fitting relationship) column and new equation with slope adjusted to 1 (second fitting relationship) column, log10(FVA_Ri) represents the logarithm to the base 10 of the fracture width FVA_Ri calculated from the i-th electrical imaging logging, i=1, 2, ..., 5; log10(FVA_C) represents the logarithm to the base 10 of the fracture width FVA_C measured in the model well; R 2 represents the correlation coefficient of the original data fitting equation (first fitting equation); R²_1 represents the fitting correlation coefficient of the new equation (second fitting equation) after adjusting the slope of the original fitting equation of the data point to 1; n_1 represents the intercept of the new equation (second fitting equation); d represents the electrical imaging logging fracture width correction coefficient of the model well calculated based on n_1.

[0147] 5) Establish a relationship between the electrical imaging logging fracture width correction coefficient d and formation contrast (Rm / Rxo) for the model well. Mud resistivity Rm is directly measured, while the flush zone resistivity Rxo represents the background resistivity in the absence of fractures. This can be directly measured by reading the shallow lateral resistivity of the tight section in the model well without fractures (see Figure 7). Specific values ​​are shown in Table 4.

[0148]

[0149] According to Table 4, the calculation formula of the fracture width correction coefficient d applicable to any well is obtained:

[0150] d=2.3464-1909.3·(Rm / Rxo)

[0151] From the above relationship, we can know that a=2.3464 and f=-1909.3.

[0152] In Figure 7, the first track indicates the depth corresponding to the rock block of the rock physics model; the second track indicates the rock block number of the rock physics model; the third track indicates the static electrical imaging image; the fourth track indicates the dynamic electrical imaging image; and the fifth to ninth tracks indicate the shallow resistivity curves and mud resistivity corresponding to the i-th electrical imaging logging, where i = 1, 2, ..., 5, and the value corresponding to each dotted line is 1187 Ω.m.

[0153] 6) Calculate the flushing zone resistivity Rxo' at the electrical imaging logging fracture in Well B-2Sa (actual well). Combine the results of a and f in step 5) and substitute these parameters into the calculation model of the electrical imaging logging fracture width correction coefficient d' of the actual well:

[0154] d'=a+f·(Rm' / Rxo') (11)

[0155] In equation (11), d' represents the actual well electrical imaging logging fracture width correction coefficient, Rm' represents the actual mud resistivity at the fracture site in the actual well, obtained directly from well logging, a and f represent the parameters obtained according to step 5), and Rxo' represents the actual well flushing zone resistivity. The specific calculated fracture parameters for Well B-2Sa are shown in Table 5.

[0156]

[0157] Table 5 explains: Depth indicates the formation depth corresponding to the fracture; Fracture width indicates the electrical imaging logging fracture width calculated according to formula (1); Flushing zone resistivity indicates the calculated flushing zone resistivity at the fracture position of the actual well; Mud resistivity indicates the mud resistivity of the actual well obtained by logging; Fracture width correction coefficient is the electrical imaging logging fracture width correction coefficient of the actual well calculated according to formula (11); Corrected fracture width is the accurate fracture width after correction by the electrical imaging logging fracture width correction coefficient of the actual well.

[0158] The calculation and correction results of various parameters and fracture width are shown in Figure 8. The first track represents the formation depth; the second track represents the electrical imaging logging image; the third track represents the fracture tadpole picked up based on the electrical imaging logging image; the fourth track represents the original fracture width calculated based on the electrical imaging logging image using formula (1); the fifth track is the resistivity and deep resistivity curves, which are the resistivity of the flushing zone at the fracture location calculated based on the electrical imaging logging image and the deep resistivity curve measured by the lateral logging, respectively; the sixth track represents the electrical imaging logging fracture width correction coefficient of the actual well calculated according to the present invention; and the seventh track represents the fracture width calculation result after correction based on the electrical imaging logging fracture width correction coefficient of the actual well.

[0159] 7) Calculate the corrected electrical imaging logging fracture width of Well B-2Sa. Combined with the electrical imaging logging fracture width correction coefficient d value of the actual well calculated in step 6) (see d value in Table 5), the corrected electrical imaging logging fracture width of Well B-2Sa is calculated using the following formula:

[0160] W = d·W0 (where d is defined in the same way as d' in formula (10))

[0161] Where d represents the electrical imaging logging fracture width correction coefficient of the actual well obtained in step 6), which corresponds to the fracture width correction coefficient column in Table 5; Wo represents the electrical imaging logging fracture width calculated according to formula (1), which corresponds to the fracture width column in Table 5; W represents the fracture width corrected according to the electrical imaging logging fracture width correction coefficient of the actual well, which corresponds to the corrected fracture width column in Table 5.

[0162] The present invention provides a method for calculating the width of electrical imaging logging fractures based on a model well. The method uses a fracture width meter and electrical imaging logging equipment to measure the fracture width of the model well, compares the changing relationship between the fracture width obtained by the fracture width meter and the fracture width obtained by electrical imaging logging, and combines formation contrast to obtain a fracture width correction coefficient for any well. The method can correct and accurately calculate the fracture width of electrical imaging logging in an actual well, providing reliable guidance for the effectiveness and productivity evaluation of fractured reservoirs.

[0163] Furthermore, the present invention provides a method for calculating the width of electrical imaging logging fractures based on a model well, which also has the following advantages: 1. The invented method for considering the influence of formation contrast on the calculation results of electrical imaging logging fracture width based on a model well can clarify the key influencing factors of the fracture width calculation results, and point out a new direction for the accurate calculation of fracture width; 2. The method of the present invention can accurately calculate the width of electrical imaging logging fractures, avoiding the use of other methods to obtain the true width of fractures under high temperature and high pressure downhole environment, which can effectively save costs, and can make the fracture widths of different oil and gas fields and different wells comparable, and has strong practicality; 3. The present invention is a method for accurately calculating the width of electrical imaging logging fractures based on a model well, which provides an effective method with strong versatility and good economy, and can be widely used in the fracture evaluation of carbonate rock and metamorphic rock buried hill fracture reservoirs, providing reliable guidance for the effectiveness and productivity evaluation of fracture reservoirs.

[0164] The following describes the model well-based electrical imaging logging fracture width calculation system provided by the present invention. The model well-based electrical imaging logging fracture width calculation system described below and the model well-based electrical imaging logging fracture width calculation method described above can be referenced to each other.

[0165] 9 , the present invention provides a model well-based electrical imaging logging fracture width calculation system, which may include:

[0166] The first data acquisition module is used to obtain the fracture width of the model well by using a fracture width meter;

[0167] The second data acquisition module is used to: prepare the mud in the model well into muds with different resistivities in batches, obtain mud resistivity data and obtain the electrical imaging logging fracture width of the model well corresponding to the muds of different batches;

[0168] The first data processing module is used to obtain different formation contrasts based on different mud resistivity data and corresponding electrical imaging logging fracture widths;

[0169] The second data processing module is used to obtain a correction coefficient for the electrical imaging logging fracture width of the model well based on the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width of the model well;

[0170] The third data processing module is used to obtain the fracture width correction coefficient of any well based on the electrical imaging logging fracture width correction coefficient and formation contrast of the model well;

[0171] The correction module is used to obtain the corrected electrical imaging logging fracture width of an actual well according to the fracture width correction coefficient of any well.

[0172] According to a model well-based electrical imaging logging fracture width calculation system provided by the present invention, the second data processing module may include:

[0173] The fitting submodule is used to obtain a fitting relationship based on the variation pattern between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well;

[0174] The model well coefficient obtaining submodule is used to obtain the electrical imaging logging fracture width correction coefficient of the model well according to the fitting relationship.

[0175] Specifically, the fitting submodule is specifically used to: obtain a first fitting relationship according to the change law between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well; perform a first optimization process on the first fitting relationship to obtain a second fitting relationship;

[0176] Among them, the first fitting relationship is:

[0177] Log10(FVA_C)=s i Log10(FVA_Ri)+n i0 ,

[0178] In the first fitting equation, FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the model well at the time i of electrical imaging logging, s i Indicates the slope of the first fitting relationship for the i-th time, n i0 represents the intercept of the first fitting relationship of the i-th time;

[0179] And, wherein the second fitting relationship is:

[0180] Log10(FVA_C)=1·Log10(FVA_Ri)+n i ,

[0181] In the second fitting equation, FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the model well at the time i of electrical imaging logging, s i It represents the slope of the first fitting relationship of the i-th time, and its value is 1, n i Represents the intercept of the second fitting relationship of the i-th order.

[0182] Specifically, the model well coefficient obtaining submodule is specifically used to: perform a second optimization process on the second fitting relationship to obtain a relationship expression between the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width, thereby obtaining the electrical imaging logging fracture width correction coefficient of the model well;

[0183] The relational expression is:

[0184] FVA_C=10 ni ·FVA_Ri,

[0185] In the relational expression, FVA_C represents the fracture width of the model well obtained by the fracture width gauge, and FVA_R i represents the width of the electrical imaging logging fracture of the model well i, 10 ni represents the correction coefficient of the fracture width of the i-th electrical imaging logging in the model well, d i =10 ni , i=1,2,...,n.

[0186] It should be noted that the expression of the fracture width correction coefficient of any well is:

[0187] d=a+f·(Rm / Rxo),

[0188] In the expression of the fracture width correction coefficient of any well, d represents the fracture width correction coefficient of the model well, Rm represents the mud resistivity of the model well, Rxo represents the resistivity of the flushing zone of the model well, Rm / Rxo represents the formation contrast of the model well, and a and f represent the parameters applicable to the fracture width correction coefficient of any well. The expression of the fracture width correction coefficient of any well is used to fit the data points of the model well to obtain a and f, where d is the fracture width correction coefficient d obtained by n electrical imaging logging. i Rm and Rxo are composed of Rm during n times of electrical imaging logging. i 、Rxo i Composition, i=1,2,...,n.

[0189] According to a model well-based electrical imaging logging fracture width calculation system provided by the present invention, the third data processing module may include:

[0190] The actual well coefficient obtaining submodule is used to obtain the fracture width correction coefficient of the actual well according to the fracture width correction coefficient of any well;

[0191] The correction submodule is used to obtain the corrected electrical imaging logging fracture width of the actual well according to the fracture width correction coefficient of the actual well;

[0192] The expression of the fracture width correction coefficient of the actual well is:

[0193] d'=a+f·(Rm' / Rxo'),

[0194] In the expression of the fracture width correction coefficient of the actual well, d' represents the fracture width correction coefficient of the actual well, a and f represent the parameters of the fracture width correction coefficient of any well, Rm' represents the mud resistivity of the actual well, and Rxo' represents the resistivity of the flushing zone of the actual well.

[0195] According to a model well-based electrical imaging logging fracture width calculation system provided by the present invention, the correction module may include:

[0196] The mud resistivity obtaining submodule is used to: obtain the mud resistivity at the actual well fracture;

[0197] The flushing zone resistivity obtaining submodule is used to obtain the flushing zone resistivity at the electrical imaging logging fracture of the actual well according to the electrical imaging logging fracture width of the actual well;

[0198] The correction calculation submodule is used to obtain the corrected electrical imaging logging fracture width of the actual well based on the mud resistivity, flushing zone resistivity, electrical imaging logging fracture width and fracture width correction coefficient of the actual well;

[0199] The expression of the corrected electrical imaging logging fracture width of the actual well is:

[0200] W = d'·W0,

[0201] In the expression of the corrected electrical imaging logging fracture width of the actual well, W represents the corrected electrical imaging logging fracture width of the actual well, W0 represents the electrical imaging logging fracture width, and d' represents the fracture width correction coefficient of the actual well.

[0202] Those skilled in the art will appreciate that the various modules or steps of the present invention described above can be implemented using a general-purpose computing device, can be centralized on a single computing device, or can be distributed across a network of multiple computing devices. Alternatively, they can be implemented using program code executable by a computing device, and thus, can be stored in a storage device and executed by the computing device. In some cases, the steps shown or described herein can be performed in a different order than that described herein, or can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module. Thus, the present invention is not limited to any particular combination of hardware and software.

[0203] FIG10 illustrates a schematic diagram of the physical structure of an electronic device. As shown in FIG10 , the electronic device may include: a processor 810, a communications interface 820, a memory 830, and a communication bus 840. The processor 810, the communications interface 820, and the memory 830 communicate with each other via the communication bus 840. The processor 810 may call logic instructions in the memory 830 to execute a method for calculating fracture width based on electrical imaging logging in a model well. The method includes:

[0204] The fracture width of the model well is obtained by using a fracture width meter;

[0205] The mud in the model well is prepared into muds with different resistivities in batches, mud resistivity data is obtained, and the electrical imaging logging fracture width of the model well corresponding to the muds of different batches is obtained;

[0206] According to the different mud resistivity data and the corresponding electrical imaging logging fracture width, the formation contrast of different times is obtained;

[0207] According to the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a correction coefficient of the fracture width of the electrical imaging logging of the model well is obtained;

[0208] The fracture width correction coefficient of any well is obtained based on the electrical imaging logging fracture width correction coefficient and formation contrast of the model well;

[0209] According to the fracture width correction coefficient of any well, the corrected electrical imaging logging fracture width of the actual well is obtained.

[0210] Furthermore, the logic instructions in the aforementioned memory 830 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product, stored in a storage medium, includes instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, mobile hard drives, read-only memories (ROMs), random access memories (RAMs), magnetic disks, or optical disks.

[0211] On the other hand, the present invention further provides a computer program product, comprising a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the method for calculating fracture width based on electrical imaging logging of a model well provided by the above methods, the method comprising:

[0212] The fracture width of the model well is obtained by using a fracture width meter;

[0213] The mud in the model well is prepared into muds with different resistivities in batches, mud resistivity data is obtained, and the electrical imaging logging fracture width of the model well corresponding to the muds of different batches is obtained;

[0214] According to the different mud resistivity data and the corresponding electrical imaging logging fracture width, the formation contrast of different times is obtained;

[0215] According to the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a correction coefficient of the fracture width of the electrical imaging logging of the model well is obtained;

[0216] The fracture width correction coefficient of any well is obtained based on the electrical imaging logging fracture width correction coefficient and formation contrast of the model well;

[0217] According to the fracture width correction coefficient of any well, the corrected electrical imaging logging fracture width of the actual well is obtained.

[0218] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the method for calculating fracture width based on electrical imaging logging of a model well provided by the above methods is implemented. The method comprises:

[0219] The fracture width of the model well is obtained by using a fracture width meter;

[0220] The mud in the model well is prepared into muds with different resistivities in batches, mud resistivity data is obtained, and the electrical imaging logging fracture width of the model well corresponding to the muds of different batches is obtained;

[0221] According to the different mud resistivity data and the corresponding electrical imaging logging fracture width, the formation contrast of different times is obtained;

[0222] According to the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a correction coefficient of the fracture width of the electrical imaging logging of the model well is obtained;

[0223] The fracture width correction coefficient of any well is obtained based on the electrical imaging logging fracture width correction coefficient and formation contrast of the model well;

[0224] According to the fracture width correction coefficient of any well, the corrected electrical imaging logging fracture width of the actual well is obtained.

[0225] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, i.e., they may be located in one location or distributed across multiple network units. Some or all of the modules may be selected based on actual needs to achieve the objectives of the present embodiment. Persons of ordinary skill in the art will be able to understand and implement the present invention without inventive effort.

[0226] Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented using software plus a necessary general-purpose hardware platform, or of course, hardware. Based on this understanding, the essence of the above technical solution, or the portion that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a magnetic disk, or an optical disk, and includes a number of instructions for causing a computer device (such as a personal computer, server, or network device) to execute the methods described in each embodiment or certain portions of the embodiments.

[0227] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

[0228] CROSS-REFERENCE TO RELATED APPLICATIONS

[0229] This application claims priority to the Chinese patent application (application number 202311422983.3) filed on October 30, 2023, the entire contents of which are incorporated herein by reference. Industrial Applicability

[0230] The present invention provides a method for calculating the width of electrical imaging logging fractures based on a model well. The method uses a fracture width meter and an electrical imaging logging instrument to measure the fracture width of the model well, compares the variation relationship between the fracture width obtained by the fracture width meter and the fracture width calculated by electrical imaging logging, and combines formation contrast to obtain a fracture width correction coefficient for any well. The method can correct and accurately calculate the electrical imaging logging fracture width of an actual well, providing reliable guidance for the effectiveness and productivity evaluation of fractured reservoirs.

Claims

1. A method for calculating fracture width based on electrical imaging logging of a model well, characterized in that: include: The fracture width of the model well is obtained by using a fracture width gauge; The mud in the model well is prepared into muds with different resistivities in batches, mud resistivity data is obtained, and the electrical imaging logging fracture width of the model well corresponding to the muds of different batches is obtained; According to the different mud resistivity data and the corresponding electrical imaging logging fracture width, the formation contrast of different times is obtained; According to the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a correction coefficient of the fracture width of the electrical imaging logging of the model well is obtained; According to the electrical imaging logging fracture width correction coefficient and formation contrast of the model well, the fracture width correction coefficient of any well is obtained; According to the fracture width correction coefficient of any well, the corrected electrical imaging logging fracture width of the actual well is obtained.

2. The method for calculating fracture width based on electrical imaging logging of a model well according to claim 1, characterized in that: The method of obtaining the electrical imaging logging fracture width correction coefficient of the model well according to the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width of the model well comprises: According to the variation law between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a fitting relationship is obtained; According to the fitting relationship, the electrical imaging logging fracture width correction coefficient of the model well is obtained.

3. The method for calculating fracture width based on electrical imaging logging of a model well according to claim 2, characterized in that: The fitting relationship is obtained according to the variation law between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, including: According to the variation law between the fracture width obtained by the fracture width meter and the fracture width of the electrical imaging logging of the model well, a first fitting relationship is obtained; Performing a first optimization process on the first fitting relationship to obtain a second fitting relationship; Among them, the first fitting relationship is: Log10(FVA_C)=s i ·Log10(FVA_Ri)+n i0 , In the first fitting equation, FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the i-th electrical imaging logging in the model well, s i represents the slope of the first fitting relationship for the ith time, n i0 represents the intercept of the first fitting relationship of the i-th time; And, wherein the second fitting relationship is: Log10(FVA_C)=1·Log10(FVA_Ri)+n i , In the second fitting equation, FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the i-th electrical imaging logging in the model well, s i It represents the slope of the first fitting relationship of the ith time, and its value is 1, n i Represents the intercept of the i-th second fitting relationship.

4. The method for calculating fracture width based on electrical imaging logging of a model well according to claim 3, characterized in that: According to the fitting relationship, the electrical imaging logging fracture width correction coefficient of the model well is obtained, which is specifically: Performing a second optimization process on the second fitting relationship to obtain a relationship expression between the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width, thereby obtaining a correction coefficient for the electrical imaging logging fracture width of the model well; The relational expression is: FVA_C=10 ni ·FVA_Ri, In the relational expression, FVA_C represents the fracture width of the model well obtained by the fracture width meter, and FVA_R i represents the fracture width of the i-th electrical imaging logging in the model well, 10 ni represents the correction coefficient of the fracture width of the i-th electrical imaging logging of the model well, d i =10 ni , i=1,2,...,n.

5. The method for calculating fracture width based on electrical imaging logging of a model well according to claim 4, characterized in that: The expression of the fracture width correction coefficient of any well is: d = a + f (Rm / Rxo), In the expression of the fracture width correction coefficient of any well, d represents the fracture width correction coefficient of the model well, Rm represents the mud resistivity of the model well, Rxo represents the resistivity of the flushing zone of the model well, Rm / Rxo represents the formation contrast of the model well, a and f represent the parameters of the fracture width correction coefficient of any well, and a and f are obtained by fitting the data points of the model well using the expression of the fracture width correction coefficient of any well, where d is obtained by the fracture width correction coefficient d of n electrical imaging logging. i Rm and Rxo are composed of Rm during n times of electrical imaging logging. i , Rxo i Composition, i=1,2,...,n.

6. The method for calculating fracture width based on electrical imaging logging of a model well according to claim 5, characterized in that: The method of obtaining the corrected electrical imaging logging fracture width of an actual well according to the fracture width correction coefficient of any well comprises: According to the fracture width correction coefficient of any well, the fracture width correction coefficient of the actual well is obtained; According to the fracture width correction coefficient of the actual well, the corrected electrical imaging logging fracture width of the actual well is obtained; Among them, the expression of the fracture width correction coefficient of the actual well is: d'=a+f·(Rm' / Rxo'), In the expression of the fracture width correction coefficient of the actual well, d' represents the fracture width correction coefficient of the actual well, a and f represent the parameters of the fracture width correction coefficient of any well, Rm' represents the mud resistivity of the actual well, and Rxo' represents the resistivity of the flushing zone of the actual well.

7. The method for calculating fracture width based on electrical imaging logging of a model well according to claim 6, characterized in that: The method of obtaining the corrected electrical imaging logging fracture width of the actual well according to the fracture width correction coefficient of the actual well comprises: Get the mud resistivity at the actual well fracture; According to the electrical imaging logging fracture width of the actual well, the resistivity of the flushing zone at the electrical imaging logging fracture of the actual well is obtained; According to the mud resistivity of the actual well, the resistivity of the flushing zone, the electrical imaging logging fracture width and the fracture width correction coefficient, the corrected electrical imaging logging fracture width of the actual well is obtained; Among them, the expression of the actual well corrected electrical imaging logging fracture width is: W = d'·W0, In the expression of the corrected electrical imaging logging fracture width of the actual well, W represents the corrected electrical imaging logging fracture width of the actual well, W0 represents the electrical imaging logging fracture width, and d' represents the fracture width correction coefficient of the actual well.

8. A model well-based electrical imaging logging fracture width calculation system, characterized in that: include: The first data acquisition module is used to obtain the fracture width of the model well through a fracture width meter; The second data acquisition module is used to: prepare the mud in the model well into muds with different resistivities in batches, obtain mud resistivity data and obtain the electrical imaging logging fracture width of the model well corresponding to the muds of different batches; The first data processing module is used to obtain different formation contrasts according to different mud resistivity data and corresponding electrical imaging logging fracture widths; The second data processing module is used to obtain a correction coefficient for the electrical imaging logging fracture width of the model well according to the fracture width obtained by the fracture width meter and the electrical imaging logging fracture width of the model well; The third data processing module is used to obtain the fracture width correction coefficient of any well according to the electrical imaging logging fracture width correction coefficient and formation contrast of the model well; The correction module is used to obtain the corrected electrical imaging logging fracture width of the actual well according to the fracture width correction coefficient of any well.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the method for calculating fracture width based on electrical imaging logging of a model well as described in any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for calculating fracture width based on electrical imaging logging of a model well as claimed in any one of claims 1 to 7 is implemented.

Images