一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data     

Zhiwei Xu  James Irving  Kyle Lindsay  John Bradford  Peimin Zhu  Klaus Holliger 《Geophysical Prospecting》2020,68(2):678-689

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.  相似文献   

A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations     

Chol Sim  Chunyou Sun  Nam Yun 《Geophysical Prospecting》2020,68(4):1253-1269

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.  相似文献   

平面不规则基础隔震结构抗扭设计研究   总被引：1，自引：0，他引：1       下载免费PDF全文

黄小宁  王威  王宁  杜永峰 《地震工程学报》2020,42(5):1283-1288

针对平面不规则结构在水平地震作用下的振动特性，通过调整隔震层隔震支座的布置，得到3种不同工况的隔震层刚心与上部结构质心、刚心相对位置关系，分别以楼层位移和层间位移为指标的扭转位移比，作为平面不规则基础隔震结构扭转响应指标，利用弹塑性时程分析方法，通过对3种不同工况的扭转指标对比分析研究，提出适用于平面不规则基础隔震结构的抗扭设计方法。结果表明：对于平面不规则结构，应在保证隔震层扭转位移比小于1.2的基础上，使隔震层的刚心和上部结构的刚心分别位于上部结构质心的两侧，可有效控制上部结构的扭转。  相似文献   

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction     

Ruth M. K. Plets  S. Louise Callard  J. Andrew G. Cooper  Joseph T. Kelley  Daniel F. Belknap  Robin J. Edwards  Antony J. Long  Rory J. Quinn  Derek W. T. Jackson 《第四纪科学杂志》2019,34(4-5):285-298

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.  相似文献   

喀喇昆仑山西北部冰川运动速度地形控制特征   总被引：2，自引：2，他引：0  

张莎莎  张震  刘时银  李鑫  黄丹妮  薛乃婷 《冰川冻土》2019,41(5):1015-1025

为了探讨地形和海拔对冰川季节和年平均运动速度的影响程度，利用2013-2018年GoLive数据与ASTER GDEM V2数据对喀喇昆仑山西北部3 389条冰川的地形（坡度、坡向、海拔）和冰川运动速度进行了综合分析。结果表明：冰川表面运动速度在物质平衡线处（3 970~4 770 m）达到最快，是冰川积极维持物质平衡的一种体现。坡度平缓地区在不同海拔下的冰川运动速度有明显的差别，但是不同坡度地区的冰川运动速度随海拔变化的趋势基本一致，均呈现先增大后减小。北坡冰川运动速度较平稳，南坡和西南坡的冰川运动速度（均为0.25 m·d^-1）最快并且变化幅度较大，最小值与最大值相差近4倍。冰川运动速度不是呈现单一的季节性变化，同时还会受到地形的控制。低海拔区域冰川运动速度在消融期（3-6月）较快，中海拔区域在消融前（11月至次年2月）较快。  相似文献   

考虑冲刷影响的筒型基础竖向极限承载力研究     

陆罗观  刘润  练继建  陈广思 《海洋工程》2019,37(3):69-77

大直径宽浅式筒型基础,阻水宽度大,在位工作期间受波浪海流作用,其周围土体易被冲刷。为研究单侧地基土体受冲刷后筒型基础的竖向极限承载力变化,通过引进冲刷率的概念,采用有限元方法研究了不同冲刷率下筒型基础的竖向极限承载力;并基于Meyerhof理论建立了计算不同冲刷率下筒型基础竖向极限承载力的极限平衡方法。研究结果表明,随着冲刷率增大,筒型基础的极限承载力出现不同程度的下降,当冲刷率为0.8时,即筒型基础单侧土体冲刷深度达6.4 m时,筒型基础的竖向极限承载力折减率为3.28%。建立的极限平衡算法可准确计算冲刷条件下筒型基础的竖向极限承载力。  相似文献   

山西省大同市新荣区石墨矿岩石矿物地球化学特征与成因     

梁钰琦  赵永  吴光存  张越  夏庆霖 《地球学报》2020,41(6):827-834

山西省大同市新荣区石墨矿带总体呈北东—南西向展布, 长约22 km, 赋矿层位为中下太古界右所堡组, 岩性为含石墨黑云斜长片麻岩。目前该矿带上已有7个勘查区, 累计查明晶质石墨矿物资源量约 5000万吨, 但对其矿床成因类型研究相对较少。本文收集整理和分析总结了区域内地质资料及化验分析结果, 采用变质岩原岩恢复的四种方法确定该矿带的成因: 从石墨矿带的矿体产状和岩石组合、岩相学和矿物学特征定性判断该矿带矿床为沉积岩区域深变质成矿类型; 矿体样品地球化学微量元素比值(Sr:Ba等)符合副变质岩特征; 稀土元素配分曲线及其研究表明矿体及围岩具有沉积岩变质特征; 以TiO2-F图解、[(al+fm)–(c+alk)]/Si图解以及(al-alk):c图解展示的岩石化学特征半定量地确认新荣区石墨矿带矿床为以黏土岩为主的区域变质成矿。新荣区石墨矿带成因类型的确定填补了该区域矿床成因类型研究的空白, 为下一步地质找矿及建模提供了重要依据; 三种图解相结合, 逐步缩小原岩范围, 锁定原岩类型, 为变质岩原岩恢复及矿床成因研究提供了新的思路和判别方法。  相似文献   

基于深度置信网络的GNSS-IR土壤湿度反演     

陈堃  沈飞  曹新运  朱逸凡 《测绘通报》2020,(9):100-105

基于大地测量型GNSS接收机获取的反射信号反演土壤湿度是GNSS领域的研究热点。为克服常规线性回归和BP神经网络算法等的缺陷，本文提出了一种基于深度置信网络的GNSS-IR土壤湿度反演方法。试验结果表明，基于该方法得到的决定系数、土壤湿度平均绝对误差和均方根误差分别为0.909 8、0.017、0.021，与线性回归和BP神经网络算法相比，与实测数据吻合度更高，可有效提高土壤湿度反演精度，证明了方法的有效性和可靠性。  相似文献   

基于噪声谱比法的内蒙古西部地区台站场地响应研究     

张珂  王鑫  王树波  张红才  宋晓燕 《大地测量与地球动力学》2020,40(6):646-651

利用内蒙古西部12个台站的地脉动噪声数据，采用噪声谱比法研究台站的场地响应情况。分析表明，台站场地响应按曲线形态可分为3类，且可能受地形地貌、局部构造和台基状况等条件影响。对比分析噪声谱比法与Moya方法的场地响应结果发现，场地响应曲线形态基本一致，只有极少数台站存在明显差异。结果表明，内蒙古西部区域大部分台站的台基状况总体较好，场地响应曲线较为平坦，无明显频率放大点。  相似文献