Single well analysis method and device for limestone reservoir pore structure type
By establishing a chart showing the relationship between pore structure type and physical properties of limestone rock samples, and correcting the interpretation results of well logging properties, the shortcomings in evaluating the production capacity of single wells in limestone oil and gas reservoirs were resolved. This enabled a scientific analysis of the pore structure type of limestone reservoirs, providing an important basis for oil and gas exploration and development.
Patent Information
- Application Number
- CN202410967757.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-07-18
- Publication Date
- 2026-01-20
AI Technical Summary
Existing technologies are insufficient to effectively evaluate the single-well production capacity of non-fractured and non-cavitary limestone oil and gas reservoirs, especially the impact of pore structure on production capacity, resulting in inadequate evaluation of the production capacity of limestone oil and gas reservoirs.
A chart was established to show the relationship between pore structure type and physical properties of limestone rock samples. By correcting the interpretation results of well logging, a cross-plot of porosity and permeability was constructed. Combining core, well logging, and test data, the pore structure type of limestone reservoir was analyzed.
It provides a scientific basis for evaluating the single-well productivity of limestone oil and gas reservoirs and for seismic prediction of favorable reservoirs in complex carbonate rocks, and improves the accuracy of evaluating the pore structure type of limestone reservoirs.
Smart Images
Figure CN121363415A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of oil and gas exploration and development, and particularly relates to a single-well analysis method and device for pore structure types of limestone reservoirs. BACKGROUND
[0002] At present, there are few basic researches on the main control factors of single-well productivity of carbonate rock oil and gas reservoirs at home and abroad. The single-well productivity evaluation of carbonate rock reservoirs in the exploration and evaluation stage plays a crucial role in the seismic prediction of complex carbonate rock favorable reservoirs, oilfield development strategies and development plan design. For non-fracture and non-cave type limestone oil and gas reservoirs, the production capacity of single wells is very different from that of single wells of clastic rock oil and gas reservoirs, and it is difficult to effectively evaluate the production capacity of single wells. Specifically, the main control factors of the production capacity of single wells of such limestone oil and gas reservoirs are not only porosity and permeability, but also have great relevance to reservoir pore structure. Therefore, it is necessary to establish a single-well evaluation method suitable for pore type limestone oil and gas reservoirs, which is very important for the production capacity evaluation of such oil and gas reservoirs and the seismic prediction of complex carbonate rock favorable reservoirs. SUMMARY
[0003] To solve the above problems, the present application aims to provide a single-well analysis method and device for pore structure types of limestone reservoirs. The main technical feature of the present application is to establish a pore structure type and physical property relationship chart of limestone rock samples, correct the physical property interpretation results of logging according to the physical property data of limestone reservoir rock samples, establish a corrected logging porosity and permeability crossplot, and apply the pore structure physical property chart of rock samples to the corrected logging physical property crossplot, thereby establishing a pore structure type columnar chart of the whole well section of the target layer of limestone. This technology can provide a scientific basis for single-well limestone productivity evaluation in the oil and gas exploration and evaluation stage and seismic prediction of complex carbonate rock favorable reservoirs. To achieve the above purpose, the present application provides the following technical solutions:
[0004] The present application provides a single-well analysis method for pore structure types of limestone reservoirs, which comprises,
[0005] Based on the porosity and permeability data of limestone rock samples and the pore microstructure type of limestone, a pore structure type and physical property relationship chart of limestone rock samples is established;
[0006] Based on the physical property curve data of logging interpretation and the physical property data of rock samples of the whole well section of the target layer of single-well limestone, a logging physical property crossplot is established;
[0007] The pore structure type and physical property relationship chart of the limestone rock samples is applied to the logging physical property crossplot, and a pore structure type result chart of the whole well section of the target layer of limestone is established.
[0008] Based on the pore structure type result map of the whole well section of the limestone target layer, a single well of the limestone reservoir pore structure type is analyzed.
[0009] Further, the pore structure type and physical property relationship chart of the limestone rock sample is established, including:
[0010] Based on the porosity and permeability data of the limestone rock sample, a physical property cross plot of the limestone rock sample is constructed;
[0011] The abnormal value in the physical property cross plot of the limestone rock sample is removed;
[0012] According to the pore microstructure type of the limestone, on the basis of statistical analysis of the physical property cross plot of the limestone rock sample after removing the abnormal value, the distribution range of different pore microstructure types is divided;
[0013] Based on the distribution range of different pore microstructure types, the corresponding pore structure type unit is obtained;
[0014] The pore structure type unit is numbered, and the pore structure type and physical property relationship chart of the limestone rock sample is established.
[0015] Further, the abnormal value includes a value less than the test instrument precision threshold.
[0016] Further, based on the physical property curve data of the single well limestone target layer well section logging interpretation and the physical property data of the rock sample, a logging physical property cross plot is established, including:
[0017] The physical property data of the single well limestone target layer well section logging interpretation is obtained;
[0018] Based on the physical property data of the rock sample, the physical property data of the single well limestone target layer well section logging interpretation is corrected, and the corrected single well logging interpretation physical property data is obtained;
[0019] Based on the corrected single well logging interpretation physical property data, a logging physical property cross plot is established.
[0020] Further, the corrected single well logging interpretation physical property data is the physical property data meeting the production application threshold standard and suitable for production application.
[0021] Further, based on the pore structure type result map of the whole well section of the limestone target layer, a single well of the limestone reservoir pore structure type is analyzed, including:
[0022] Based on the core data, logging data and test data, a limestone physical property, pore structure type and productivity contribution expert knowledge database is established;
[0023] The pore structure type result map of the whole well section of the limestone target layer is combined with the limestone physical property, pore structure type and productivity contribution expert knowledge database to analyze the single well of the limestone reservoir pore structure type.
[0024] Further, the limestone physical property, pore structure type and productivity contribution expert knowledge database is established based on the core data, logging data and test data, and the database comprises:
[0025] The pore structure type of the limestone is determined based on the core data;
[0026] The porosity and permeability data of the limestone reservoir are obtained according to the logging data;
[0027] The productivity contribution data are obtained according to the test data; the test data comprise downhole cable formation test and oil test data;
[0028] The limestone reservoir physical property and productivity contribution corresponding pore structure type relationship database is established based on the pore structure type of the limestone, the porosity and permeability data of the limestone reservoir and the productivity contribution data.
[0029] Further, the physical property data comprise porosity and permeability.
[0030] The application further provides a single well analysis device for the pore structure type of the limestone reservoir, and the device comprises,
[0031] A first establishing module is configured to establish a pore structure type and physical property relationship chart of the limestone rock sample based on the porosity and permeability data of the limestone rock sample and the pore microstructure type of the limestone;
[0032] A second establishing module is configured to establish a logging physical property crossplot based on the physical property curve data of the well logging interpretation of the whole well section of the limestone target layer and the physical property data of the rock sample;
[0033] A third establishing module is configured to apply the pore structure type and physical property relationship chart of the limestone rock sample to the logging physical property crossplot to establish a pore structure type result map of the whole well section of the limestone target layer;
[0034] An analysis module is configured to analyze the single well of the limestone reservoir pore structure type based on the pore structure type result map of the whole well section of the limestone target layer.
[0035] Further, the first establishing module performs the following steps:
[0036] The physical property crossplot of the limestone rock sample is constructed based on the porosity and permeability data of the limestone rock sample;
[0037] Eliminate the abnormal value in the physical property cross plot of the limestone rock sample;
[0038] According to the pore microstructure type of the limestone, on the basis of the statistical analysis of the physical property cross plot of the limestone rock sample after eliminating the abnormal value, the distribution range of different pore microstructure types is divided;
[0039] Based on the distribution range of different pore microstructure types, the corresponding pore structure type unit is obtained;
[0040] The pore structure type unit is numbered, and the pore structure type and physical property relationship chart of the limestone rock sample is established.
[0041] The technical effects and advantages of the present application are as follows:
[0042] The present application utilizes core, logging and test data to establish the physical property, pore structure type and productivity contribution database of the limestone rock sample; constructs the porosity-permeability physical property cross plot of the limestone rock sample; the core one is to establish the pore structure type and physical property relationship chart of the limestone reservoir rock sample on this basis; further, the logging physical property interpretation result is corrected according to the physical property data of the limestone rock sample; the core two is to establish the corrected logging porosity and permeability cross plot, and finally, the pore structure physical property chart established by the rock sample is applied to the corrected logging physical property cross plot, and the pore structure type result map of the whole well section of the limestone target layer is established, the single well productivity evaluation can be carried out, and the geological basis for the complex carbonate rock favorable reservoir seismic prediction is provided. Through research and literature search, no method example of constructing the pore structure type chart based on the core physical property to evaluate the pore structure type of the limestone reservoir of a single well has been found in China.
[0043] Other features and advantages of the present application will be set forth in the following description, and in part will become apparent to those skilled in the art from the description, or can be learned by practice of the present application. The objects and other advantages of the present application can be achieved and obtained by the structure indicated in the specification and drawings. BRIEF DESCRIPTION OF DRAWINGS
[0044] Figure 1 is the technical solution of the present application;
[0045] Figure 2 is the single well limestone porosity-permeability cross plot provided by the embodiment of the present application;
[0046] Figure 3 is the single well limestone pore structure type zoning map provided by the embodiment of the present application;
[0047] Figure 4 is the limestone pore structure type zoning chart provided by the embodiment of the present application;
[0048] Figure 5 is a single-well limestone original logging interpretation result graph provided by the embodiment of the present application;
[0049] Figure 6 is a single-well limestone logging interpretation result graph after rock sample data correction provided by the embodiment of the present application;
[0050] Figure 7 is a single-well limestone pore structure type classification result graph provided by the embodiment of the present application. DETAILED DESCRIPTION
[0051] 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 a part of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the protection scope of the present application.
[0052] To solve the problems in the prior art, the present application discloses a single-well analysis method for pore structure types of limestone reservoirs, as shown in Figure 1 The method comprises,
[0053] Step 1: establishing a pore structure type and physical property relationship graph of limestone rock samples based on porosity and permeability data of the limestone rock samples and pore microstructure types of the limestone;
[0054] Step 2: establishing a logging physical property crossplot based on physical property curve data of single-well limestone target layer full well section logging interpretation and physical property data of rock samples;
[0055] Step 3: applying the pore structure type and physical property relationship graph of the limestone rock samples to the logging physical property crossplot to establish a pore structure type result graph of the limestone target layer full well section;
[0056] Step 4: analyzing single-well pore structure types of limestone reservoirs based on the pore structure type result graph of the limestone target layer full well section.
[0057] In one specific embodiment of the present application, Step 1: the establishment of the pore structure type and physical property relationship graph of the limestone rock samples comprises:
[0058] Step 101: constructing a physical property crossplot of limestone rock samples based on porosity and permeability data of the limestone rock samples;
[0059] Step 102: based on the property cross plot, data quality control is carried out, and outliers in the property cross plot of the limestone rock sample are removed; wherein the outliers include values less than the test instrument precision threshold. Specifically, the data value range is checked, and the values less than the test instrument precision are removed; for example, in the embodiment, the values less than 0.001 mD (the precision of the test instrument used is 0.001 mD) are removed.
[0060] Step 103: according to the pore microstructure type of the limestone, on the basis of statistical analysis of the property cross plot of the limestone rock sample after removing the outliers, the distribution range of different pore microstructure types is divided; the pore microstructure type of the limestone includes secondary pore residual intergranular pore, dissolution pore, primary pore, etc.
[0061] Step 104: based on the distribution range of different pore microstructure types, the corresponding pore structure type unit is obtained; specifically, according to the actual production given property threshold value (such as the lower limit value of porosity and permeability given by logging interpretation), and combining the pore structure characteristics to give the division standard.
[0062] Step 105: the pore structure type unit is numbered, and a pore structure type and property relationship chart of the limestone rock sample is established.
[0063] In one specific embodiment of the present application, step 2: based on the property curve data of the single-well limestone target layer full well section logging interpretation and the property data of the rock sample, a logging property cross plot is established, which includes:
[0064] Step 201: according to logging operation, the property data of the single-well limestone target layer full well section logging interpretation is obtained;
[0065] Step 202: based on the property data of the rock sample, the property data of the single-well limestone target layer full well section logging interpretation is corrected, and the corrected single-well logging interpretation property data is obtained; specifically, the property data of the rock sample and the property data of the single-well limestone target layer full well section logging interpretation are fitted according to the actual data to obtain the best relationship, so that the correlation coefficient of the relationship of the two reaches the production application threshold value.
[0066] The core property data and the logging interpretation property data always have a certain deviation. The core property data is used as the standard, the logging interpretation property data is corrected according to the best relationship obtained above, such as a neural network method, to obtain new corrected logging interpretation property data.
[0067] The corrected single-well logging interpretation property data is the property data that meets the production application threshold value standard and is suitable for production application (such as pore structure analysis and productivity evaluation).
[0068] Step 203: based on the corrected single well logging interpretation of physical property data, a logging physical property crossplot is established.
[0069] In one specific embodiment of the present application, step 3: the pore structure type of the limestone rock sample is applied to the logging physical property crossplot to establish a pore structure type result map of the whole well section of the limestone target layer, which can determine the pore structure type of the whole well section.
[0070] Rock pore structure refers to the geometric shape, size, distribution, mutual connection and configuration relationship of pore throat, which is the basis for the study of oil and gas enrichment regularity, reservoir evaluation and development plan compilation, and has very important research value for carbonate reservoirs with pores as the main reservoir space. The carbonate rock deposition process is jointly controlled by biological, chemical and mechanical effects, and after deposition, various diagenesis is carried out, resulting in more complex pore structure and stronger heterogeneity than clastic rock, and especially the complex pore permeability relationship restricted by pore structure often becomes one of the factors restricting the oil and gas exploitation of carbonate reservoirs. Physical properties mainly include porosity and permeability.
[0071] In one specific embodiment of the present application, step 4: based on the pore structure type result map of the whole well section of the limestone target layer, a single well of the limestone reservoir pore structure type is analyzed, including:
[0072] Step 401: based on core data, logging data and test data, a limestone physical property, pore structure type and productivity contribution expert knowledge database is established; including:
[0073] Step 4011: based on core data, the pore structure type of the limestone is determined;
[0074] Step 4012: according to the logging data, the porosity and permeability data of the limestone reservoir are obtained;
[0075] Step 4013: according to the test data, the productivity contribution data are obtained; the test data include downhole cable formation test and oil test data;
[0076] Step 4014: based on the pore structure type of the limestone, the porosity and permeability data of the limestone reservoir and the productivity contribution data, a limestone reservoir physical property, productivity contribution corresponding pore structure type relationship database is established.
[0077] Step 402: the pore structure type result map of the whole well section of the limestone target layer is combined with the limestone physical property, pore structure type and productivity contribution expert knowledge database to analyze a single well of the limestone reservoir pore structure type.
[0078] The application further provides a single-well analysis device for pore structure types of a limestone reservoir, the device comprising,
[0079] a first establishing module configured to establish a pore structure type and physical property relationship chart of the limestone rock sample based on the porosity and permeability data of the limestone rock sample and the pore microstructure type of the limestone;
[0080] a second establishing module configured to establish a logging physical property crossplot based on the physical property curve data of the single-well limestone target layer full well section logging interpretation and the physical property data of the rock sample;
[0081] a third establishing module configured to apply the pore structure type and physical property relationship chart of the limestone rock sample to the logging physical property crossplot to establish a pore structure type result chart of the limestone target layer full well section;
[0082] an analysis module configured to analyze the single well of the limestone reservoir pore structure type based on the pore structure type result chart of the limestone target layer full well section.
[0083] Further, the first establishing module performs the following steps:
[0084] constructing a physical property crossplot of the limestone rock sample based on the porosity and permeability data of the limestone rock sample;
[0085] eliminating abnormal values in the physical property crossplot of the limestone rock sample;
[0086] dividing the distribution range of different pore microstructure types based on the statistical analysis of the physical property crossplot of the limestone rock sample after eliminating abnormal values according to the pore microstructure type of the limestone;
[0087] acquiring corresponding pore structure type units based on the distribution range of different pore microstructure types;
[0088] numbering the pore structure type units to establish the pore structure type and physical property relationship chart of the limestone rock sample.
[0089] As to the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment related to the method, and thus will not be described in detail here.
[0090] The following will be further described with reference to specific embodiments.
[0091] This time, taking a certain limestone oil and gas reservoir as an example, a probe well is selected, which has complete coring data, completes core physical property analysis test, and obtains downhole cable formation test and oil test data, as well as logging data collection.
[0092] (1) Establish a limestone physical property, pore structure type and productivity contribution database, including core thin section, scanning electron microscope, mercury injection method to determine the pore structure type of carbonate rock, rock physical property analysis to obtain porosity and permeability data of limestone reservoir, downhole wireline formation test and oil test data to obtain productivity contribution data. Establish a relationship database of limestone reservoir physical property and productivity contribution corresponding to pore structure type.
[0093] (2) Build a porosity-permeability property intersection graph of limestone rock samples. According to the porosity and permeability data of limestone rock samples, create a porosity-permeability property intersection graph of single well Figure 2 ) and carry out data quality control to eliminate outliers.
[0094] (3) Establish a pore structure type and physical property relationship chart of limestone rock samples. According to the pore microstructure type of limestone, divide the distribution range of different pore structure types based on statistical analysis of the property intersection graph Figure 3 ), the horizontal coordinate is effective porosity (unit: %), the vertical coordinate is permeability, unit: mD; and according to the relationship between physical property and pore structure of the basic database, the range of each pore structure type is divided, and the pore structure type is marked as I, II, III and IV; the divided pore structure type unit is numbered to form a pore structure type identification chart Figure 4 ), which can be used for well logging property data plotting to determine the pore structure type identification based on well logging property data in the layer section without core data.
[0095] (4) Correct the well logging property result graph according to the limestone rock sample physical property data. Obtain the porosity and permeability property curve data of the well logging interpretation of the whole well section of the limestone target layer of single well Figure 5 ); establish the corresponding relationship formula of the well logging interpreted physical property data and the physical property data of the rock sample by neural network method, and correct the well logging interpreted physical property data of single well Figure 6 ) according to this, the correction standard: the correlation coefficient of the relationship formula reaches the limestone oil and gas reservoir production application threshold value of 85% and above.
[0096] (5) Generate single well limestone pore structure type result graph based on pore structure type chart. Establish the corrected well logging porosity and permeability intersection graph; apply the established pore structure type and physical property relationship chart to the corrected well logging porosity and permeability intersection graph; generate the pore structure type result graph of the whole well section of the limestone target layer Figure 7 ), and combined with the basic database, the single well productivity can be evaluated (Table 1), and the productivity contribution ability of each section can be determined.
[0097] Table 1 Corresponding relationship chart of pore structure type and productivity contribution ability
[0098] Number Physical property characteristics Pore structure type Contributing ability to productivity 1 Medium-high porosity and low permeability II Lower 2 Low porosity and medium-low permeability III Higher 3 High porosity and medium-high permeability I High 4 Low-medium-high porosity and low permeability IV Low or none
[0099] According to Table 1, the contribution ability of each section to productivity is clear. In combination with the porosity and permeability of each section, the contribution ability of each section to productivity can be determined. Figure 7 It can be seen that, Figure 7 The section with the pore structure type I in the rightmost column is of high porosity and medium-high permeability, indicating high contribution ability to productivity; the section with the pore structure type II is of medium-high porosity and low permeability, indicating lower contribution ability to productivity; the section with the pore structure type III is of low porosity and medium-low permeability, indicating higher contribution ability to productivity; and the section with the pore structure type IV is of low-high porosity and low permeability, indicating low contribution ability to productivity or no contribution to productivity. Thus, the contribution ability of the whole section of a single well to productivity can be obtained.
[0100] Finally, it should be noted that the above only describes the preferred embodiments of the present application and is not intended to limit the present application. Although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions recorded in the foregoing embodiments or equivalently replace some technical features thereof, and any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims
1. A method of single well analysis of the type of pore structure of a limestone reservoir, characterized in that, The method comprises, based on the porosity, permeability data of the limestone rock sample and the pore microstructure type of the limestone, a pore structure type and property relationship chart of the limestone rock sample is established; based on the property curve data interpreted from the full well section logging of the single well limestone target layer and the property data of the rock sample, a logging property crossplot is established; the pore structure type and property relationship chart of the limestone rock sample is applied to the logging property crossplot to establish a pore structure type result map of the full well section of the limestone target layer; based on the pore structure type result map of the full well section of the limestone target layer, a single well of the limestone reservoir pore structure type is analyzed.
2. The method for single well analysis of pore structure type of limestone reservoir according to claim 1, characterized in that, The establishment of the pore structure type and property relationship chart of the limestone rock sample comprises: based on the porosity, permeability data of the limestone rock sample, a property crossplot of the limestone rock sample is constructed; abnormal values in the property crossplot of the limestone rock sample are eliminated; according to the pore microstructure type of the limestone, on the basis of statistical analysis of the property crossplot of the limestone rock sample after the elimination of abnormal values, the distribution range of different pore microstructure types is divided; based on the distribution range of different pore microstructure types, the corresponding pore structure type unit is obtained; the pore structure type unit is numbered to establish the pore structure type and property relationship chart of the limestone rock sample.
3. The method for single well analysis of pore structure types of limestone reservoirs according to claim 2, characterized in that, The abnormal values include values less than the threshold value of the testing instrument accuracy.
4. The method for single well analysis of pore structure type of limestone reservoir according to claim 1, characterized in that, The establishment of the logging property crossplot based on the property curve data interpreted from the full well section logging of the single well limestone target layer and the property data of the rock sample comprises: the property data interpreted from the full well section logging of the single well limestone target layer is obtained; the property data interpreted from the full well section logging of the single well limestone target layer is corrected based on the property data of the rock sample to obtain corrected property data interpreted from the single well logging; the logging property crossplot is established based on the corrected property data interpreted from the single well logging.
5. The single well analysis method of the limestone reservoir pore structure type according to claim 4, wherein the corrected property data interpreted from the single well logging is property data meeting the production application threshold value standard and suitable for production application.
6. The method for single well analysis of pore structure types of limestone reservoirs according to claim 1, characterized in that, Based on the pore structure type result map of the full well section of the limestone target layer, the single well of the limestone reservoir pore structure type is analyzed, comprising: based on the core data, logging data and test data, an expert knowledge database of limestone properties, pore structure types and productivity contribution is established; the pore structure type result map of the full well section of the limestone target layer is combined with the expert knowledge database of limestone properties, pore structure types and productivity contribution to analyze the single well of the limestone reservoir pore structure type.
7. The method for single well analysis of pore structure types of limestone reservoirs according to claim 6, characterized in that, The establishment of the expert knowledge database of limestone properties, pore structure types and productivity contribution based on the core data, logging data and test data comprises: based on the core data, the pore structure type of the limestone is determined; according to the logging data, the porosity, permeability data of the limestone reservoir are obtained; according to the test data, productivity contribution data are obtained; the test data include downhole cable formation testing and oil testing data; The base stone is limestone, the pore structure type of the limestone, the porosity and permeability data of the limestone reservoir, and the productivity contribution data are used to establish the relationship database of the physical property of the limestone reservoir and the pore structure type corresponding to the productivity contribution.
8. The single-well analysis method for the pore structure type of a limestone reservoir according to any one of claims 1-7, characterized in that, The physical property data includes porosity and permeability.
9. A single well analysis device for a limestone reservoir pore structure type, characterized by, The device comprises, A first establishing module is configured to establish a pore structure type and physical property relationship chart of the limestone rock sample based on the porosity and permeability data of the limestone rock sample and the pore microstructure type of the limestone. A second establishing module is configured to establish a logging physical property crossplot based on the physical property curve data interpreted from the full well section logging of the single-well limestone target layer and the physical property data of the rock sample. A third establishing module is configured to apply the pore structure type and physical property relationship chart of the limestone rock sample to the logging physical property crossplot to establish a pore structure type result map of the full well section of the limestone target layer. An analyzing module is configured to analyze the single-well limestone reservoir pore structure type based on the pore structure type result map of the full well section of the limestone target layer.
10. The apparatus for single well analysis of pore structure types of limestone reservoirs according to claim 9, characterized in that, The first establishing module performs the following steps: Based on the porosity and permeability data of the limestone rock sample, a physical property crossplot of the limestone rock sample is constructed. Anomalous values in the physical property crossplot of the limestone rock sample are removed. According to the pore microstructure type of the limestone, the distribution range of different pore microstructure types is divided based on the statistical analysis of the physical property crossplot of the limestone rock sample after removing the anomalous values. Based on the distribution range of different pore microstructure types, the corresponding pore structure type unit is obtained. The pore structure type unit is numbered to establish the pore structure type and physical property relationship chart of the limestone rock sample.