重磁异常的人机联作校正－迭代反演

本文提出了三度体重磁异常的人机联作校正－迭代反演方法。该方法用二度半组合多边形棱柱体来逼近三度体，从而把三度体重磁异常反演问题转化为二度半体的反演问题；为了消除组合体迭加场的影响，该方法采用了一种校正－迭代技术。理论模型反演计算表明，该方法实际可行。     

岩石节理粗糙度系数的分形特征

岩石节理粗糙度系数JRC是估算节理抗剪强度和变形指标最重要的参数。通过对简易纵剖面仪获取的节理表面轮廓曲线的分形研究,讨论了节理表面轮廓曲线的自相似性和JRC的自相似性,并根据实测统计资料的分析,指出了分形理论研究JRC的适用条件和有效的使用方法。由实测统计资料的JRC尺寸效应自相似性分析,认为JRC尺寸效应具分形结构。本文介绍了一种确定JRC尺寸效应分维数D的方法,由此确定的分维数D具有明确的物理意义。     

Study of the characteristics of earthquake precursory anomalies of He and H2 Beijing No. 2 Cotton Mill well and their relation with earthquakes

This paper systematically summarized the relation of the abnormal change of He and H₂ dissolved gasses in deep hot-water well at Beijing No. 2 Cotton Mill and their relation with seismicity of the region. An empirical formula for the relation of the precursory anomaly of He and H₂ gasses with the magnitude and occurrence time of earthquake has been derived. Two events of magnitude greater than 5 in North China had been successfully predicted by using this formula. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 490–497, 1993.     

大型抽水试验最优化设计研究

抽水试验最优化设计,是指针对一个具体的地下水系统优选出一个抽水试验设计,它使抽水试验的费用最小化,而试验所获资料能满足识别参数所要求的可靠性,它可以概化为一个非线性的混合整数规划问题,其决策变量是确定抽水井和观测井的位置及井数、抽水量和观测频率等,采用搜索算法求解了这一优化问题,得到了理论最优解,并进行了实例演算。     

重力资料的三维密度反演

分别论述了在空间域与频率域中由重力资料反演三维密度分布的理论及其实现的方法,并提出了一种简便、实用的频率域三维密度反演算法。     

3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用

3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换（NHOPRT）地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.     

Bayes判别分析法与地震短临预测

建立了一套完整的Bayes判别分析方法,并提出了明确的衡量综合预测结果所冒风险大小的风险代价比Kdn该方法可用于不同时间尺度的地震综合预测。选取云南地区1994～1996年的水氡观测资料,采用X2统计检验法识别前兆异常,利用Bayes判别分析方法对云南地区的强震活动性进行了内符检验和外推预测。在风险代价比Kdn取4的情况下,内符检验的R值可达0．54,外推时空占有率0．07,获得了满意的效果。     

韧性剪切作用下千枚岩中SiO2的迁移模拟实验—以金山金矿为例

实验较好地模拟了韧性剪切作用过程中ＳｉＯ２“释放－局部迁移”的迁移过程。实验结果表明：４００℃条件下,绢云母,绿泥石等矿物在韧性剪切作用下分解,解释出大量ＳｉＯ２,在成矿流体的参与下发生局部迁移,聚集。要使ＳｉＯ２发生大规模的迁移,富集,还需要一个相对开放的环境（如脆性剪切变形）。     

Estimation and improvement in the geolocation accuracy of rational polynomial coefficients with minimum GCPs using KOMPSAT-3A

ABSTRACT

In this paper, we propose a method to regenerate Rational Polynomial Coefficients (RPCs) using KOMPSAT-3A imagery and to reduce the geolocation error using minimum ground control points (GCPs). To estimate the new RPCs, the physical sensor model fitted to KOMPSAT-3A imagery was utilized and virtual GCPs over the study area were created. The size of the virtual grid used was 20x20x20. To remove the sensor-related errors in physical sensor model, three different image correction models (image coordinate translation model, shift and drift model, and affine transformation model) were additionally applied. We evaluated our proposed method in two areas within Korea, one in urban (Seoul) and one in rural (Goheung) areas. The results showed that there was a significant improvement after applying the suggested approach in the two areas. The image coordinate translation model is suggested in terms of GCP requirement and expected errors estimated from the error propagation analysis using Gauss–Markov Model (GMM).     

Multi-scale remote sensing to support insurance policies in agriculture: from mid-term to instantaneous deductions

ABSTRACT

Climate change is today one of the biggest issues for farmers. The increasing number of natural disasters and change of seasonal trends is making insurance companies more interested in new technologies that can somehow support them in quantifying and mapping risks. Remotely sensed data, with special focus on free ones, can certainly provide the most of information they need, making possible to better calibrate insurance fees in space and time. In this work, a prototype of service based on free remotely sensed data is proposed with the aim of supporting insurance companies’ strategies. The service is thought to calibrate annual insurance rates, longing for their reduction at such level that new customers could be attracted. The study moves from the entire Piemonte region (NW Italy), to specifically focus onto the Cuneo province (Southern Piemonte), which is mainly devoted to agriculture. MODIS MOD13Q1-v6 and Sentinel-2 L2A image time series were jointly used. NDVI maps from MODIS data were useful to describe the midterm phenological trends of main crops at regional level in the period 2000–2018; differently, Sentinel-2 data permitted to map local crop differences at field level in 2016 and 2017 years. With reference to MODIS data, the average phenological behavior of main crop classes in the area, obtained from the CORINE Land Cover map Level 3, was considered using a time series decomposition approach. Trend analyses showed that the most of the crop classes alternated three phases (about 7 years) suggesting that, presently, this is probably the time horizon to be considered to tune mid-term algorithms for risk estimates in the agricultural context. Crop classes trends were consequently split into three phases and each of them modeled by a first-order polynomial function used to update correspondent insurance risk rate. Sentinel-2 data were used to map phenological anomalies at field level for the 2016 and 2017 growing seasons; shifts from class average behavior were considered to locally and temporarily tune insurance premium around its average trend as described at the previous step. Synthesizing, one can say that this approach, integrating MODIS and Sentnel-2 data, makes possible to locally and temporarily calibrate premiums of indexed insurance policies by describing the average trends of crop performance (NDVI) at regional level by MODIS data and refining it at field and specific crop level by Sentinel-2 data.