首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 328 毫秒
1.
区间B样条小波有限元GPR模拟双相随机混凝土介质   总被引:2,自引:2,他引:0       下载免费PDF全文
冯德山  王珣 《地球物理学报》2016,59(8):3098-3109
基于可分离小波理论,由一维区间B样条小波尺度函数的张量积构造二维B样条小波基,并将它作为GPR波动方程求解的插值函数,通过引入转换矩阵,实现小波系数空间与雷达电磁场之间的转换.应用Galerkin算法,推导了二维区间B样条小波有限元GPR波动方程离散格式,求出了2阶1尺度与2阶2尺度BSWI尺度函数的积分值及联系系数,给出了该算法的详细求解过程.编制了BSWI的Matlab模拟程序,应用该程序对两个典型实例进行了正演,结果表明:BSWI能采用较少的单元达到与FEM相似的精度,而BSWI算法尺度提升能提高解的精度,但耗时会急剧增加.最后,将BSWI算法应用于双相随机混凝土模型,说明随机介质模型理论能灵活、有效地描述实际混凝土介质的分布,正演剖面与实测剖面特征更相符,能更真实地模拟雷达波的传播过程,可为提高GPR的探测效果和解释准确性提供理论基础.  相似文献   

2.
Based on the measurement of the arrival time of maxima magnitude from band-pass filtering signals which were determined using a new Morlet wavelet multiple-filter method,we develop a method for measuring intrinsic and attenuative dispersion of the first cycle direct P-wave.We determine relative group delays of spectral components of direct P-waves for 984 ray paths from SML and ALS stations of the Taiwan Central Weather Bureau Seismic Network(CWBSN).Using continuous relaxation model,we deduce a new transfer...  相似文献   

3.
The gravitational potential of the Earth is usually modeled by means of a series expansion in terms of spherical harmonics. However, the computation of the series coefficients requires preferably homogeneous distributed global data sets. Since one of the most important features of wavelet functions is the ability to localize both in the spatial and in the frequency domain, regional and local structures may be modeled by means of a spherical wavelet expansion. In general, applying wavelet theory a given input data set is decomposed into a certain number of frequency-dependent detail signals, which can be interpreted as the building blocks of a multi-resolution representation. On the other hand, there is no doubt that the low-frequency part of the geopotential can be modeled appropriately by means of spherical harmonics. Hence, the main idea of this paper is to derive a combined model consisting of an expansion in spherical harmonics for the low-frequency part and an expansion in spherical wavelets for the remaining medium and high-frequency parts of the gravity field. Furthermore, an appropriate parameter estimation procedure is outlined to solve for the unknown model coefficients.  相似文献   

4.
小波尺度域含气储层地震波衰减特征   总被引:22,自引:4,他引:18       下载免费PDF全文
黏弹性衰减因子Q的可靠估计可通过Q反褶积来提高地震资料的分辨率并有助于振幅分析. 本文从小波理论出发,结合地震波在黏弹性介质中的传播方程,推导出小波尺度域地震波能量衰减公式. 能量衰减公式具有下列性质:(1)Q值越大,能量衰减得越慢;Q值越小,能量衰减越严重;(2)尺度越小,信号中保留的能量越少;(3)对于脉冲源来说在理想的无衰减介质(即Q趋近于∞)中传播时,信号在不同尺度内的能量相同. 利用尺度能量公式,可从反射地震资料中直接估计品质因子Q(即衰减因子),也可以提取不同尺度的能量衰减剖面作为储层描述的属性参数,用来进行岩性识别和指示气藏,与经典的谱比法相比,避免了谱比法所面临的双时窗问题以及进行谱估计的窗选择问题. 理论模型试验表明了本文方法的正确性和有效性.  相似文献   

5.
Rainfall–runoff processes appear to be highly nonlinear in Bayinbluk watersheds of the northwestern China. In this study, the time‐scale wavelet transform has been used for the analysis of this nonstationary system. The Haar and Morlet wavelet transform were used to analyse the rainfall–runoff conversion relationship. Wavelet power spectrum and change point methods are also employed to analyse rainfall rates and runoffs measured at daily to half‐hourly sampling rate. The four experimental sites (Luoto, Haer, Kuce and Shengl) are located in the Tianshan Mountains (Xinjiang province, China). Correlation analysis and wavelet transform are first applied to runoff process in different underlying surfaces. Wavelet analyses of rainfall rates and runoffs also give meaningful information on the temporal variability of the rainfall–runoff relationship. Change point and wavelet power spectrum analysis provide simple interpretation of energy distribution between different scales. The results indicate that wavelet transform is a good method for analysing the nonlinear relationship of temporal–spatial responses between rainfall and runoff. This method allowed quantification of the processes affecting runoff and provided an insight into their implications in surface water management. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

6.
消除探地雷达数据的子波衰减和频散的反滤波方法   总被引:1,自引:1,他引:0       下载免费PDF全文
消除探地雷达数据的子波衰减和频散可以很好地提高探地雷达的勘探深度和勘探分辨率.常用的消除探地雷达数据的子波衰减和频散方法为反Q滤波方法.该方法需要利用地下介质的Q参数,但是正确求取地下介质的Q参数很困难.针对这一问题,本文提出了一种消除探地雷达数据的子波衰减和频散的反滤波方法.该方法以地下介质反射系数是随机数为前提,利用地下介质等效滤波器具有最小相位这个特性,通过求取等效滤波器的振幅谱来求取等效滤波器的反滤波器.最后,利用该反滤波器对探地雷达数据进行反滤波,实现消除探地雷达数据的子波衰减和频散.  相似文献   

7.
通过阐述地面沉陷机理,分析地下岩体内应力状态和地下开挖空间的形成过程.引入数学积分模型建立地面沉陷理论模型,重点说明单元下沉盆地和特征点下沉,推导出地面沉陷预报公式.并利用有理分式解法求解地面沉陷预报数值.此法的计算精度高,适用面广,有效地解决了常规解法工作量大、计算精度低的缺点.为一种有效、实用的计算地面沉陷预报值的方法.  相似文献   

8.
本文将小波有限元法引入到流体饱和多孔隙介质二维波动方程的正演模拟中,以二维Daubechies小波的尺度函数代替多项式函数作为插值函数,构造二维张量积小波单元.引入一类特征函数解决了Daubechies小波没有显式解析表达式所带来的基函数积分值计算问题,并推导出计算分数节点上Daubechies小波函数值的递推公式,从而构造出由小波系数空间到波场位移空间的快速小波变换.数值模拟结果表明该方法是有效的.  相似文献   

9.
不依赖源子波的跨孔雷达时间域波形反演   总被引:1,自引:0,他引:1       下载免费PDF全文
刘四新  孟旭  傅磊 《地球物理学报》2016,59(12):4473-4482
波形反演是近年来较热门的反演方法,其分辨率可以达到亚波长级别.在波形反演的实际应用中,源子波的估计十分重要.传统方法使用反褶积来估计源子波并随着反演过程更新,该方法在合成数据波形反演中效果较好,但在实际数据反演过程中存在一系列的问题.由于实际数据信噪比较低,在源子波估计过程中需要大量的人为干涉,且结果并不一定可靠.本文使用一种基于褶积波场的新型目标函数,令反演过程不再依赖源子波.详细推导了针对跨孔雷达波形反演的梯度及步长公式,实现介电常数和电导率的同步反演.针对一个合成数据模型同时反演介电常数和电导率,结果表明该方法能够反演出亚波长尺寸异常体的形状和位置.接着,将该方法应用到两组实际数据中,并与基于估计源子波的时间域波形反演结果进行比较.结果表明不依赖源子波的时间域波形反演结果分辨率更高,也更准确.  相似文献   

10.
利用具有良好时频局部化性质的小波基函数时频分析真实地震动,提出一种可行的基于小波分析法调整地震动的方法。实际地震波(El Centro波和TH2TG055波)经过调幅、小波变换和标准反应谱拟合得到调整后的地震波,将其输入到拟建结构模型中进行地震反应分析,对比层间剪力、层间位移及层间位移角等几方面的计算结果,表明选取小波分析后的地震波作为结构时程分析输入能够得到比较准确的数值,满足结构抗震设计的需要,在工程应用方面具有一定的应用参考价值。  相似文献   

11.
The non‐stationary rocking response of liquid storage tanks under seismic base excitations including soil interaction has been developed based on the wavelet domain random vibration theory. The ground motion has been characterized through statistical functionals of wavelet coefficients of the ground acceleration history. The tank–liquid–foundation system is modelled as a multi‐degree‐of‐freedom (MDOF) system with both lateral and rocking motions of vibration of the foundation. The impulsive and convective modes of vibration of the liquid in the tank have been considered. The wavelet domain coupled dynamic equations are formulated and then solved to get the expressions of instantaneous power spectral density function (PSDF) in terms of functionals of input wavelet coefficients. The moments of the instantaneous PSDF are used to obtain the stochastic responses of the tank in the form of coefficients of hydrodynamic pressure, base shear and overturning base moment for the largest expected peak responses. Parametric variations are carried out to study the effects of various governing parameters like height of liquid in the tank, height–radius ratio of the tank, ratio of total liquid mass to mass of foundation, and shear wave velocity in the soil medium, on the responses of the tank. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

12.
逆散射级数法预测层间多次波的算法改进及其策略   总被引:12,自引:5,他引:7       下载免费PDF全文
波动方程预测层间多次波有两类算法.一是由Berkhout和Verschuur(1998)提出的基于CFP(共聚焦点)延拓的算法,另一是由Weglein(1997)等人提出的基于逆散射级数法的算法.ISS(逆散射级数)算法具有不依赖速度模型的优点,但是其计算成本很高.本文改进了Weglein提出的1-D预测公式,提高了计算效率.从理论方面,改进算法的计算速度提高了大约12倍.实际结果表明其计算速度更快,提高近80倍左右,主要原因在于后者的空间复杂度小.此外,由于时间-空间域比伪深度-波数域的层间多次波预测有更多的优点,本文推导了1.5-D时间-空间域ISS层间多次波预测算法,该算法有计算速度快、预测噪音小、适应观测系统能力强和不依赖地震子波等优点.在实际应用中,计算成本高是ISS预测层间多次波算法的主要障碍.因此,本文还讨论了预测层间多次波的策略.结合本文提出的改进算法和应用策略,能够在实际地震资料处理中较好地应用ISS预测层间多次波.  相似文献   

13.
The spectral-profile and wavelet analyses are applied for estimating the periodicity in the depth variations of the elastic wave velocities in the detailed seismoacoustic model of the Kola Superdeep Borehole. The geological confinement of the revealed multiscale periodicity and its relationship with the geodynamical processes in the Earth??s crust are considered. The hierarchical sequence of the sizes of the blocks identified in the section is determined; this hierarchy is consistent with the geophysical structure of the medium.  相似文献   

14.
Some new imaging formulas for seismic reflection wave and their theoretical basis are given. Phenomena of wave propagation should be characterized by instantaneous spectrum and expressed by complex function of three variables (time, space and frequency) in mathematics. Various physical parameters of medium are also complex functions of two variables (space and frequency). The relationship between reflection coefficient of medium and spectrum of reflected wave is given. Multi-reflection and filter of formations are considered in inversion formulas. Problems in classical convolution model and wave equation are illustrated. All these inversion formulas can be used to image underground medium by wavelet transform and method of “3-basic colors”. Different colors mean different media.  相似文献   

15.
16.
One way to investigate the magnetosphere–ionosphere coupling is through the simultaneous observation of different parameters measured at different locations of the geospace environment and try to determine some relationships among them. The main objective of this work is to examine how the solar energetic particles and the interplanetary medium conditions may affect the space and time configuration of the ring current at low-latitudes and also to get a better understanding on how these particles interfere with the lower ionosphere in the South Atlantic Magnetic Anomaly region (SAMA). To accomplish this, the cosmic noise absorption (CNA) and the horizontal component of the Earth's magnetic field data measured from sites located in the SAMA region were compared with the proton and electron fluxes, interplanetary medium conditions (solar wind and the north–south component of the interplanetary magnetic field measured on board satellites), the SYM-H index and magnetometer data from Kakioka (KAK-Japan), located significantly outside the SAMA region. The time series analyzed correspond to the geomagnetic disturbance that occurred on August 25–30, 1998. The analysis was performed by implementing wavelet techniques, with particular attention to singularities detection, which highlights the presence of transient signals. The results are discussed in terms of the first three wavelet decomposition levels of the parameters. The magnitude of wavelet coefficients of the solar wind and proton flux at the two energy ranges analyzed is timely well correlated, indicating that these two signals are energetically linked. The larger wavelet coefficient amplitude of KAK and VSS magnetograms shows time delays that are compatible with an asymmetric configuration of the ring current, considering that at the storm time, VSS was at the dawn sector of the magnetosphere and KAK at the dusk side. The wavelet analysis of CNA signals reveals that the signal may be sensitive to the ionization produced by energetic electrons and protons as well. The time delays observed in wavelet coefficients may give an indication of the different accelerating process to which the particles are submitted when traveling along the magnetic field lines, from higher to lower latitudes, and the likely contribution of these particles to the ionization measured as an absorption of the cosmic noise in the lower ionosphere.  相似文献   

17.

The spectral-profile and wavelet analyses are applied for estimating the periodicity in the depth variations of the elastic wave velocities in the detailed seismoacoustic model of the Kola Superdeep Borehole. The geological confinement of the revealed multiscale periodicity and its relationship with the geodynamical processes in the Earth’s crust are considered. The hierarchical sequence of the sizes of the blocks identified in the section is determined; this hierarchy is consistent with the geophysical structure of the medium.

  相似文献   

18.
利用零偏移距VSP资料在小波域计算介质Q值   总被引:8,自引:3,他引:5       下载免费PDF全文
赵伟  葛艳 《地球物理学报》2008,51(4):1202-1208
讨论了地震子波为一般零相位子波情况下,利用零偏移距垂直地震剖面(VSP)资料在时-频域计算Q值方法.以Morlet小波为基本小波,基于单程波传播理论, 在地震子波为一般的零相位子波情况下,推导出了在时-频域计算Q值的公式.通过对横坐标作恰当的非线性变换后,该公式在形式上与地震子波为脉冲时的计算公式类似.文中采用VSP正演记录,验证了本文提出的方法的正确性.文中还讨论了在时-频域恰当地选择有效频带,降低反射波对Q值计算结果影响的方法.  相似文献   

19.
Jan F. Adamowski   《Journal of Hydrology》2008,353(3-4):247-266
In this study, a new method of stand-alone short-term spring snowmelt river flood forecasting was developed based on wavelet and cross-wavelet analysis. Wavelet and cross-wavelet analysis were used to decompose flow and meteorological time series data and to develop wavelet based constituent components which were then used to forecast floods 1, 2, and 6 days ahead. The newly developed wavelet forecasting method (WT) was compared to multiple linear regression analysis (MLR), autoregressive integrated moving average analysis (ARIMA), and artificial neural network analysis (ANN) for forecasting daily stream flows with lead-times equal to 1, 2, and 6 days. This comparison was done using data from the Rideau River watershed in Ontario, Canada. Numerical analysis was performed on daily maximum stream flow data from the Rideau River station and on meteorological data (rainfall, snowfall, and snow on ground) from the Ottawa Airport weather station. Data from 1970 to 1997 were used to train the models while data from 1998 to 2001 were used to test the models. The most significant finding of this research was that it was demonstrated that the proposed wavelet based forecasting method can be used with great accuracy as a stand-alone forecasting method for 1 and 2 days lead-time river flood forecasting, assuming that there are no significant trends in the amplitude for the same Julian day year-to-year, and that there is a relatively stable phase shift between the flow and meteorological time series. The best forecasting model for 1 day lead-time was a wavelet analysis model. In testing, it had the lowest RMSE value (13.8229), the highest R2 value (0.9753), and the highest EI value (0.9744). The best forecasting model for 2 days lead-time was also a wavelet analysis model. In testing, it had the lowest RMSE value (31.7985), the highest R2 value (0.8461), and the second highest EI value (0.8410). It was also shown that the proposed wavelet based forecasting method is not particularly accurate for longer lead-time forecasting such as 6 days, with the ANN method providing more accurate results. The best forecasting model for 6 days lead-time was an ANN model, with the wavelet model not performing as well. In testing, the wavelet model had an RMSE of 57.6917, an R2 of 0.4835, and an EI of 0.4366.  相似文献   

20.
本文推导了基于位移激励计算单自由度体系拟速度谱公式,通过构造的脉冲位移时程对公式精度进行了验证;之后利用小波变换去除强震记录噪声而保留地震动永久位移,再基于去趋势项方法和滤波方法去除永久位移后,计算拟速度谱。算例结果表明:短周期段内,不保留永久位移的位移激励拟速度谱值与保留永久位移的位移激励拟速度谱值相差很小;中长周期段内,不保留永久位移的位移激励拟速度谱值总体上小于保留永久位移的位移激励拟速度谱值,且不保留永久位移时,滤波方法引起的拟速度谱降幅大于去趋势项方法所引起的拟速度谱降幅。因此,基于位移激励计算中长周期结构的地震反应时,应保留地震动永久位移,或基于去趋势项方法去除永久位移。   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号