CN103105629A - Equal apparent resistivity surface diagram electromagnetic prospecting data interpretation method - Google Patents
Equal apparent resistivity surface diagram electromagnetic prospecting data interpretation method Download PDFInfo
- Publication number
- CN103105629A CN103105629A CN 201310038374 CN201310038374A CN103105629A CN 103105629 A CN103105629 A CN 103105629A CN 201310038374 CN201310038374 CN 201310038374 CN 201310038374 A CN201310038374 A CN 201310038374A CN 103105629 A CN103105629 A CN 103105629A
- Authority
- CN
- China
- Prior art keywords
- elevation
- apparent resistivity
- curved surface
- autocad
- equal
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000010586 diagram Methods 0.000 title abstract 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 239000003245 coal Substances 0.000 claims description 4
- JTJMJGYZQZDUJJ-UHFFFAOYSA-N phencyclidine Chemical compound C1CCCCN1C1(C=2C=CC=CC=2)CCCCC1 JTJMJGYZQZDUJJ-UHFFFAOYSA-N 0.000 claims description 3
- 239000004215 Carbon black (E152) Substances 0.000 claims description 2
- 229930195733 hydrocarbon Natural products 0.000 claims description 2
- 150000002430 hydrocarbons Chemical class 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 230000007704 transition Effects 0.000 abstract 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
Images
Abstract
The invention discloses an equal apparent resistivity surface diagram electromagnetic prospecting data interpretation method which utilizes an equal apparent resistivity surface diagram to represent electrical characteristics of the stratigraphic position or a specified elevation curved surface. The method comprises the following steps: the stratigraphic position and the specified elevation curved surface of the equal apparent resistivity surface diagram are drawn according to requirements, and an electric map with a triangulation network is generated; transition of a file format of a map is conducted, and then curved surface elevation of an electric file is read according to measure point coordinates; and apparent resistivity of each measure point elevation is captured in apparent resistivity-depth data by utilizing obtained curved surface elevation, and apparent resistivity surface diagram data are formed. The equal apparent resistivity surface diagram electromagnetic prospecting data interpretation method has the advantages of adapting to three-dimensional property of a real geologic structure, avoiding position interpreting errors by unusual judgment, improving correctness of a geology conclusion and being applied to electromagnetic prospecting data interpretation.
Description
Technical field
The invention belongs to the geophysical survey field, be specifically related to a kind of electricity and electromagnetic prospecting interpretation methods.
Background technology
At CSAMT(Controlled Source Audio-frequency Magneto-Telluric), MT(Magneto-Telluric), TEM(Transient Electro-Magnetic) etc. in the electromagnetic prospecting engineering, the equal electrical resistivity planimetric map of making certain degree of depth is that data is processed the basic map of explaining
[1-3], be the important component part of exploration results report.Three-dimensional nature due to true geologic structure, process in the electromagnetic survey data and explain that midplane figure can not show the electrical of specific stratum well, being rich in water or adopting electrical property features such as resistivity distribution on the elevation curved surface of sky situation, the electrical situation of oil storage gas bearing reservoir, appointment as the coal seam, cause an abnormal layer misconstruction that judges, affect the correctness of geology conclusion.
Documents
[1] State Economic and Trade Commission. controlled source audio-frequency magnetotelluric method (CSAMT) exploration engineering rules [SY-T5772-2002 the People's Republic of China (PRC) oil and gas industry standard] .2002-05-28 issue, 2002-08-01 implements
[2] Minitry of Geology and Mineral Resources of the People's Republic of China (PRC). magnetotelluric sounding method technical regulation [DZ/T0173-1997 the People's Republic of China (PRC) geological and mineral industry standard] .1997-03-19 issue, 1997-11-01 implements
[3] Minitry of Geology and Mineral Resources of the People's Republic of China (PRC). ground transient electromagnetic method technical regulation [DZ/T0187-1997 the People's Republic of China (PRC) geological and mineral industry standard] .1997-07-1 issue, 1998-01-15 implements
Summary of the invention
Affect the defective of the correctness of geology conclusion for the layer bit-errors that overcomes abnormal judgement in prior art, the invention provides a kind of equal electrical resistivity surface chart electromagnetic survey interpretation methods, to improve the accuracy of electromagnetic survey data interpretation.
In order to solve above technical matters, the technical solution adopted in the present invention is as follows:
A kind of equal electrical resistivity surface chart electromagnetic survey interpretation methods, the electrical property feature on the elevation curved surface that it is characterized in that utilizing the equal electrical resistivity surface chart to represent stratigraphic horizon or appointment specifically comprises the following steps:
Step 1 is drawn the stratigraphic horizons such as coal seam, hydrocarbon zone of equal electrical resistivity surface chart, the elevation curved surface of appointment as required, prepares corresponding AutoCAD electronics map, keeps the triangulation network in AutoCAD digital map process;
Step 2 utilizes the Conversion of Spatial Data disposal system FME10 of Canadian Safe Software company exploitation to carry out the conversion of AutoCAD map file format, then reads the curved surface elevation of AutoCAD e-file according to the measuring point coordinate;
Step 3, the apparent resistivity that utilizes the curved surface elevation that obtains to intercept each measuring point elevation in apparent resistivity-depth data forms apparent resistivity surface chart data, the intercept method of intercepting surface chart data is for replacing the degree of depth of planimetric map with the elevation of each measuring point, obtain the apparent resistivity value of certain degree of depth by linear interpolation, then utilize the mapping softwares such as AutoCAD or Surfer to generate the equal electrical resistivity surface chart.
The present invention has beneficial effect.The present invention utilizes the equal electrical resistivity surface chart to represent electrical property feature on the elevation curved surface of stratigraphic horizon or appointment, the three-dimensional nature that has adapted to true geologic structure, can show well the electrical of specific stratum, being rich in water or adopting electrical property features such as resistivity distribution on the elevation curved surface of sky situation, the electrical situation of oil storage gas bearing reservoir, appointment as the coal seam, avoid a layer misconstruction of abnormal judgement, improved the correctness of geology conclusion.
Description of drawings
Fig. 1 is seat earth equal electrical resistivity surface chart, and wherein (a) surveys district's seat earth elevation surface chart, (b) is the projection of equal electrical resistivity on the elevation curved surface.
Embodiment
Below in conjunction with accompanying drawing, technical scheme of the present invention is described in further detail.
With in transient electromagnetic method exploration, specify seat earth the apparent resistivity surface chart be plotted as example.
The seat earth altimetric data of district's boring is surveyed in step 1, utilization, measuring point coordinate listed according to table 1, and generating the AutoCAD coal-seam floor contour map is elevation map, as shown in (a) in Fig. 1;
Table 1 measuring point longitude and latitude Excel file layout
| ? | A | B | C |
| 1 | X | Y | Z |
| 2 | 952800 | 97000 | ? |
| 3 | 952850 | 97000 | ? |
| 4 | 952900 | 97000 | ? |
| 5 | 952950 | 97000 | ? |
| 6 | 953000 | 97000 | ? |
| 7 | 952800 | 97050 | ? |
| 8 | 952850 | 97050 | ? |
| 9 | 952900 | 97050 | ? |
| 10 | 952950 | 97050 | ? |
| 11 | 953000 | 97050 | ? |
| 12 | 952800 | 97100 | ? |
| 13 | 952850 | 97100 | ? |
| 14 | 952900 | 97100 | ? |
| 15 | 952950 | 97100 | ? |
| 16 | 953000 | 97100 | ? |
| 17 | 952800 | 97150 | ? |
| 18 | 952850 | 97150 | ? |
| 19 | 952900 | 97150 | ? |
| 20 | 952950 | 97150 | ? |
| 21 | 953000 | 97150 | ? |
| 22 | 952800 | 97200 | ? |
| 23 | 952850 | 97200 | ? |
| 24 | 952900 | 97200 | ? |
| 25 | 952950 | 97200 | ? |
| 26 | 953000 | 97200 | ? |
* in showing, X and Y represent respectively the longitude and latitude of measuring point.
Step 2, utilize the Conversion of Spatial Data disposal system FME10(Feature Manipulate Engineering of Canadian Safe Software company exploitation), carry out the conversion of AutoCAD map file format according to prompting, the measuring point value of contour that reads is put under the Z sign of table 1 C row, separately made table 2;
Table 2 reads the Excel file after the measuring point elevation
| ? | A | B | C |
| 1 | X | Y | Z |
| 2 | 952800 | 97000 | 1181.81 |
| 3 | 952850 | 97000 | 1175.6 |
| 4 | 952900 | 97000 | 1142.76 |
| 5 | 952950 | 97000 | 1128.98 |
| 6 | 953000 | 97000 | 1127.56 |
| 7 | 952800 | 97050 | 1123.18 |
| 8 | 952850 | 97050 | 1123.69 |
| 9 | 952900 | 97050 | 1149.66 |
| 10 | 952950 | 97050 | 1179.87 |
| 11 | 953000 | 97050 | 1186.42 |
| 12 | 952800 | 97100 | 1155.76 |
| 13 | 952850 | 97100 | 1153.97 |
| 14 | 952900 | 97100 | 1129.25 |
| 15 | 952950 | 97100 | 1116.15 |
| 16 | 953000 | 97100 | 1115.36 |
| 17 | 952800 | 97150 | 1131.99 |
| 18 | 952850 | 97150 | 1164.01 |
| 19 | 952900 | 97150 | 1173.67 |
| 20 | 952950 | 97150 | 1181.83 |
| 21 | 953000 | 97150 | 1188.71 |
| 22 | 952800 | 97200 | 1116.05 |
| 23 | 952850 | 97200 | 1113.79 |
| 24 | 952900 | 97200 | 1109.98 |
| 25 | 952950 | 97200 | 1109.43 |
| 26 | 953000 | 97200 | 1139.44 |
* in showing, X, Y and Z represent respectively longitude and latitude and the elevation of measuring point.
The stratigraphic horizon curved surface elevation that obtains in step 3, use table 2, intercepting surface chart data in apparent resistivity-depth data; Intercept method is identical with the apparent resistivity planimetric map, namely replace the degree of depth of planimetric map with the elevation of each measuring point, obtain the apparent resistivity value of certain degree of depth by linear interpolation, then utilize the Surfer mapping software to generate the equal electrical resistivity surface chart, as shown in (b) in Fig. 1, expression be equal electrical resistivity projective distribution on (a) curved surface in Fig. 1.
Claims (2)
1. equal electrical resistivity surface chart electromagnetic survey interpretation methods, the electrical property feature on the elevation curved surface that it is characterized in that utilizing the equal electrical resistivity surface chart to represent stratigraphic horizon or appointment specifically comprises the following steps:
Step 1 is drawn the stratigraphic horizons such as coal seam, hydrocarbon zone of equal electrical resistivity surface chart, the elevation curved surface of appointment as required, prepares corresponding AutoCAD electronics map, keeps the triangulation network in AutoCAD digital map process;
Step 2 utilizes the Conversion of Spatial Data disposal system FME10 of Canadian Safe Software company exploitation to carry out the conversion of AutoCAD map file format, then reads the curved surface elevation of AutoCAD e-file according to the measuring point coordinate;
Step 3, the apparent resistivity that utilizes the curved surface elevation that obtains to intercept each measuring point elevation in apparent resistivity-depth data forms apparent resistivity surface chart data.
2. a kind of equal electrical resistivity surface chart electromagnetic survey interpretation methods as claimed in claim 1, it is characterized in that the intercept method of intercepting surface chart data in described step 3 is for replacing the degree of depth of planimetric map with the elevation of each measuring point, obtain the apparent resistivity value of certain degree of depth by linear interpolation, then utilize the mapping softwares such as AutoCAD or Surfer to generate the equal electrical resistivity surface chart.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310038374.8A CN103105629B (en) | 2013-01-31 | 2013-01-31 | A kind of equal electrical resistivity surface chart electromagnetic survey interpretation methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310038374.8A CN103105629B (en) | 2013-01-31 | 2013-01-31 | A kind of equal electrical resistivity surface chart electromagnetic survey interpretation methods |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103105629A true CN103105629A (en) | 2013-05-15 |
| CN103105629B CN103105629B (en) | 2015-08-19 |
Family
ID=48313623
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310038374.8A Expired - Fee Related CN103105629B (en) | 2013-01-31 | 2013-01-31 | A kind of equal electrical resistivity surface chart electromagnetic survey interpretation methods |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103105629B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104965229A (en) * | 2015-06-10 | 2015-10-07 | 山东科技大学 | Method used for determining range of old mined-out area of coal mine |
| CN107944103A (en) * | 2017-11-14 | 2018-04-20 | 中国煤炭地质总局水文地质局 | A kind of seat earth elevation isopleth Height attribute extracting method |
| CN113960684A (en) * | 2021-09-29 | 2022-01-21 | 江苏大学 | Apparent resistivity-depth profile generation method for short offset electromagnetic exploration |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO329836B1 (en) * | 2008-07-07 | 2011-01-03 | Advanced Hydrocarbon Mapping As | Method for Transforming and Imaging Electromagnetic Exploration Data for Submarine Hydrocarbon Reservoirs |
| CN102508311B (en) * | 2011-11-30 | 2013-10-16 | 河北煤炭科学研究院 | Tunnel advanced detection data multi-parameter spatial mapping method |
-
2013
- 2013-01-31 CN CN201310038374.8A patent/CN103105629B/en not_active Expired - Fee Related
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104965229A (en) * | 2015-06-10 | 2015-10-07 | 山东科技大学 | Method used for determining range of old mined-out area of coal mine |
| CN107944103A (en) * | 2017-11-14 | 2018-04-20 | 中国煤炭地质总局水文地质局 | A kind of seat earth elevation isopleth Height attribute extracting method |
| CN113960684A (en) * | 2021-09-29 | 2022-01-21 | 江苏大学 | Apparent resistivity-depth profile generation method for short offset electromagnetic exploration |
| CN113960684B (en) * | 2021-09-29 | 2024-03-19 | 江苏大学 | Apparent resistivity-depth profile generation method for short offset electromagnetic exploration |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103105629B (en) | 2015-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zehner et al. | Workflows for generating tetrahedral meshes for finite element simulations on complex geological structures | |
| Yang et al. | Three-dimensional inversion of airborne time-domain electromagnetic data with applications to a porphyry deposit | |
| Oyedele et al. | Integration of 3D seismic and well log data in the optimal reservoir characterisation of emi field, offshore Niger Delta oil province, Nigeria | |
| Xue et al. | Extracting the virtual reflected wavelet from TEM data based on regularizing method | |
| Bujakowski et al. | Integrated seismic and magnetotelluric exploration of the Skierniewice, Poland, geothermal test site | |
| CN108121009B (en) | Development late stage complex fault block velocity-variant depth mapping method | |
| Aweto et al. | Estimating aquifer parameters with geoelectric soundings: case study from the shallow Benin formation at Orerokpe, Western Niger Delta, Nigeria | |
| Lu et al. | 3D electromagnetic modeling of graphitic faults in the Athabasca Basin using a finite-volume time-domain approach with unstructured grids | |
| CN103105629B (en) | A kind of equal electrical resistivity surface chart electromagnetic survey interpretation methods | |
| Ebong et al. | 3D structural modelling and fluid identification in parts of Niger Delta Basin, southern Nigeria | |
| Zhen et al. | Integration of engineering geological investigation data and construction of a 3D geological structure model in the Xiong’an new area | |
| KR101157792B1 (en) | A method for 3-d geological structure analysis by using structure index | |
| Tonini et al. | 3D geo-mapping based on surface data for preliminary study of underground works: A case study in Val Topina (Central Italy) | |
| Özgüler et al. | Geophysical Investigations in Denizli Geothermal Fields, Turkey | |
| Zhang et al. | Method and application of urban 3D rapid modeling of geology based on CAD borehole logs | |
| Zhang et al. | The launch of a large regional gravity information system in China | |
| Pollock et al. | 3D exploratory analysis of descriptive lithology records using regular expressions | |
| Lei et al. | Research on the channel wave field characters of goaf in coal mine and its application | |
| Jordan et al. | Measuring and modeling fault density for CO2 storage plume-fault encounter probability estimation | |
| Tianfu et al. | Borehole databases and 3D geological model of Jurassic-Cretaceous strata in Dongsheng area, North Odors Basin | |
| Shafie et al. | 2D Gravity Inversion Technique in The Study of Cheshire Basin | |
| Rabeh | Using 3-D magnetic modeling to evaluate subsurface structures of the Siwa Oasis, Western Desert, Egypt | |
| Wood et al. | Model creation based on digital borehole records and interpreted geological cross‐sections | |
| Amoka et al. | Application of Geographic Information System for mineral exploration in Nigeria | |
| Soyer et al. | Geophysics over high temperature geothermal fields: Lessons from magnetotelluric and gravity joint inversions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yan Shu Inventor after: Chen Bin Inventor after: Xue Guoqiang Inventor after: Qiu Weizhong Inventor after: Zhu Die Inventor after: Di Qingyun Inventor after: Li Gang Inventor before: Yan Shu Inventor before: Chen Bin Inventor before: Zhu Die |
|
| COR | Change of bibliographic data | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150819 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |