分形分维理论在地质灾害发育及空间分布规律中的应用——以长安区滑坡、崩塌地质灾害为例

长安区区内地层岩性、地质构造多样,加之人类经济工程活动、降雨、地震等诱发影响,秦岭北坡及黄土台塬边缘滑坡崩塌地质灾害发育,是陕西省地质灾害较严重的县(区)之一。本文以长安区为研究区域,在对长安区自然地理及区域地质环境条件等资料充分收集的基础上,结合长安区滑坡崩塌灾害点的实际调查资料,详细分析了长安区的滑坡崩塌灾害发育特征及空间分布规律,并采用网络覆盖法,对长安区滑坡崩塌地质灾害的空间分布进行了分形分维计算。求得地质灾害的空间分布分维值为0.876 6,反应出长安区滑坡崩塌地质灾害空间分布特性的复杂性。研究结果表明,分形分维理论能合理评价地质灾害空间分布特征,对类似地区的地质灾害空间分布特征分析具有借鉴意义。     

Dasymetric Spatiotemporal Interpolation

This research applies the principles of dasymetric mapping to spatiotemporal interpolation by extending the spatial concepts of zone and area to their temporal analogs of interval and duration, respectively. An example application of dasymetric spatiotemporal interpolation using crime event data is presented. Results indicate that dasymetric spatiotemporal interpolation significantly improves the accuracy of estimates over areal or duration weighting. In addition, even when dasymetric interpolation in either the spatial or temporal dimension is relatively weak, combining dasymetric estimation in both space and time dimensions simultaneously has the potential to amplify the accuracy of the overall dasymetric estimation.     

MORPHOMETRIC ANALYSIS OF KARST DEPRESSIONS ON A MEDITERRANEAN KARST MASSIF

Dolines are the most typical landform in karst landscapes. They play an important role in karst hydrogeology; they are effective sediment traps that provide information about past climates and karst evolution, and can also be indicators of tectonic activity. The morphometic analysis of those depressions can help in their study and interpretation. However, detailed identification and delineation of dolines becomes increasingly difficult when the study area is very large, of limited access or covered by vegetation. Here, we have used a method developed by the authors for automatic identification and delineation of karst depressions based on the use of a pit removal algorithm. Next a complete morphometric analysis of the karst depressions has been carried out. The main advantage of the methodology is that by delineating and identifying the depth of each doline from its rim, many morphometric statistics can be obtained, including the relation of doline depth and area with respect to altitude, the main doline orientations, doline density, the fractal dimension of doline area and frequency, and point field analysis, among others. The results are interpreted from a geomorphological point of view and the mapped dolines can be easily integrated into geomorphological mapping. A case study is presented from the large Sierra Gorda karst massif in southern Spain, where 3100 depressions have been identified. A strong structural control of doline shape and orientation is demonstrated; the main areas of doline density are highlighted and doline size is shown to follow a fractal law.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

基于非均质地层参数表征的颗粒流模型构建

非均质地层参数的获取与表征方法是浅层地质体改造的重要依据。基于非均质地层的地质建模,表征具有随机不确定性和模糊不确定性的土体表观参数和物理力学参数。利用 Weibull 分布统计描述局部地层信息,并应用Diamond-Square分形插值方法,结合差分盒子维计算非均质地层分形维数,进一步合理演算整体地层参数。采用PFC3D数值模拟软件,结合马尔科夫链蒙特卡洛方法采样数据,对非均质地层块石颗粒及土颗粒进行建模。通过非均质地层模拟结果,发现其破坏过程与应力应变曲线要比均质土体复杂得多,由于块石大颗粒的含量、分布特性对土样力学性质影响显著。本文建模方法有助于非均质地层参数表征研究。     

2011~2020年勘查地球化学数据处理研究进展

我国建立了包含海量数据的高质量的勘查地球化学数据库,为矿产勘查、环境评价和地质调查等提供了重要的数据支撑。如何高效处理勘查地球化学数据,并从中发掘和识别深层次信息一直是勘查地球化学学科研究的热点和前沿领域。本文在系统调研国内外学者过去十年发表的论著基础上,对勘查地球化学数据处理方法进行分析与对比,从勘查地球化学数据库建设、地球化学异常识别及其不确定性评价等方面概述了我国近十年来在该领域取得的主要研究进展,包括:(1)分形与多重分形模型由于考虑了地球化学空间模式的复杂性和尺度不变性,在全球范围内得到极大的发展和推广,我国学者引领了基于分形与多重分形的勘查地球化学数据处理;(2)机器学习和大数据思维开始在该领域启蒙,并迅速得到关注,正在成为研究热点和前沿领域,我国学者率先开展基于机器学习算法的勘查地球化学大数据挖掘研究;(3)我国学者需要进一步加强勘查地球化学数据缺失值处理以及成分数据闭合效应研究。今后该领域应进一步加强对弱缓地球化学异常识别、异常不确定性评价以及异常识别与其形成机理相结合等方面的研究。     

Brittleness index predictions from Lower Barnett Shale well-log data applying an optimized data matching algorithm at various sampling densities

The capability of accurately predicting mineralogical brittleness index(BI)from basic suites of well logs is desir-able as it provides a useful indicator of the fracability of tight formations.Measuring mineralogical components in rocks is expensive and time consuming.However,the basic well log curves are not well correlated with BI so correlation-based,machine-learning methods are not able to derive highly accurate BI predictions using such data.A correlation-free,optimized data-matching algorithm is configured to predict BI on a supervised basis from well log and core data available from two published wells in the Lower Barnett Shale Formation(Texas).This transparent open box(TOB)algorithm matches data records by calculating the sum of squared errors be-tween their variables and selecting the best matches as those with the minimum squared errors.It then applies optimizers to adjust weights applied to individual variable errors to minimize the root mean square error(RMSE)between calculated and predicted(BI).The prediction accuracy achieved by TOB using just five well logs(Gr,pb,Ns,Rs,Dt)to predict BI is dependent on the density of data records sampled.At a sampling density of about one sample per 0.5 ft BI is predicted with RMSE～0.056 and R2～0.790.At a sampling density of about one sample per 0.1 ft BI is predicted with RMSE～0.008 and R2～0.995.Adding a stratigraphic height index as an additional(sixth)input variable method improves BI prediction accuracy to RMSE～0.003 and R2～0.999 for the two wells with only 1 record in 10,000 yielding a BI prediction error of＞±0.1.The model has the potential to be applied in an unsupervised basis to predict BI from basic well log data in surrounding wells lacking mineralogical measure-ments but with similar lithofacies and burial histories.The method could also be extended to predict elastic rock properties in and seismic attributes from wells and seismic data to improve the precision of brittleness index and fracability mapping spatially.     

梅江河流域年均降雨量空间插值方法研究

分析了降雨空间分布的影响要素及降雨空间插值的主要方法,并在此基础上,利用地理信息系统(GIS)的空间分析功能,综合考虑高程、坡度、坡向等影响降雨空间分布的多种要素,应用协克里金方法,探讨流域降雨量的空间插值问题及其应用,并与反距离权重法、克里金方法等单要素方法进行了降雨空间插值的效果比较.结果表明:应用协克里金方法,在综合考虑高程、坡度、坡向等影响要素的条件下,对降雨进行空间插值,其结果具有比较好的精度.     

Time domain method for calculating free field motion of a layered half-space subjected to obliquely incident body waves

In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.     

基于粒径分维值检验的松散介质渗透系数空间估值——以河北平原为例

为了克服松散含水介质非均质性各向异性对介质渗透系数的求取所带来的困难,以水动力条件为线索,利用互相关法寻找同一时代地层及同一时代地表水的主、次流方向.在主流方向上,建立了石家庄-束鹿-衡水段介质颗粒粒径随距离递变的数学模型,并以此为基础推求无钻孔地段的粒径.同时,运用分形理论,计算了中、上更新统地层主流方向平均粒径空间分布的分维值,并将其作为检验粒径估值正确与否的判据.通过不同岩性的渗透系数-深度模型,实现了用平均粒径预测渗透系数的全过程.