成都国际会展中心输气管线府河穿越隧道涌水量预测

成都拟建国际会展中心供气工程府河输气隧道在府河河底以下40m,采用“竖井＋隧道＋竖井”方案。隧址位于府河及其Ⅰ级阶地,隧道涌水问题严重。为保证施工、运营的安全,在分析隧址区水文地质条件的基础上,对隧道涌水做出初步预测评价,为施工提供依据.     

Bayesian inversion of joint SH seismic and seismoelectric data to infer glacier system properties

In glacial studies, properties such as glacier thickness and the basement permeability and porosity are key to understand the hydrological and mechanical behaviour of the system. The seismoelectric method could potentially be used to determine key properties of glacial environments. Here we analytically model the generation of seismic and seismoelectric signals by means of a shear horizontal seismic wave source on top of a glacier overlying a porous basement. Considering a one-dimensional setting, we compute the seismic waves and the electrokinetically induced electric field. We then analyse the sensitivity of the seismic and electromagnetic data to relevant model parameters, namely depth of the glacier bottom, porosity, permeability, shear modulus and saturating water salinity of the glacier basement. Moreover, we study the possibility of inferring these key parameters from a set of very low noise synthetic data, adopting a Bayesian framework to pay particular attention to the uncertainty of the model parameters mentioned above. We tackle the resolution of the probabilistic inverse problem with two strategies: (1) we compute the marginal posterior distributions of each model parameter solving multidimensional integrals numerically and (2) we use a Markov chain Monte Carlo algorithm to retrieve a collection of model parameters that follows the posterior probability density function of the model parameters, given the synthetic data set. Both methodologies are able to obtain the marginal distributions of the parameters and estimate their mean and standard deviation. The Markov chain Monte Carlo algorithm performs better in terms of numerical stability and number of iterations needed to characterize the distributions. The inversion of seismic data alone is not able to constrain the values of porosity and permeability further than the prior distribution. In turn, the inversion of the electric data alone, and the joint inversion of seismic and electric data are useful to constrain these parameters as well as other glacial system properties. Furthermore, the joint inversion reduces the uncertainty of the model parameters estimates and provides more accurate results.     

Confidence regions for GPS baselines by Bayesian statistics

 The global positioning system (GPS) model is distinctive in the way that the unknown parameters are not only real-valued, the baseline coordinates, but also integers, the phase ambiguities. The GPS model therefore leads to a mixed integer–real-valued estimation problem. Common solutions are the float solution, which ignores the ambiguities being integers, or the fixed solution, where the ambiguities are estimated as integers and then are fixed. Confidence regions, so-called HPD (highest posterior density) regions, for the GPS baselines are derived by Bayesian statistics. They take care of the integer character of the phase ambiguities but still consider them as unknown parameters. Estimating these confidence regions leads to a numerical integration problem which is solved by Monte Carlo methods. This is computationally expensive so that approximations of the confidence regions are also developed. In an example it is shown that for a high confidence level the confidence region consists of more than one region. Received: 1 February 2001 / Accepted: 18 July 2001     

基于AIS的海上渔船捕捞活动聚集区提取方法研究

船舶自动识别系统(Automatic Identification System, AIS)不仅是海上交通监管的有效工具,也为研究海上交通运输及其相关产业活动特征提供了一种良好的数据源。基于海上渔船AIS数据,本研究利用高斯混合模型(Gaussian Mixed Model,GMM)识别渔船捕捞活动状态,提出一种将核密度估计(Kernel Density Estimation,KDE)与热点分析(Hot Spot Analysis, HSA)相融合用于渔船捕捞活动聚集区提取的方法。结果显示：与单一KDE或HSA方法相比,二者相融合的方法将KDE的距离衰减效应与HSA统计指数相结合,在渔船捕捞活动聚集区提取中的应用效果较好、效率较高;采用该融合方法,基于2018年9—12月AIS数据,实现对渤海海峡周边海域渔船捕捞活动聚集区的提取,发现不同月份,渔船捕捞活动聚集区的分布范围和空间形态特征具有一定差异性,烟威近岸海域和渤海海峡是主要的捕捞活动聚集区,其结果可为该海域捕捞活动管理和海洋生态保护提供技术方法和决策支持。     

西藏地区复杂地形下的降水空间分布估算模型

本文提供了一个描述西藏地区年、季降水量空间分布的估算模型.利用卫星遥测数字化地形高程资料和西藏地区仅有的27个常规气象站的多年平均降水整编资料,根据地形坡向站点分为三类.再采用多元逐步回归方法,建立西藏地区的年、季降水量和经度、纬度、海拔高度、坡度、坡向、遮蔽度等6个地理、地形因子之间的关系模型,估算西藏地区降水量的空间分布.结果表明,此方法建立的关于西藏地区降水量与诸因子之间方程的相关性显著,平均绝对误差、相对误差分别为0.93mm和1.16%,对估算模型进行F检验,均通过置信度为0.95的相关检验,回归效果较显著.分析表明估算降水能够定量、定性地再现西藏地区的实际降水分布.     

Estimation of a continuous spatio-temporal population model

In this paper we propose a continuous spatio-temporal model that describes population change in a region in terms of population growth, migration drift towards regions with better economic or climate conditions, and population diffusion from more populated to less populated areas. Finite-differences are used to approximate the space and time derivatives. The model is estimated by using population data from the US census corresponding to the period 1790–1910. People tend to migrate from east to west, and to relocate towards regions with lower precipitation levels and more abundant coal and iron resources. Also population growth tends to be larger in zones with higher precipitation levels and higher temperatures.     

粤港澳大湾区海岛空间分布特征探析

海岛是海洋资源开发利用价值的核心内容之一。为合理保护及开发利用粤港澳大湾区海岛资源,促进海岛产业布局优化,从海岛资源空间分布角度出发,探析海岛在粤港澳大湾区发展中重点平台的作用。以粤港澳大湾区纳入中国海域海岛标准名录的850个海岛为研究对象,在对其地理坐标、陆域面积等基础数据进行数理统计的基础上,通过ArcGIS软件平台,采用最邻近距离法、核密度估计法、洛伦兹曲线和基尼系数方法,对粤港澳大湾区海岛的空间分布特征进行了分析与研究。结果表明:粤港澳大湾区海岛在空间上趋于集聚分布,呈现总体带状分布,局部组团分布的特点,形成了湾区西部海岛分布带和大鹏湾-大亚湾区、九州列岛区、万山-佳蓬担杆列岛区、川山群岛区4个明显的集聚区;粤港澳大湾区海岛以沿岸岛为主,面积分布呈现不均衡性的特点,基尼系数达0.58。     

唐山与通海两大地震的强余震的差别及其原因探讨

1970年通海地震序列与1976年唐山地震序列的对比表明,两次地震的余震衰减,余震强度和频度,余震分布区的形状等有区别.对此提出了一种解释：若余震分布区域的长轴长度a与短轴长度b之比a/b越大,则震源体越接近于一个面,在一个面上进行的调整过程应当比较简单,余震的衰减比较快,强余震比较少而小;反之则震源体越表现出一个体的形态,在一个体内进行的调整过程应当比较复杂,余震衰减比较慢,强余震比较多而大.对于那些在主震后能够很快进行定位并勾画出余震区域轮廓的情况,可以尝试利用a/b快速判断后面强余震的多少以及余震衰减的快慢.     

模型参数估计的优效性

本文根据多参数的R-C不等式和优效计量定理,得到了函数模型和随机模型参数估计的方差下界公式。函数模型的参数估计是优效估计量,但随机模型的参数估计即方差分量估计不是优效估计量;作为特例,单位权方差的估计是优效估计量,同时证明了函数模型和随机模型参数估计的不相关性。     

分析海浪方向谱的扩展本征矢方法Ⅰ.方法的导出

基于交叉谱矩阵可以按本征值划分为信号和含有噪音部分的思想，提出了一种自直接测量数据估计海浪方向谱的方法。该方法称为扩展本征矢方法（EEV），可应用于单点测量系统、仪器阵列以及由二者构成的复合阵列。现有的某些估计方法（如最大似然方法及其扩展形式等）仅是此方法的某种特例。