基于GPS和水准数据的微粒群算法联合反演祁连山北缘断裂三维滑动速率     

刘杰 《大地测量与地球动力学》2019,39(9):906-909

将微粒群算法与位错理论模型相结合，采用中国地壳运动观测网络提供的青藏高原地区2001~2004年GPS测量数据和2000~2006年水准测量数据，通过常规定权和附有相对权比的方法对祁连山北缘断裂的三维滑动速率进行联合反演，并与蚁群算法反演结果进行对比。结果表明，微粒群算法收敛速度快、稳定性高，结合经典位错理论模型，是一种可以有效求解断层三维滑动速率反演问题的优化算法，在大地测量反演领域极具应用潜力。  相似文献   

Gully erosion spatial modelling: Role of machine learning algorithms in selection of the best controlling factors and modelling process     

《地学前缘(英文版)》2020,11(6):2207-2219

This investigation assessed the efficacy of 10 widely used machine learning algorithms (MLA) comprising the least absolute shrinkage and selection operator (LASSO), generalized linear model (GLM), stepwise generalized linear model (SGLM), elastic net (ENET), partial least square (PLS), ridge regression, support vector machine (SVM), classification and regression trees (CART), bagged CART, and random forest (RF) for gully erosion susceptibility mapping (GESM) in Iran. The location of 462 previously existing gully erosion sites were mapped through widespread field investigations, of which 70% (323) and 30% (139) of observations were arbitrarily divided for algorithm calibration and validation. Twelve controlling factors for gully erosion, namely, soil texture, annual mean rainfall, digital elevation model (DEM), drainage density, slope, lithology, topographic wetness index (TWI), distance from rivers, aspect, distance from roads, plan curvature, and profile curvature were ranked in terms of their importance using each MLA. The MLA were compared using a training dataset for gully erosion and statistical measures such as RMSE (root mean square error), MAE (mean absolute error), and R-squared. Based on the comparisons among MLA, the RF algorithm exhibited the minimum RMSE and MAE and the maximum value of R-squared, and was therefore selected as the best model. The variable importance evaluation using the RF model revealed that distance from rivers had the highest significance in influencing the occurrence of gully erosion whereas plan curvature had the least importance. According to the GESM generated using RF, most of the study area is predicted to have a low (53.72%) or moderate (29.65%) susceptibility to gully erosion, whereas only a small area is identified to have a high (12.56%) or very high (4.07%) susceptibility. The outcome generated by RF model is validated using the ROC (Receiver Operating Characteristics) curve approach, which returned an area under the curve (AUC) of 0.985, proving the excellent forecasting ability of the model. The GESM prepared using the RF algorithm can aid decision-makers in targeting remedial actions for minimizing the damage caused by gully erosion.  相似文献   

埋深对恐龙化石的影响效果研究          下载免费PDF全文

陈诚  贾超  张尚坤  杜圣贤  罗文强  田京祥 《山东国土资源》2019,35(5)

已发现的恐龙化石大多数埋藏在不同深度的地层中,不同埋深化石所受风化破坏程度不尽相同。为深入研究恐龙化石的风化机理,探究恐龙化石的保护措施和方法,重点分析埋深产生的侧向压力对恐龙化石强度和破坏特性的影响。通过有限差分软件FLAC3D进行数值模拟,揭示埋深因素对恐龙化石保存的影响程度。试验结果表明:围岩对恐龙化石施加的侧向压力是影响恐龙化石变形和强度特性的一个重要因素。在弹性变形阶段,恐龙化石的初始强度、峰值强度随侧向压力的增大而不断增大。当应力超过了恐龙化石极限强度后进入塑性变形,其初始强度、峰值强度逐渐减小,最终达到残余强度。  相似文献   

基于最大球算法定量表征土体空间孔隙网络     

解朝阳  张巍  姚东方  王坚  朱雨辰 《工程地质学报》2020,28(1):60-68

孔隙网络控制着土体渗流、排水固结与基质吸力等重要工程性质。本文介绍了多孔介质孔隙网络最大球建模基本原理与算法。以显微CT扫描振捣干法生成的南京粉砂试样为例,采用最大球算法建立了试样三维重构模型表征单元体(REV)的空间孔隙网络球棍模型,计算得到了样品REV尺度的孔隙网络参数,统计发现,孔隙半径、喉道半径、孔隙配位数、孔隙截面形状因子、喉道截面形状因子与喉道长度等孔隙参数均近似服从正态分布,孔隙体积近似服从衰减型指数分布。孔隙与喉道半径分别分布在100μm与65μm以内,两者数学期望分别为40. 0μm与18. 0μm;配位数分布在25以内,数学期望为5. 1;孔隙与喉道截面形状因子分别分布在0. 01~0. 04与0. 01~0. 05的区间内,两者数学期望分别为0. 019与0. 033;喉道长度分布在100~800μm以内,数学期望为292. 22μm。同时发现,样品中体积小于1. 5×10^7μm^3的小孔隙数量超过90%。本方法可应用于土体细观孔隙结构的定量表征。  相似文献   

基于CUDA的地震相干体并行算法   总被引：1，自引：0，他引：1  

张全 《地质与勘探》2020,56(1):147-153

相干体技术在地震勘探资料解释方面得到了广泛的应用,由于相干体技术处理的对象是三维地震数据体,所以算法运算时间较长。为了缩短解释周期,本文充分发挥GPU并行计算优势,对C3相干体算法进行并行化分析。从硬盘读取数据到GPU上计算相干值并写入硬盘的整个过程进行分析,剔除了冗余数据的读取,完成了C3相干体算法的并行化设计与实现。最后分别对串行算法与并行算法进行性能测试,结果表明本文设计的并行算法在保证精度的前提下达到了16倍左右的加速比,对加快地震资料解释具有重要意义。  相似文献   

珠江三角洲海相沉积软土压缩和剪切变形特性试验研究   总被引：3，自引：0，他引：3       下载免费PDF全文

张岩  耿济世  毛磊  刘芙荣 《地震工程学报》2018,40(4):745-751

为了研究近代海相沉积软土的变形特征,本文对珠江三角洲河口地区的全新世海相沉积的软土进行了室内单轴压缩与不同围压下作用下的三轴固结不排水(CU)、不固结不排水(UU)的剪切试验。试验结果表明:原状土和重塑土样的孔隙比随压力的增大而下降,且与加压方式有密切关系,加压幅度越小,土的扰动破坏越小;原状土样固结系数随着压力增加而递减,整体近似符合指数函数分布,重塑土样固结系数随着压力增加逐渐增大,整体近似符合线性函数分布。在相同的固结压力下,原状土样的强度低于重塑土样的强度。土样的应力应变曲线为应变硬化型,且围压对重塑土的影响小于原状土。  相似文献   

基于分区间泛化Kriging近似模型的船舶局部结构形状优化     

张干锋  王德禹 《海洋工程》2018,36(6):46-57

为提高Kriging近似模型在船舶结构性能多维度响应预测方面的适用性,对常规Kriging近似模型进行分区间泛化改进:一是对设计样本点的各个维度(分量)进行划区,并在每一个分区间内采用最优拉丁超立方(OLhd)取样;二是引入比例系数w1组合高斯与指数型相关函数提高模型对数据的泛化能力,从而在每个划分的区间内建立泛化的Kriging近似模型。通过SCH测试函数,验证了构造的分区间泛化Kriging模型有效性。结合参数化建模和改进Kriging近似模型对某过渡肘板和舱口角隅边界进行形状优化,优化第一阶段由改进的近似模型通过多岛遗传算法得到全局初步的最优解,第二阶段在初步解的基础上缩小优化变量范围,由少量的FEM计算即可在小范围内搜寻到精确的最优设计变量。结果表明:分区间泛化Kriging近似模型在预测多维度响应时较常规Kriging模型预测精度更高;分阶段的形状优化流程在保证极小误差和缩小计算成本的情况下能够得到理想的应力分布和重量优化结果,有助于船舶结构的轻量化研究。  相似文献   

基于贝叶斯框架的各向异性页岩储层岩石物理反演技术   总被引：1，自引：0，他引：1       下载免费PDF全文

刘财  符伟  郭智奇  刘喜武  刘宇巍 《地球物理学报》2018,61(6):2589-2600

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究，针对页岩黏土含量高、层间微裂缝发育等特点，利用Backus平均理论描述页岩黏土矿物弹性参数，利用Chapman理论计算与水平微裂缝有关的VTI各向异性，并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法，并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息，并以测井纵、横波速度作为约束条件建立反演的目标函数，同时利用粒子群算法进行最优化搜索.计算结果表明，基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高，并且在一定范围内变化，这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关，同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关，而与横波各向异性参数γ相关性较小.另外，页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.  相似文献   

基于F-范数的不确定性平差模型的解算方法     

鲁铁定  朱国红 《大地测量与地球动力学》2018,38(6):557-561

为提高基于F-范数的不确定性平差模型的解算效率，给出直接迭代算法进行参数估计。该算法无需SVD，解算过程简单且易于编程计算，同时给出迭代不收敛时的SVD-解方程算法。二元线性拟合及沉降观测AR模型的算例结果表明，这2种算法正确可行，与SVD-迭代算法具有等价性。当迭代收敛时，宜使用直接迭代算法，收敛速度更快，解算效率更高；当迭代不收敛时，可釆用SVD-解方程算法。  相似文献   

Full moment tensor inversion constrained by double-couple focal mechanism for induced seismicity*          下载免费PDF全文

Yuyang Tan  Haijiang Zhang  Junlun Li  Chen Yin  Furong Wu 《地震科学（英文版）》2020,33(4):177-193

In this study, we propose a new method to determine full moment tensor solution for induced seismicity. This method generalizes the full waveform matching algorithm we have developed to determine the double-couple (DC) focal mechanism based on the neighbourhood algorithm. One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching. Through synthetic tests using a real passive seismic survey geometry, the results show the new constraint can help better recover the DC components of inverted moment tensors. We further investigate how errors in the velocity model and source location affect the moment tensor solution. The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components. We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests. While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints, both solutions are able to fit the observed waveforms at similar levels. The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.  相似文献