Coal seam seismic exploration data analysis method and system

Abstract

The present application is suitable for coal seam exploration technical field, provide a kind of coal seam seismic exploration data analysis method and system, the method includes: the historical seismic exploration data of micro logging and historical record data are analyzed, obtain the interference data related to micro logging and obtain the excitation preparation data of the micro logging, the interference data are used to characterize the interference of first target to seismic exploration, the first target includes predictable interference target and / or dynamic interference target;Based on the interference data and the excitation position, determine the first time period of the dynamic interference target through the interference range, in the second time period, based on the interference data, determine the future interference period of the micro logging by the predictable interference target, the whole interpretation process greatly reduces the difficulty of data interpretation and analysis, so as to provide strong geological basis for the safety production of large dip angle coal seam and other mining areas.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of coal seam exploration, and particularly relates to a coal seam seismic exploration data analysis method and system. BACKGROUND

[0002] Seismic exploration can provide technical support for formulating coal mining design. Seismic exploration is to artificially excite (or excite in the well) seismic waves to the underground, and record the reflected seismic waves from the underground on the ground surface or stratum through instrument equipment. The obtained seismic record is processed through a series of processes, and finally the seismic profile is obtained.

[0003] It is particularly important to obtain the geological information of the coal seam through seismic exploration for the mining of the mine, especially for the coal seam with large dip angle. In the related technology, the coal seam geological exploration data is analyzed to analyze the distribution of the coal seam. However, the analysis ability of the seismic exploration data in the related technology is limited, and the ideal analysis result cannot be obtained, thereby affecting the quality of the geological data. SUMMARY

[0004] The purpose of the embodiment of the present application is to provide a coal seam seismic exploration data analysis method and system, which aims to solve the problems in the background technology.

[0005] The embodiment of the present application is realized as follows. In a first aspect, a coal seam seismic exploration data analysis method comprises:

[0006] Analyzing the historical seismic exploration data and the historical record data of the micro logging, obtaining the interference data related to the micro logging and the excitation preparation data of the micro logging, wherein the excitation preparation data comprises an excitation position, an excitation time period and an excitation amount of explosive, and the interference data is used to represent the interference of a first target on the seismic exploration, and the first target comprises a predictable interference target and / or a dynamic interference target;

[0007] Based on the interference data and the excitation position, a first time period in which the dynamic interference target passes through the interference range is determined, and in a second time period, based on the interference data, a future interference time period of the micro logging by the predictable interference target is determined, wherein the interference range is determined based on the micro logging, and the second time period covers the excitation time period and is at least ahead of a preset time length;

[0008] According to the first time period, the future interference time period and the historical seismic exploration data, a target item of the excitation preparation data is updated to obtain first excitation data, and first seismic exploration data collected according to the first excitation data is obtained, wherein the target item comprises the excitation time period and the excitation amount of explosive;

[0009] The first seismic exploration data is analyzed by using a preset denoising scheme, wherein the preset denoising scheme is used to remove the interference caused by the first target.

[0010] In a second aspect, the present application further provides a coal seam seismic exploration data analysis system, which applies the coal seam seismic exploration data analysis method described above, and comprises: a first analysis module, configured to analyze historical seismic exploration data and historical record data of micro logging, obtain interference data related to the micro logging, and obtain firing preparation data of the micro logging, wherein the firing preparation data comprises a firing position, a firing time period and a firing explosive amount, the interference data is used to represent the interference of the first target on seismic exploration, and the first target comprises a predictable interference target and / or a dynamic interference target;

[0011] A time period determination module is configured to determine a first time period during which the dynamic interference target passes through an interference range based on the interference data and the firing position, and determine a future interference time period of the predictable interference target on the micro logging based on the interference data during a second time period, wherein the interference range is determined based on the micro logging, and the second time period covers the firing time period and is at least a preset length of time in advance;

[0012] A firing data updating module is configured to update a target item of the firing preparation data according to the first time period, the future interference time period and the historical seismic exploration data, obtain first firing data, and obtain first seismic exploration data collected according to the first firing data, wherein the target item comprises the firing time period and the firing explosive amount;

[0013] A second analysis module is configured to analyze the first seismic exploration data by using a preset denoising scheme, wherein the preset denoising scheme is used to remove the interference caused by the first target.

[0014] The coal seam seismic exploration data analysis method and system provided by the present application obtain interference data and determine related time information from the analysis of historical seismic exploration data, update firing preparation data in combination with historical seismic exploration data, thereby obtaining first firing data, and thus theoretically obtain first seismic exploration data with high quality, and process the first seismic exploration data by using means such as near trace cutting in the preset denoising scheme, thereby further improving the analysis quality of coal seam seismic data, obtaining geological information with high precision and high signal-to-noise ratio, greatly reducing the difficulty of data analysis in the entire analysis process, and thus providing strong geological basis for the safety production of large-dip-angle coal seams and other mining areas, and providing technical support for coal mining design. BRIEF DESCRIPTION OF DRAWINGS

[0015] Figure 1It is a main flow chart of a coal seam seismic exploration data analysis method.

[0016] Figure 2 It is a flow chart of obtaining interference data related to micro logging in the coal seam seismic exploration data analysis method.

[0017] Figure 3 It is a flow chart of determining a first time period in the coal seam seismic exploration data analysis method.

[0018] Figure 4 It is a type and development condition chart of interference waves in the coal seam seismic exploration data analysis method.

[0019] Figure 5 It is a main structure chart of a coal seam seismic exploration data analysis system. DETAILED DESCRIPTION

[0020] In order to make the purpose, technical scheme and advantages of the present application clearer, the present application is further described in detail below in combination with the drawings and examples. It should be understood that the specific examples described herein are only used to explain the present application and do not limit the present application.

[0021] The specific implementation of the present application is described in detail below in combination with specific examples.

[0022] The coal seam seismic exploration data analysis method and system provided by the present application solve the technical problems in the background art.

[0023] As shown in Figure 1 FIG. 1 is a main flow chart of a coal seam seismic exploration data analysis method provided by an embodiment of the present application, and the coal seam seismic exploration data analysis method includes:

[0024] Step S10: analyze historical seismic exploration data and historical record data of micro logging, obtain interference data related to micro logging and acquire excitation preparation data of the micro logging, wherein the excitation preparation data includes excitation position, excitation time period and excitation explosive quantity, the interference data is used to represent the interference of a first target on seismic exploration, and the first target includes a predictable interference target and / or a dynamic interference target;

[0025] Specifically, the historical seismic exploration data records the characteristic information of the interference target, such as the predictable interference target: having a global array response to all the receivers in the coverage; the historical data records the interference process, based on the analysis of both, the first target, that is, the predictable interference target and / or the dynamic interference target, can be determined, the predictable interference target generally refers to the interference source with a fixed rule or can be predicted in advance, the occurrence time, frequency, intensity or spatial position of which can be predicted based on the historical data or external information, such as meteorological events such as rainfall, strong wind, thunderstorm, etc., vehicle vibration on the road around the mining area (such as the daily 7:00-9:00 commuting period), with seasonal / periodic regularity; while the dynamic interference target refers to the interference source that appears randomly and cannot be predicted in advance, which has no fixed rule in time, space or frequency spectrum characteristics, and needs to be monitored in real time or analyzed and identified afterwards, such as random operation of excavators, drilling machines, etc., which produce wideband impact vibration, etc., both of which can interfere with seismic exploration; after the first target is determined, the interference data can be obtained based on the characteristics of the first target, and the excitation preparation data of the microlog is obtained, the interference rule is analyzed based on the historical data, which is helpful to accurately and quickly obtain the interference data, and is also helpful to the effective updating and optimization of the excitation preparation data in the subsequent.

[0026] Step S11: based on the interference data and the excitation position, determining a first time period in which the dynamic interference target passes through the interference range, and based on the interference data, determining a future interference period of the predictable interference target to the microlog in a second time period, wherein the interference range is determined based on the microlog, and the second time period covers the excitation period and at least a preset time in advance;

[0027] Specifically, the interference data is used to record the interference process or interference trend of the first target to the microlog, based on the interference data of the dynamic interference target and the excitation position, the time when the dynamic interference target enters the interference range, that is, the first time period, can be predicted, and the dynamic interference target appears in a relatively far position and needs to be monitored; for the predictable target, there is a certain rule, which can be monitored starting from a preset time in advance of the excitation period, and the preset time is at least the preparation time of the excitation, and based on the interference data of the predictable interference target, the future interference period of the predictable interference target to the microlog can be determined, the determination of the first period and the future interference period can be directly used for updating and optimizing the excitation preparation data in the subsequent.

[0028] Step S12: updating the target item of the excitation preparation data according to the first time period, the future interference period and the historical seismic exploration data, obtaining first excitation data, and acquiring first seismic exploration data collected according to the first excitation data, wherein the target item includes the excitation period and the excitation amount;

[0029] Specifically, the first target can interfere with the excitation process of the microlog under the excitation preparation data. By the first time period, the future interference period and the historical seismic exploration data, the target term can be updated, such as the excitation period, to avoid or reduce the interference of the first target, to realize the preset avoidance strategy combined with the historical seismic exploration data, and if necessary, to combine other micrologs under the same or similar conditions with good interpretation effect of the seismic exploration data to change the excitation amount in the target microlog, thereby updating the excitation preparation data combined with the target term, to obtain reasonable first excitation data for execution, and then to obtain the first seismic exploration data collected according to the first excitation data, which helps to improve the quality of coal seam seismic data.

[0030] Step S13: analyzing the first seismic exploration data by a preset denoising scheme, wherein the preset denoising scheme is used to remove the interference caused by the first target.

[0031] Specifically, the foregoing is based on the interpretation of historical seismic exploration data to obtain interference data and determine related time information, and then update the excitation preparation data combined with the historical seismic exploration data to obtain the first excitation data, so that theoretically, the first seismic exploration data of high quality can be obtained. Further, the near trace cut-off method in the preset denoising scheme can be used to analyze the latest first seismic exploration data, thereby further improving the interpretation quality of coal seam seismic data to obtain high-precision and high-signal-to-noise ratio geological information, greatly reducing the difficulty of data interpretation and analysis, and providing technical support for coal mining design.

[0032] In application, the embodiment obtains interference data related to the micro logging and obtains the firing preparation data of the micro logging by analyzing historical seismic exploration data and historical record data of the micro logging, analyzes interference rules through historical data, helps to accurately and quickly obtain interference data, and helps to effectively update and optimize the firing preparation data subsequently; based on the interference data and the firing position, a first time period of the dynamic interference target passing through the interference range is determined, and in a second time period, based on the interference data, a future interference time period of the predictive interference target to the micro logging is determined, the determination of the first time period and the future interference time period can be directly used for the update and optimization of the firing preparation data subsequently; the target item of the firing preparation data is updated according to the first time period, the future interference time period and the historical seismic exploration data, first firing data is obtained, and first seismic exploration data collected according to the first firing data is obtained, which can realize the preset avoidance strategy in combination with the historical seismic exploration data, and other micro loggings under the same or similar conditions with good interpretation effect are combined as needed, and then the target item is combined to update the firing preparation data, so as to obtain reasonable first firing data to execute to obtain the first seismic exploration data collected according to the first firing data, which helps to improve the quality of coal seam seismic data; the first seismic exploration data is analyzed through the preset denoising scheme, which can suppress limited noise interference while maximizing the retention of coal seam effective signals to ensure the reliability of structure interpretation; the present application obtains interference data and determines related time information from the interpretation of historical seismic exploration data, and then updates the firing preparation data in combination with the historical seismic exploration data to obtain the first firing data, so that theoretically, first seismic exploration data with high quality can be obtained, and the first seismic exploration data is processed by means such as near trace cutting in the preset denoising scheme, so as to further improve the interpretation quality of coal seam seismic data, so as to obtain geological information with high precision and high signal-to-noise ratio, the whole interpretation process greatly reduces the difficulty of data interpretation and analysis, thereby providing strong geological basis for safety production in large-dip-angle coal seams and other mining areas, and providing technical support for coal mining design.

[0033] As Figure 2 shown, as a preferred embodiment of the present application, the analysis of the historical seismic exploration data and the historical record data of the micro logging to obtain the interference data related to the micro logging comprises:

[0034] Step S101: performing spectrum feature recognition and geophone array analysis on the historical seismic exploration data to determine a suspicious interference time period caused by suspicious factors, the suspicious factors including environmental vibration and weather change;

[0035] Specifically, the dynamic interference (such as vehicles, engineering machinery), the occurrence time and position are unpredictable (such as temporary construction vehicles, drilling rig moving), low-frequency vibration (vehicle engine), high-frequency impact (mechanical collision), the spectrum is expressed as a burst peak or wideband noise, the interference energy decays with distance, and can only affect local geophones; it is expressed as a burst spectrum sudden increase + local array response.

[0036] The predictable interference (such as weather), periodicity or predictability: such as thunderstorm (wideband electromagnetic pulse), strong wind (low-frequency environmental vibration), rainfall (surface wave attenuation), thunderstorm: instantaneous high-frequency electromagnetic interference (spectrum energy sudden increase); strong wind: continuous low-frequency environmental noise (0.1-10Hz); rainfall: the change of surface medium causes the distortion of reflected signal, global influence: weather interference usually covers the entire exploration area, and all geophones can be simultaneously affected; the characteristic of predictable interference is expressed as: periodic spectrum characteristic + global array response, when the corresponding combination characteristics are identified based on historical seismic exploration data, the corresponding suspicious interference period is determined.

[0037] Step S102: identifying based on historical record data, determining a first target in the suspicious interference period and obtaining a feature identifier of the first target;

[0038] In the case where the future interference period has been determined, it is not necessary to make a general identification on the historical record data, but an identification on the historical record data in the suspicious interference period can be directly made, for the dynamic target, the historical record data includes image data, and the feature identifier includes a standard state photo of the dynamic interference target, numbering information, equipment code and the like; for the predictable interference target, the historical record data is report information recorded by a professional department, such as actual meteorological information, and the feature identifier includes a WMO code given by the World Meteorological Organization or other defined code.

[0039] Step S103: obtaining the interference data of the first target according to the feature identifier.

[0040] Specifically, the interference data of the dynamic interference target is derived from the latest image data near the microlog and the fixed-point intelligence data set outside the interference range, and the interference data of the predictable interference target is the local recorded prediction data or the data from a professional department, such as the environmental interference data from the meteorological department system.

[0041] It can be understood that, by performing spectrum characteristic identification and geophone array analysis on historical seismic exploration data, a suspicious interference period caused by a suspicious factor is determined, then an identification based on historical record data is performed, a first target in the suspicious interference period is determined, and a feature identifier of the first target is obtained, which improves the efficiency of data identification and ensures that interference data can be obtained accordingly.

[0042] As Figure 3 shown, as a preferred embodiment of the present application, based on the interference data and the excitation position, determining the first time period of the dynamic interference target passing through the interference range, in the second time period, based on the interference data, determining the future interference period of the predictive interference target to the microlog includes:

[0043] Step S111: according to the feature identification of the dynamic interference target, determining the position of the dynamic interference target and the first interference degree;

[0044] Specifically, the feature identification includes the standard state photo, number information, equipment code and the like of the dynamic interference target, according to the feature identification, the single frame image meeting the preset similarity threshold between the standard state photo can be determined in the image data, the position of the dynamic interference target can be determined according to the shooting position in the single frame image, or the positioning information of the dynamic interference target can be directly determined according to the equipment code, the dynamic interference target includes trucks, engineering vehicles, agricultural machinery and the like, the first interference degree can be determined by quantitative and qualitative analysis, the quantitative can be realized by fixed-point monitoring, such as vehicle low-frequency (5-15Hz) energy proportion, >40% is significant interference; for engineering machinery, harmonic distortion degree (THD), THD>15% needs to be warned, etc.

[0045] Step S112: according to the first interference degree and the excitation position, reversely determining the interference range centered on the microlog;

[0046] Specifically, according to the first interference degree, the farthest distance that can just propagate to the microlog position and ignore the noise is reversely calculated, the distance is converted into the straight line distance on the ground, then taking the microlog as the center and the distance as the radius, the interference range of the dynamic interference target can be determined, for example, the interference source is excavator, according to the following steps: Wherein P(r): noise amplitude at distance r, P0 is the initial amplitude of the interference source, α represents the stratum absorption coefficient (related to frequency, unit: dB / m), reverse calculation, when P(r) drops to threshold value Pth, rmax value, then according to d ground = r max ·cosθ can calculate the ground projection distance, assuming uniform medium, θ is the propagation angle; interference source: excavator (initial amplitude initial amplitude P0=1000μPa, main frequency 30Hz), absorption coefficient α=0.05==0.05dB / m (medium hard clay), noise threshold: Pth=50Pth=50μPa (corresponding to SNR=10dB), solving rmax≈120m, ground projection: dground=120×cos45°≈85m, that is, the interference range is a circular area with the microlog as the center and the radius of 85 meters.

[0047] Step S113: According to the position and moving information of the dynamic interference target, determine the moving direction and moving speed of the dynamic interference target;

[0048] Specifically, the average speed is calculated according to the moving path and time length (or the instantaneous speed is directly taken as the moving speed when the difference between the average speed and the instantaneous speed is too large for a long time), and the moving direction is determined according to the path.

[0049] Step S114: According to the moving direction, the first position of the dynamic interference target and the interference range, determine the distance information from the first position to the nearest edge interference point and the farthest edge interference point in the interference range, and based on the distance information and the moving speed, determine the first time period of the dynamic interference target passing through the interference range.

[0050] Specifically, the nearest edge interference point and the farthest edge interference point are the two intersection points of the moving direction and the interference range, for example, the two intersection points of a straight line and a circular area, the first position (for example, the current position or the updated position) of the dynamic interference target is known, the distance to the two intersection points is known, and the moving speed is known, so the time taken from the first position of the dynamic interference target to the nearest edge interference point and the farthest edge interference point can be calculated, and the current time during the calculation is superimposed to determine the first time period, which is updated.

[0051] The above can accurately determine the interference of the microlog based on the moving and interference information of the dynamic interference target, thereby facilitating the updating of the target item of the subsequent excitation preparation data.

[0052] As a preferred embodiment of the present application, based on the interference data and the excitation position, the first time period of the dynamic interference target passing through the interference range is determined, and in the second time period, based on the interference data, the future interference time period of the predictive interference target on the microlog is determined, which includes:

[0053] Step: In the second time period, the future interference time period is determined according to the interference data, and the future interference time period is divided according to the second interference degree, wherein the interference data of the predictive interference target includes external prediction data. Since the second time period is limited, the second interference degree at this time is the interference degree of the predictive interference target.

[0054] It should be understood that the interference data of the predictive interference target is derived from external prediction data, such as environmental interference data given by a meteorological department system, which generally has an interference degree level, such as the level of wind and rain and the amount of rainfall.

[0055] As a preferred embodiment of the present application, the target item of the shot preparation data is updated according to the first time period, the future interference time period and the historical seismic exploration data, to obtain first shot data, and first seismic exploration data collected according to the first shot data is obtained, wherein the target item includes the shot time period and the shot dose includes:

[0056] Step: judging whether the future interference time period overlaps in time with the shot time period;

[0057] Step: if it is determined that the future interference time period overlaps in time with the shot time period, a non-future interference time period is taken as a new candidate time period of the shot time period, a new shot time period is determined based on the candidate time period, it is judged whether the new shot time period overlaps in time with the first time period, and if yes, the process jumps to the pause instruction step;

[0058] Specifically, the first time period is an activity time period of a dynamic interference target (such as a vehicle or a machine), the future interference time period is a high-noise time period caused by environmental vibration (such as wind and rain), and the shot time period is a time window in which seismic exploration is planned to be performed. For example, the shot time period is 09:00-9:50. Half an hour before the shot, it is learned that the predictable interference target changes (the prediction is inaccurate), at this time, the environmental future interference time period (heavy rain) is 09:35-10:00 (overlapping). Thus, the new shot time period is determined to be 15:00-16:00, and the dynamic future interference time period (excavator operation or vehicle driving in) is 11:30-12:00 (not overlapping with the new shot time period). The processing process is to perform the shot process according to the new shot time period. If it is determined that the new shot time period is 11:50-12:30 and 14:20-15:20 (overlapping with the new shot time period), a stop instruction can be sent to the excavator, and at least the excavator is required to evacuate (out of the interference range) or stop a period of time before 15:00. The above can avoid data pollution caused by environmental vibration and ensure the acquisition quality. Alternatively, when the candidate time periods are few and the future interference time period exceeds a preset time length, the future interference time period is divided according to a second interference degree. The time when the second interference degree meets a preset slight level and avoids the dynamic future interference time period is selected as the new shot time. The period of time is avoided to be delayed for too long to affect the term of the new shot time. The rainfall of the preset slight level is small, a geophone with a higher natural frequency can be used to reduce low-frequency interference, and / or the shot capacity is updated in the subsequent process, and the shot energy is appropriately increased (to compensate for the rainwater absorption loss).

[0059] Step: if it is determined that the future interference time period does not overlap in time with the shot time period, the original shot time period is kept, and it is judged whether the first time period overlaps in time with the shot time period, and if yes, the process jumps to the pause instruction step;

[0060] Specifically, when the environmental vibration does not affect the firing time, the original firing period is directly maintained, and if the future interference period of the dynamic interference target overlaps with the firing time, the pause indication step should be performed, such as sending a stop command to the excavator, requiring early evacuation (out of the interference range) or stopping.

[0061] The pause indication step includes: pausing the dynamic interference target, that is, requiring early evacuation (out of the interference range) or stopping, and making a time-limited movement indication based on the new firing period or the firing period.

[0062] Specifically, the pause requires full takeover of the dynamic interference target, while the time-limited movement requires completion before the preparation time of the firing period.

[0063] It should be understood that the above can accurately avoid environmental interference (such as weather), actively control dynamic interference (such as vehicles), minimize the risk of data pollution, improve construction efficiency, and be suitable for complex environments such as mining areas and cities, so that the quality of coal seam seismic data can be significantly improved before formal interpretation, and the reliability of subsequent structural interpretation can be ensured.

[0064] As a preferred embodiment of the present application, the target item of the firing preparation data is updated according to the first time period, the future interference period, and the historical seismic exploration data, to obtain first firing data, and first seismic exploration data collected according to the first firing data is obtained, wherein the target item includes the firing period and the firing charge amount.

[0065] Under the new firing period, or under the condition of maintaining the original firing period, the first historical seismic exploration sub-data that meets the preset similarity condition of the firing position is determined in the historical seismic exploration data of the microlog.

[0066] Specifically, under the new firing period, or maintaining the original firing period, a new preparation parameter condition is represented, one of the two cases is taken, and the preset similarity condition includes historical seismic exploration data (sub-data set) in the same coal seam inclination region, and ensures that the parameters such as the observation system (geophone arrangement, sampling rate) of the historical data are consistent with the current microlog (generally from other micrologs).

[0067] The firing charge amount is modified according to the signal-to-noise ratio of the first historical seismic exploration sub-data.

[0068] Optionally, the modification of the firing charge amount according to the signal-to-noise ratio of the first historical seismic exploration sub-data includes:

[0069] Compare the change of signal-to-noise ratio of the first historical seismic exploration sub-data, and determine whether the surface wave changes, and determine the change amount of the corresponding excitation amount of the explosive.

[0070] Specifically, historical seismic exploration data (sub-data set) in the same coal seam inclination region as the target micro logging is extracted, and the parameters such as the observation system (geophone arrangement, sampling rate) of the historical data are ensured to be consistent with the current micro logging, the signal-to-noise ratio SNR is calculated based on the common shot point gather (CSP) or common receiver point gather (CRP) of the historical seismic exploration data, SNR = 10lg{effective signal energy (coal seam reflection time window) / noise energy (first arrival front or reflection rear time window)}; the excitation amount and well depth of each shot data are labeled, and a three-dimensional scatter plot (well depth, explosive amount, SNR) is drawn, a trend surface (such as polynomial regression) is fitted, and then the micro logging explosive amount is determined according to the fitted trend surface; or the historical seismic exploration data is input into a neural network model that has been trained, and finally the micro logging explosive amount can also be obtained, and the calculation results are as follows: when the SNR increases with the increase of the explosive amount, the explosive amount should be appropriately increased (such as increasing by a first preset increment, for example, +10%-20%), according to the following: the increase of the explosive amount can enhance the seismic wave energy and effectively suppress random noise, thereby improving the signal-to-noise ratio SNR; when the SNR increases with the increase of the explosive amount but the surface wave is enhanced, the explosive amount is increased by a small amplitude (such as increasing by a second preset increment, the second preset increment is lower than the first preset increment, for example, +5%-10%) and is matched with FK filtering; according to the following: a large explosive amount can excite stronger low-frequency surface waves, and FK filtering is needed to suppress the interference of surface waves while improving the signal energy, so as to balance the energy and noise, when the SNR is saturated or decreases at a certain explosive amount, the historical best explosive amount is not adjusted; when the SNR is always low (lower than a preset signal-to-noise ratio), the well depth is insufficient, which leads to energy loss, at this time, the well depth should be increased instead of the explosive amount; the above can optimize the explosive amount through historical data driving, for example, improve the quality of the micro logging seismic data.

[0071] As shown in Figure 4 The types and development of interference waves are shown, so that effective measures can be taken to suppress various interference waves to the greatest extent. The main interference waves in the study area are surface waves, refracted waves, microseisms and acoustic waves; in view of the problem, a detailed embodiment is given, which avoids the corresponding high-frequency random noise and low-frequency waves and has directional low-frequency vibration characteristics on the basis of avoiding the first target interference, further, the first seismic exploration data is analyzed through a preset denoising scheme, wherein the preset denoising scheme is used to remove the interference caused by the first target.

[0072] Based on the preset order in the preset denoising scheme, the domain transformation and physical cutting of the first seismic exploration data are performed.

[0073] Specifically, the first seismic exploration data is processed in the order of acoustic wave, refracted wave and surface wave. For the acoustic wave, mute is first used to directly cut off the data of near offset (e.g. <300m), and then FK filtering is used to cut off the low-velocity energy (velocity <500m / s) in the FK domain. For the refracted wave, first arrival cut-off is used to remove the refracted wave within the first arrival time window (e.g. 0-300ms), and then Radon transform is used to separate linear events (refracted wave) and hyperbolic reflection wave. For the surface wave, FK filtering is first used to cut off the low-frequency sector energy (velocity <1000m / s, frequency <15Hz), and then wavelet transform is used to threshold process the low-frequency wavelet coefficients. The above can suppress limited noise interference (the operation in the foregoing embodiment has avoided a lot of unnecessary interference) and maximize the retention of coal seam effective signals, thereby ensuring the reliability of structural interpretation.

[0074] As shown in Figure 5 as another preferred embodiment of the present application, in another aspect, a coal seam seismic exploration data analysis system 200 comprises:

[0075] A first analysis module 210 is configured to analyze historical seismic exploration data and historical record data of micro logging, obtain interference data related to the micro logging, and obtain firing preparation data of the micro logging, wherein the firing preparation data comprises firing position, firing time period and firing explosive amount, and the interference data is used to represent the interference of a first target on seismic exploration, and the first target comprises a predictable interference target and / or a dynamic interference target.

[0076] A time period determination module 220 is configured to determine a first time period during which the dynamic interference target passes through an interference range based on the interference data and the firing position, and determine a future interference time period of the predictable interference target on the micro logging based on the interference data within a second time period, wherein the interference range is determined based on the micro logging, and the second time period covers the firing time period and is at least ahead of a preset time length.

[0077] A firing data updating module 230 is configured to update a target item of the firing preparation data according to the first time period, the future interference time period and the historical seismic exploration data, to obtain first firing data, and to obtain first seismic exploration data collected according to the first firing data, wherein the target item comprises the firing time period and the firing explosive amount.

[0078] A second analysis module 240 is configured to analyze the first seismic exploration data by using a preset denoising scheme, wherein the preset denoising scheme is used to remove the interference caused by the first target.

[0079] As a preferred embodiment of the present application, the first analysis module 210 is specifically used for:

[0080] Spectrum feature recognition and geophone array analysis are performed on historical seismic exploration data to determine a suspicious interference period caused by suspicious factors, including environmental vibration and weather change;

[0081] Identification is performed based on historical record data to determine a first target in the suspicious interference period and obtain a feature identifier of the first target;

[0082] The interference data of the first target is obtained according to the feature identifier.

[0083] It should be noted that the implementation methods and corresponding method steps of the above two embodiments completely correspond to the specific implementation description of the coal seam seismic exploration data analysis method in the foregoing embodiments, and will not be described here.

[0084] In order to be able to load the above-mentioned method and system to run smoothly, the system can include more or less components than described above, or combine certain components, or different components, for example, it can include input and output devices, network access devices, buses, processors and memories, etc.

[0085] The processor can be a central processing unit (CPU), and can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc. The above processor is the control center of the system, and is connected with various parts by various interfaces and lines.

[0086] The above only describes the preferred embodiments of the present application and is not used to limit the present application, and any modification, equivalent replacement and improvement within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims

1. A method for analyzing coal seam seismic exploration data, characterized in that, include: Analyze historical seismic exploration data and historical records of micrologging wells to obtain interference data related to the micrologging wells and acquire excitation preparation data of the micrologging wells. The excitation preparation data includes excitation location, excitation time period and excitation charge. The interference data is used to characterize the interference of a first target on seismic exploration. The first target includes predictive interference targets and / or dynamic interference targets. Based on the interference data and the excitation location, a first time period in which the dynamic interference target passes through the interference range is determined. In a second time period, based on the interference data, the future interference period of the predictive interference target on the micro-logging well is determined. The interference range is determined based on the micro-logging well, and the second time period covers the excitation period and is at least preset in advance. Based on the first time period, the future interference period, and the historical seismic exploration data, the target items of the excitation preparation data are updated to obtain the first excitation data, and the first seismic exploration data collected based on the first excitation data is obtained, wherein the target items include the excitation period and the excitation charge. The first seismic exploration data is analyzed by a preset denoising scheme, wherein the preset denoising scheme is used to remove interference caused by non-first targets.

2. The method for analyzing coal seam seismic exploration data according to claim 1, characterized in that, The analysis of historical seismic exploration data and historical records of micrologging was used to obtain interference data related to micrologging, including: The historical seismic exploration data is subjected to spectral feature identification and detector array analysis to determine the period of suspicious interference caused by suspicious factors, including environmental vibration and weather changes. Based on the historical data, the first target within the suspected interference period is identified and its feature identifier is obtained. The interference data of the first target is obtained based on the feature identifier.

3. The method for analyzing coal seam seismic exploration data according to claim 2, characterized in that, The process of determining a first time period in which the dynamic interference target passes through the interference range based on the interference data and the excitation location, and determining the future interference period of the predictive interference target on the micro-logging well based on the interference data during a second time period, includes: Based on the feature identifier of the dynamic interference target, the location of the dynamic interference target and the first degree of interference are determined; Based on the first level of interference and the excitation location, the interference range centered on the micro-logging well is determined in reverse. Based on the position and movement information of the dynamic interference target, determine the movement direction and movement speed of the dynamic interference target; Based on the direction of movement, the first position of the dynamic interference target, and the interference range, the distance information from the first position of the dynamic interference target to the nearest edge interference point and the farthest edge interference point in the interference range is determined. Based on the distance information and the movement speed, the first time period during which the dynamic interference target passes through the interference range is determined.

4. The method for analyzing coal seam seismic exploration data according to any one of claims 1-3, characterized in that, The process of determining a first time period in which the dynamic interference target passes through the interference range based on the interference data and the excitation location, and determining the future interference period of the predictive interference target on the micro-logging well based on the interference data during a second time period, includes: During the second time period, the future interference period is determined based on the interference data. The future interference period is distinguished according to the second interference level. The interference data of the predictive interference target includes external prediction data.

5. The method for analyzing coal seam seismic exploration data according to claim 2 or 3, characterized in that, The target items of the excitation preparation data are updated based on the first time period, the future interference period, and the historical seismic exploration data to obtain the first excitation data, and the first seismic exploration data collected based on the first excitation data is obtained. The target items include the excitation period and the excitation charge amount, which include: Determine whether the future interference period overlaps with the excitation period; If it is determined that the future interference period and the excitation period overlap, then the non-future interference period is used as a candidate period for the new excitation period. The new excitation period is determined based on the candidate period. It is then determined whether the new excitation period overlaps with the first period. If it does, the process proceeds to the pause instruction step. If it is determined that the future interference period and the excitation period do not overlap, the original excitation period is maintained, and it is determined whether the first time period and the excitation period overlap. If they do, the process jumps to the pause instruction step. The pause instruction step includes: issuing a pause instruction to the dynamic interference target, and issuing a time-limited movement instruction based on the new excitation period or the excitation period.

6. The method for analyzing coal seam seismic exploration data according to claim 5, characterized in that, The target items of the excitation preparation data are updated based on the first time period, the future interference period, and the historical seismic exploration data to obtain the first excitation data, and the first seismic exploration data collected based on the first excitation data is obtained. The target items include the excitation period and the excitation charge amount, which include: Under the new excitation period, or under the condition of maintaining the original excitation period, determine the first historical seismic exploration sub-data that meets the preset similarity conditions of the excitation location in the historical seismic exploration data of the micro-logging. The amount of excitation charge is modified based on the signal-to-noise ratio of the first historical seismic exploration data.

7. The method for analyzing coal seam seismic exploration data according to claim 6, characterized in that, The step of modifying the excitation charge based on the signal-to-noise ratio of the first historical seismic exploration data includes: Compare the changes in the signal-to-noise ratio of the first historical seismic exploration data and determine whether the surface wave has changed, and determine the corresponding change in the amount of excitation charge.

8. The method for analyzing coal seam seismic exploration data according to any one of claims 1-3, characterized in that, The first seismic exploration data is analyzed using a preset denoising scheme, wherein the preset denoising scheme is used to remove interference from non-target sources, including: Based on the preset order in the preset denoising scheme, the first seismic survey data is subjected to domain transformation and physical cut-off.

9. A coal seam seismic exploration data analysis system, characterized in that, The method for analyzing coal seam seismic exploration data as described in any one of claims 1-8 includes: The first analysis module is used to analyze historical seismic exploration data and historical records of micro-logging wells, obtain interference data related to micro-logging wells, and acquire the excitation preparation data of the micro-logging wells. The excitation preparation data includes excitation location, excitation time period, and excitation charge. The interference data is used to characterize the interference of a first target on seismic exploration. The first target includes predictive interference targets and / or dynamic interference targets. The time period determination module is used to determine, based on the interference data and the excitation location, a first time period in which the dynamic interference target passes through the interference range, and within a second time period, based on the interference data, to determine the future interference period of the predictable interference target on the micro-logging well, wherein the interference range is determined based on the micro-logging well, and the second time period covers the excitation period and is at least preset in advance. The excitation data update module is used to update the target items of the excitation preparation data according to the first time period, the future interference period, and the historical seismic exploration data to obtain the first excitation data and acquire the first seismic exploration data collected according to the first excitation data, wherein the target items include the excitation period and the excitation charge. The second parsing module is used to parse the first seismic exploration data using a preset denoising scheme, wherein the preset denoising scheme is used to remove interference caused by non-first targets.

10. The coal seam seismic exploration data analysis system according to claim 9, characterized in that, The first parsing module is specifically used for: Spectral feature identification and detector array analysis are performed on historical seismic exploration data to determine the period of suspicious interference caused by suspicious factors, including environmental vibrations and weather changes. Based on historical data, the first target within the suspected interference period is identified and its feature identifier is obtained. The interference data of the first target is obtained based on the feature identifier.

Images