CN105785435B - Seismic physical model experimental system and its analog data acquisition method and apparatus - Google Patents
Seismic physical model experimental system and its analog data acquisition method and apparatus Download PDFInfo
- Publication number
- CN105785435B CN105785435B CN201610187426.1A CN201610187426A CN105785435B CN 105785435 B CN105785435 B CN 105785435B CN 201610187426 A CN201610187426 A CN 201610187426A CN 105785435 B CN105785435 B CN 105785435B
- Authority
- CN
- China
- Prior art keywords
- physical model
- liquid
- discrete
- time series
- seismic physical
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 239000007788 liquid Substances 0.000 claims abstract description 123
- 238000002474 experimental method Methods 0.000 claims abstract description 105
- 230000005284 excitation Effects 0.000 claims abstract description 74
- 238000004088 simulation Methods 0.000 claims abstract description 40
- 238000005070 sampling Methods 0.000 claims description 42
- 238000012937 correction Methods 0.000 claims description 28
- 238000002604 ultrasonography Methods 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 239000012530 fluid Substances 0.000 claims 8
- 230000002463 transducing effect Effects 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 13
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000013480 data collection Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/22—Transmitting seismic signals to recording or processing apparatus
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (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
本发明提供了一种地震物理模型实验系统及其模拟数据采集方法和装置,涉及地震物理模型采集技术领域。该系统包括:地震物理模型、模拟采集控制系统、激发换能器、接收换能器、液体槽装置,地震物理模型的底面为平面,且其置于填充有液体的液体槽装置中;液体槽装置中的液体平面高于地震物理模型;激发换能器和接收换能器与液体平面完全接触,且接触平面与液体平面高度相同;激发换能器和接收换能器均与模拟采集控制系统连接。本发明可以解决当前传统的地震物理模拟实验超声波能量衰减,不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低,以及实验过程需要大量人工干预,地震物理模型数据采集时间长,采集效率低的问题。
The invention provides an earthquake physics model experiment system and its simulation data acquisition method and device, and relates to the technical field of earthquake physics model acquisition. The system includes: an earthquake physical model, an analog acquisition control system, an exciting transducer, a receiving transducer, and a liquid tank device. The bottom surface of the seismic physical model is a plane, and it is placed in a liquid tank device filled with liquid; the liquid tank The liquid plane in the device is higher than the seismic physical model; the excitation transducer and the receiving transducer are in full contact with the liquid plane, and the contact plane is at the same height as the liquid plane; both the excitation transducer and the receiving transducer are connected to the analog acquisition control system connect. The present invention can solve the problem of ultrasonic energy attenuation in the current traditional earthquake physics simulation experiment, which is not enough to detect the internal conditions of the earthquake physics model, and the signal-to-noise ratio of the collected data is low, and the experimental process requires a lot of manual intervention, and the data acquisition time of the earthquake physics model is long. long and low collection efficiency.
Description
技术领域technical field
本发明涉及地震物理模型采集技术领域,尤其涉及一种地震物理模型实验系统及其模拟数据采集方法和装置。The invention relates to the technical field of earthquake physics model acquisition, in particular to an earthquake physics model experiment system and a simulation data acquisition method and device thereof.
背景技术Background technique
当前,随着油气勘探技术的发展,我国在现阶段的油气勘探开发的重点目标区几乎都涉及起伏地表地形,如柴达木盆地、四川盆地和塔里木盆地等处。在起伏地表区域,由于地形的高程变化较大,使得野外地震勘探采集到的地震资料信噪比低、波场复杂,增加了后期地震资料处理与解释的周期与难度,因此增加了在起伏地表区域油气勘探的难度。At present, with the development of oil and gas exploration technology, the key target areas of oil and gas exploration and development in my country at the present stage almost all involve undulating surface terrain, such as Qaidam Basin, Sichuan Basin and Tarim Basin. In the undulating surface area, due to the large elevation change of the terrain, the signal-to-noise ratio of the seismic data collected by field seismic exploration is low and the wave field is complex, which increases the period and difficulty of later seismic data processing and interpretation. The difficulty of regional oil and gas exploration.
通过地震物理模型来进行地震物理模拟是地震勘探技术中一种有效的理论与实践相互验证的手段,对于研究起伏地表区域的地震波场传播规律,解决油气勘探中急需研究的复杂地表地震勘探技术问题等都具有十分重要的意义。当前传统的通过地震物理模型来进行地震物理模拟实验的方式如图1所示,通过计算机控制激发换能器101和接收换能器102分别放置于地震物理模型103表面某点,然后控制激发换能器101发射超声波(模拟野外地震勘探中的炸药震源激发),通过接收换能器102接收并记录从物理模型内部反射回来的超声波(模拟野外检波器接收地下地层界面反射回来的地震波场),从而在实验室就可以完成与野外实际勘探中相类似的采集过程。Seismic physical simulation through seismic physical model is an effective means of mutual verification of theory and practice in seismic exploration technology. It is useful for studying the propagation law of seismic wave field in undulating surface areas and solving complex surface seismic exploration technology problems that urgently need to be studied in oil and gas exploration. etc. are of great significance. The current traditional way of carrying out earthquake physics simulation experiments through earthquake physics models is shown in Figure 1. The excitation transducer 101 and the receiving transducer 102 are placed on a certain point on the surface of the earthquake physics model 103 by computer control, and then the excitation transducer is controlled to Transducer 101 emits ultrasonic waves (simulating the excitation of explosive source in the field seismic exploration), receives and records the ultrasonic waves reflected from the inside of the physical model through the receiving transducer 102 (simulating the field geophone receiving the seismic wave field reflected back from the interface of the underground formation), Therefore, the collection process similar to that in actual field exploration can be completed in the laboratory.
当前的激发换能器101和接收换能器102均是直接与地震物理模型103的表面接触的。而对于表面有起伏的地震物理模型,激发换能器和接收换能器与地震物理模型表面是硬接触的耦合方式,容易导致激发换能器和接收换能器与地震物理模型表面接触不良,造成超声波能量穿过接触面时发生衰减,超声波能量不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低。另外,由于地震物理模型表面起伏不平,激发换能器和接收换能器难以放置于地震物理模型表面,实验过程可能还需要大量的人工干预,造成地震物理模型数据采集时间长,采集效率低。Both the excitation transducer 101 and the reception transducer 102 are currently in direct contact with the surface of the seismic physics model 103 . For the seismic physical model with undulating surface, the excitation transducer and the receiving transducer are coupled with the surface of the seismic physical model in a hard contact coupling mode, which may easily lead to poor contact between the excitation transducer and the receiving transducer and the surface of the seismic physical model. The ultrasonic energy is attenuated when it passes through the contact surface, and the ultrasonic energy is not enough to detect the internal conditions of the seismic physical model, and the signal-to-noise ratio of the collected data is low. In addition, due to the uneven surface of the seismic physical model, it is difficult to place the excitation transducer and the receiving transducer on the surface of the seismic physical model. The experimental process may require a lot of manual intervention, resulting in long data acquisition time and low acquisition efficiency of the seismic physical model.
发明内容Contents of the invention
本发明的实施例提供一种地震物理模型实验系统及其模拟数据采集方法和装置,以解决当前传统的地震物理模拟实验超声波能量衰减,不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低,以及实验过程需要大量人工干预,地震物理模型数据采集时间长,采集效率低的问题。Embodiments of the present invention provide an earthquake physics model experiment system and its simulation data acquisition method and device, to solve the problem that the ultrasonic energy attenuation in the current traditional earthquake physics simulation experiment is not enough to detect the internal conditions of the earthquake physics model, and the collected data The signal-to-noise ratio is low, and the experimental process requires a lot of manual intervention, the seismic physical model data acquisition time is long, and the acquisition efficiency is low.
为达到上述目的,本发明采用如下技术方案:To achieve the above object, the present invention adopts the following technical solutions:
一种地震物理模型实验系统,包括:地震物理模型、模拟采集控制系统、激发换能器、接收换能器、液体槽装置,所述地震物理模型的底面为平面,且所述地震物理模型置于填充有液体的所述液体槽装置中;所述液体槽装置中的液体平面高于所述地震物理模型;所述激发换能器和接收换能器与所述液体平面完全接触,且接触平面与所述液体平面高度相同;所述激发换能器和接收换能器均与所述模拟采集控制系统连接;An earthquake physics model experimental system, comprising: an earthquake physics model, an analog acquisition control system, an exciting transducer, a receiving transducer, and a liquid tank device, the bottom surface of the earthquake physics model is a plane, and the earthquake physics model is placed In the liquid tank device filled with liquid; the liquid level in the liquid tank device is higher than the seismic physics model; the excitation transducer and the receiving transducer are in full contact with the liquid plane, and contact The plane is at the same height as the liquid plane; both the exciting transducer and the receiving transducer are connected to the analog acquisition control system;
所述激发换能器,用于向所述地震物理模型发射第一超声波信号;the excitation transducer is configured to transmit a first ultrasonic signal to the seismic physics model;
所述接收换能器,用于接收所述地震物理模型反射的第二超声波信号;The receiving transducer is used to receive the second ultrasonic signal reflected by the seismic physics model;
所述模拟采集控制系统,用于控制激发换能器的超声波发射端处于预先设置的激发点坐标,并控制接收换能器的超声波接收端处于预先设置的接收点坐标,并从所述第二超声波信号中采集时间信息,形成第一离散时间序列。The analog acquisition control system is used to control the ultrasonic transmitting end of the excitation transducer to be at the preset excitation point coordinates, and control the ultrasonic receiving end of the receiving transducer to be at the preset receiving point coordinates, and from the second Time information is collected from the ultrasonic signal to form a first discrete time sequence.
具体的,所述液体槽装置中的液体为一预设恒定温度的水。Specifically, the liquid in the liquid tank device is water at a preset constant temperature.
一种地震物理模型实验系统的模拟数据采集方法,应用于上述的地震物理模型实验系统,所述方法包括:A simulation data acquisition method of an earthquake physics model experiment system, applied to the above-mentioned earthquake physics model experiment system, said method comprising:
在水平方向上建立平面直角坐标系;Establish a plane Cartesian coordinate system in the horizontal direction;
接收一次实验的实验坐标信息,所述实验坐标信息包括在液体平面上的激发点坐标和接收点坐标;receiving experimental coordinate information of an experiment, the experimental coordinate information including excitation point coordinates and receiving point coordinates on the liquid plane;
控制激发换能器的超声波发射端处于所述激发点坐标,并控制接收换能器的超声波接收端处于所述接收点坐标;controlling the ultrasonic transmitting end of the excitation transducer to be at the coordinates of the excitation point, and controlling the ultrasonic receiving end of the receiving transducer to be at the coordinates of the receiving point;
控制所述激发换能器向地震物理模型发射第一超声波信号;controlling the excitation transducer to transmit a first ultrasonic signal to the seismic physics model;
控制接收换能器接收地震物理模型反射的第二超声波信号;controlling the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model;
以预先设置的采样间隔和采样点数从所述第二超声波信号中采集时间信息,形成第一离散时间序列;collecting time information from the second ultrasonic signal with a preset sampling interval and number of sampling points to form a first discrete time sequence;
对所述第一离散时间序列进行校正,生成第二离散时间序列;correcting the first discrete time series to generate a second discrete time series;
根据多次实验的所述第二离散时间序列生成地震物理模型实验模拟数据。The earthquake physical model experiment simulation data is generated according to the second discrete time series of multiple experiments.
进一步的,该地震物理模型实验系统的模拟数据采集方法,还包括:Further, the simulation data acquisition method of the earthquake physics model experiment system also includes:
确定所述激发点坐标处的地震物理模型表面的第一高程数据,并确定所述接收点坐标处的地震物理模型表面的第二高程数据;determining the first elevation data of the surface of the physical earthquake model at the coordinates of the excitation point, and determining the second elevation data of the surface of the physical earthquake model at the coordinates of the receiving point;
确定当前液体平面的第三高程数据。Determine the third elevation data of the current liquid level.
进一步的,该地震物理模型实验系统的模拟数据采集方法,还包括:Further, the simulation data acquisition method of the earthquake physics model experiment system also includes:
确定当前超声波信号在液体槽装置的液体中的传播速度。Determines the current propagation speed of the ultrasonic signal in the liquid of the liquid tank assembly.
此外,该地震物理模型实验系统的模拟数据采集方法,还包括:In addition, the simulation data acquisition method of the earthquake physics model experiment system also includes:
根据所述第一高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的激发换能器的超声波发射端到地震物理模型表面之间的液体层引起的第一超声波旅行时间差;According to the first elevation data, the third elevation data and the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device, determine the first wave caused by the liquid layer between the ultrasonic transmitting end of the current excitation transducer and the surface of the seismic physics model - ultrasonic travel time difference;
所述第一超声波旅行时间差为:其中,h0为所述第三高程数据;v0为超声波信号在液体槽装置的液体中的传播速度;(sxi,syi)为当前第i次实验的激发点坐标;为第i次实验的第一高程数据。The first ultrasonic travel time difference is: Wherein, h 0 is the third elevation data; v 0 is the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device; (sx i , sy i ) is the excitation point coordinates of the current i-th experiment; is the first elevation data of the i-th experiment.
此外,该地震物理模型实验系统的模拟数据采集方法,还包括:In addition, the simulation data acquisition method of the earthquake physics model experiment system also includes:
根据所述第二高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的接收换能器的超声波接收端到地震物理模型表面之间的液体层引起的第二超声波旅行时间差;According to the second elevation data, the third elevation data and the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device, determine the first wave caused by the liquid layer between the ultrasonic receiving end of the current receiving transducer and the surface of the seismic physics model Two ultrasound travel time difference;
所述第二超声波旅行时间差为:其中,(rxi,ryi)为当前第i次实验的接收点坐标;为第i次实验的所述第二高程数据。The second ultrasonic travel time difference is: Among them, (rx i , ry i ) is the receiving point coordinates of the current i-th experiment; is the second elevation data of the i-th experiment.
具体的,对所述第一离散时间序列进行校正,生成第二离散时间序列,包括:Specifically, correcting the first discrete time series to generate a second discrete time series, including:
根据所述第一超声波旅行时间差、第二超声波旅行时间差和所述采样间隔确定一校正参数;所述校正参数为:其中,dt为所述采样间隔;表示的值采用四舍五入方式取整数;A correction parameter is determined according to the first ultrasonic travel time difference, the second ultrasonic travel time difference and the sampling interval; the correction parameter is: Wherein, dt is described sampling interval; express The value of is rounded to an integer;
根据所述第一离散时间序列与所述校正参数,生成所述第二离散时间序列;所述第二离散时间序列为:Generate the second discrete time series according to the first discrete time series and the correction parameters; the second discrete time series is:
其中,si′(j)为第i次实验的第二离散时间序列的第j个采样点的离散时间值;si(j+k)为第i次实验的第一离散时间序列的第j+k个采样点的离散时间值;M为所述采样点数。Among them, s i ′(j) is the discrete time value of the jth sampling point of the second discrete time series of the i-th experiment; s i (j+k) is the first discrete time series of the i-th experiment Discrete time values of j+k sampling points; M is the number of sampling points.
一种地震物理模型实验系统的模拟数据采集装置,应用于上述的地震物理模型实验系统,所述装置包括:An analog data acquisition device of an earthquake physics model experiment system, applied to the above earthquake physics model experiment system, said device comprising:
坐标系建立单元,用于在水平方向上建立平面直角坐标系;The coordinate system establishment unit is used to establish a plane Cartesian coordinate system in the horizontal direction;
实验坐标信息接收单元,用于接收一次实验的实验坐标信息,所述实验坐标信息包括在液体平面上的激发点坐标和接收点坐标;The experimental coordinate information receiving unit is used to receive the experimental coordinate information of an experiment, and the experimental coordinate information includes the coordinates of the excitation point and the coordinates of the receiving point on the liquid plane;
换能器位置控制单元,用于控制激发换能器的超声波发射端处于所述激发点坐标,并控制接收换能器的超声波接收端处于所述接收点坐标;The transducer position control unit is used to control the ultrasonic transmitting end of the exciting transducer to be at the excitation point coordinates, and control the ultrasonic receiving end of the receiving transducer to be at the receiving point coordinates;
超声波信号发射控制单元,用于控制所述激发换能器向地震物理模型发射第一超声波信号;An ultrasonic signal transmitting control unit, configured to control the excitation transducer to transmit a first ultrasonic signal to the seismic physical model;
超声波信号接收控制单元,用于控制接收换能器接收地震物理模型反射的第二超声波信号;The ultrasonic signal receiving control unit is used to control the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model;
第一离散时间序列生成单元,用于以预先设置的采样间隔和采样点数从所述第二超声波信号中采集时间信息,形成第一离散时间序列;A first discrete time series generating unit, configured to collect time information from the second ultrasonic signal with a preset sampling interval and number of sampling points to form a first discrete time series;
校正单元,用于对所述第一离散时间序列进行校正,生成第二离散时间序列;a correction unit, configured to correct the first discrete time series to generate a second discrete time series;
地震物理模型实验模拟数据生成单元,用于根据多次实验的所述第二离散时间序列生成地震物理模型实验模拟数据。The earthquake physics model experiment simulation data generating unit is configured to generate the earthquake physics model experiment simulation data according to the second discrete time series of multiple experiments.
进一步的,该地震物理模型实验系统的模拟数据采集装置,还包括:Further, the analog data acquisition device of the earthquake physics model experiment system also includes:
高程数据确定单元,用于确定所述激发点坐标处的地震物理模型表面的第一高程数据,确定所述接收点坐标处的地震物理模型表面的第二高程数据,以及确定当前液体平面的第三高程数据。An elevation data determining unit, configured to determine the first elevation data of the surface of the seismic physics model at the coordinates of the excitation point, determine the second elevation data of the surface of the earthquake physics model at the coordinates of the receiving point, and determine the first elevation data of the current liquid level Three elevation data.
进一步的,该地震物理模型实验系统的模拟数据采集装置,还包括:Further, the analog data acquisition device of the earthquake physics model experiment system also includes:
传播速度确定单元,用于确定当前超声波信号在液体槽装置的液体中的传播速度。The propagation velocity determining unit is used for determining the propagation velocity of the current ultrasonic signal in the liquid of the liquid tank device.
进一步的,该地震物理模型实验系统的模拟数据采集装置,还包括:Further, the analog data acquisition device of the earthquake physics model experiment system also includes:
第一超声波旅行时间差确定单元,用于根据所述第一高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的激发换能器的超声波发射端到地震物理模型表面之间的液体层引起的第一超声波旅行时间差;The first ultrasonic travel time difference determination unit is used to determine the current ultrasonic transmitting end of the excitation transducer to the seismic physics according to the first elevation data, the third elevation data and the propagation speed of the ultrasonic signal in the liquid of the liquid tank device The first ultrasonic travel time difference caused by the liquid layer between the model surfaces;
所述第一超声波旅行时间差为:其中,h0为所述第三高程数据;v0为超声波信号在液体槽装置的液体中的传播速度;(sxi,syi)为当前第i次实验的激发点坐标;为第i次实验的第一高程数据。The first ultrasonic travel time difference is: Wherein, h 0 is the third elevation data; v 0 is the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device; (sx i , sy i ) is the excitation point coordinates of the current i-th experiment; is the first elevation data of the i-th experiment.
进一步的,该地震物理模型实验系统的模拟数据采集装置,还包括:Further, the analog data acquisition device of the earthquake physics model experiment system also includes:
第二超声波旅行时间差确定单元,用于根据所述第二高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的接收换能器的超声波接收端到地震物理模型表面之间的液体层引起的第二超声波旅行时间差;The second ultrasonic travel time difference determining unit is used to determine the current ultrasonic receiving end of the receiving transducer to the seismic physics The second ultrasonic travel time difference caused by the liquid layer between the model surfaces;
所述第二超声波旅行时间差为:其中,(rxi,ryi)为当前第i次实验的接收点坐标;为第i次实验的所述第二高程数据。The second ultrasonic travel time difference is: Among them, (rx i , ry i ) is the receiving point coordinates of the current i-th experiment; is the second elevation data of the i-th experiment.
此外,该校正单元,包括:In addition, the calibration unit includes:
校正参数确定模块,用于根据所述第一超声波旅行时间差、第二超声波旅行时间差和所述采样间隔确定一校正参数;所述校正参数为:其中,dt为所述采样间隔;表示的值采用四舍五入方式取整数;A correction parameter determination module, configured to determine a correction parameter according to the first ultrasonic travel time difference, the second ultrasonic travel time difference and the sampling interval; the correction parameter is: Wherein, dt is described sampling interval; express The value of is rounded to an integer;
第二离散时间序列生成模块,用于根据所述第一离散时间序列与所述校正参数,生成所述第二离散时间序列;所述第二离散时间序列为:The second discrete time series generation module is used to generate the second discrete time series according to the first discrete time series and the correction parameters; the second discrete time series is:
其中,si′(j)为第i次实验的第二离散时间序列的第j个采样点的离散时间值;si(j+k)为第i次实验的第一离散时间序列的第j+k个采样点的离散时间值;M为所述采样点数。Among them, s i ′(j) is the discrete time value of the jth sampling point of the second discrete time series of the i-th experiment; s i (j+k) is the first discrete time series of the i-th experiment Discrete time values of j+k sampling points; M is the number of sampling points.
本发明实施例提供的一种地震物理模型实验系统及其模拟数据采集方法和装置,所应用的激发换能器和接收换能器并非与地震物理模型采用硬接触的方式,而是通过中间设置的液体层来与地震物理模型间隔。这样,该地震物理模型实验系统最终形成的地震物理模型实验模拟数据信噪比较高,结果较为准确。而在地震物理模型实验系统的模拟数据采集方法中,通过控制激发换能器向地震物理模型发射第一超声波信号,并控制接收换能器接收地震物理模型反射的第二超声波信号;根据第二超声波信号采集的时间信息形成第一离散时间序列;进而还需要将述第一离散时间序列进行校正,以克服液体层对超声波传播的影响,生成第二离散时间序列;之后根据多次实验的第二离散时间序列即可生成地震物理模型实验模拟数据。可见,本发明能够避免当前传统的地震物理模拟实验超声波能量由于硬接触方式的接触不良而衰减,不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低,以及实验过程需要大量人工干预,地震物理模型数据采集时间长,采集效率低的问题。The embodiment of the present invention provides an earthquake physics model experiment system and its simulation data acquisition method and device. The liquid layer is separated from the seismic physics model. In this way, the simulated data of the earthquake physics model experiment finally formed by the earthquake physics model experiment system has a high signal-to-noise ratio and relatively accurate results. In the analog data acquisition method of the earthquake physics model experiment system, the excitation transducer is used to transmit the first ultrasonic signal to the earthquake physics model, and the receiving transducer is controlled to receive the second ultrasonic signal reflected by the earthquake physics model; according to the second The time information collected by the ultrasonic signal forms a first discrete time series; furthermore, the first discrete time series needs to be corrected to overcome the influence of the liquid layer on ultrasonic propagation, and generate a second discrete time series; then according to the first discrete time series of multiple experiments Two discrete time series can generate the experimental simulation data of the earthquake physical model. It can be seen that the present invention can avoid the attenuation of the ultrasonic energy in the current traditional earthquake physics simulation experiment due to the poor contact of the hard contact mode, which is not enough to detect the internal situation of the earthquake physics model, and the signal-to-noise ratio of the collected data is low, and the experimental process requires A large number of manual interventions, long data acquisition time for seismic physical models, and low acquisition efficiency.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only These are some embodiments of the present invention. For those skilled in the art, other drawings can also be obtained according to these drawings without any creative effort.
图1为现有技术中的传统的通过地震物理模型来进行地震物理模拟实验部署示意图;FIG. 1 is a schematic diagram of the deployment of a traditional earthquake physics simulation experiment using an earthquake physics model in the prior art;
图2为本发明实施例提供的一种地震物理模型实验系统的结构示意图;Fig. 2 is the structural representation of a kind of earthquake physics model experiment system provided by the embodiment of the present invention;
图3为本发明实施例提供的一种地震物理模型实验系统的模拟数据采集方法的流程图一;Fig. 3 is the flow chart one of the simulation data acquisition method of a kind of earthquake physics model experiment system provided by the embodiment of the present invention;
图4为本发明实施例提供的一种地震物理模型实验系统的模拟数据采集方法的流程图二;Fig. 4 is the flow chart two of the simulation data acquisition method of a kind of earthquake physics model experiment system provided by the embodiment of the present invention;
图5为本发明实施例中的地震物理模型实验结果数据和现有技术的实验结果数据的示意图;Fig. 5 is the schematic diagram of the experimental result data of earthquake physics model in the embodiment of the present invention and the experimental result data of prior art;
图6为本发明实施例提供的一种地震物理模型实验系统的模拟数据采集装置的结构示意图一;Fig. 6 is a structural schematic diagram 1 of an analog data acquisition device of an earthquake physics model experiment system provided by an embodiment of the present invention;
图7为本发明实施例提供的一种地震物理模型实验系统的模拟数据采集装置的结构示意图二。Fig. 7 is a second structural schematic diagram of an analog data acquisition device of an earthquake physics model experiment system provided by an embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.
如图2所示,本发明实施例提供一种地震物理模型实验系统20,包括:地震物理模型201、模拟采集控制系统202、激发换能器203、接收换能器204、液体槽装置205。该地震物理模型201的底面为平面,且该地震物理模型201放置于填充有液体206,且放置于水平面上的液体槽装置205中;液体槽装置205中的液体平面高于地震物理模型201;激发换能器203和接收换能器204与液体平面完全接触,且接触平面与液体平面高度相同;激发换能器203和接收换能器204均与模拟采集控制系统202连接。As shown in FIG. 2 , an embodiment of the present invention provides an earthquake physics model experiment system 20 , including: an earthquake physics model 201 , a simulation acquisition control system 202 , an exciting transducer 203 , a receiving transducer 204 , and a liquid tank device 205 . The bottom surface of the earthquake physics model 201 is a plane, and the earthquake physics model 201 is placed in a liquid tank device 205 filled with liquid 206 and placed on a horizontal plane; the liquid level in the liquid tank device 205 is higher than the earthquake physics model 201; The exciting transducer 203 and the receiving transducer 204 are in full contact with the liquid plane, and the contact plane is at the same height as the liquid plane; both the exciting transducer 203 and the receiving transducer 204 are connected to the analog acquisition control system 202 .
该激发换能器203,可以向地震物理模型201发射第一超声波信号。The excitation transducer 203 can transmit the first ultrasonic signal to the seismic physics model 201 .
该接收换能器204,可以接收地震物理模型201反射的第二超声波信号。The receiving transducer 204 can receive the second ultrasonic signal reflected by the seismic physics model 201 .
模拟采集控制系统202,可以控制激发换能器203的超声波发射端处于预先设置的激发点坐标,并控制接收换能器204的超声波接收端处于预先设置的接收点坐标,并从所述第二超声波信号中采集时间信息,形成第一离散时间序列。The analog acquisition control system 202 can control the ultrasonic transmitting end of the excitation transducer 203 to be at the preset excitation point coordinates, and control the ultrasonic receiving end of the receiving transducer 204 to be at the preset receiving point coordinates, and from the second Time information is collected from the ultrasonic signal to form a first discrete time sequence.
此处,激发换能器203和接收换能器204与液体平面完全接触,且接触平面与液体平面高度相同,是为了保证激发换能器203和接收换能器204均可以与液体平面良好耦合,且不侵入液体平面,利于高程数据的准确。Here, the exciting transducer 203 and the receiving transducer 204 are in full contact with the liquid plane, and the contact plane is at the same height as the liquid plane, in order to ensure that both the exciting transducer 203 and the receiving transducer 204 can be well coupled with the liquid plane , and does not invade the liquid plane, which is conducive to the accuracy of elevation data.
本发明实施例所涉及的地震物理模型201可以参照野外实际的地质体结构与实际测量得到的地表起伏高程数据,其与实际的地址体结构呈一定比例,如1:5000或1:10000等,但不仅局限于此。The earthquake physical model 201 involved in the embodiment of the present invention can refer to the actual geological body structure in the field and the surface relief elevation data obtained by actual measurement, which is in a certain ratio to the actual address body structure, such as 1:5000 or 1:10000, etc. But it's not limited to that.
值得说明的是,液体槽装置205中的液体206可以为一预设恒定温度的水。该预设恒定温度大于0℃,以保证水不结冰。另外,该液体206还可以采用粘稠度较低的其他液体,此处不再赘述。It should be noted that the liquid 206 in the liquid tank device 205 can be water at a preset constant temperature. The preset constant temperature is greater than 0° C. to ensure that the water does not freeze. In addition, the liquid 206 can also use other liquids with lower viscosity, which will not be repeated here.
本发明实施例提供的一种地震物理模型实验系统,其中所应用的激发换能器和接收换能器并非与地震物理模型采用硬接触的方式,而是通过中间设置的液体层来与地震物理模型间隔。这样,该地震物理模型实验系统最终形成的地震物理模型实验模拟数据信噪比较高,结果较为准确。避免当前传统的地震物理模拟实验超声波能量由于硬接触方式的接触不良而衰减,不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低,以及实验过程需要大量人工干预,地震物理模型数据采集时间长,采集效率低的问题。An embodiment of the present invention provides an experimental system for an earthquake physics model, wherein the applied excitation transducer and receiving transducer are not in hard contact with the earthquake physics model, but are in contact with the earthquake physics model through a liquid layer in the middle. Model interval. In this way, the simulated data of the earthquake physics model experiment finally formed by the earthquake physics model experiment system has a high signal-to-noise ratio and relatively accurate results. To avoid the attenuation of ultrasonic energy in the current traditional earthquake physics simulation experiment due to poor contact in the hard contact method, it is not enough to detect the internal conditions of the earthquake physics model, and the signal-to-noise ratio of the collected data is low, and the experimental process requires a lot of manual intervention. The collection time of physical model data is long and the collection efficiency is low.
此外,如图3所示,本发明实施例提供一种地震物理模型实验系统的模拟数据采集方法,应用于上述图2所示的地震物理模型实验系统,该方法包括:In addition, as shown in FIG. 3, an embodiment of the present invention provides a simulated data acquisition method for an earthquake physics model experiment system, which is applied to the earthquake physics model experiment system shown in FIG. 2 above. The method includes:
步骤301、在水平方向上建立平面直角坐标系。Step 301, establishing a plane Cartesian coordinate system in the horizontal direction.
步骤302、接收一次实验的实验坐标信息,实验坐标信息包括在液体平面上的激发点坐标和接收点坐标。Step 302: Receive experimental coordinate information of an experiment, where the experimental coordinate information includes the coordinates of the excitation point and the coordinates of the receiving point on the liquid plane.
步骤303、控制激发换能器的超声波发射端处于激发点坐标,并控制接收换能器的超声波接收端处于接收点坐标。Step 303 , controlling the ultrasonic transmitting end of the exciting transducer to be at the excitation point coordinates, and controlling the ultrasonic receiving end of the receiving transducer to be at the receiving point coordinates.
步骤304、控制激发换能器向地震物理模型发射第一超声波信号。Step 304, controlling the excitation transducer to transmit the first ultrasonic signal to the seismic physical model.
步骤305、控制接收换能器接收地震物理模型反射的第二超声波信号。Step 305: Control the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model.
步骤306、以预先设置的采样间隔和采样点数从第二超声波信号中采集时间信息,形成第一离散时间序列。Step 306: Collect time information from the second ultrasonic signal with a preset sampling interval and number of sampling points to form a first discrete time sequence.
步骤307、对第一离散时间序列进行校正,生成第二离散时间序列。Step 307. Correct the first discrete time series to generate a second discrete time series.
步骤308、根据多次实验的第二离散时间序列生成地震物理模型实验模拟数据。Step 308 , generating earthquake physical model experimental simulation data according to the second discrete time series of multiple experiments.
本发明实施例提供的一种地震物理模型实验系统的模拟数据采集方法,所应用的激发换能器和接收换能器并非与地震物理模型采用硬接触的方式,而是通过中间设置的液体层来与地震物理模型间隔。这样,通过控制激发换能器向地震物理模型发射第一超声波信号,并控制接收换能器接收地震物理模型反射的第二超声波信号;根据第二超声波信号采集的时间信息形成第一离散时间序列;进而还需要将述第一离散时间序列进行校正,以克服液体层对超声波传播的影响,生成第二离散时间序列;之后根据多次实验的第二离散时间序列即可生成地震物理模型实验模拟数据。可见,本发明能够避免当前传统的地震物理模拟实验超声波能量由于硬接触方式的接触不良而衰减,不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低,以及实验过程需要大量人工干预,地震物理模型数据采集时间长,采集效率低的问题。The embodiment of the present invention provides a simulated data acquisition method for an earthquake physics model experiment system. The applied excitation transducer and receiving transducer are not in hard contact with the earthquake physics model, but through a liquid layer set in the middle to be spaced from the earthquake physics model. In this way, by controlling the excitation transducer to transmit the first ultrasonic signal to the seismic physical model, and controlling the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model; the first discrete time series is formed according to the time information collected by the second ultrasonic signal ; and furthermore, the first discrete time series needs to be corrected to overcome the influence of the liquid layer on ultrasonic propagation, and generate a second discrete time series; then, the earthquake physical model experimental simulation can be generated according to the second discrete time series of multiple experiments data. It can be seen that the present invention can avoid the attenuation of the ultrasonic energy in the current traditional earthquake physics simulation experiment due to the poor contact of the hard contact mode, which is not enough to detect the internal situation of the earthquake physics model, and the signal-to-noise ratio of the collected data is low, and the experimental process requires A large number of manual interventions, long data acquisition time for seismic physical models, and low acquisition efficiency.
为了使本领域的技术人员更好的了解本发明,下面列举一个更为详细的实施例,如图4所示,本发明实施例提供一种地震物理模型实验系统的模拟数据采集方法,应用于上述图2所示的地震物理模型实验系统,该方法包括:In order to enable those skilled in the art to better understand the present invention, a more detailed embodiment is enumerated below, as shown in Figure 4, the embodiment of the present invention provides a kind of simulation data collection method of earthquake physics model experiment system, is applied to The earthquake physics model experiment system shown in above-mentioned Fig. 2, this method comprises:
步骤401、在水平方向上建立平面直角坐标系。Step 401, establishing a plane Cartesian coordinate system in the horizontal direction.
此处的平面直角坐标系的原点可以在液体槽装置中任意选取。The origin of the plane Cartesian coordinate system here can be arbitrarily selected in the liquid tank device.
步骤402、接收一次实验的实验坐标信息,实验坐标信息包括在液体平面上的激发点坐标和接收点坐标。Step 402: Receive experimental coordinate information of an experiment, where the experimental coordinate information includes the coordinates of the excitation point and the coordinates of the receiving point on the liquid plane.
具体的,每进行一次实验,既需要获取一次实验坐标信息,以使得后续通过实验坐标信息来控制激发换能器和接收换能器的位置。Specifically, every time an experiment is performed, it is necessary to obtain the experimental coordinate information once, so that the positions of the exciting transducer and the receiving transducer can be controlled subsequently by using the experimental coordinate information.
步骤403、确定激发点坐标处的地震物理模型表面的第一高程数据,并确定接收点坐标处的地震物理模型表面的第二高程数据,确定当前液体平面的第三高程数据。Step 403: Determine the first elevation data of the surface of the seismic physics model at the coordinates of the excitation point, determine the second elevation data of the surface of the seismic physics model at the coordinates of the receiving point, and determine the third elevation data of the current liquid level.
步骤404、确定当前超声波信号在液体槽装置的液体中的传播速度。Step 404. Determine the propagation velocity of the current ultrasonic signal in the liquid in the liquid tank device.
步骤405、根据第一高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的激发换能器的超声波发射端到地震物理模型表面之间的液体层引起的第一超声波旅行时间差。Step 405, according to the first elevation data, the third elevation data and the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device, determine the current wave caused by the liquid layer between the ultrasonic transmitting end of the excitation transducer and the surface of the seismic physics model First ultrasound travel time difference.
此处,第一超声波旅行时间差为:其中,h0为第三高程数据;v0为超声波信号在液体槽装置的液体中的传播速度;(sxi,syi)为当前第i次实验的激发点坐标;为第i次实验的第一高程数据。Here, the first ultrasonic travel time difference is: Among them, h 0 is the third elevation data; v 0 is the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device; (sx i , sy i ) is the excitation point coordinates of the current i-th experiment; is the first elevation data of the i-th experiment.
步骤406、根据第二高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的接收换能器的超声波接收端到地震物理模型表面之间的液体层引起的第二超声波旅行时间差。Step 406, according to the second elevation data, the third elevation data and the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device, determine the current wave caused by the liquid layer between the ultrasonic receiving end of the receiving transducer and the surface of the seismic physics model Second ultrasound travel time difference.
此处,第二超声波旅行时间差为:其中,(rxi,ryi)为当前第i次实验的接收点坐标;为第i次实验的第二高程数据。Here, the second ultrasonic travel time difference is: Among them, (rx i , ry i ) is the receiving point coordinates of the current i-th experiment; is the second elevation data of the i-th experiment.
步骤407、控制激发换能器的超声波发射端处于激发点坐标,并控制接收换能器的超声波接收端处于接收点坐标。Step 407 , controlling the ultrasonic transmitting end of the exciting transducer to be at the excitation point coordinates, and controlling the ultrasonic receiving end of the receiving transducer to be at the receiving point coordinates.
步骤408、控制激发换能器向地震物理模型发射第一超声波信号。Step 408, controlling the excitation transducer to transmit the first ultrasonic signal to the seismic physical model.
步骤409、控制接收换能器接收地震物理模型反射的第二超声波信号。Step 409, controlling the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model.
步骤410、以预先设置的采样间隔和采样点数从第二超声波信号中采集时间信息,形成第一离散时间序列。Step 410: Collect time information from the second ultrasonic signal with a preset sampling interval and number of sampling points to form a first discrete time series.
此处的第一离散时间序列可以表示为:{Si(j),j=1,2,...,M;i=1,2,...,N};其中,i表示第i次实验,j表示第j次采样,M为上述的采样点数,N为实验总次数。The first discrete time series here can be expressed as: {S i (j),j=1,2,...,M; i=1,2,...,N}; where i represents the i-th is the number of experiments, j represents the jth sampling, M is the number of sampling points mentioned above, and N is the total number of experiments.
步骤411、根据第一超声波旅行时间差、第二超声波旅行时间差和采样间隔确定一校正参数。Step 411. Determine a correction parameter according to the first ultrasonic travel time difference, the second ultrasonic travel time difference and the sampling interval.
其中,该校正参数为:其中,dt为采样间隔;表示的值采用四舍五入方式取整数。Among them, the correction parameters are: Among them, dt is the sampling interval; express The value of is rounded to an integer.
步骤412、根据第一离散时间序列与校正参数,生成第二离散时间序列。Step 412: Generate a second discrete time series according to the first discrete time series and the correction parameters.
其中,该第二离散时间序列为:Wherein, the second discrete time series is:
其中,si′(j)为第i次实验的第二离散时间序列的第j个采样点的离散时间值;si(j+k)为第i次实验的第一离散时间序列的第j+k个采样点的离散时间值;M为采样点数。Among them, s i ′(j) is the discrete time value of the jth sampling point of the second discrete time series of the i-th experiment; s i (j+k) is the first discrete time series of the i-th experiment Discrete time values of j+k sampling points; M is the number of sampling points.
在步骤412之后返回执行步骤402以进行下一次实验,若实验次数满足预先设置的次数阈值,则继续执行步骤413。After step 412, return to step 402 for the next experiment, and if the number of experiments meets the preset number threshold, continue to execute step 413.
步骤413、根据多次实验的第二离散时间序列生成地震物理模型实验模拟数据。Step 413 , generating earthquake physical model experimental simulation data according to the second discrete time series of multiple experiments.
为了说明本发明实施例提供的一种地震物理模型实验系统的模拟数据采集方法的效果,可以参见图5所示,其中图5的(a)部分为现有技术中的传统的地震物理模拟实验所采集得到的数据示意图;图5的(b)部分为采用本发明实施例提供的地震物理模型实验系统的模拟数据采集方法得到的数据示意图;图5的(c)部分为采用本发明实施例提供的地震物理模型实验系统的模拟数据采集方法,但不进行液体层影响校正过程,得到的最终的数据示意图。图5中纵坐标为时间,横坐标为地震数据道号。比较图5中椭圆形框51和矩形框52所指示的相同位置处的地震波场,可以看出:地震物理模型内部结构所对应的地震反射信号可以通过本发明提供的方法获得。将图5中椭圆形框51中的信息进行比较可知,图5的(b)部分得到的地震记录中的地震信号连续性更强;而由图5中矩形框52中信息的比较可知,图5的(b)部分得到的地震记录垂向分辨率更强。图5中的箭头指出了图5的(a)部分和(b)部分两种采集方法得到的不同地震信息,其中,图5中的一号箭头53指出了传统模拟采集得到的地震记录中的直达波干扰。图5中二号箭头54指出了本发明的模拟采集得到的地震记录中由于水层的存在而发生的多次波,在分析数据中需要注意。图5中三号箭头55指示处是未经后期校正处理的地震波场信号。可见,未经后期校正处理的地震波场信号随着道号而偏移。In order to illustrate the effect of the simulation data acquisition method of a kind of earthquake physics model experiment system that the embodiment of the present invention provides, can refer to as shown in Fig. 5, wherein the part (a) of Fig. 5 is the traditional earthquake physics simulation experiment in the prior art The data schematic diagram that is gathered; The (b) part of Fig. 5 is the data schematic diagram that adopts the simulated data acquisition method of the earthquake physics model experiment system that the embodiment of the present invention provides; The (c) part of Fig. 5 is that adopts the embodiment of the present invention The simulation data acquisition method of the earthquake physics model experiment system is provided, but the liquid layer influence correction process is not performed, and the final data schematic diagram is obtained. In Fig. 5, the ordinate is time, and the abscissa is the seismic data track number. Comparing the seismic wavefields at the same positions indicated by the oval box 51 and the rectangular box 52 in Fig. 5, it can be seen that the seismic reflection signal corresponding to the internal structure of the seismic physical model can be obtained by the method provided by the present invention. Comparing the information in the oval frame 51 in Fig. 5 shows that the seismic signal continuity in the seismic record obtained in part (b) of Fig. 5 is stronger; and by comparing the information in the rectangular frame 52 in Fig. 5, it can be seen that The vertical resolution of seismic records obtained in part (b) of 5 is stronger. The arrows in Fig. 5 have pointed out the different seismic information obtained by the two acquisition methods of part (a) and part (b) of Fig. 5, wherein, arrow 53 in Fig. 5 has pointed out the difference in the seismic records obtained by traditional analog acquisition Direct wave interference. The second arrow 54 in Fig. 5 points out the multiple waves generated by the existence of the water layer in the seismic records obtained by the simulated acquisition of the present invention, which should be paid attention to when analyzing the data. The position indicated by the third arrow 55 in Fig. 5 is the seismic wave field signal without post-correction processing. It can be seen that the seismic wave field signal without post-correction processing shifts with the track number.
本发明实施例提供的一种地震物理模型实验系统的模拟数据采集方法,所应用的激发换能器和接收换能器并非与地震物理模型采用硬接触的方式,而是通过中间设置的液体层来与地震物理模型间隔。这样,通过控制激发换能器向地震物理模型发射第一超声波信号,并控制接收换能器接收地震物理模型反射的第二超声波信号;根据第二超声波信号采集的时间信息形成第一离散时间序列;进而还需要将述第一离散时间序列进行校正,以克服液体层对超声波传播的影响,生成第二离散时间序列;之后根据多次实验的第二离散时间序列即可生成地震物理模型实验模拟数据。可见,本发明能够避免当前传统的地震物理模拟实验超声波能量由于硬接触方式的接触不良而衰减,不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低,以及实验过程需要大量人工干预,地震物理模型数据采集时间长,采集效率低的问题。The embodiment of the present invention provides a simulated data acquisition method for an earthquake physics model experiment system. The applied excitation transducer and receiving transducer are not in hard contact with the earthquake physics model, but through a liquid layer set in the middle to be spaced from the earthquake physics model. In this way, by controlling the excitation transducer to transmit the first ultrasonic signal to the seismic physical model, and controlling the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model; the first discrete time series is formed according to the time information collected by the second ultrasonic signal ; and furthermore, the first discrete time series needs to be corrected to overcome the influence of the liquid layer on ultrasonic propagation, and generate a second discrete time series; then, the earthquake physical model experimental simulation can be generated according to the second discrete time series of multiple experiments data. It can be seen that the present invention can avoid the attenuation of the ultrasonic energy in the current traditional earthquake physics simulation experiment due to the poor contact of the hard contact mode, which is not enough to detect the internal situation of the earthquake physics model, and the signal-to-noise ratio of the collected data is low, and the experimental process requires A large number of manual interventions, long data collection time for seismic physical models, and low collection efficiency.
对应于上述图3、图4所示的方法实施例,如图6所示,本发明实施例提供一种地震物理模型实验系统的模拟数据采集装置,应用于上述的地震物理模型实验系统,该装置包括:Corresponding to the above method embodiments shown in Fig. 3 and Fig. 4, as shown in Fig. 6, an embodiment of the present invention provides an analog data acquisition device for an earthquake physics model experiment system, which is applied to the above earthquake physics model experiment system, the Devices include:
坐标系建立单元61,用于在水平方向上建立平面直角坐标系。The coordinate system establishing unit 61 is configured to establish a plane Cartesian coordinate system in the horizontal direction.
实验坐标信息接收单元62,用于接收一次实验的实验坐标信息,实验坐标信息包括在液体平面上的激发点坐标和接收点坐标。The experimental coordinate information receiving unit 62 is configured to receive the experimental coordinate information of an experiment, the experimental coordinate information includes the excitation point coordinates and the receiving point coordinates on the liquid plane.
换能器位置控制单元63,用于控制激发换能器的超声波发射端处于激发点坐标,并控制接收换能器的超声波接收端处于接收点坐标。The transducer position control unit 63 is used to control the ultrasonic transmitting end of the exciting transducer to be at the excitation point coordinates, and control the ultrasonic receiving end of the receiving transducer to be at the receiving point coordinates.
超声波信号发射控制单元64,用于控制激发换能器向地震物理模型发射第一超声波信号。The ultrasonic signal transmission control unit 64 is configured to control the excitation transducer to transmit the first ultrasonic signal to the seismic physical model.
超声波信号接收控制单元65,用于控制接收换能器接收地震物理模型反射的第二超声波信号。The ultrasonic signal receiving control unit 65 is configured to control the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model.
第一离散时间序列生成单元66,用于以预先设置的采样间隔和采样点数从第二超声波信号中采集时间信息,形成第一离散时间序列。The first discrete time sequence generating unit 66 is configured to collect time information from the second ultrasonic signal at a preset sampling interval and number of sampling points to form a first discrete time sequence.
校正单元67,用于对第一离散时间序列进行校正,生成第二离散时间序列。A correction unit 67, configured to correct the first discrete time series to generate a second discrete time series.
地震物理模型实验模拟数据生成单元68,用于根据多次实验的第二离散时间序列生成地震物理模型实验模拟数据。The earthquake physics model experiment simulation data generation unit 68 is configured to generate the earthquake physics model experiment simulation data according to the second discrete time series of multiple experiments.
进一步的,如图7所示,该地震物理模型实验系统的模拟数据采集装置,还可以包括:Further, as shown in Figure 7, the analog data acquisition device of the earthquake physics model experiment system may also include:
高程数据确定单元69,可以确定激发点坐标处的地震物理模型表面的第一高程数据,确定接收点坐标处的地震物理模型表面的第二高程数据,以及确定当前液体平面的第三高程数据。The elevation data determination unit 69 can determine the first elevation data of the surface of the seismic physical model at the coordinates of the excitation point, the second elevation data of the surface of the seismic physical model at the coordinates of the receiving point, and the third elevation data of the current liquid level.
进一步的,如图7所示,该地震物理模型实验系统的模拟数据采集装置,还可以包括:Further, as shown in Figure 7, the analog data acquisition device of the earthquake physics model experiment system may also include:
传播速度确定单元70,可以确定当前超声波信号在液体槽装置的液体中的传播速度。The propagation velocity determining unit 70 can determine the propagation velocity of the current ultrasonic signal in the liquid of the liquid tank device.
进一步的,如图7所示,该地震物理模型实验系统的模拟数据采集装置,还包括:Further, as shown in Figure 7, the analog data acquisition device of the earthquake physics model experiment system also includes:
第一超声波旅行时间差确定单元71,可以根据第一高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的激发换能器的超声波发射端到地震物理模型表面之间的液体层引起的第一超声波旅行时间差。The first ultrasonic travel time difference determination unit 71 can determine the current ultrasonic transmitting end of the excitation transducer to the surface of the seismic physics model according to the first elevation data, the third elevation data and the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device The difference in travel time of the first ultrasonic wave is caused by the liquid layer between.
第一超声波旅行时间差为:其中,h0为第三高程数据;v0为超声波信号在液体槽装置的液体中的传播速度;(sxi,syi)为当前第i次实验的激发点坐标;为第i次实验的第一高程数据。The first ultrasonic travel time difference is: Among them, h 0 is the third elevation data; v 0 is the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device; (sx i , sy i ) is the excitation point coordinates of the current i-th experiment; is the first elevation data of the i-th experiment.
第二超声波旅行时间差确定单元72,可以根据第二高程数据、第三高程数据和超声波信号在液体槽装置的液体中的传播速度,确定当前的接收换能器的超声波接收端到地震物理模型表面之间的液体层引起的第二超声波旅行时间差。The second ultrasonic travel time difference determination unit 72 can determine the current ultrasonic receiving end of the receiving transducer to the surface of the seismic physics model according to the second elevation data, the third elevation data and the propagation velocity of the ultrasonic signal in the liquid of the liquid tank device The second ultrasonic travel time difference is caused by the liquid layer in between.
第二超声波旅行时间差为:其中,(rxi,ryi)为当前第i次实验的接收点坐标;为第i次实验的第二高程数据。The second ultrasonic travel time difference is: Among them, (rx i , ry i ) is the receiving point coordinates of the current i-th experiment; is the second elevation data of the i-th experiment.
此外,如图7所示,该校正单元67,可以包括:In addition, as shown in FIG. 7, the correction unit 67 may include:
校正参数确定模块671,用于根据第一超声波旅行时间差、第二超声波旅行时间差和采样间隔确定一校正参数;校正参数为:其中,dt为采样间隔;表示的值采用四舍五入方式取整数。The correction parameter determination module 671 is used to determine a correction parameter according to the first ultrasonic travel time difference, the second ultrasonic travel time difference and the sampling interval; the correction parameter is: Among them, dt is the sampling interval; express The value of is rounded to an integer.
第二离散时间序列生成模块672,用于根据第一离散时间序列与校正参数,生成第二离散时间序列;第二离散时间序列为:The second discrete time series generation module 672 is used to generate a second discrete time series according to the first discrete time series and correction parameters; the second discrete time series is:
其中,si′(j)为第i次实验的第二离散时间序列的第j个采样点的离散时间值;si(j+k)为第i次实验的第一离散时间序列的第j+k个采样点的离散时间值;M为采样点数。Among them, s i ′(j) is the discrete time value of the jth sampling point of the second discrete time series of the i-th experiment; s i (j+k) is the first discrete time series of the i-th experiment Discrete time values of j+k sampling points; M is the number of sampling points.
值得说明的是,本发明实施例提供的一种地震物理模型实验系统的模拟数据采集装置的具体实现方式可以参见上述方法实施例,此处不再赘述。It is worth noting that, for the specific implementation manner of the analog data acquisition device of the earthquake physics model experiment system provided by the embodiment of the present invention, reference may be made to the above-mentioned method embodiment, which will not be repeated here.
本发明实施例提供的一种地震物理模型实验系统的模拟数据采集装置,所应用的激发换能器和接收换能器并非与地震物理模型采用硬接触的方式,而是通过中间设置的液体层来与地震物理模型间隔。这样,通过控制激发换能器向地震物理模型发射第一超声波信号,并控制接收换能器接收地震物理模型反射的第二超声波信号;根据第二超声波信号采集的时间信息形成第一离散时间序列;进而还需要将述第一离散时间序列进行校正,以克服液体层对超声波传播的影响,生成第二离散时间序列;之后根据多次实验的第二离散时间序列即可生成地震物理模型实验模拟数据。可见,本发明能够避免当前传统的地震物理模拟实验超声波能量由于硬接触方式的接触不良而衰减,不足以探测地震物理模型的内部情况,且所采集的数据信噪比较低,以及实验过程需要大量人工干预,地震物理模型数据采集时间长,采集效率低的问题。The embodiment of the present invention provides an analog data acquisition device for an earthquake physics model experiment system. The applied excitation transducer and receiving transducer are not in hard contact with the earthquake physics model, but through a liquid layer arranged in the middle to be spaced from the earthquake physics model. In this way, by controlling the excitation transducer to transmit the first ultrasonic signal to the seismic physical model, and controlling the receiving transducer to receive the second ultrasonic signal reflected by the seismic physical model; the first discrete time series is formed according to the time information collected by the second ultrasonic signal ; and furthermore, the first discrete time series needs to be corrected to overcome the influence of the liquid layer on ultrasonic propagation, and generate a second discrete time series; then, the earthquake physical model experimental simulation can be generated according to the second discrete time series of multiple experiments data. It can be seen that the present invention can avoid the attenuation of the ultrasonic energy in the current traditional earthquake physics simulation experiment due to the poor contact of the hard contact mode, which is not enough to detect the internal situation of the earthquake physics model, and the signal-to-noise ratio of the collected data is low, and the experimental process requires A large number of manual interventions, long data acquisition time for seismic physical models, and low acquisition efficiency.
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art should understand that the embodiments of the present invention may be provided as methods, systems, or computer program products. Accordingly, the present invention can take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It should be understood that each procedure and/or block in the flowchart and/or block diagram, and a combination of procedures and/or blocks in the flowchart and/or block diagram can be realized by computer program instructions. These computer program instructions may be provided to a general purpose computer, special purpose computer, embedded processor, or processor of other programmable data processing equipment to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing equipment produce a An apparatus for realizing the functions specified in one or more procedures of the flowchart and/or one or more blocks of the block diagram.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing apparatus to operate in a specific manner, such that the instructions stored in the computer-readable memory produce an article of manufacture comprising instruction means, the instructions The device realizes the function specified in one or more procedures of the flowchart and/or one or more blocks of the block diagram.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded onto a computer or other programmable data processing device, causing a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process, thereby The instructions provide steps for implementing the functions specified in the flow chart or blocks of the flowchart and/or the block or blocks of the block diagrams.
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。In the present invention, specific examples have been applied to explain the principles and implementation methods of the present invention, and the descriptions of the above examples are only used to help understand the method of the present invention and its core idea; meanwhile, for those of ordinary skill in the art, according to this The idea of the invention will have changes in the specific implementation and scope of application. To sum up, the contents of this specification should not be construed as limiting the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610187426.1A CN105785435B (en) | 2016-03-29 | 2016-03-29 | Seismic physical model experimental system and its analog data acquisition method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610187426.1A CN105785435B (en) | 2016-03-29 | 2016-03-29 | Seismic physical model experimental system and its analog data acquisition method and apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105785435A CN105785435A (en) | 2016-07-20 |
CN105785435B true CN105785435B (en) | 2018-08-03 |
Family
ID=56392298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610187426.1A Active CN105785435B (en) | 2016-03-29 | 2016-03-29 | Seismic physical model experimental system and its analog data acquisition method and apparatus |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105785435B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108061752A (en) * | 2016-11-07 | 2018-05-22 | 中国石油化工股份有限公司 | A kind of ultrasound test system of curved surface geological model and ultrasonic acquisition method |
CN106772581A (en) * | 2016-12-13 | 2017-05-31 | 西京学院 | A kind of three-dimensional relief surface physical analogy acquisition method based on reconfiguration technique |
CN112014886B (en) * | 2020-09-08 | 2023-10-27 | 中煤科工集团西安研究院有限公司 | Excitation source array hybrid acquisition method and system for earthquake physical simulation |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1299064A (en) * | 2000-11-14 | 2001-06-13 | 石油大学(北京) | Multi-channel collecting method and device for physical analog of earthquake |
US6681887B1 (en) * | 1998-12-18 | 2004-01-27 | Westerngeco | Method for the determination of local wave heights and acoustic sensor in marine seismic signals |
CN101246640A (en) * | 2007-02-15 | 2008-08-20 | 中国石油化工股份有限公司 | Analog ultrasonic wave earthquake signal physical excitation and receiving system and method thereof |
CN103558632A (en) * | 2013-11-20 | 2014-02-05 | 中国石油天然气集团公司 | Efficient method and device for acquiring single-trace data of seismic physical model |
CN203606507U (en) * | 2013-11-19 | 2014-05-21 | 中国石油天然气集团公司 | Multi-trace seismic physical model data acquisition device |
CN105277967A (en) * | 2014-07-22 | 2016-01-27 | 中国石油化工股份有限公司 | Water tank physical model ultrasonic automatic detection system and method |
-
2016
- 2016-03-29 CN CN201610187426.1A patent/CN105785435B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6681887B1 (en) * | 1998-12-18 | 2004-01-27 | Westerngeco | Method for the determination of local wave heights and acoustic sensor in marine seismic signals |
CN1299064A (en) * | 2000-11-14 | 2001-06-13 | 石油大学(北京) | Multi-channel collecting method and device for physical analog of earthquake |
CN101246640A (en) * | 2007-02-15 | 2008-08-20 | 中国石油化工股份有限公司 | Analog ultrasonic wave earthquake signal physical excitation and receiving system and method thereof |
CN203606507U (en) * | 2013-11-19 | 2014-05-21 | 中国石油天然气集团公司 | Multi-trace seismic physical model data acquisition device |
CN103558632A (en) * | 2013-11-20 | 2014-02-05 | 中国石油天然气集团公司 | Efficient method and device for acquiring single-trace data of seismic physical model |
CN105277967A (en) * | 2014-07-22 | 2016-01-27 | 中国石油化工股份有限公司 | Water tank physical model ultrasonic automatic detection system and method |
Non-Patent Citations (4)
Title |
---|
三维地震物理模型的研究;魏建新 等;《石油地球物理勘探》;20021231;第37卷(第6期);第556-561页 * |
地震物理模型技术的应用与发展;郝守玲 等;《勘探地球物理进展》;20020430;第25卷(第2期);第34-42页 * |
基于物理模型的煤田地震属性响应特征的关键技术研究;戴世鑫;《中国博士学位论文全文数据库 基础科学辑》;20130515(第05期);第17页第1段、第18页第2.6节、第21-24页第2.7节 * |
孔洞物理模型数据的地震响应特征分析;季敏 等;《石油地球物理勘探》;20090430;第44卷(第2期);第196-200页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105785435A (en) | 2016-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100545679C (en) | The earthquake poststack forward modeling method of relief surface | |
MX2011010913A (en) | Interferometric seismic data processing. | |
US7460437B2 (en) | Seismic data processing method and system for migration of seismic signals incorporating azimuthal variations in the velocity | |
CN102540250A (en) | Azimuth fidelity angle domain imaging-based fractured oil and gas reservoir seismic exploration method | |
CN103323821B (en) | Water surface target radiation noise simulation method | |
BR102015021137B1 (en) | PROCESS FOR SURVEYING A REAL UNDERGROUND FORMATION, COMPUTER SYSTEM TO GENERATE AN IMAGE OF A REAL UNDERGROUND FORMATION AND NON TRANSIENT COMPUTER READIBLE | |
CN116520431B (en) | Method, device and medium for constructing broadband layered sound velocity structure of shallow seabed sediments | |
CN110058309A (en) | A three-dimensional time-lapse seismic observation system and method | |
CN105785435B (en) | Seismic physical model experimental system and its analog data acquisition method and apparatus | |
CN106324665A (en) | Method and system of inverting fracture density | |
EA037479B1 (en) | Diving wave illumination using migration gathers | |
CN111443384B (en) | Determination of node depth and water column velocity | |
BR102013004900A2 (en) | METHODS AND SYSTEMS FOR DEMYSTICATING MARINE SEISMIC WAVE FIELDS USING FUNCTIONAL COST MINIMIZATION | |
CN104199088A (en) | Incident angle gather extraction method and system | |
CN109709602A (en) | Remote detection sound wave offset imaging method, device and system | |
US10317543B2 (en) | Estimation of a far field signature in a second direction from a far field signature in a first direction | |
CN109738944B (en) | Wide-angle reflection-based seismic acquisition parameter determination method and device | |
US20220026593A1 (en) | Implicit property modeling | |
CN102967884B (en) | Reliability of wave impedance inversion data evaluation method and device | |
CN115421195B (en) | Method, device, equipment and storage medium for generating velocity field in seismic exploration | |
CN110579798A (en) | Method and system for seismic acquisition and observation with equal reflection angle intervals | |
US9829593B2 (en) | Determination of an impulse response at a subsurface image level | |
Cheng et al. | Locating leaking buried pipes based on ground microseismic records in 3D space | |
CN115755177A (en) | Broadband seismic constraint inversion imaging method for seabed sulfide ore body structure | |
CN108802822B (en) | The direct prestack time migration method of guarantor's width and device in direction anisotropy medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |