CN112651547B - Method for predicting size of sandstone reservoir of curtreta - Google Patents

Method for predicting size of sandstone reservoir of curtreta Download PDF

Info

Publication number
CN112651547B
CN112651547B CN202011460495.8A CN202011460495A CN112651547B CN 112651547 B CN112651547 B CN 112651547B CN 202011460495 A CN202011460495 A CN 202011460495A CN 112651547 B CN112651547 B CN 112651547B
Authority
CN
China
Prior art keywords
river
sand
concave
sand body
drilling
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011460495.8A
Other languages
Chinese (zh)
Other versions
CN112651547A (en
Inventor
蒙晓灵
兰义飞
杨勇
刘海锋
王金成
王准备
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202011460495.8A priority Critical patent/CN112651547B/en
Publication of CN112651547A publication Critical patent/CN112651547A/en
Application granted granted Critical
Publication of CN112651547B publication Critical patent/CN112651547B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Husbandry (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Game Theory and Decision Science (AREA)
  • Agronomy & Crop Science (AREA)
  • Development Economics (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method for predicting the scale of a sandstone reservoir of a curtret river, which comprises the following steps: s001: selecting and sampling each drilling well in a curved flow reach of river field, and counting the interval between each drilling well and the thickness of the sand body of each drilling well; s002: calculating the sand body transfer rate of two adjacent drilling wells in sequence from west to east, and analyzing the value, wherein when the thickness of the sand body is thickened, the sand body transfer rate is positive value, the sediment of the concave bank of the curved flow river can be reflected, and when the thickness of the sand body is thinned, the sand body transfer rate is negative value, the sediment of the convex bank of the curved flow river can be reflected; s003: judging the cone-shaped river concave bank run sand cone-shaped gradient rate and the cone-shaped river convex bank run sand cone-shaped gradient rate according to the sand cone gradient rate of each well drilling; s004: respectively calculating the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river; s005: and calculating the river scale according to the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river.

Description

Method for predicting size of sandstone reservoir of curtreta
Technical Field
The invention belongs to the technical field of geological research of oil and gas field exploitation, and particularly relates to a method for predicting the scale of a sandstone reservoir of a curtreta river.
Background
The scale prediction of the curvelet sand body is mainly based on the prediction of the seismic data at home and abroad, and the quality of the seismic data is relatively good in the northern desert area of the Erdos basin, so that a good resolution effect can be obtained. However, for loess tableland features in the south of the basin, the seismic data have been difficult to collect and analyze. The loess plateau landform seismic data has poor data quality due to thicker loess layers, and is an important technical attack direction in the aspects of accuracy and precision of sand body identification. Especially for deep burial (> 4500 m) of a similar Qingyang gas field, the sandstone reservoir of a thin reservoir (4-8 m) greatly increases the difficulty of predicting the reservoir scale by seismic data, and new sandstone reservoir scale prediction methods and means are needed.
Disclosure of Invention
The present invention provides a method of curvedly sandstone reservoir scale prediction for overcoming, or at least partially solving or alleviating, the above problems.
Therefore, the invention provides a method for predicting the scale of a sandstone reservoir of a curtrey river, which is characterized by comprising the following steps of:
S001: selecting and sampling each drilling well in a curved flow reach of river field, and counting the interval between each drilling well and the thickness of the sand body encountered by the drilling well;
S002: calculating the sand body transfer rate of two adjacent drilling wells in sequence from west to east, and analyzing the value, wherein when the thickness of the sand body is thickened and the sand body transfer rate is positive, the sediment of the concave bank of the curved stream river is reflected, and when the thickness of the sand body is thinned and the sand body transfer rate is negative, the sediment of the convex bank of the curved stream river is reflected;
s003: judging the cone-shaped river concave bank run sand cone-shaped gradient rate and the cone-shaped river convex bank run sand cone-shaped gradient rate according to the sand cone gradient rate of each well drilling;
s004: respectively calculating the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river;
S005: and calculating the river scale according to the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river.
In step S001, a sample point (finished well) W 1......Wn is screened using the prescribed well spacing 0.6Km as a basis discussion.
In step S002, the sand gradient K H refers to the change in sand thickness per unit distance; the formula of K H is: k H=(Hwn-Hwn-1)/S;
Wherein K H is the thickness gradient of the sand body, and the unit is m/km; wn is the sample point (well completion); hw n is the thickness of the sand body of a certain well point drill meeting mountain, and the unit is m; s is the distance between two adjacent wells.
In step S002, the side erosion of the concave side is the main part of the river channel sand body of the cursive river, the development of the convex side is the beach deposition, the deposited sand body has two asymmetric characteristics, on the concave side of the cursive river, K H is more than 0, on the convex side, K H is less than 0, namely, the river channel is a complete river channel of one period, the river sand body of the cursive river is the superposition of the side of the multi-period river channel, and the K H value shows regular change.
In step S003, the curvelet sand body is deposited by lateral overlapping of a multi-stage river channel, and on the basis of identifying each stage river channel, the gradient K H Concave recess of the curvelet sand body in the area is determined to be between 8m/km and 23m/km according to analysis of a plurality of concave bank deposition sample points of the Qingyang gas field 100 and a real drilling data large database.
In step S004, a K H Concave recess value is determined according to a land sediment sample point of the gas field in the target area and the maximum value of the thickness of the solid drilling sand body, and the maximum extension distance of the land of the curved flow river is calculated according to a K H formula and a K H Concave recess value range.
In step S004, an average value of K H Convex part values is determined according to all land deposition sample points of the gas field in the target area, and the maximum extension distance of the lands of the curved stream is calculated according to the formula K H and the average value of K H Convex part values.
In step S005, the size of the river channel is the sum of the maximum extension distance of the concave side of the river and the maximum extension distance of the convex side of the river.
The invention aims to solve the problems, analyzes the sediment characteristics of the sediment sandstone reservoir of the area by introducing the concept of the ' sand thickness gradient rate ', and establishes a set of curved river sand scale prediction modes suitable for the terrains of eastern and southern loess terrains of the basin according to the ' asymmetry of the two sides of the sediment of the curved river by perfecting the sand prediction method, thereby improving the prediction precision of the scale of the block sand in the area and becoming one of powerful supports for guiding the effective deployment and scale development of gas fields.
Drawings
FIG. 1 is a partial distribution diagram of sample points in the present invention;
FIG. 2 is a partial identification view of a primary channel in which a portion of the sample points in FIG. 1 are located;
FIG. 3 is a schematic diagram of a complete river at which the sample points of FIG. 1 are located;
FIG. 4 is a schematic cross-sectional view of a multiple-stage sand deposit on a side beach of a curveflow;
Fig. 5 is a plan view of the curveflow river side beach sand deposit in fig. 4.
Detailed Description
The invention is further illustrated by a specific example.
Example 1
In the method, a main geological target mountain section of the Qingyang gas field is a deposition area of a curvelan, and the following steps are adopted:
S001: selecting and sampling each drilling well in a curved flow reach of river field, and counting the interval between each drilling well and the thickness of the sand body encountered by the drilling well;
S002: calculating the sand body transfer rate of two adjacent drilling wells in sequence from west to east, and analyzing the value, wherein when the thickness of the sand body is thickened and the sand body transfer rate is positive, the sediment of the concave bank of the curved stream river is reflected, and when the thickness of the sand body is thinned and the sand body transfer rate is negative, the sediment of the convex bank of the curved stream river is reflected;
s003: judging the cone-shaped river concave bank run sand cone-shaped gradient rate and the cone-shaped river convex bank run sand cone-shaped gradient rate according to the sand cone gradient rate of each well drilling;
s004: respectively calculating the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river;
S005: and calculating the river scale according to the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river.
Example 2
S001: selecting and sampling each drilling well in a curved flow reach of river field, and counting the interval between each drilling well and the thickness of the sand body encountered by the drilling well;
In step S001, the drill is completed by analyzing the data of the previous evaluation block and using the well spacing of 0.6Km as the basis for the ancient world on the erdos basin to screen the sample points, i.e. the drill W 1......Wn.
Example 3
S002: and calculating the sand body transition rate of two adjacent wells in sequence from west to east, analyzing the value, and reflecting the concave sides of the curved flow river when the thickness of the sand body is thickened and reflecting the convex sides of the curved flow river when the thickness of the sand body is thinned.
In step S002, the sand gradient K H refers to the change in sand thickness per unit distance; the formula of K H is: k H=(Hwn-Hwn-1)/S;
Wherein K H is the thickness gradient of the sand body, and the unit is m/km; wn is the sample point, and the well is completed; hw n is the thickness of the sand body of a certain well point drill meeting mountain, and the unit is m; s is the distance between two adjacent wells.
As shown in fig. 1, statistics of intervals between the sample points and thicknesses of sand bodies encountered by each well drilling are completed according to the selected well completion data, and the sand body gradient rate K H of two adjacent points is calculated from west to east in sequence according to a sand body gradient rate formula, and is analyzed, and the analysis results are shown in table 1.
TABLE 1 Sand thickness transition rate at actual drilling points
Example 4
In step S002, the side erosion of the concave side is the main part of the river channel sand body of the cursive river, the development of the convex side is the beach deposition, the deposited sand body has two asymmetric characteristics, on the concave side of the cursive river, K H is more than 0, on the convex side, K H is less than 0, namely, the river channel is a complete river channel of one period, the river sand body of the cursive river is the superposition of the side of the multi-period river channel, and the K H value shows regular change.
Referring to fig. 2 and 3, k H >0, the thickness of the sand body is thickened, which can reflect the concave sides of the curved stream river; k H is smaller than 0, and the thickness of the sand body is reduced, so that the convex bank of the curved stream river can be reflected. When K H is more than 0 and K H is less than 0, the sand body is thickened and then thinned, and the real drilling data show that the process describes a first-stage complete river channel.
Example 5
In step S003, the curvelet sand body is deposited by lateral overlapping of a multi-stage river channel, and on the basis of identifying each stage river channel, the gradient K H Concave recess of the curvelet sand body in the area is determined to be between 8m/km and 23m/km according to analysis of a plurality of concave bank deposition sample points of the Qingyang gas field 100 and a real drilling data large database.
Referring to fig. 1, 4 and 5, since the concave beach sand deposit of the curveflow river is divided into a plurality of periods, the thickness of the sand deposited from the first period of the beach sand deposit to the last period of the beach sand deposit is deeper and deeper, and the value of the gradient rate K H of the sand is larger and larger, when the value of K H becomes negative, this is indicated as the convex side of the curveflow river, and when the value of K H becomes positive again, this is indicated as the concave side of the previous period, but on the whole is also on the convex side of the curveflow river. In fig. 4, K H Convex part 1、KH Convex part 2、KH Convex part 3 and K H Convex part 4 refer to the rate of the body transition of the asperity bank run sand at different times, respectively.
And determining the integral sand thickness gradient rate K H Concave recess of the concave bank of the curvata according to data of more than 100 curvata sediment sample points of the Qingyang gas field and experience.
Example 6
S004: respectively calculating the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river;
In step S004, a K H Concave recess value is determined according to a land sediment sample point of the gas field in the target area and the maximum value of the thickness of the solid drilling sand body, and the maximum extension distance of the land of the curved flow river is calculated according to a K H formula and a K H Concave recess value range.
According to actual measurement, obtaining that the maximum thickness of the curved stream river sand body is 13.5m (14 m), namely Hw n≤14m,Hwn-1 approaches to or more than 0, substituting the limit value of Hw n and the limit value of Hw n-1 into K H=(Hwn-Hwn-1)/S, namely K H Concave recess (14 m-0 m)/0.6 km, obtaining the K H Concave recess value which is less than or equal to 23m/km, and substituting K H Concave recess into K H=(Hwn-Hwn-1)/S when K H Concave recess E [8, 23] is determined according to real-drilling data analysis, wherein S is the extension distance of the concave bank of the curved stream river; s epsilon 0,1.75 is obtained reversely, namely the maximum extension distance of the concave bank is 1.75km.
Example 7
In step S004, an average value of the K H Convex part values is determined according to the land deposition sample points of all the periods of the gas field in the target area, and the maximum extension distance of the land of the curved stream river is calculated according to the K H formula and the average value of the K H Convex part values.
As shown in fig. 1, according to the real drilling data of a plurality of deposited sample points of the curved stream, the sample points cover at least two periods of complete river channels, the average value of K H Convex part is calculated according to the data in table 1, the average value of K H Convex part is (-0.5 m/km-5.2 m/km-1.6/km)/3= -2.43m/km, and the average value of K H Convex part and Hw n are substituted into K H=(Hwn-Hwn-1)/S, wherein S is the extension distance of the convex bank of the curved stream; i.e., -2.43m/km = (14 m-0 m)/S, S = -5.76km is obtained, the minus sign represents the direction of the curved stream river land direction, i.e., the maximum distance of extension of the curved stream river land is 5.76km.
Example 8
In step S005, the scale of the river channel of the curved stream is the sum of the maximum extension distance of the concave side of the curved stream and the maximum extension distance of the convex side of the curved stream, i.e. 1.75km+5.76km is about 7.5km, the eastern and western maximum width of the curved stream in this period is 7.5km, and the well position deployment is guided according to the maximum width of the curved stream, so that the risk is reduced.
The present invention is not limited to the above-mentioned embodiments, and any changes or substitutions that can be easily understood by those skilled in the art within the technical scope of the present invention are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims. The components and structures not specifically described in this embodiment are well known in the art and are not described in detail herein.

Claims (2)

1. The method for predicting the size of the sandstone reservoir of the curtret river is characterized by comprising the following steps of:
S001: selecting and sampling each drilling well in a curved flow reach of river field, and counting the interval between each drilling well and the thickness of sand encountered by the drilling well;
S002: calculating the sand body transfer rate of two adjacent drilling wells in sequence from west to east, and analyzing the value, wherein when the thickness of the sand body is thickened and the sand body transfer rate is positive, the sediment of the concave bank of the curved stream river is reflected, and when the thickness of the sand body is thinned and the sand body transfer rate is negative, the sediment of the convex bank of the curved stream river is reflected;
s003: judging the cone-shaped river concave bank run sand cone-shaped gradient rate and the cone-shaped river convex bank run sand cone-shaped gradient rate according to the sand cone gradient rate of each well drilling;
s004: respectively calculating the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river;
s005: calculating the river scale according to the maximum extension distance of the concave sides of the curveflow river and the maximum extension distance of the convex sides of the curveflow river;
in step S001, screening a sample point, namely a drilling W 1 ......Wn, by adopting a specified well spacing of 0.6 Km;
In step S002, the sand gradient K H refers to the change in sand thickness per unit distance; the formula of K H is: k H=(Hwn-Hwn-1)/S;
Wherein K H is the thickness gradient of the sand body, and the unit is m/km; w n is the drilling sample point; hw n is the thickness of the sand body in the drilling process, and the unit is m; s is the distance between two adjacent wells;
In the step S002, the side erosion of the concave side is the main part, the development side beach of the convex side is deposited, the deposited sand has two-side asymmetric characteristics, K H is more than 0 on the concave side of the curved river, K H is less than 0 on the convex side, namely, a first-stage complete river is represented, the curved river sand is overlapped on the side of the multi-stage river, and the K H value shows regular change;
In the step S003, the curvelet sand body is laterally overlapped and deposited in a multi-stage river channel, and on the basis of identifying each stage of river channel, the transfer rate K H Concave recess of the curvelet sand body in the area is determined according to analysis of a plurality of concave bank deposition sample points of the Qingyang gas field 100 and a real drilling data large database;
in step S004, determining a K H Concave recess value according to a concave bank deposition sample point of the gas field in the target area and an actually measured drilling and sand body thickness maximum value, and calculating a maximum extension distance of the concave bank of the curved flow river according to a K H formula and a K H Concave recess value range;
In step S004, a K H Convex part value average value is determined according to the land sediment sample point of the gas field in the target area, and the maximum extension distance of the land of the curved stream is calculated according to a K H formula and the K H Convex part value average value.
2. The method according to claim 1, wherein in step S005, the size of the river channel of the curveway is a sum of a maximum extension distance of the concave side of the curveway and a maximum extension distance of the convex side of the curveway.
CN202011460495.8A 2020-12-11 2020-12-11 Method for predicting size of sandstone reservoir of curtreta Active CN112651547B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011460495.8A CN112651547B (en) 2020-12-11 2020-12-11 Method for predicting size of sandstone reservoir of curtreta

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011460495.8A CN112651547B (en) 2020-12-11 2020-12-11 Method for predicting size of sandstone reservoir of curtreta

Publications (2)

Publication Number Publication Date
CN112651547A CN112651547A (en) 2021-04-13
CN112651547B true CN112651547B (en) 2024-05-28

Family

ID=75354819

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011460495.8A Active CN112651547B (en) 2020-12-11 2020-12-11 Method for predicting size of sandstone reservoir of curtreta

Country Status (1)

Country Link
CN (1) CN112651547B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104453877A (en) * 2014-11-05 2015-03-25 长江大学 Underground deep-buried meander point dam sand body historical reconstruction method
CN105607146A (en) * 2015-09-10 2016-05-25 中国海洋石油总公司 Quantitative characterization method of meandering river sand body scale
CN107918150A (en) * 2016-10-10 2018-04-17 中国石油化工股份有限公司 Based on single river recognition methods under offshore oilfield wide spacing
CN109736773A (en) * 2018-11-23 2019-05-10 中国石油天然气股份有限公司 Track tracking method for river sand horizontal well
CN111663932A (en) * 2019-03-06 2020-09-15 中国石油天然气股份有限公司 Method and apparatus for determining subsurface reservoir structure

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104453877A (en) * 2014-11-05 2015-03-25 长江大学 Underground deep-buried meander point dam sand body historical reconstruction method
CN105607146A (en) * 2015-09-10 2016-05-25 中国海洋石油总公司 Quantitative characterization method of meandering river sand body scale
CN107918150A (en) * 2016-10-10 2018-04-17 中国石油化工股份有限公司 Based on single river recognition methods under offshore oilfield wide spacing
CN109736773A (en) * 2018-11-23 2019-05-10 中国石油天然气股份有限公司 Track tracking method for river sand horizontal well
CN111663932A (en) * 2019-03-06 2020-09-15 中国石油天然气股份有限公司 Method and apparatus for determining subsurface reservoir structure

Also Published As

Publication number Publication date
CN112651547A (en) 2021-04-13

Similar Documents

Publication Publication Date Title
CN108010103B (en) Rapid and fine generation method of complex river terrain
CN101697229B (en) Method for extracting region of interest of medical image
CN103616731B (en) Method and device for determining altered volcanic rock effective reservoir in oil and gas exploration
CN103399346B (en) A kind of well shake associating impedance initial value modeling method
CN104237964A (en) Paleochannel type sandstone-type uranium ore deposit positioning method
CN111080021B (en) Sand body configuration CMM neural network prediction method based on geological information base
CN111399056B (en) Method for predicting crack strength based on divided azimuth filtering
CN113608278B (en) Sandstone-type uranium ore positioning method in red heterolayer of sedimentary basin
CN107193046A (en) A kind of sand thickness Forecasting Methodology and system based on spectrum inverting
CN108490491A (en) A kind of beach body prediction technique indicating inverting based on waveform
CN112651547B (en) Method for predicting size of sandstone reservoir of curtreta
CN110428497B (en) Braided river training image generation method
CN109239783B (en) Method for accepting or rejecting earthquake information in well earthquake combination
CN112782760B (en) Method for dissecting braided river reservoir structure by using discontinuous boundaries of seismic reservoir
CN104502962A (en) Method and device for designing shot points
WO2023004530A1 (en) Method for locating sandstone-type uranium deposits in red variegated layer in a sedimentary basin
CN111159869B (en) Automatic calculation method for quantitative parameters between sandstone uranium deposit and ancient river course
CN111997604B (en) Single-stage river boundary identification method
CN117518257A (en) Method for determining top and bottom surfaces of reservoir by applying seismic attributes
CN112198554B (en) Method for establishing high-precision inversion initial model driven by seismic waveform
CN110795513A (en) Method for predicting distribution of river facies source storage ectopic type compact oil gas dessert area
CN114624775B (en) Comprehensive quantitative restoration method for ancient water depth of sedimentary lake basin
CN114594518B (en) Fine stratum contrast method for complex fault blocks in later development period based on well-seismic alternation
CN114660658A (en) Prediction method for constructing soft coal distribution area
CN110320572B (en) Method and system for identifying sedimentary facies

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant