A method and device for logging quantitative prediction of carbonate rock corrosion intensity

Abstract

The application provides a logging quantitative prediction method for the dissolution intensity of carbonate rocks, which selects a standard sample of dissolution affected by dissolution more obviously as a basis for establishing quantitative evaluation of the dissolution intensity, and further verifies the standard sample of dissolution. Then, a target dissolution intensity development index parameter reflecting the dissolution intensity is obtained, and the dissolution intensity of the pore-type bioclastic limestone reservoir is described. The evaluation method fully combines the basic geological data, and constructs a quantitative evaluation index of the dissolution intensity according to the target dissolution intensity development index parameter, so that the evaluation of the development scale difference of the pore-type bioclastic limestone reservoir in the target interval is more in line with the actual situation of the oilfield. The method can provide a specific and quantitative method for the description of the dissolution intensity of the pore-type bioclastic limestone reservoir, make up for the lack of diagenetic transformation between the carbonate rock facies and the reservoir, accurately describe the dissolution intensity of the carbonate rock, and lay a foundation for the deterministic description of the reservoir unit.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological exploration and oil and gas exploration, and particularly relates to a method and device for logging quantitative prediction of dissolution intensity of carbonate rocks. BACKGROUND

[0002] Carbonate rock is the main reservoir rock type of carbonate rock reservoirs, and is developed on a large scale in some areas, and is the main business field of current oil and gas exploration and development. Unlike clastic rocks, carbonate rocks are easily affected by diagenesis, resulting in complex and irregular pore-permeability relationships in reservoirs, and may have large span of physical property distribution range of the same lithofacies or overlapping of physical property distribution of different lithofacies.

[0003] At present, the research on diagenesis of carbonate rocks is mainly qualitative description. Among the many diageneses of carbonate rocks, dissolution, as the diagenesis with the largest influence range and degree, plays an important role in reservoir reconstruction.

[0004] Therefore, how to accurately describe the dissolution intensity of carbonate rocks has become a technical problem to be solved in the field. SUMMARY

[0005] The embodiment of the present application provides a method for logging quantitative prediction of dissolution intensity of carbonate rocks, to solve the problem of how to accurately describe the dissolution intensity of carbonate rocks in the prior art.

[0006] The embodiment of the present application provides a method for logging quantitative prediction of dissolution intensity of carbonate rocks, comprising:

[0007] obtaining a dissolution standard sample associated with a target interval and intensity grade data corresponding to the dissolution standard sample;

[0008] obtaining a measured intersection line for the target interval according to the dissolution standard sample, and determining whether the measured intersection line deviates from a preset limestone lithology line;

[0009] If deviated, the obtained dissolution intensity development index parameter is correlated and screened according to the dissolution standard sample, to obtain a target dissolution intensity development index parameter;

[0010] constructing a quantitative evaluation index of dissolution intensity according to the target dissolution intensity development index parameter, and quantitatively predicting the dissolution intensity of the target interval to be measured according to the quantitative evaluation index of dissolution intensity.

[0011] Optionally, obtaining a dissolution standard sample associated with a target interval and intensity grade data corresponding to the dissolution standard sample comprises:

[0012] determine a lithofacies division strategy and porosity-permeability crossplot data for the target interval according to the obtained basic geological data; wherein the lithofacies division strategy is a lithofacies division strategy obtained by core description based on sedimentary characteristics;

[0013] obtain a dissolution candidate sample according to the lithofacies division strategy and the porosity-permeability crossplot data;

[0014] select a dissolution standard sample from the dissolution candidate sample according to multi-dimensional characteristic information of the target interval obtained, and intensity grade data corresponding to the dissolution standard sample.

[0015] Optionally, the determining the lithofacies division strategy and the porosity-permeability crossplot data for the target interval according to the obtained basic geological data comprises:

[0016] obtain single-well lithofacies research data and coring well basic data in the basic geological data;

[0017] obtain a lithofacies division strategy and porosity-permeability data samples of the target interval according to the single-well lithofacies research data;

[0018] make porosity-permeability crossplot of the porosity-permeability data samples of the target interval according to the coring well basic data to obtain porosity-permeability crossplot data.

[0019] Optionally, the obtaining the dissolution candidate sample according to the lithofacies division strategy and the porosity-permeability crossplot data comprises:

[0020] analyze the porosity-permeability crossplot data of the target interval according to the lithofacies division strategy to obtain a first result that physical properties of the target interval cannot be distinguished;

[0021] determine that an influencing factor of the target interval is a dissolution factor according to the porosity-permeability crossplot data of the target interval and the basic geological data;

[0022] eliminate data samples with complex and chaotic porosity-permeability crossplot from the porosity-permeability crossplot data of the target interval according to the first result and the dissolution factor to obtain the dissolution candidate sample.

[0023] Optionally, the multi-dimensional characteristic information of the target interval at least comprises diagenetic phenomenon characteristic information, reservoir-permeability space type characteristic information and reservoir-permeability space development and distribution characteristic information;

[0024] Correspondingly, the selecting the dissolution standard sample from the dissolution candidate sample according to the multi-dimensional characteristic information of the target interval obtained comprises:

[0025] According to the diagenetic phenomenon characteristic information, the reservoir space type characteristic information and the reservoir space development and distribution characteristic information, a dissolution standard sample is selected from the dissolution candidate samples.

[0026] Optionally, the strength grade data of the standard sample at least includes no development strength grade, weak development strength grade, medium development strength grade and development strength grade.

[0027] Optionally, according to the dissolution standard sample, a measurement intersection line for the target interval is obtained, and whether the measurement intersection line deviates from a preset limestone lithology line is judged.

[0028] According to the dissolution standard sample, a standard sample logging value is obtained.

[0029] According to the standard sample logging value, a first intersection line of a neutron porosity logging value and a density logging value and a second intersection line of a sonic time difference logging value and the density logging value are generated.

[0030] According to the comparison between the first intersection line and the second intersection line and the preset limestone lithology line respectively, whether the first intersection line and the second intersection line respectively deviate from the preset limestone lithology line is judged.

[0031] If deviating, it is determined that the target interval has dissolution characteristics under dissolution.

[0032] Optionally, according to the dissolution standard sample, a correlation screening is performed on a dissolution strength development index parameter obtained, to obtain a target dissolution strength development index parameter.

[0033] A candidate dissolution strength development index parameter associated with the dissolution standard sample is obtained.

[0034] According to a preset analysis strategy, the candidate dissolution strength development index parameter is screened to obtain the target dissolution strength development index parameter.

[0035] Optionally, the target dissolution strength development index parameter at least includes a density logging value, a neutron porosity logging value, a deep resistivity logging value, a microsphere resistivity logging value and a sonic time difference logging value.

[0036] Optionally, according to the target dissolution strength development index parameter, a quantitative evaluation index of dissolution strength is constructed.

[0037] According to the target dissolution strength development index parameter, a multivariate linear regression strategy is adopted to construct a dissolution strength development index calculation formula.

[0038] According to the dissolution intensity development index calculation formula, a calculation result for the dissolution intensity development index is obtained, and the calculation result is compared with the intensity grade data of the associated dissolution standard sample;

[0039] If the comparison result is consistent, the dissolution intensity development index calculation formula is determined as a quantitative evaluation index of the dissolution intensity.

[0040] Optionally, the method further comprises:

[0041] According to the dissolution intensity development index calculation formula, a dissolution intensity quantitative evaluation curve is obtained;

[0042] According to the intensity grade data of the standard sample, a dissolution intensity division standard reference curve is obtained;

[0043] The dissolution intensity quantitative evaluation curve and the dissolution intensity division standard reference curve are compared;

[0044] If the comparison result is consistent, the dissolution intensity quantitative evaluation curve is determined as a quantitative evaluation index of the dissolution intensity.

[0045] Technical effects and advantages of the present application:

[0046] The present application provides a logging quantitative prediction method for carbonate rock dissolution intensity. The logging quantitative prediction method for carbonate rock dissolution intensity provided by the present application considers the sedimentary background and diagenesis at the same time. On the basis of basic lithofacies division, dissolution standard samples which are significantly affected by dissolution are selected as the basis for establishing quantitative evaluation of dissolution intensity, and the dissolution standard samples are further verified. Then, the principal component analysis method and the multiple linear regression method are used to obtain a target dissolution intensity development index parameter which can reflect the dissolution intensity, so as to realize the characterization of the dissolution intensity of the pore-type bioclastic limestone reservoir. The evaluation method fully combines basic geological data, especially constructs a quantitative evaluation index of dissolution intensity according to the target dissolution intensity development index parameter, and is more in line with the actual situation of the oilfield in the evaluation of the development scale difference of the pore-type bioclastic limestone reservoir in the target interval. The method is more in line with the production needs of the oilfield site. The research results can provide a specific quantitative method for the characterization of the dissolution intensity of the pore-type bioclastic limestone reservoir, make up for the lack of diagenetic transformation between carbonate rock lithofacies and reservoirs, and accurately describe the dissolution intensity of carbonate rocks, thereby laying a foundation for deterministic characterization of reservoir units. BRIEF DESCRIPTION OF DRAWINGS

[0047] Figure 1 is a flowchart of the logging quantitative prediction method for carbonate rock dissolution intensity provided by the first embodiment of the present application;

[0048] Figure 2is a porosity-permeability crossplot provided by the first embodiment of the present application.

[0049] Figure 3 is a slice photo schematic diagram corresponding to the dissolution intensity grade in the first embodiment of the present application.

[0050] Figure 4 is a first crossplot of the neutron porosity logging value and the density logging value provided by the first embodiment of the present application.

[0051] Figure 5 is a second crossplot of the acoustic traveltime logging value and the density logging value provided by the first embodiment of the present application.

[0052] Figure 6 is a principal component analysis schematic diagram provided by the first embodiment of the present application.

[0053] Figure 7 is a multiple linear regression analysis schematic diagram provided by the first embodiment of the present application.

[0054] Figure 8 is a dissolution intensity index single-well columnar chart of a test well provided by the first embodiment of the present application.

[0055] Figure 9 is a schematic diagram of a logging quantitative prediction device for carbonate rock dissolution intensity provided by the second embodiment of the present application.

[0056] Figure 10 is a schematic diagram of an electronic device provided by the third embodiment of the present application. DETAILED DESCRIPTION

[0057] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the scope of protection of the present application.

[0058] At present, the research on carbonate rock diagenesis is mainly qualitative description. Among the numerous diagenesis of carbonate rock, dissolution, as the diagenesis with the largest influence range and influence degree, plays an important role in reservoir reconstruction. To this end, the embodiments of the present application provide a method to solve the technical problem of how to accurately describe the dissolution intensity of carbonate rock in the prior art.

[0059] The present application will be described in detail through multiple embodiments and drawings.

[0060] First Embodiment

[0061] The first embodiment of the present application provides a logging quantitative prediction method for carbonate rock dissolution intensity, and the following will be described in combination with Figure 1 The logging quantitative prediction method for carbonate rock dissolution intensity will be described in detail.

[0062] Step S101: Obtain a dissolution standard sample associated with a target interval and intensity grade data corresponding to the dissolution standard sample.

[0063] Step S102: Obtain a measured intersection line for the target interval according to the dissolution standard sample, and determine whether the measured intersection line deviates from a preset limestone lithology line.

[0064] Step S103: If deviated, perform correlation screening on the obtained dissolution intensity development index parameter according to the dissolution standard sample to obtain a target dissolution intensity development index parameter.

[0065] Step S104: Construct a quantitative evaluation index of dissolution intensity according to the target dissolution intensity development index parameter, and perform quantitative prediction on the dissolution intensity of the target interval to be measured according to the quantitative evaluation index of dissolution intensity.

[0066] In an embodiment, the step S101 of obtaining a dissolution standard sample associated with a target interval and intensity grade data corresponding to the dissolution standard sample includes the following contents:

[0067] Specifically, first, determine a lithofacies division strategy and porosity-permeability intersection data for the target interval according to the obtained basic geological data. In this step, the basic geological data is the pre-measured geological data, and the lithofacies division strategy is a basic lithofacies division strategy obtained by core description based on sedimentary characteristics. This step specifically includes: obtaining single-well lithofacies research data and coring well basic data in the basic geological data, and obtaining a lithofacies division strategy and porosity-permeability data sample for the target interval according to the single-well lithofacies research data, wherein the single-well lithofacies research data at least includes coarse-grained particle limestone facies data, fine-grained particle limestone data, particle-dominated mud particle limestone facies data, mud-dominated particle limestone facies data, grain-mud limestone facies data, and micritic limestone facies data, and the corresponding lithofacies division strategy is obtained according to different single-well lithofacies research data. Then, according to the coring well basic data, the porosity-permeability data sample of the target interval is intersected to obtain porosity-permeability intersection data. See Figure 2 , Figure 2 is a basic lithofacies division scheme porosity-permeability intersection graph provided by the embodiment of the present application.

[0068] After obtaining the lithofacies division strategy and the porosity-permeability crossplot data of the target interval, a dissolution candidate sample is obtained according to the lithofacies division strategy and the porosity-permeability crossplot data. Specifically, the porosity-permeability crossplot data of the target interval is analyzed according to the lithofacies division strategy to obtain a first result that the physical property of the target interval cannot be distinguished, or a second result that the physical property of the target interval can be distinguished, for example, the porosity-permeability crossplot data of the target interval is analyzed according to the lithofacies division strategy to obtain the second result that the physical property of the target interval can be distinguished, and the second result is directly used for strength grade division. Figure 2 It can be seen that the grain limestone and the micritic limestone in the area have good zoning, and the second result of the physical property distinction can be directly used for strength grade division. The porosity-permeability distribution of the siltstone and the granular mudstone in the area is more complex and chaotic. Then, the influencing factor of the target interval is determined to be the dissolution factor according to the porosity-permeability crossplot data and the basic geological data of the target interval, and the data sample with complex and chaotic porosity-permeability crossplot is removed from the porosity-permeability crossplot data of the target interval according to the first result and the dissolution factor to obtain the dissolution candidate sample.

[0069] After obtaining the dissolution candidate sample, a dissolution standard sample is selected from the dissolution candidate sample according to the obtained multi-dimensional characteristic information of the target interval, and strength grade data corresponding to the dissolution standard sample is obtained. Specifically, in an embodiment, the multi-dimensional characteristic information of the target interval at least includes: diagenetic phenomenon characteristic information, reservoir-permeability space type characteristic information and reservoir-permeability space development and distribution characteristic information. The diagenetic phenomenon characteristic information is used to determine that the influencing factor of the target interval is the dissolution factor, and the reservoir-permeability space type characteristic information and the reservoir-permeability space development and distribution characteristic information at least include: porosity type characteristic information, secondary porosity development degree characteristic information and secondary porosity proportion characteristic information.

[0070] Correspondingly, the dissolution standard sample is selected from the dissolution candidate sample according to the obtained multi-dimensional characteristic information of the target interval, which includes: the dissolution standard sample is selected from the dissolution candidate sample according to the diagenetic phenomenon characteristic information, the reservoir-permeability space type characteristic information and the reservoir-permeability space development and distribution characteristic information. In an embodiment, as shown in Table 1 and Table 2, Figure 3 as shown in Table 1 and Table 2, Figure 3 which are schematic diagrams of the thin section photos corresponding to the dissolution strength grade division in the embodiments of the present application. According to the dissolution strength control secondary porosity development degree, the dissolution strength grading standard can be established, according to the high and low of the development degree, it is determined that the dissolution strength can be divided into non-development strength grade, weak development strength grade, medium development strength grade and development strength grade from strong to weak. At the same time, the characteristics of different grades of dissolution strength are determined. In addition, the above four grades of strength grade are valued, that is, the non-development strength grade is set to 1, the weak development strength grade is set to 2, the medium development strength grade is set to 3 and the development strength grade is set to 4.

[0071] Table 1

[0072]

[0073] In addition, in an embodiment, according to the intensity level data of the standard sample, a dissolution intensity division standard reference curve is obtained, which is used for subsequent verification of the accuracy of the dissolution intensity value calculated by logging.

[0074] In an embodiment, the step S102 of obtaining the measurement intersection line for the target interval according to the dissolution standard sample and judging whether the measurement intersection line deviates from the preset limestone lithology line specifically includes the following contents:

[0075] Specifically, first, the standard sample logging value is obtained according to the dissolution standard sample, wherein the standard sample logging value is the logging value of the dissolution standard sample. Then, the first intersection line of the neutron porosity logging value and the density logging value and the second intersection line of the acoustic time difference logging value and the density logging value are generated according to the standard sample logging value. The neutron porosity logging value and the acoustic time difference logging value increase with the decrease of the density logging value and have a good linear relationship. After the dissolution of the dissolution interval of the target interval, the porosity is improved, and the corresponding neutron porosity value and acoustic time difference logging value are higher at the same density logging value. That is, after the dissolution, the dissolution interval is more biased towards the direction of good physical property distribution in the neutron porosity logging value-density logging value crossplot and the acoustic time difference logging-density logging value crossplot. Finally, the first intersection line and the preset limestone lithology line are compared, and the second intersection line and the preset limestone lithology line are compared, and it is judged whether the first intersection line and the second intersection line deviate from the preset limestone lithology line, respectively. Specifically, refer to the first intersection line and the second intersection line shown in Figure 4 and Figure 5 , Figure 4 the first intersection line of the neutron porosity logging value and the density logging value provided by the embodiments of the present application. Figure 5 the second intersection line of the acoustic time difference logging value and the density logging value provided by the embodiments of the present application. From Figure 4 and Figure 5 It can be seen that the neutron, density, and acoustic logging curves of the dissolution interval deviate towards the direction of better physical property, especially the neutron-density crossplot, which deviates obviously from the limestone lithology line, which conforms to the general rule of the reconstruction of the reservoir by the dissolution interval. That is, if the measurement intersection line deviates from the preset limestone lithology line, it is determined that the target interval has the dissolution characteristics under the dissolution.

[0076] After determining that the target interval has the dissolution characteristics under the dissolution, in an example, step S103 performs correlation screening on the obtained dissolution intensity development index parameters according to the dissolution standard sample to obtain target dissolution intensity development index parameters, which specifically include the following contents: Specifically, first, candidate dissolution intensity development index parameters associated with the dissolution standard sample are obtained, wherein, in an embodiment, the candidate dissolution intensity development index parameters at least include density logging values, neutron porosity logging values, deep resistivity logging values, microsphere resistivity logging values, and acoustic time difference logging values, and natural gamma logging values. Then, according to a preset analysis strategy, the candidate dissolution intensity development index parameters are screened to obtain target dissolution intensity development index parameters. In an embodiment, the preset analysis strategy is principal component analysis, and the target dissolution intensity development index parameters at least include density logging values, neutron porosity logging values, deep resistivity logging values, microsphere resistivity logging values, and acoustic time difference logging values. See Figure 6 and Table 2 below, Figure 6 is a principal component analysis diagram provided by the embodiment of the present application. Table 2 is a principal component analysis correlation coefficient chart provided by the embodiment of the present application.

[0077] Table 2

[0078]

[0079] After obtaining the target dissolution intensity development index parameters, a quantitative evaluation index of the dissolution intensity is constructed according to the target dissolution intensity development index parameters, in an example, step S103 of constructing the quantitative evaluation index of the dissolution intensity according to the target dissolution intensity development index parameters includes the following contents:

[0080] Specifically, first, according to the target dissolution intensity development index parameters, a multivariate linear regression strategy is used to construct a dissolution intensity development index calculation formula. In an example, the dissolution intensity development index calculation formula is as follows: Figure 7 is a multivariate linear regression analysis diagram provided by the embodiment of the present application. After linear regression, the dissolution intensity development index calculation formula is as follows:

[0081]

[0082] wherein, LOG(K_CORR) represents the dissolution intensity development index; DEN represents the density logging value, CNL represents the neutron porosity logging value, LLD represents the deep resistivity logging value, MSFL represents the microsphere resistivity logging value, and AC represents the acoustic time difference logging value.

[0083] Then, a calculation result for the dissolution intensity development index is obtained according to the dissolution intensity development index calculation formula, and the calculation result is compared with the intensity grade data of the standard sample associated with the dissolution effect, and if the comparison result is consistent, the dissolution intensity development index calculation formula is determined as the quantitative evaluation index of the dissolution intensity.

[0084] In an embodiment, the method further comprises: obtaining a dissolution intensity quantitative evaluation curve according to the dissolution intensity development index calculation formula, obtaining a dissolution intensity division standard reference curve according to the intensity grade data of the standard sample, comparing the dissolution intensity quantitative evaluation curve with the dissolution intensity division standard reference curve, and if the comparison result is consistent, determining the dissolution intensity quantitative evaluation curve as the quantitative evaluation index of the dissolution intensity. Figure 8 , Figure 8 A dissolution intensity index columnar chart of a test well is provided in the embodiment. After the dissolution intensity dissolution index is divided according to the range, the dissolution intensity dissolution index is compared with the dissolution intensity, and it is found that there is a good matching relationship between the two, so that the dissolution intensity logging curve can be popularized in the oilfield.

[0085] After the quantitative evaluation index of the dissolution intensity is constructed according to the target dissolution intensity development index parameter, the dissolution intensity of the target layer to be measured can be quantitatively predicted according to the quantitative evaluation index of the dissolution intensity.

[0086] The logging quantitative prediction method of the dissolution intensity of the carbonate rock provided in the application considers the sedimentary background and diagenesis at the same time, and on the basis of the basic lithofacies division, the dissolution effect standard sample which is significantly affected by the dissolution effect is selected as the basis for establishing the dissolution intensity quantitative evaluation, and the dissolution effect standard sample is further verified. Then, the principal component analysis method and the multiple linear regression method are used to obtain the target dissolution intensity development index parameter which can reflect the dissolution intensity, so as to realize the characterization of the dissolution intensity of the pore-type bioclastic limestone reservoir. The evaluation method fully combines the basic geological data, especially the quantitative evaluation index of the dissolution intensity constructed according to the target dissolution intensity development index parameter, and the evaluation of the development scale difference of the pore-type bioclastic limestone reservoir in the target layer is more consistent with the actual situation of the oilfield, and the method is more consistent with the production needs of the oilfield site. The research results can provide a specific quantitative method for the characterization of the dissolution intensity of the pore-type bioclastic limestone reservoir, make up for the lack of diagenetic transformation between the carbonate rock lithofacies and the reservoir, and accurately describe the dissolution intensity of the carbonate rock, thereby laying a foundation for the deterministic characterization of the reservoir unit.

[0087] Second embodiment

[0088] In the first embodiment, a method for logging quantitative prediction of carbonate rock dissolution intensity is provided. Correspondingly, a device for logging quantitative prediction of carbonate rock dissolution intensity is provided in the third embodiment of the present application. Since the device embodiment is basically similar to the first embodiment of the method, it is described simply, and the relevant part is described in the method embodiment. The device embodiment described below is only illustrative.

[0089] Please refer to Figure 9 A schematic diagram of a device for logging quantitative prediction of carbonate rock dissolution intensity is provided in the second embodiment of the present application.

[0090] The device for logging quantitative prediction of carbonate rock dissolution intensity comprises: a dissolution standard sample obtaining unit 901 configured to obtain a dissolution standard sample associated with a target interval and intensity grade data corresponding to the dissolution standard sample; a measurement intersection line processing unit 902 configured to obtain a measurement intersection line for the target interval according to the dissolution standard sample and determine whether the measurement intersection line deviates from a preset limestone lithology line; a target dissolution intensity development index parameter obtaining unit 903 configured to, if the measurement intersection line deviates from the preset limestone lithology line, perform correlation screening on a dissolution intensity development index parameter obtained according to the dissolution standard sample to obtain a target dissolution intensity development index parameter; and a quantitative evaluation index constructing unit 904 configured to construct a quantitative evaluation index of dissolution intensity according to the target dissolution intensity development index parameter and perform quantitative prediction of the dissolution intensity of the target interval to be measured according to the quantitative evaluation index of dissolution intensity.

[0091] Optionally, the dissolution standard sample obtaining unit comprises:

[0092] A data processing unit configured to determine a facies division strategy and porosity-permeability intersection data for the target interval according to the obtained basic geological data; wherein the facies division strategy is a basic facies division strategy obtained by core description based on sedimentary characteristics;

[0093] A dissolution candidate sample obtaining unit configured to obtain a dissolution candidate sample according to the facies division strategy and the porosity-permeability intersection data;

[0094] A dissolution standard sample selecting unit configured to select a dissolution standard sample and intensity grade data corresponding to the dissolution standard sample from the dissolution candidate sample according to the obtained multi-dimensional feature information of the target interval.

[0095] Optionally, the data processing unit comprises:

[0096] A data processing first sub-unit configured to obtain single-well facies research data and coring well basic data in the basic geological data;

[0097] The data processing second sub-unit is configured to obtain a lithofacies division strategy and a porosity-permeability data sample of the target interval according to the single-well lithofacies study data.

[0098] The data processing third sub-unit is configured to make a porosity-permeability crossplot of the porosity-permeability data sample of the target interval according to the coring well basic data to obtain porosity-permeability crossplot data.

[0099] Optionally, the dissolution candidate sample obtaining unit comprises:

[0100] The first result obtaining unit is configured to analyze the porosity-permeability crossplot data of the target interval according to the lithofacies division strategy to obtain a first result of undistinguishable physical properties of the target interval.

[0101] The dissolution factor determining unit is configured to determine, according to the porosity-permeability crossplot data of the target interval and the basic geological data, that the influencing factor of the target interval is a dissolution factor.

[0102] The data excluding unit is configured to exclude, according to the first result and the dissolution factor, a data sample with complex and chaotic porosity-permeability crossplot from the porosity-permeability crossplot data of the target interval to obtain a dissolution candidate sample.

[0103] Optionally, the multi-dimension feature information of the target interval at least comprises diagenetic phenomenon feature information, reservoir-permeability space type feature information and reservoir-permeability space development and distribution feature information.

[0104] Correspondingly, the dissolution standard sample selecting unit is configured to select a dissolution standard sample from the dissolution candidate sample according to the diagenetic phenomenon feature information, the reservoir-permeability space type feature information and the reservoir-permeability space development and distribution feature information.

[0105] Optionally, the strength grade data of the standard sample at least comprises no development strength grade, weak development strength grade, medium development strength grade and development strength grade.

[0106] Optionally, the dissolution standard sample selecting unit comprises:

[0107] The standard sample logging value obtaining unit is configured to obtain a standard sample logging value according to the dissolution standard sample.

[0108] The crossplot line generating unit is configured to generate a first crossplot line of the neutron porosity logging value and the density logging value and a second crossplot line of the acoustic time difference logging value and the density logging value according to the standard sample logging value.

[0109] A judging unit is configured to compare the first intersection line and the second intersection line with the preset limestone lithology line respectively, and determine whether the first intersection line and the second intersection line deviate from the preset limestone lithology line respectively.

[0110] A determining unit is configured to determine that the target layer section has a dissolution feature under dissolution if the deviation exists.

[0111] Optionally, the target dissolution intensity development index parameter obtaining unit comprises:

[0112] A candidate dissolution intensity development index parameter obtaining unit is configured to obtain a candidate dissolution intensity development index parameter associated with the dissolution standard sample.

[0113] A target dissolution intensity development index parameter screening unit is configured to screen the candidate dissolution intensity development index parameter according to a preset analysis strategy to obtain a target dissolution intensity development index parameter.

[0114] Optionally, the target dissolution intensity development index parameter at least comprises a density logging value, a neutron porosity logging value, a deep resistivity, a microsphere resistivity, and an acoustic time difference logging value.

[0115] Optionally, the quantitative evaluation index constructing unit comprises:

[0116] A dissolution intensity development index calculation formula constructing unit is configured to construct a dissolution intensity development index calculation formula according to the target dissolution intensity development index parameter by using a multiple linear regression strategy.

[0117] A first comparison unit is configured to obtain a calculation result of the dissolution intensity development index according to the dissolution intensity development index calculation formula, and compare the calculation result with the intensity grade data of the associated dissolution standard sample.

[0118] A quantitative evaluation index determining unit is configured to determine the dissolution intensity development index calculation formula as a quantitative evaluation index of dissolution intensity if the comparison result is consistent.

[0119] Optionally, the quantitative evaluation index constructing unit further comprises:

[0120] A dissolution intensity quantitative evaluation curve obtaining unit is configured to obtain a dissolution intensity quantitative evaluation curve according to the dissolution intensity development index calculation formula.

[0121] A dissolution intensity division standard reference curve obtaining unit is configured to obtain a dissolution intensity division standard reference curve according to the intensity grade data of the standard sample.

[0122] The second comparison unit is configured to compare the dissolution strength quantitative evaluation curve with the dissolution strength division standard reference curve.

[0123] The quantitative evaluation index determination unit is configured to determine the dissolution strength quantitative evaluation curve as a quantitative evaluation index of the dissolution strength if the comparison result is consistent.

[0124] The third embodiment

[0125] Corresponding to the above-mentioned method embodiments, the third embodiment of the present application further provides an electronic device. Figure 10 As shown in the figure, Figure 10 is a schematic diagram of an electronic device provided in the third embodiment of the present application. The electronic device comprises at least one processor 1001, at least one communication interface 1002, at least one memory 1003 and at least one communication bus 1004; optionally, the communication interface 1002 can be an interface of a communication module, such as an interface of a GSM module; the processor 1001 can be a processor CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present application. The memory 1003 can include a high-speed RAM memory, and can also include a non-volatile memory, for example, at least one disk memory. Among them, the memory 1003 stores a program, and the processor 1001 calls the program stored in the memory 1003 to execute the method provided in the above-mentioned embodiments of the present application.

[0126] The fourth embodiment

[0127] Corresponding to the above-mentioned method, the fourth embodiment of the present application further provides a computer storage medium. The computer storage medium stores a computer program, and the computer program is run by a processor to execute the method provided in the above-mentioned embodiments of the present application.

[0128] Although the present application is disclosed with the preferred embodiments as above, it is not intended to limit the present application, and any person skilled in the art can make possible changes and modifications without departing from the spirit and scope of the present application, therefore the protection scope of the present application should be subject to the scope defined by the claims of the present application.

[0129] In a typical configuration, the computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0130] Memory can include non-persistent memory in computer-readable media, random access memory (RAM), and / or non-volatile memory such as read-only memory (ROM) or flash memory (flash RAM). Memory is an example of computer-readable media.

[0131] 1. Computer-readable media includes permanent and non-permanent, removable and non-removable media, which can be implemented by any method or technology to store information. Information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device. According to the definition herein, computer-readable media does not include non-transitory computer-readable media, such as modulated data signals and carriers.

[0132] 2. Those skilled in the art should understand that the embodiments of the present application can be provided as a method, system or computer program product. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0133] It should be noted that the embodiments of the present application can involve the use of user data. In actual application, user-specific personal data can be used in the schemes described herein in accordance with the requirements of applicable laws and regulations of the country (for example, user's explicit consent, user's actual notification, etc.), within the scope permitted by applicable laws and regulations.

[0134] It should be noted that the user information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties, and the collection, use and processing of related data need to comply with relevant laws, regulations and standards of relevant countries and regions, and provide corresponding operation portal for user to choose authorization or refusal.

Claims

1. A method for logging quantitative prediction of carbonate rock dissolution intensity, characterized in that, The method comprises the following steps: obtaining a dissolution standard sample associated with a target interval and intensity grade data corresponding to the dissolution standard sample, comprising: determining a lithofacies division strategy and porosity-permeability crossplot data for the target interval according to the obtained basic geological data; wherein the lithofacies division strategy is a lithofacies division strategy obtained by core description based on sedimentary characteristics; obtaining a dissolution candidate sample according to the lithofacies division strategy and the porosity-permeability crossplot data; selecting a dissolution standard sample from the dissolution candidate sample according to the obtained multi-dimensional feature information of the target interval, and intensity grade data corresponding to the dissolution standard sample; obtaining a measurement crossplot for the target interval according to the dissolution standard sample, and determining whether the measurement crossplot deviates from a preset limestone lithology line, comprising: obtaining a standard sample logging value according to the dissolution standard sample; generating a first crossplot of neutron porosity logging value and density logging value and a second crossplot of acoustic time difference logging value and density logging value according to the standard sample logging value; comparing the first crossplot and the second crossplot with the preset limestone lithology line respectively, and determining whether the first crossplot and the second crossplot deviate from the preset limestone lithology line respectively; if deviating, determining that the target interval has dissolution characteristics under dissolution; if deviating, performing correlation screening on the obtained dissolution intensity development index parameter according to the dissolution standard sample to obtain a target dissolution intensity development index parameter; constructing a quantitative evaluation index of dissolution intensity according to the target dissolution intensity development index parameter, and quantitatively predicting the dissolution intensity of the target interval to be measured according to the quantitative evaluation index of dissolution intensity, wherein the construction of the quantitative evaluation index of dissolution intensity according to the target dissolution intensity development index parameter comprises: adopting a multiple linear regression strategy to construct a dissolution intensity development index calculation formula according to the target dissolution intensity development index parameter; comparing the calculation result of the dissolution intensity development index calculation formula with the intensity grade data of the associated dissolution standard sample; if the comparison result is consistent, determining the dissolution intensity development index calculation formula as the quantitative evaluation index of dissolution intensity.

2. The method for logging quantitative prediction of carbonate rock dissolution intensity according to claim 1, characterized in that, The determination of the lithofacies division strategy and the porosity-permeability crossplot data for the target interval according to the obtained basic geological data comprises: obtaining single-well lithofacies research data and coring well basic data in the basic geological data; obtaining a lithofacies division strategy and porosity-permeability data samples of the target interval according to the single-well lithofacies research data; performing porosity-permeability crossplot on the porosity-permeability data samples of the target interval according to the coring well basic data to obtain porosity-permeability crossplot data.

3. The method for logging quantitatively predicting the dissolution intensity of carbonate rock according to claim 1, characterized in that, The obtaining of the dissolution candidate sample according to the lithofacies division strategy and the porosity-permeability crossplot data comprises: analyzing the porosity-permeability crossplot data of the target interval according to the lithofacies division strategy to obtain a first result that physical properties cannot be distinguished for the target interval; determine, according to the poroperm crossplot data and the basic geological data of the target interval, that the influencing factor of the target interval is a dissolution factor; remove, according to the first result and the dissolution factor, a complex and chaotic data sample from the poroperm crossplot data of the target interval to obtain a dissolution candidate sample.

4. The method of quantitative prediction of carbonate rock dissolution intensity from well logging data according to claim 1, characterized in that, The multi-dimensional characteristic information of the target interval at least includes: diagenetic phenomenon characteristic information, reservoir and permeation space type characteristic information, and reservoir and permeation space development and distribution characteristic information; Correspondingly, the selecting of the dissolution standard sample from the dissolution candidate sample according to the obtained multi-dimensional characteristic information of the target interval includes: selecting the dissolution standard sample from the dissolution candidate sample according to the diagenetic phenomenon characteristic information, the reservoir and permeation space type characteristic information, and the reservoir and permeation space development and distribution characteristic information.

5. The method for logging quantitatively predicting the karst corrosion intensity of carbonate rock according to claim 1, characterized in that, The strength grade data of the standard sample at least includes: no development strength grade, weak development strength grade, medium development strength grade, and development strength grade.

6. The method of quantitative prediction of carbonate rock dissolution intensity from well logging data according to claim 1, characterized in that, The correlation screening of the obtained dissolution strength development index parameter according to the dissolution standard sample to obtain a target dissolution strength development index parameter includes: obtaining a candidate dissolution strength development index parameter associated with the dissolution standard sample; screening the candidate dissolution strength development index parameter according to a preset analysis strategy to obtain a target dissolution strength development index parameter.

7. The method for logging quantitatively predicting the karst corrosion intensity of the carbonate rock according to claim 6, characterized in that, The target dissolution strength development index parameter at least includes: density logging value, neutron porosity logging value, deep resistivity logging value, microsphere resistivity logging value, and acoustic time difference logging value.

8. The method for logging quantitatively predicting the karst corrosion intensity of carbonate rock according to claim 1, characterized in that, Further includes: obtaining a dissolution strength quantitative evaluation curve according to the dissolution strength development index calculation formula; obtaining a dissolution strength division standard reference curve according to the strength grade data of the standard sample; comparing the dissolution strength quantitative evaluation curve and the dissolution strength division standard reference curve; if the comparison result is consistent, determining the dissolution strength quantitative evaluation curve as a quantitative evaluation index of the dissolution strength.

9. A well logging quantitative prediction device for the intensity of carbonate rock karst erosion, characterized in that, includes: a dissolution standard sample obtaining unit for obtaining a dissolution standard sample associated with the target interval and strength grade data corresponding to the dissolution standard sample, including: determining, according to the obtained basic geological data, a facies division strategy and poroperm crossplot data for the target interval; wherein the facies division strategy is a basic facies division strategy obtained by core description based on sedimentary characteristics; obtaining a dissolution candidate sample according to the facies division strategy and the poroperm crossplot data; selecting a dissolution standard sample from the dissolution candidate sample according to the obtained multi-dimensional characteristic information of the target interval, and strength grade data corresponding to the dissolution standard sample; a measurement intersection line processing unit for obtaining a measurement intersection line for the target interval according to the dissolution standard sample, and judging whether the measurement intersection line deviates from a preset limestone lithology line, including: obtaining a standard sample logging value according to the dissolution standard sample; generating a first crossplot of the neutron porosity logging value and the density logging value and a second crossplot of the acoustic time difference logging value and the density logging value according to the standard sample logging value; comparing the first crossplot and the second crossplot with the preset limestone lithology line respectively, and judging whether the first crossplot and the second crossplot deviate from the preset limestone lithology line respectively; if deviating, determining that the target layer section has the dissolution characteristics under the dissolution effect; a target dissolution intensity development index parameter obtaining unit, configured to, if deviating, perform correlation screening on the obtained dissolution intensity development index parameter according to the dissolution effect standard sample, to obtain a target dissolution intensity development index parameter; a quantitative evaluation index constructing unit, configured to construct a quantitative evaluation index of the dissolution intensity according to the target dissolution intensity development index parameter, and perform quantitative prediction on the dissolution intensity of the target layer section to be measured according to the quantitative evaluation index of the dissolution intensity, wherein the construction of the quantitative evaluation index of the dissolution intensity according to the target dissolution intensity development index parameter comprises: adopting a multiple linear regression strategy to construct a dissolution intensity development index calculation formula according to the target dissolution intensity development index parameter; obtaining a calculation result of the dissolution intensity development index according to the dissolution intensity development index calculation formula, and comparing the calculation result with the intensity grade data of the associated dissolution effect standard sample; if the comparison result is consistent, determining the dissolution intensity development index calculation formula as the quantitative evaluation index of the dissolution intensity.

Images