CN103235345B - Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect - Google Patents

Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect Download PDF

Info

Publication number
CN103235345B
CN103235345B CN201310170292.9A CN201310170292A CN103235345B CN 103235345 B CN103235345 B CN 103235345B CN 201310170292 A CN201310170292 A CN 201310170292A CN 103235345 B CN103235345 B CN 103235345B
Authority
CN
China
Prior art keywords
wire frame
slope
time period
transient electromagnetic
attenuation curve
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
CN201310170292.9A
Other languages
Chinese (zh)
Other versions
CN103235345A (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.)
Xian Research Institute Co Ltd of CCTEG
Original Assignee
Xian Research Institute Co Ltd of CCTEG
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 Xian Research Institute Co Ltd of CCTEG filed Critical Xian Research Institute Co Ltd of CCTEG
Priority to CN201310170292.9A priority Critical patent/CN103235345B/en
Publication of CN103235345A publication Critical patent/CN103235345A/en
Application granted granted Critical
Publication of CN103235345B publication Critical patent/CN103235345B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention relates to technical field of geophysical exploration, be specifically related to a kind of method utilizing attenuation curve slope to correct little wire frame transient electromagnetic inductive effect.The inductive effect that its small multi-turn loop correcting the use of little wire frame transient electromagnetic brings image data, eliminates the impact that in existing little wire frame transient electromagnetic processing procedure, inductive effect is brought.The process step that the present invention adopts is: (1) starts; (2) ground large wire frame attenuation curve of secondary field morphological feature is added up; (3) mutational range in little wire frame TEM measured data is selected; (4) slope that in computer capacity, each time period is corresponding; (5) starting point is selected; (6) new slope value is directly composed for the time period being later than starting point; (7) the new slope value of weight coefficient backstepping early than the time period of starting point is introduced; (8) whether, be then terminate, no if being satisfied with according to each measuring point slope of a curve variance yields size (as being less than 0.3) artificial judgment calibration result, then return step (3).

Description

Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect
One, technical field:
The present invention relates to technical field of geophysical exploration, be specifically related to a kind of method utilizing attenuation curve slope to correct little wire frame transient electromagnetic inductive effect.
Two, background technology:
Little wire frame transient electromagnetic method be widely used in that underground coal mine is visited in advance in recent years, the detection of roof and floor water damage, the field such as lateral wall and the inner water damage detection of workplace a kind of electromagnetism class method of exploration.The method is due to underground construction space constraint, the large transmitting wire frame (the general length of side is greater than 100m) that ground routine can not be adopted to use, the length of side is often used to be less than the small multi-turn loop of 2m, this can bring the wire frame self-induction and mutual inductance impact that can not ignore in original sampled data, actual measurement little wire frame transient electromagnetic attenuation curve of secondary field and ground large wire frame transient electromagnetic attenuation curve of secondary field is made to have bigger difference, a lot of more ripe process and inversion method cannot be applied, affect data processing Explanation Accuracy, be difficult to the needs meeting Practical Project.
From data, also do not remove the technology of inductive effect impact at present, common practices is directly deleted at early stage influenced data, do not participate in data processing, make explanation results occur superficial part blind area.
Three, summary of the invention:
The present invention is in order to solve the weak point in above-mentioned background technology, a kind of method utilizing attenuation curve slope to correct little wire frame transient electromagnetic inductive effect is provided, the inductive effect that its small multi-turn loop correcting the use of little wire frame transient electromagnetic brings image data, eliminates the impact that in existing little wire frame transient electromagnetic processing procedure, inductive effect is brought.
For achieving the above object, the technical solution used in the present invention is: a kind of method utilizing attenuation curve slope to correct little wire frame transient electromagnetic inductive effect, is characterized in that: described method flow step is:
(1) start;
(2) ground large wire frame attenuation curve of secondary field morphological feature is added up;
(3) mutational range in little wire frame TEM measured data is selected;
(4) slope that in computer capacity, each time period is corresponding;
(5) starting point is selected;
(6) new slope value is directly composed for the time period being later than starting point;
(7) the new slope value of weight coefficient backstepping early than the time period of starting point is introduced;
(8) whether, be then terminate, no if being satisfied with according to each measuring point slope of a curve variance yields size (as being less than 0.3) artificial judgment calibration result, then return step (3).
Described concrete operation steps is:
Described step (2): ground large wire frame transient electromagnetic attenuation curve of secondary field morphological feature under statistics log-log coordinate system;
Described step (3): under log-log coordinate system, selects little wire frame transient electromagnetic secondary
Field die-away curve does not meet the mutational range of ground large wire frame transient electromagnetic attenuation curve of secondary field form statistical law;
Described step (4): during according to each in the selected scope of following formulae discovery step (2)
Between the slope of section
；
Described step (5): select slope maximum a bit as starting point;
Described step (6): make the calibration slopes of the time period compared with initial point evening directly equal it
Original slope;
Described step (7): from starting point according to the time of following formula backstepping compared with initial point morning
The calibration slopes of section
In formula represent slope after correcting, , represent weight coefficient, weight coefficient can according to such as
Lower formula is determined
Until backstepping is to the time point the earliest of actual measurement, complete whole flow process.
Compared with prior art, the advantage that the present invention has and effect as follows: The present invention gives the rapid and convenient method effectively removing inductive effect impact, add the usable range of raw data, greatly reduce and explain blind area, be conducive to the safety in production ensureing mine development end work.
Four, accompanying drawing illustrates:
Fig. 1 is the general flow chart utilizing attenuation curve slope to correct little wire frame transient electromagnetic inductive effect embodiment of the method;
Fig. 2 is ground large wire frame transient electromagnetic attenuation curve of secondary field aspect graph;
Fig. 3 is mine multiturn little wire frame transient electromagnetic attenuation curve of secondary field aspect graph;
Fig. 4 is primary curve and correction 1 curve comparison figure;
Fig. 5 is correction 1 curve and correction 2 curve comparison figure.
Five, embodiment:
Below in conjunction with embodiment and Figure of description, the present invention is described in detail.
See Fig. 1, a kind of method utilizing attenuation curve slope to correct little wire frame transient electromagnetic inductive effect, process step is:
(1) start;
(2) ground large wire frame attenuation curve of secondary field morphological feature is added up;
(3) mutational range in little wire frame TEM measured data is selected;
(4) slope that in computer capacity, each time period is corresponding;
(5) starting point is selected;
(6) new slope value is directly composed for the time period being later than starting point;
(7) the new slope value of weight coefficient backstepping early than the time period of starting point is introduced;
(8) whether artificial judgment calibration result is satisfied with, and is, then terminate, no, then return step (3).
Realize especially by following steps:
(1), under log-log coordinate system, the mutational range little wire frame transient electromagnetic attenuation curve of secondary field not meeting ground large wire frame transient electromagnetic attenuation curve of secondary field form statistical law is selected;
(2), according to following formula
Calculate each time period slope in (2) selected scope;
(3), select slope maximum a bit as starting point;
(4), the calibration slopes of the time period compared with initial point evening is made directly to equal its original slope;
(5), from starting point according to the calibration slopes of following formula backstepping compared with the initial point time period early
In formula represent slope after correcting, , represent weight coefficient.Weight coefficient can be determined according to following formula
Until backstepping is to the time point the earliest of actual measurement, complete whole flow process.
As the curve after finishing a flow process is still unsatisfied with, thinks and do not eliminate inductive impact completely, above process can be repeated, till satisfied to curve after correction.
For measured data, effect of the present invention is described below:
As can be seen from Figures 2 and 3, under log-log coordinate system, it is close that ground large wire frame TEM die-away curve part shows as global slopes, similar straight line, and the slope of multiturn little wire frame TEM die-away curve early part and late period part slope notable difference larger, similar two straight lines, connecting portion has turning point to exist, and the reason that early part numerical value is larger is exactly because be subject to the impact of inductive effect.
Fig. 4 and Fig. 5 provides a little wire frame TEM measured data, can see that the ground large wire frame TEM die-away curve form of the tracing pattern after correction 1 time and Fig. 2 is similar, eliminate the impact of partial inductance effect, but do not obtain optimal effect yet, after carrying out 1 correction again, achieve more preferably inductive effect eradicating efficacy.
Embodiment:
t U U’
150.917 240955 4384.16 4.32765 1.07239 0.50567 0.33937
160.922 182513 4092.52 5.35115 1.33705 0.50739 0.34464
171.895 128231 3747.03 5.43462 1.63555 0.50792 0.35059
183.902 88845 3355.27 10.2375 2.68895 0.51462 0.35897
195.908 46499.9 2830.57 14.0128 3.30917 0.52182 0.37561
208.883 18931.3 2289.34 10.4566 2.5966 0.51784 0.3979
222.89 9603.6 1934.27 7.15877 1.99382 0.51302 0.41612
237.866 6029.18 1699.06 7.14171 1.96625 0.51373 0.43142
253.874 3786.59 1494.83 5.51399 1.65738 0.51107 0.44676
270.851 2649.95 1342.76 5.45994 1.66574 0.5114 0.46034
288.86 1864.64 1206.23 5.58124 1.59963 0.51207 0.47554
307.838 1307.25 1089.49 3.32158 1.06466 0.50775 0.48989
328.816 1050.17 1015.65 1.50474 0.98893 0.50353 0.49647
350.828 952.607 952.607 2.44434 2.44434
373.807 815.764 815.764 3.10131 3.10131
398.788 667.474 667.474 2.46399 2.46399
425.77 568.036 568.036 2.38937 2.38937
454.753 485.331 485.331 2.07317 2.07317
485.737 423.343 423.343 1.67432 1.67432
518.722 379.244 379.244 1.79753 1.79753
553.708 337.26 337.26 1.88176 1.88176
590.695 298.621 298.621 1.64686 1.64686
629.682 268.786 268.786 1.46141 1.46141
671.64 244.605 244.605 1.43853 1.43853
716.631 222.822 222.822 1.39596 1.39596
764.591 203.556 203.556 1.35389 1.35389
815.586 186.518 186.518 1.34834 1.34834
869.549 171.082 171.082 1.3288 1.3288
927.515 157.023 157.023 1.33224 1.33224
989.483 144.06 144.06 1.29856 1.29856
1055.45 132.478 132.478 1.277 1.277
1125.42 122.052 122.052 1.28784 1.28784
1200.37 112.328 112.328 1.30175 1.30175
1280.34 103.282 103.282 1.25838 1.25838
1365.29 95.2611 95.2611 1.21047 1.21047
1456.24 88.1074 88.1074 1.20275 1.20275
1553.2 81.5346 81.5346 1.17779 1.17779
1657.12 75.5466 75.5466 1.18887 1.18887
1768.08 69.944 69.944 1.21227 1.21227
1886.01 64.6779 64.6779 1.2084 1.2084
2011.95 59.8181 59.8181 1.21951 1.21951
2145.89 55.2966 55.2966 1.17175 1.17175
2288.8 51.273 51.273 1.11535 1.11535
2441.72 47.7047 47.7047 1.09917 1.09917
2604.65 44.435 44.435 1.06664 1.06664
2778.54 41.4751 41.4751 1.08985 1.08985
2964.45 38.6485 38.6485 1.12507 1.12507
3162.36 35.9381 35.9381 1.11969 1.11969
3373.24 33.4321 33.4321 1.12895 1.12895
3598.13 31.0828 31.0828 1.13059 1.13059
3838 28.8956 28.8956 1.12953 1.12953
4093.87 26.8641 26.8641 1.14591 1.14591
4366.72 24.9495 24.9495 1.15781 1.15781
4658.55 23.149 23.149 1.06102 1.06102
4970.39 21.6111 21.6111 0.96118 0.96118
5302.24 20.3094 20.3094 0.91134 0.91134
5656.04 19.1483 19.1483 0.8303 0.8303
6033.85 18.1474 18.1474 0.85748 0.85748
6436.64 17.1692 17.1692 0.89228 0.89228
6866.41 16.207 16.207 0.84884 0.84884
7325.17 15.3412 15.3412 0.83072 0.83072
7814.91 14.5382 14.5382 0.71446 0.71446
8337.63 13.881 13.881 0.5814 0.5814
8895.34 13.3682 13.3682 0.51406 0.51406
9490.04 12.9308 12.9308 0.41743 0.41743
10123.7 12.5866 12.5866 0.43028 0.43028
10799.4 12.2415 12.2415 0.45518 0.45518
11520 11.8868 11.8868 0.40088 0.40088
12288.6 11.583 11.583 0.37454 0.37454
13109.1 11.3059 11.3059 0.34916 0.34916
13984.7 11.0536 11.0536 0.31866 0.31866
14918.2 10.8283 10.8283 0.31105 0.31105
15913.7 10.6129 10.6129 0.304 0.304
16976.1 10.4064 10.4064 0.29177 0.29177
18109.6 10.212 10.212 0.28568 0.28568
19318.9 10.0251 10.0251 0.28513 0.28513
20609.2 9.84201 9.84201 0.28729 0.28729
21985.5 9.66092 9.66092 0.30726 0.30726
23453.8 9.4709 9.4709 0.33248 0.33248
25019.9 9.26953 9.26953 0.3541 0.3541
26691.1 9.05972 9.05972 0.38351 0.38351
28474.1 8.8378 8.8378 0.32993 0.32993
30376.1 8.65125 8.65125 0.27148 0.27148
32405.1 8.50072 8.50072 0.2611 0.2611
34569.9 8.35839 8.35839 0.23003 0.23003
36878.7 8.23501 8.23501 0.28276 0.28276
39341.4 8.08585 8.08585 0.34545 0.34545
41969.1 7.90725 7.90725 0.35126 0.35126
44771.6 7.72974 7.72974 0.37946 0.37946
47761 7.54247 7.54247 0.23755 0.23755
50950.4 7.42753 7.42753 0.07545 0.07545
54352.6 7.39139 7.39139 0.01227 0.01227
57982.7 7.38553 7.38553 -0.09894 -0.09894
61855.7 7.43293 7.43293 0.14912 0.14912
65986.5 7.36162 7.36162 0.4336 0.4336
70393.2 7.15814 7.15814 0.5054 0.5054
75094.8 6.92801 6.92801 0.67481 0.67481
80110.1 6.63226 6.63226 0.42399 0.42399
85460.3 6.45293 6.45293 0.12752 0.12752
91167.4 6.39996 6.39996 0.12752 0.12752

Claims (1)

1. utilize attenuation curve slope to correct a method for little wire frame transient electromagnetic inductive effect, it is characterized in that: described method flow step is:
(1) start;
(2) ground large wire frame attenuation curve of secondary field morphological feature is added up;
(3) mutational range in little wire frame TEM measured data is selected;
(4) slope that in computer capacity, each time period is corresponding;
(5) starting point is selected;
(6) new slope value is directly composed for the time period being later than starting point;
(7) the new slope value of weight coefficient backstepping early than the time period of starting point is introduced;
(8) whether, be then terminate, no if being satisfied with according to each measuring point slope of a curve variance yields size artificial judgment calibration result, then return step (3);
Described step (2): ground large wire frame transient electromagnetic attenuation curve of secondary field morphological feature under statistics log-log coordinate system;
Described step (3): under log-log coordinate system, selects the mutational range little wire frame transient electromagnetic attenuation curve of secondary field not meeting ground large wire frame transient electromagnetic attenuation curve of secondary field form statistical law;
Described step (4): the slope selecting each time period in scope according to following formulae discovery ;
Described step (5): select slope maximum a bit as starting point;
Described step (6): make the calibration slopes of the time period compared with initial point evening directly equal its original slope;
Described step (7): from starting point according to the calibration slopes of following formula backstepping compared with the initial point time period early ; In formula represent slope after correcting, , represent weight coefficient, weight coefficient can be determined according to following formula
Until backstepping is to the time point the earliest of actual measurement, complete whole flow process.
CN201310170292.9A 2013-05-10 2013-05-10 Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect Active CN103235345B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310170292.9A CN103235345B (en) 2013-05-10 2013-05-10 Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310170292.9A CN103235345B (en) 2013-05-10 2013-05-10 Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect

Publications (2)

Publication Number Publication Date
CN103235345A CN103235345A (en) 2013-08-07
CN103235345B true CN103235345B (en) 2015-11-18

Family

ID=48883396

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310170292.9A Active CN103235345B (en) 2013-05-10 2013-05-10 Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect

Country Status (1)

Country Link
CN (1) CN103235345B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105068129A (en) * 2015-07-06 2015-11-18 中煤科工集团西安研究院有限公司 Smooth processing method for wavelet data of transient electromagnetic data
CN107765319B (en) * 2017-10-09 2019-04-30 中国矿业大学(北京) A kind of Mine Transient Electromagnetic Method resistivity correction method
CN108776357B (en) * 2018-04-16 2019-10-01 中煤科工集团西安研究院有限公司 The bearing calibration of sedimentary formation transient electromagnetic method electromagnetic interference and device

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202075419U (en) * 2011-03-03 2011-12-14 西安理工大学 Transient electromagnetic method shallow-layer coincident loop multi-coil difference receiving device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2009219487B2 (en) * 2008-02-25 2011-11-17 Shell Internationale Research Maatschappij B.V. Method of determining a transient electromagnetic response of a formation

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202075419U (en) * 2011-03-03 2011-12-14 西安理工大学 Transient electromagnetic method shallow-layer coincident loop multi-coil difference receiving device

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
卡尔曼滤波压制瞬变电磁干扰技术及应用;范涛 等;《煤田地质与勘探》;20130430;第41卷(第2期);全文 *
洞内瞬变电磁（TEM）探测技术研究;杜庆丰 等;《物探化探计算技术》;20071031;第29卷;摘要，第166页第2栏倒数第1-3段，第167页第7段 *
瞬变电磁数据与处理方法探讨;杜庆丰 等;《物探与化探》;20060228;第30卷(第1期);全文 *
瞬变电磁法大电流、小回线装置在矿井巷道内探测的应用及数据精细处理;崔健;《中国优秀硕士学位论文全文库 基础科学辑》;20110815;全文 *

Also Published As

Publication number Publication date
CN103235345A (en) 2013-08-07

Similar Documents

Publication Publication Date Title
CN103235345B (en) Attenuation curve slope is utilized to correct the method for little wire frame transient electromagnetic inductive effect
CN102609687B (en) Subway construction drawing and engineering parameter automatic identification method
CN105719259A (en) Pavement crack image detection method
CN102779231B (en) Based on contiguous coordinate transformation parameter computing method
CN103761405A (en) Three-dimensional underground integrated pipelines auto-generating system with error self-checking function
CN103745487B (en) The mixed compression sensing method of Bayes's EO-1 hyperion solution based on structural sparse priori
CN109934915B (en) Three-dimensional modeling method for deep alteration minerals
WO2021243967A1 (en) Three-dimensional resistivity tomography method and system
CN105095540B (en) The normal plane scanning analytic method of relative position relation between a kind of well
CN107588878A (en) A kind of easy accurate interpretation method of the geostress survey based on rock kasier point features
CN105547166A (en) 2D laser scanning based rapid detection method for heading window of heading machine
CN109188481A (en) A kind of new method being fitted GPS height anomaly
CN104778678A (en) Pore throat recognition method with consideration of pore throat tail end
CN103874198A (en) Method for precisely positioning wireless mobile terminal in mine environment
CN204881535U (en) Pipe diameter measuring device in pipeline of harbour
CN103267985A (en) Method for rectifying mine transient electromagnetic inductive effect by utilizing emission current
WO2024078134A1 (en) Excavation tunnel full-waveform inversion method based on multi-parameter constraint and structure correction
CN103489197B (en) A kind of urban aerial image corner feature matching process
CN104715160A (en) Soft measurement modeling data outlier detecting method based on KMDB
CN104504697A (en) Automatic extracting method for neck-edge line of denture repair preparation mesh model
CN105572733B (en) A kind of seismic velocity spectrum automatic pick method
CN102903084B (en) Wavelet field image noise variance method of estimation under a kind of α stable model
CN114236624B (en) Method and system for estimating fracturing modification space volume based on electromagnetic method
CN203949607U (en) A kind of reinforcing bar is measured slide calliper rule
CN105824042A (en) Maximum longitudinal distance design method based on optimal lighting energy

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB02 Change of applicant information

Address after: 710075 Xi'an high tech Industrial Development Zone, Shaanxi, Jin Jin Road, No. 82

Applicant after: Xi'an Research Institute, CCTEG

Address before: 710075 Xi'an high tech Industrial Development Zone, Shaanxi, Jin Jin Road, No. 82

Applicant before: China Coal Technology & Engineering Group Xi'an Research Institute

COR Change of bibliographic data

Free format text: CORRECT: APPLICANT; FROM: CHINA COAL TECHNOLOGY + ENGINEERING GROUP XI AN RESEARCH INSTITUTE TO: XI AN RESEARCH INSTITUTE OF CHINA COAL TECHNOLOGY + ENGINEERING GROUP CO., LTD.

C14 Grant of patent or utility model
GR01 Patent grant