CN108227006A - Applied to the time shift controlled-source audiomagnetotelluric detection method during cbm development - Google Patents

Applied to the time shift controlled-source audiomagnetotelluric detection method during cbm development Download PDF

Info

Publication number
CN108227006A
CN108227006A CN201710055199.1A CN201710055199A CN108227006A CN 108227006 A CN108227006 A CN 108227006A CN 201710055199 A CN201710055199 A CN 201710055199A CN 108227006 A CN108227006 A CN 108227006A
Authority
CN
China
Prior art keywords
apparent resistivity
eyeball
time shift
data
cbm development
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.)
Pending
Application number
CN201710055199.1A
Other languages
Chinese (zh)
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.)
China Petroleum and Chemical Corp
Original Assignee
China Petroleum and Chemical Corp
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 China Petroleum and Chemical Corp filed Critical China Petroleum and Chemical Corp
Priority to CN201710055199.1A priority Critical patent/CN108227006A/en
Publication of CN108227006A publication Critical patent/CN108227006A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/083Controlled source electromagnetic [CSEM] surveying

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

A kind of time shift controlled-source audiomagnetotelluric detection method applied to during cbm development, it includes the following steps:Before the cbm development of target area and in development process, according to same observation procedure, same target area or trial zone are repeated several times with requirement according to intervals and six and obtain several eyeball apparent resistivity curves;It to aforementioned data processing, obtains along layer difference apparent resistivity data body, the range that coal seam reservoirs water system dynamic change and the cbm development process during coal bed gas well mining of obtaining involve achievees the purpose that Quantitative Monitoring cbm development.The method of the present invention has been knitted in Guizhou to be applied, and achieve successful experience in golden block, Yanchuan South block cbm development.

Description

Applied to the time shift controlled-source audiomagnetotelluric detection during cbm development Method
Technical field
The present invention relates to the applications of time shift controlled-source audiomagnetotelluric (CSAMT) technology in cbm development field. In cbm development target area, controlled-source audiomagnetotelluric (CSAMT) technology is applied under being constrained based on Seismic-geologic Models, Change in apparent resistivity is observed to obtain coal seam reservoirs and upper and lower stratum hydrodynamic force change information by time shift, reaches monitoring coal bed gas The purpose of development effectiveness.
Background technology
The theoretical research of dynamic detection is carried out during cbm development and discussing at home and abroad has article to deliver.In State's mining industry university coal bed gas resource and key lab of the Filling process Ministry of Education and Xinjiang University's geology and mineral engineering institute (Ge Yanyan, Fu Xuehai etc., 2013 years) pass through Transient Electromagnetic Technique (Transient Electromagnetic in somewhere Method, TEM), utilize coal measures change in resistance during transient electromagnetic method detection achievement inverting coal bed gas well mining, monitoring Coal seam reservoirs water system dynamic change characterization during coal bed gas well mining.
Fig. 1 can be seen that shallow earth's surface formation variation on electric section is little, and what changing occurred in apparent resistivity is to enter After coal measures, on electric section, there is variation in resistivity, after mining half a year detection apparent resistivity numerical value with respect to before mining Increase, declaratives coal seam reservoirs water is discharged, apparent resistivity increase.
Detection of coal seam bedding cutting apparent resistivity is substantially uniform before Fig. 2 can be seen that mining, without apparent strong and weak exceptions area Domain, reflection underground medium maintain the electrical uniform state of relative equilibrium;Detection of coal seam bedding cutting apparent resistivity after mining half a year There is strong and weak anomalous variation in numerical value.In figure red line area marking be apparent resistivity riser region, blue line tab area be regarding Resistivity reduces region.
The dynamic change of apparent resistivity shows the otherness before and after the anisotropism on stratum, reservoir reconstruction, coal seam reservoirs and top The hydraulic connection of bottom plate reflects the direction of propagation of underground water during coal bed gas well mining.
The dynamic change of apparent resistivity not only reflects the otherness before and after reservoir reconstruction, but also indicates that coal bed gas well is arranged The migration direction of underground water during adopting, result of the test provide new approaches for the Groundwater movement monitoring during cbm development And new technology.
Transient Electromagnetic Technique (TEM) provides new approaches although being monitored for the Groundwater movement during cbm development, There are following deficiencies:
1) constraint of Seismic-geologic Models is not based on, the change of divergence situation of the apparent resistivity on section is a kind of comprehensive Close effect, it is impossible to single reflection coal seam reservoirs change in apparent resistivity situation;
2) cbm development area more mature at present is generally possible to carry out cbm development mainly in North China and south Single coal seam thickness there was only 2~10 meters.The variation zone reflected on TEM electric sections is in 200~550 meters of ranges of absolute altitude, thus The sea level changes along coal seam bedding cutting arrived cannot be simply considered that it is range that coal bed gas extraction is involved;
3) above method application does not obviously account for seasonality to electrically observing influence, and six months difference is exactly season On the contrary.Because rainy season and dry season are observed same target area, electrical value difference is different to be increased;
4) Fig. 3 is carefully analyzed, the section left side shallow-layer variation zone observed after half a year is also just deep layer variation zone, corresponding Absolute altitude it is relatively low, surface layer, which electrically changes, can directly affect deeper formation change in apparent resistivity.
Invention content
During being directed to cbm development, how economic, applicable, effective ways are carried out to monitor Coal bed gas extraction effect proposes a kind of time shift controlled-source audiomagnetotelluric (CSAMT) applied to during cbm development Detection method.
The technical scheme is that:
A kind of time shift controlled-source audiomagnetotelluric detection method applied to during cbm development, it include with Lower step:
(1) several eyeballs are chosen according to cbm development target area to be detected, before cbm development, according to Certain observation procedure obtains several eyeball apparent resistivity curves in target area for the first time;In target area CBM Drainage After adopting, when drilling well or the experiment hydraulic fracturing of small well group carry out, with corresponding in second of acquisition target area of identical observation procedure Eyeball apparent resistivity curve;It is more later at regular intervals with identical observation procedure with cbm development process It is secondary to obtain corresponding eyeball apparent resistivity curve and eyeball distribution characteristics data are formed with the apparent resistivity curve of each eyeball Body;
(2) quantitative inversion calculating is carried out to each eyeball apparent resistivity curve obtained each time, obtains achievement apparent resistance Rate data volume, achievement apparent resistivity data body is arranged in chronological order, and achievement apparent resistivity is obtained by mathematic interpolation The edge layer difference apparent resistivity data body of data volume;
(3) the difference boundary of coal seam reservoirs apparent resistivity data body is drawn a circle to approve in the plane, temporally observes the variation of these boundaries, Obtain the range involved during coal seam reservoirs water system dynamic change and cbm development during coal bed gas well mining.
Observation eyeball sets resistivity curve instrument as controlled-source audiomagnetotelluric in the step (1) of the present invention Instrument further, data is carried out using the V8CSAMT instruments of Canadian phoenix company production or the instrument of model suitable with its Acquisition.
The regional extent of eyeball should cover and comprehensively more than cbm development target area in the step (1) of the present invention;It is real The measurement and positioning of measuring point uses high-precision GPS positioning measurement instrument, performs GB/T18314-2009《Global positioning system (GPS) Specifications of surveys》;Observation system, receiving point away from, receive the indexs of correlation such as line-spacing, the distance of controllable source and reception area, perform SY/T 5772-2012《Controlled source audio-frequency magnetotelluric method (CSAMT) exploration engineering regulation》.
The time interval for obtaining eyeball apparent resistivity curve in the step (1) of the present invention every time is 2 years, first number Do not carry out in the area of observation coverage there are no drilling well or when small well group is tested according to acquisition and carry out.
The step (1) of the present invention is since being observed second, it is desirable that accomplishes the observation system with observation for the first time every time System, excitaton source, receiving point, construction parameter, acquisition instrument, construction seasons are identical.
The step (2) of the present invention is specially:
(a) Geologic modeling is carried out using the drilling well in target area and seismic data, the geological layering for obtaining each eyeball is deep Degrees of data and formation lithology demarcate each formation lithology change in apparent resistivity range;
(b) under the constraint of above-mentioned geological model, each eyeball apparent resistivity curve that different time sections obtain is carried out Quantitative inversion calculates;The interference of variable, systematic error, geography information in gatherer process is removed, is obtained corresponding to each period Achievement apparent resistivity data body, the difference of each achievement apparent resistivity data body be only during coal-bed gas exploitation different times draw Caused by the coal seam reservoirs SEA LEVEL VARIATION single factors risen.
Dependent variable causes the method for interference in removal gatherer process:The electrical method measuring point acquired twice compares and analyzes, right In there are the point of notable difference, if due to changing caused by acquisition environmental change or humanity activities, then removal or again Secondary acquisition.
Because systematic error causes the method for interference in removal gatherer process:Unminding area's progress, electrical method data are adopted twice Collection, makes the difference quantitative inversion resistivity twice, can be quantified in target series of strata and obtain the system difference of gathered data twice.
Because geography information variation causes the method for interference in removal gatherer process:The geographic basis of each collection point twice It compares and analyzes, noise spot is generated to electrical method acquisition condition, and (such as nearby there is newly-built wire facility interference in second of collection point Source), remove interference source according to actual conditions.
(c) each achievement apparent resistivity data body is subjected to spacial migration according to observation time, obtains time shift apparent resistivity sky Between data volume, time shift apparent resistivity spatial data body is recorded as:Time shift controlled-source audiomagnetotelluric apparent resistivity data Body;
(d) in step (a) under the constraint of geological model, each time shift apparent resistivity spatial data that step (c) is obtained Body is divided into according to geological layering depth data and formation lithology along layer time shift apparent resistivity spatial data body, is recorded as:Along layer Time shift controlled-source audiomagnetotelluric apparent resistivity data body;
(e) difference is asked for by adjacent time to the edge layer time shift apparent resistivity data body that step (d) obtains, obtains along layer Difference apparent resistivity data body, is recorded as:Along layer time shift controlled-source audiomagnetotelluric apparent resistivity variance data body.
The present invention step (e) be additionally included in along layer time shift apparent resistivity data body by adjacent time ask for difference it Before, different depth correction and flattening technical finesse are carried out to the edge layer apparent resistivity spatial data body that step (d) obtains.
The step (3) of the present invention is in data exception point, data handling procedure caused by possibility in data acquisition It is studied exceptions area caused by possible;The reason of using difference is caused along layer difference apparent resistivity data body profile analysis, knows The range that other coal seam reservoirs water system dynamic change situation and cbm development process pressure drop involve.
The identification of data exception point:It is determined according to industry standard, is normal in the error range for filling perhaps in specification Error fills error perhaps more than specification, is exactly data exception point.The value that the size of data exception point should be abnormal point subtracts just Constant error is poor.Multiple data exception points are linked to be piece, are exactly exceptions area.
There are many reason of causing data exception point.In coal-bed gas exploitation area, only exclude what claim 6 (b) described Dependent variable in gatherer process, systematic error, geography information interference exceptions area, it is because of coal-bed gas exploitation area that could react Water and the coal bed gas content changed region during cbm development.
Using along layer difference apparent resistivity data body section, above-mentioned exceptions area boundary point more accurately can be analyzed and drawn a circle to approve With abnormal point value size, corresponding isogram and delineation variation boundary model on these spot projections to plan view, will be sketched out It encloses, so as to the model for obtaining coal-bed gas exploitation area coal seam reservoirs water system dynamic change situation and cbm development process pressure drop involves It encloses.
Beneficial effects of the present invention:
The present invention has carried out exploratory application in Rocks in Zhijin, Guizhou and Yanchuan South block, in practice it has proved that:Cbm development mistake Detect coal bed gas extraction dynamic effect using time shift controlled-source audiomagnetotelluric method in journey, Method And Principle be it is feasible, It is and economic and practical.By first Application, preferable geological effect is achieved knitting the golden block coal bed gas extraction area of observation coverage.Thus It is as follows to generate social benefit:
1) it is for the first time cbm development dynamic detection and decision, searches out the method and technology of " economy is applicable in, is effective ";
2) controlled-source audiomagnetotelluric (CSAMT) method application field has been expanded, first Application is led in cbm development Domain belongs to innovation;
3) application conditions and timing are that time shift controlled-source audiomagnetotelluric (CSAMT) method is applied to coal bed gas Development behavior detection is successfully crucial;
4) popularization and application of the method will be that coal gas industry (clean energy resource) plays a role in promoting, and generate huge society Benefit.
Description of the drawings
Fig. 1:Apparent resistivity intends the comparison of section before the A line minings of somewhere, after mining half a year.
Fig. 2:Apparent resistivity bedding cutting compares before the coal seam mining of somewhere 3, after mining half a year.
Fig. 3:Three dimensional seismic data well, shake joint interpretation section.
Fig. 4:Guizhou Rocks in Zhijin, Guizhou CSAMT method NW106 line the old and new Data Inversion comparison section.
Fig. 5:Guizhou Rocks in Zhijin, Guizhou three-dimensional apparent resistivity data body Lungtan Formation bottom bedding cutting comparison plan view.
Specific embodiment
During cbm development, by certain time interval (generally 2 years), target observation area is used can Control source audio magneto-electrotelluric exploration (CSAMT) method carries out " six is same " data acquisition, is carried out under Seismic-geologic Models constraint " five unifications " data processing, by analysis and research, it is practical with reference to coal bed gas extraction scene, by coal seam regard electric rate change association as Coal seam reservoirs hydrological variation situation, so as to reach coal bed gas extraction detection result.
Time shift controlled-source audiomagnetotelluric exploration (CSAMT) method is first Application to coal bed gas extraction detection field, And preferable geological effect is obtained, is the good method of one " economy is applicable in, is effective ", is worthy of popularization.
Guizhou knits golden coal bed gas and surveys area, has carried out 3-d seismic exploration, small well group drilling well, real to underground coal bed gas stratum Show well, shake data combine explanation (see Fig. 3), establish Seismic-geologic Models.Under the constraint of above-mentioned Seismic-geologic Models, this Controlled-source audiomagnetotelluric observation has been carried out before surveying the small well group drilling well in area;Small well group mining has carried out controllable again after 2 years Source audio magnetotelluric (CSAMT) is observed, and realizes time shift controlled-source audiomagnetotelluric observation technology.
According to physical difference, the stratum that Guizhou can be knitted to golden coal bed gas survey area is divided into 4 physical property floor from top to bottom:Three is folded System flies celestial pass group (T1f);Permian System Changxing Group (P2c);Permian Longtan group (P2l) and Permian System EMEISHAN BASALTS (P2 β).The compositions such as outwash, lam, sandy loam are predominantly remained in 4th system (Q), and resistivity is extremely low resistance between 10-50, By the 4th it is relatively thin (brill encounter the 4th be that thickness is both less than 5m) according to drilling well thickness statistical result, in conjunction with electrical method inverting, Surface layer superficial part does not occur the low-resistance stratum of certain scale, therefore the layer is not divided in practice interpretation, and merging is included into three In folded system;Triassic System Feixianguan group (T1f), lithology is based on limestone, siliceous limestone, carbonate rock, resistivity 1500- 3000 Ω .m are high resistant;Permian System Changxing Group (P2c) principal lithologic presss from both sides siltstone, based on argillaceous limestone for limestone, resistivity It is middle high resistant between 600-1250 Ω .m;Permian Longtan group (P2l) coal measure strata resistivity 100-750 Ω .m are low-resistance; Permian System EMEISHAN BASALTS (P2 β) resistivity is more than 3500 Ω .m, is high resistant.
Guizhou knit golden coal bed gas survey the acquisition of area's controlled-source audiomagnetotelluric (CSAMT) data respectively at 2013 years and What 2015 annual winters carried out, accomplish that (excitation, reception parameter want one for homologous (excitation), same to point (reception), same parameter substantially Cause), same to instrument (or observation instrument model consistent or same class), same to mode (observation system is consistent), will with season etc. " six is same " It asks;Data processing is completed in March, 2016, is under based on above-mentioned Seismic-geologic Models constraints, with same treatment system " five the is identical " requirement pair such as method, same model constraint inverting is corrected in system, same treatment flow, same treatment parameter, identical near field Data is handled twice;Analysis and research are fully to combine coal bed gas extraction well site actual conditions and mining process, right Exceptions area caused by possibility is studied in data exception point, data handling procedure caused by possibility in data acquisition;Profit With differential section or data volume, the reason of causing difference is analyzed;In the case where seismic data establishes geological model constraint, time shift is analyzed Changes in resistance rate section, so as to identify, coal stores up during Permian Longtan group (P2l) coal measure strata resistivity coal bed gas well mining Layer water system dynamic change situation, the range that monitoring cbm development process pressure drop involves, reaches dynamic detection cbm development Effect (see Fig. 4, Fig. 5).
Fig. 4 is Guizhou Rocks in Zhijin, Guizhou CSAMT method NW106 line the old and new Data Inversion comparison section.In Seismic-geologic Models about Inversion result shows on 2015 annual inverting apparent resistance inverting sections between 4480-4840 points under beam, height above sea level 1200- Region between 1050m, compared with 2013 annual data inverting resistivity profiles, in the low resistance region of " down warping region " form, apparent resistance Rate is significantly lower, and the resistivity of the corresponding low resistance region both sides, which has, slightly to be increased.With reference to Zhi Jince areas coal bed gas extraction scene Situation, which, which is in, knits 2U1P wells and knits near 1 well of examination, thus it is speculated that should be this two mouthfuls of well hydraulic fracturings and mining pressure drop The range involved, so as to form low resistance region.
Fig. 5 is Guizhou Rocks in Zhijin, Guizhou three-dimensional apparent resistivity data body Lungtan Formation bottom bedding cutting comparison plan view.It will be each Two dimensional inversion section is combined into three-dimensional data, carries out 3D grid processing, obtains intending three-dimensional apparent resistivity data body.Along earthquake Lungtan Formation bottom in Geologic modeling along coal seam attributes extraction, obtain new, old two groups of coal seam reservoirs bedding cuttings and difference is cut Piece.By analysis and research, the larger region of difference value in figure can be considered and knit 2U1P wells and knit 1 well hydraulic fracturing of examination and mining The range involved.
It promotes and applies requirement and is applicable in following 6 conditions:
1) coal seam buried depth should be less than controlled-source audiomagnetotelluric exploration (CSAMT) method in itself best depth of exploration it is (small In 1000 meters), preferably between 300~800 meters;
2) lithology of area of observation coverage earth's surface exposure is preferably single, to reduce when data acquire because excitation, condition of acceptance variation are drawn Play data exception;
3) data collection of different time requires to accomplish " six is same ":I.e. same observation system, same source excitation receive with point, are same Construction parameter is constructed with instrument acquisition, with season, utmostly reduces influence of the non-geologic(al) factor to time shift gathered data;
4) with reference to SY/T 5772-2012《Controlled source audio-frequency magnetotelluric method (CSAMT) exploration engineering regulation》With new treatment technology With algorithm, under the geological model constraint established in drilling well and seismic data, each eyeball apparent resistance for being obtained to different time sections Rate curve is uniformly processed, and is accomplished " five unifications ".I.e.:Unified processing system is uniformly processed flow, parameter is uniformly processed, system Method is corrected in one near field, unified model constraint inverting, utmostly reduces in data processing procedure human factor or during System level gray correlation pair Move the influence of controlled source audio-frequency magnetotelluric data processing achievement;
5) in analysis and research, it is desirable that with reference to coal bed gas extraction well site actual conditions and mining process, data are acquired Exceptions area caused by possibility is studied in data exception point, data handling procedure caused by possibility in the process;It is cutd open using difference The reason of face or data volume, analysis causes difference;Under the geological model constraint established in drilling well and seismic data, time shift is analyzed Changes in resistance rate section so as to identify coal seam reservoirs water system dynamic change situation during coal bed gas well mining, monitors coal bed gas The range that development process involves reaches dynamic detection cbm development effect;
6) opportunity of data acquisition is accurately held.There are no drilling well or not in the first time data acquisition preferably area of observation coverage It is carried out, and using this apparent resistivity as the background of the area of observation coverage when carrying out small well group experiment;Second of data acquisition is preferably being seen It surveys in area during drilling well or small well group test well carry out and is carried out after hydraulic fracturing, coal bed gas extraction, be conducive in this way by with first The acquisition of secondary data is compared, detect well water force split, the range that coal bed gas extraction is involved, so as to more accurate Judgement well water force split, coal bed gas extraction geological effect;Later by primary progress repeated data acquisition in 2 years, so as to reach Coal bed gas extraction dynamic detection effect.
Here it makes additional remarks to time shift controlled-source audiomagnetotelluric (CSAMT) data acquisition " six is same " as follows:
" same to observation system ":Cbm development target area or trial zone are faced, first has to clear and definite observation scope.In face of seeing Survey range will on the whole, long-range ground overall situation set out, the observation system of design " economic, be applicable in, be effective ", and each data are adopted Same observation system will be used during collection;
" same source excitation is received with point ":Refer to from first time data acquire when identified controllable source excitation position and every A receiving point position requires during each data acquisition repeated with it afterwards.Using GPS positioning and repeated fill level measurement, Precision reaches meter level.As the cbm development frequent activity of target area, topographical surface feature increase, small part observation point can be caused It cannot accomplish that " same source excitation is received with point " happens.Faced with this situation, it should ensure controlled-source audiomagnetotelluric (CSAMT) on the basis of accuracy of data acquisition and confidence level, as possible using GPS measurements carry out a position translation and close to original position into Row data acquire, and are corrected when follow-up data is handled;
" same to construction parameter ":When referring to data collection each time, controllable source shooting parameter (source electrode away from, direction, power, Frequency, duration etc.), receiving point parameter (pole span, buried depth, direction, duration etc.) requirement it is consistent;
" being acquired with instrument ":Refer to that data acquisition request is with same (or type) instrument each time.Due to instrument Upgrading or update lead to the promotions such as precision, the sensitivity of data acquisition so that primary before the mass ratio of latter secondary data acquisition It is good, equivalent to increase a variable.It is therefore desirable to when first time data acquire, as far as possible using state-of-the-art instrument at that time Device is acquired, and can ensure to keep the relative stability of acquisition instrument during data redundancy acquisition afterwards in this way.
" constructing with season ":Refer to that data acquisition request was carried out in the same season of different year each time.This is Because controlled-source audiomagnetotelluric (CSAMT) data are affected when acquiring by earth's surface water content.
It makes additional remarks to time shift controlled-source audiomagnetotelluric (CSAMT) data processing " five unifications " as follows:
" unified processing system is uniformly processed flow, parameter is uniformly processed, unified near field correction method, unified model constraint After inverting " refers to data acquisition new each time, with new theoretical, advanced processing system, it is desirable that previous all data The data of acquisition carries out " five unifications " data processing, to eliminate the influence brought by data processing.
2013 year and 2015 year, knitted respectively in Guizhou gold, Shanxi Yanchuan Nan Meicengqikaicaiqu has carried out 2 times can Control source audio magneto-electrotelluric exploration (CSAMT) " six is same " data acquisition.By being based under Seismic-geologic Models constraint " five is identical " Data processing, the result shows that, it knits golden block and obtains preferable geological effect, and Yanchuan south is poor.Analysis main cause be There are following 4 points differences in data acquisition link:
1) golden block area of observation coverage coal seam buried depth is knitted less than 300~800 meters, is just in controlled-source audiomagnetotelluric It explores in (CSAMT) optimum range;And Yanchuan South block area of observation coverage coal seam buried depth is in controllable source sound at 1000~1400 meters Outside frequency magneto-electrotelluric exploration (CSAMT) optimum range;
2) it is Triassic system limestone to knit golden block area of observation coverage exposure lithology, and lithology is single, and earth's surface height difference is smaller;And Yanchuan is southern It is that huge thick loess is covered on the basement rock of exposure that block area of observation coverage earth's surface, which is the 4th, and lithology cross directional variations are larger, and earth's surface is cut Seriously, height difference is larger;
3) knit the golden block area of observation coverage be located at " day without three days it is fine, put down without three rulers " Guizhou Province is domestic, earth's surface water-bearing layer is by season Section influences smaller;And the Yanchuan South block area of observation coverage is in that Shanxi Province is domestic, and earth's surface water-bearing layer is larger by seasonal effect, 2013 years Data acquisition is that in winter (dry season), the acquisition of 2015 annual datas is in summer (rainy season);
4) when knitting golden block area of observation coverage first time data acquisition, only 1 mouthful of drilling well in the area of observation coverage, small well group experiment does not have also There is development, second of data acquisition forces in small well group well water and split, carried out after 1 year and a half after coal bed gas extraction, so inspection It surveys with obvious effects;And during the acquisition of Yanchuan south area of observation coverage first time data, carry out small well group experiment, second number in the area of observation coverage Production is built in business development, scale according to acquisition to carry out in the process, drilling well number be continuously increased, surface pipeline network and station construction It is more etc. artificial disturbing factor, it handles and explains and is relative complex.

Claims (10)

1. a kind of time shift controlled-source audiomagnetotelluric detection method applied to during cbm development, it is characterized in that it is wrapped Include following steps:
(1) several eyeballs are chosen according to cbm development target area to be detected, before cbm development, according to certain Observation procedure obtains several eyeball apparent resistivity curves in target area for the first time;After the coal bed gas extraction of target area, bore When well or the experiment hydraulic fracturing of small well group carry out, regarded with corresponding to eyeball in second of acquisition target area of identical observation procedure Resistivity curve;With cbm development process, repeatedly obtain and correspond to identical observation procedure at regular intervals later Eyeball apparent resistivity curve is simultaneously formed eyeball distribution characteristics data volume with the apparent resistivity curve of each eyeball;
(2) quantitative inversion calculating is carried out to each eyeball apparent resistivity curve obtained each time, obtains achievement apparent resistivity number According to body, achievement apparent resistivity data body is arranged in chronological order, achievement apparent resistivity data is obtained by mathematic interpolation The edge layer difference apparent resistivity data body of body;
(3) the difference boundary of coal seam reservoirs apparent resistivity data body is drawn a circle to approve in the plane, is temporally observed the variation of these boundaries, is obtained The range that coal seam reservoirs water system dynamic change and cbm development involve in the process during coal bed gas well mining.
2. the time shift controlled-source audiomagnetotelluric detection side according to claim 1 applied to during cbm development Method, it is characterized in that observation eyeball sets resistivity curve instrument as controlled-source audiomagnetotelluric instrument in step (1), into One step carries out data collection using the V8CSAMT instruments of Canadian phoenix company production or the instrument of model suitable with its.
3. the time shift controlled-source audiomagnetotelluric detection side according to claim 1 applied to during cbm development Method, it is characterized in that the regional extent of eyeball should cover and more than cbm development target area in step (1) comprehensively;Eyeball Measurement and positioning uses high-precision GPS positioning measurement instrument, performs GB/T18314-2009《Global positioning system (GPS) measuring gage Model》;Observation system, receiving point away from, receive the indexs of correlation such as line-spacing, the distance of controllable source and reception area, perform SY/T 5772- 2012《Controlled source audio-frequency magnetotelluric method (CSAMT) exploration engineering regulation》.
4. the time shift controlled-source audiomagnetotelluric detection side according to claim 1 applied to during cbm development Method, it is characterized in that the time interval for obtaining eyeball apparent resistivity curve in step (1) every time is 2 years, the acquisition of first time data Do not carry out when small well group is tested there are no drilling well or in the area of observation coverage and carry out.
5. the time shift controlled-source audiomagnetotelluric detection side according to claim 1 applied to during cbm development Method, it is characterized in that step (1) is since being observed second, it is desirable that accomplish every time and the observation system of observation for the first time, excitation Source, receiving point, construction parameter, acquisition instrument, construction seasons are identical.
6. the time shift controlled-source audiomagnetotelluric detection side according to claim 1 applied to during cbm development Method, it is characterized in that step (2) is specially:
(a) Geologic modeling is carried out using the drilling well in target area and seismic data, obtains the geological layering depth number of each eyeball According to and formation lithology, each formation lithology change in apparent resistivity range is demarcated;
(b) under the constraint of above-mentioned geological model, each eyeball apparent resistivity curve obtained to different time sections quantifies Inversion Calculation;The interference of variable, systematic error, geography information in gatherer process is removed, obtains the achievement corresponding to each period Apparent resistivity data body, the difference of each achievement apparent resistivity data body should be regarded as during coal-bed gas exploitation caused by different times Caused by coal seam reservoirs SEA LEVEL VARIATION single factors;
(c) each achievement apparent resistivity data body is subjected to spacial migration according to observation time, obtains time shift apparent resistivity space number According to body, time shift apparent resistivity spatial data body is recorded as:Time shift controlled-source audiomagnetotelluric apparent resistivity data body;
(d) in step (a) under the constraint of geological model, by each time shift apparent resistivity spatial data body that step (c) obtains by It is divided into along layer time shift apparent resistivity spatial data body, is recorded as according to geological layering depth data and formation lithology:Along layer time shift Controlled-source audiomagnetotelluric apparent resistivity data body;
(e) difference is asked for by adjacent time to the edge layer time shift apparent resistivity data body that step (d) obtains, obtains and regarded along layer difference Resistivity data body, is recorded as:Along layer time shift controlled-source audiomagnetotelluric apparent resistivity variance data body.
7. the time shift controlled-source audiomagnetotelluric detection side according to claim 6 applied to during cbm development Method, it is characterized in that step (e) is additionally included in before asking for difference by adjacent time along layer time shift apparent resistivity data body, to step Suddenly the edge layer apparent resistivity spatial data body that (d) is obtained carries out different depth correction and flattening technical finesse.
8. the time shift controlled-source audiomagnetotelluric detection side according to claim 1 applied to during cbm development Method, it is characterized in that further comprising the steps of in step (3):Exclude the data exception point in data acquisition and exceptions area;It will It is projected on plan view along layer difference apparent resistivity data body section, obtains corresponding isogram, draw a circle to approve coal seam reservoirs water system Change the range on boundary, obtain coal-bed gas exploitation area coal seam reservoirs water system dynamic change situation and cbm development process pressure drop wave And range.
9. the time shift controlled-source audiomagnetotelluric detection side according to claim 6 applied to during cbm development Method, it is characterized in that the method for the interference of variable is in removal gatherer process described in step (b):The eyeball that will be acquired twice Eyeball on apparent resistivity curve is compared, and to be acquired again to there are the eyeballs of notable difference;The removal acquisition The method of the interference of systematic error is in the process:While production zone acquires, to unminding region using identical method Carry out the acquisition of eyeball apparent resistivity twice, it is poor that the eyeball apparent resistivity acquired twice is made, obtain systematic error and The interference of removal system error in quantitative inversion calculating process;The method of the interference of geography information is in removal gatherer process:It is right The each geographic basis of eyeball twice is compared, and excludes to generate acquisition eyeball condition the corresponding geographic basis of interference.
10. the time shift controlled-source audiomagnetotelluric detection according to claim 8 applied to during cbm development Method, it is characterized in that the data exception point is the data point of error range allowed more than industry standard, the exception Area is the region that multiple data exception points are connected.
CN201710055199.1A 2017-01-24 2017-01-24 Applied to the time shift controlled-source audiomagnetotelluric detection method during cbm development Pending CN108227006A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710055199.1A CN108227006A (en) 2017-01-24 2017-01-24 Applied to the time shift controlled-source audiomagnetotelluric detection method during cbm development

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710055199.1A CN108227006A (en) 2017-01-24 2017-01-24 Applied to the time shift controlled-source audiomagnetotelluric detection method during cbm development

Publications (1)

Publication Number Publication Date
CN108227006A true CN108227006A (en) 2018-06-29

Family

ID=62656461

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710055199.1A Pending CN108227006A (en) 2017-01-24 2017-01-24 Applied to the time shift controlled-source audiomagnetotelluric detection method during cbm development

Country Status (1)

Country Link
CN (1) CN108227006A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111965712A (en) * 2020-10-21 2020-11-20 国网江西省电力有限公司电力科学研究院 Method for correcting static effect of controllable source audio magnetotelluric method
CN111983689A (en) * 2020-08-27 2020-11-24 青海省第四地质勘查院 Near-source electromagnetic seismoelectric joint GR inversion method
CN112379449A (en) * 2020-10-30 2021-02-19 中国石油天然气集团有限公司 Method and device for processing electromagnetic data of controllable source

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103492909A (en) * 2011-02-25 2014-01-01 弗朗西斯科·奎德拉·卡纳莱斯 Method of dynamic seismic prospecting
CN104122592A (en) * 2014-07-31 2014-10-29 中国地质大学(武汉) Time shift magnetotelluric signal acquisition and inversion method
US20150034307A1 (en) * 2013-08-01 2015-02-05 Conocophillips Company Dynamic in-situ measurement of reservoir wettability
CN105717547A (en) * 2015-12-22 2016-06-29 吉林大学 Anisotropy medium magnetotelluric meshless value simulation method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103492909A (en) * 2011-02-25 2014-01-01 弗朗西斯科·奎德拉·卡纳莱斯 Method of dynamic seismic prospecting
US20150034307A1 (en) * 2013-08-01 2015-02-05 Conocophillips Company Dynamic in-situ measurement of reservoir wettability
CN104122592A (en) * 2014-07-31 2014-10-29 中国地质大学(武汉) Time shift magnetotelluric signal acquisition and inversion method
CN105717547A (en) * 2015-12-22 2016-06-29 吉林大学 Anisotropy medium magnetotelluric meshless value simulation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡祖志等: "涩北气藏时移大地电磁监测技术可行性研究", 《石油地球物理勘探》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111983689A (en) * 2020-08-27 2020-11-24 青海省第四地质勘查院 Near-source electromagnetic seismoelectric joint GR inversion method
CN111965712A (en) * 2020-10-21 2020-11-20 国网江西省电力有限公司电力科学研究院 Method for correcting static effect of controllable source audio magnetotelluric method
CN111965712B (en) * 2020-10-21 2021-03-02 国网江西省电力有限公司电力科学研究院 Method for correcting static effect of controllable source audio magnetotelluric method
CN112379449A (en) * 2020-10-30 2021-02-19 中国石油天然气集团有限公司 Method and device for processing electromagnetic data of controllable source
CN112379449B (en) * 2020-10-30 2023-05-26 中国石油天然气集团有限公司 Processing method and device for electromagnetic data of controllable source

Similar Documents

Publication Publication Date Title
CN105651676B (en) A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern
Norden et al. Geological modelling of the Triassic Stuttgart Formation at the Ketzin CO2 storage site, Germany
CN101408624B (en) Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer
CN110795878B (en) Tunnel water inflow prediction method
US9857500B2 (en) Method to constrain a basin model with curie depth
Coccia et al. Application of Refraction Microtremor (ReMi) technique for determination of 1-D shear wave velocity in a landslide area
Imposa et al. 3D subsoil model of the San Biagio ‘Salinelle’mud volcanoes (Belpasso, SICILY) derived from geophysical surveys
Jianjun et al. An application of the high-density electrical resistivity method for detecting slide zones in deep-seated landslides in limestone areas
CN102590874A (en) Method for detecting ground surface crack of upland coal-mining subsidence paddy field
CN108227006A (en) Applied to the time shift controlled-source audiomagnetotelluric detection method during cbm development
CN115097534A (en) Method and system for identifying space morphology of dry hot rock fracturing transformation
Baroncini‐Turricchia et al. Integrating MRS data with hydrologic model‐Carrizal Catchment (Spain)
CN107272080A (en) Site standardization analysis method for underground water leakage channel investigation
Wilson et al. Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site
CN107945271A (en) Three-dimensional pressure field modeling method based on geological mass tracking
MacLennan et al. Using depth to surface resistivity for mapping CO2
Tilbury et al. Pluto 4D—Australia’s first 4D over a gas field is an outstanding success
Trichandi et al. Shear-wave velocity imaging of weathered granite in La Campana (Chile) from Bayesian inversion of micro-tremor H/V spectral ratios
Ekine et al. Delineation of hydrocarbon bearing reservoirs from surface seismic and well log data (Nembe Creek) in Niger Delta oil field
Iovenitti et al. EGS exploration methodology project using the Dixie Valley geothermal system, Nevada, status update
Thabit et al. Predicting aquifer characteristics to enhance hydrogeophysical model in fractured/karstified rocks of Dammam Formation at Bahr Al-Najaf Basin, Iraq
Chen et al. Description technology of fractured-vuggy carbonate reservoir in halahatang oil field, tarim basin—take the ha 7 test area as an example
Mariita Application of geophysics to geothermal energy exploration and monitoring of its exploitation
Levon et al. A multi-disciplinary modeling approach to determine economic viability of infill horizontal well development in the Green River Basin
Abdrakhimova Improving groundwater flow model parametrization techniques

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180629