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.  相似文献   

Evaluating soil moisture estimation from ground-penetrating radar hyperbola fitting with respect to a systematic time-domain reflectometry data collection in a boreal podzolic agricultural field     

Chameera Illawathure  Gary Parkin  Sébastien Lambot  Lakshman Galagedara 《水文研究》2020,34(6):1428-1445

The upper 30 cm of the soil profile, which hosts the majority of the root biomass, can be considered as the shallow agricultural root zone of most temperate crops. The electromagnetic wave velocity in the soil obtained from reflection hyperbolas in ground-penetrating radar (GPR) data can be used to estimate soil moisture (SM). Finding shallow hyperbolas in a radargram and minimizing the subjective error associated with the hyperbola fitting are the main challenges in this approach. Nevertheless, we were motivated by the recent improvements of hyperbola fitting algorithms, which can reduce the subjective error and processing time. To overcome the difficulty of finding very shallow hyperbolas, we applied the hyperbola fitting method to reflections ranging from 27- to 50-cm depth using a 500-MHz centre-frequency GPR and compared the estimated moisture with vertically installed, 30-cm-long time-domain reflectometry (TDR) sensors. We also compared TDR and GPR sample areas in a 2-D plane using different GPR survey types and different hyperbola depths. SM measured with TDR and GPR were not significantly different according to Mann–Whitney's test. Our analyses showed that a root mean square error of 0.03 m³ m⁻³ was found between the two methods. In conclusion, the proposed method might be suitable to estimate SM with an acceptable accuracy within the root zone if the soil profile is fairly uniform within the application depth range.  相似文献   

An extended finite volume method and fixed-stress approach for modeling fluid injection–induced tensile opening in fractured reservoirs     

Rajdeep Deb  Patrick Jenny 《国际地质力学数值与分析法杂志》2020,44(8):1128-1144

Fluid injection–induced tensile opening is modeled using an extended finite volume method (XFVM). An embedded fracture strategy is used for the flow problem, that is, the fractures are discretized using finite volume segments without resolving the grid around them. Further, the discontinuities across fractures are modeled using special basis functions. The fracture openings due to enhanced fluid pressure and the associated shear slip due to traction free boundary condition on the fracture segments are both modeled using these special discontinuity basis functions. Mass transfer between fractures and matrix is modeled using the pressure difference. The enhancement of fracture storativity due to tensile opening leads to stronger coupling between flow and mechanics. An iterative scheme relying on the fixed-stress approach for fractures, which conserves the stress dependent terms over each iteration of the flow problem, has been introduced. Tensile opening has been simulated for single fractures embedded in two- and three-dimensional matrices. The convergence criterion for sequentially implicit fixed-stress scheme for fractures embedded in elastic media is established and has been validated numerically. Further, for 2D simulations, the effect of the matrix permeability for fracture propagation due to tensile opening has been studied.  相似文献   

频率域海洋可控源电磁垂直各向异性三维反演   总被引：1，自引：0，他引：1       下载免费PDF全文

彭荣华  胡祥云  李建慧  刘亚军 《地球物理学报》2019,62(6):2165-2175

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性，若仅对海洋可控源电磁（MCSEM）数据进行常规各向同性反演，有可能无法获得准确的反演解释结果，从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性（VTI）条件下频率域海洋可控源电磁数据三维反演算法.其中，三维正演采用基于二次场控制方程的交错网格有限体积法，并利用直接矩阵分解技术来求解离散所得的大型线性方程组，有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后，对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演，结果表明：各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释，而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布，展现出对海底电性各向异性结构更为优良的反演解释能力.  相似文献   

PBO钻孔应变同震响应分析     

陈燚飞  吕品姬 《大地测量与地球动力学》2019,39(7):765-770

以2017-09-08墨西哥8.1 级地震为例，分析14个PBO综合钻孔应变仪和地震仪同址记录的地震波信号，结果表明，应变仪记录的地震波符合射线传播规律。对不同震中距的站点进行基于S变换的应变仪与地震仪的同震响应特征分析，二者同震响应一致性较好，应变仪频段主要在0.25 Hz以下。取震后1 h应变数据进行分析发现，同震阶段应变主方向指向震中。  相似文献   

不同类型洞体应变仪的观测差异性分析          下载免费PDF全文

孙伶俐  蒋玲霞  吕琳  张辉  邓建刚 《地震工程学报》2020,42(1):129-135

运用SPSS中聚类及PCA方法对湖北省7个台站三类基线(石英管、铟钢棒、铟钢管)洞体应变观测记录进行相关统计分析,归纳总结其潮汐参数变化特征差异,并针对观测中自然环境因素影响进行定量分析。结果显示:(1)三类基线中,铟钢棒类M2波潮汐因子、中误差及M1波相对噪声水平较低,铟钢棒管年均漂移量最小;(2)自然环境因素(气压、降水、雷雨、温度)对湖北地区洞体应变观测影响较为显著,等量气压影响下应变量最小的为铟钢管型,等量温度影响下应变量最小的为铟钢棒型。  相似文献   

基于GPS应变、地震应变率与震源应力场对帕米尔高原现今构造变形特征的分析          下载免费PDF全文

潘正洋  邵志刚  周云  王洵  解孟雨  艾萨·伊斯马伊力 《地球物理学报》2020,63(3):1143-1154

研究帕米尔高原的构造变形特征对于理解印度板块向北推挤过程中的应变分配方式以及应力转换模式具有重要的意义.本文利用区域GPS应变场、地震应变场与震源应力场分析帕米尔高原的构造形变特征.主要结论为:(1)该区域变形主要以NNW-SSE或近N-S向的挤压为主,在高原内部伴有明显的近ENE-WSW或E-W向拉张,应力方向在帕米尔高原与塔吉克盆地区域呈现逆时针旋转的趋势,而在塔里木盆地则显示几乎与帕米尔高原的一致的应力状态,这可能与两侧盆地块体的强度差异有关.(2)安德森断层参数A∅显示帕米尔高原北缘与西侧区域为逆断层应力状态,在高原内部为正断层应力状态,这与GPS应变的结果显示的应变主要集中在主帕米尔断裂与阿莱谷地附近而在高原内部应变较低是一致的,另外应力在喀喇昆仑断裂北段的方向基本平行于断层走向,以及断层北端较低的滑动速率,这说明了地壳挤压缩短可能是帕米尔高原主要的的构造变形特征,并不支持由于边界走滑断裂导致的应变分异或者块体挤出的模式.(3)综合考虑地震应变方向与SHmax从帕米尔北部NNW-SSE方向到天山北部的近N-S方向的转换,GPS应变方向在帕米尔高原两侧盆地都存在不同程度的旋转,应力场安德森参数A∅显示的应力状态以及SKS的结果显示的近ENE-WSW方向,我们认为印度板块向北推挤与天山造山带碰撞导致帕米尔高原不对称的径向逆冲是帕米尔高原现今构造变形的主要成因与构造模式.  相似文献   

Mapping 3D visibility in an urban street environment from mobile LiDAR point clouds     

Yi Zhao  Jianping Wu  Song Shu  Handong Liang  Min Liu 《地理信息系统科学与遥感》2020,57(6):797-812

ABSTRACT

Visibility determination is a key requirement in a wide range of national and urban applications, such as national security, landscape management, and urban design. Mobile LiDAR point clouds can depict the urban built environment with a high level of details and accuracy. However, few three-dimensional visibility approaches have been developed for the street-level point-cloud data. Accordingly, an approach based on mobile LiDAR point clouds has been developed to map the three-dimensional visibility at the street level. The method consists of five steps: voxelization of point-cloud data, construction of lines-of-sight, construction of sectors of sight, construction of three-dimensional visible space, and calculation of volume index. The proposed approach is able to automatically measure the volume of visible space and openness at any viewpoint along a street. This approach has been applied to three study areas. The results indicated that the proposed approach enables accurate simulation of visible space as well as high-resolution (1 m × 1 m) mapping of the visible volume index. The proposed approach can make a contribution to the improvement of urban planning and design processes that aim at developing more sustainable built environments.  相似文献   

多层构造油藏纵向细分计算单元对储量参数及结果的影响     

颜冠山  刘宗宾  宋洪亮  韩雪芳  王欣然 《地质与资源》2020,29(4):342-350

为满足油田中后期精细开发及后续调整挖潜的需要，亟需进行更精细的储量研究，其中最重要的一环就是纵向细分计算单元.对于多层构造油藏，纵向细分主要根据油层组内隔夹层分布特征、小层物性特征，细分到以小层或者分布特征和物性相近且纵向上连续分布的小层组合为计算单元.根据纵向细分计算单元前后储量参数选值的理论推导，结合储层的沉积展布特征，认为含油面积范围内无储层尖灭时，仅平面油层厚度差异较大的油藏平均有效厚度变小.含油面积范围内储层尖灭时，正常三角洲沉积储层的油藏各参数一般都变小；而辫状河三角洲沉积储层的油藏一般平均有效厚度变小，平均有效孔隙度和平均含油饱和度变大.这一结论可以有效指导储量评估过程中纵向细分计算单元方案划分，为同类油田的精细储量研究提供技术支撑.  相似文献   

Dynamic characteristics of marine soft clay under variable phase difference and initial static shear stress     

Wenqian Dong  Xiuqing Hu  Hongtao Fu 《Marine Georesources & Geotechnology》2020,38(7):770-785

Abstract

Under seismic loading, the soil layer is subjected to multidirectional cyclic shear stress with different amplitudes and frequencies because of the coupling of multiple shear waves and the soil element within a slope or behind a retaining wall is subjected to initial static shear stress before subjected to cyclic loading. Due to the complexity of seismic loading propagation, a phase difference exists between the initial static shear stress and cyclic shear stress. To investigate the influence of the phase difference and initial static shear stress on cyclic shear strain, cyclic modulus, and cyclic strength, a series of laboratory tests are performed on Wenzhou marine soft clay by multi-directional simple shear system, which can simulate the actual state better by controlling the horizontal cyclic stress in the x and y directions simultaneously. As the phase difference varies from 0° to 90°, the dynamic shear modulus increases and cyclic strain accumulation decreases with an increasing number of cycles. The shear strain increases with the initial shear stress.  相似文献