地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。     

雄安新区地热地质模型探究：来自地球物理的证据

雄安新区内地热资源丰富，区内有牛驼镇地热田、容城地热田和高阳地热田，地热资源开发利用较早，但是对其深部热源机制仍未形成统一观点。为了研究雄安新区内地热田深部热源机制，在新区及外围进行了深反射地震和长周期大地电磁探测，对取得的同剖面的深反射地震和大地电磁数据进行处理和综合解释，探明了研究区从地表至莫霍面范围内地质构造和电性结构。下地壳结构在深反射地震剖面与大地电磁剖面上有很好的对应关系。电阻率低值区对应着在深反射地震剖面上存在一系列反射同相轴，且同相轴可以延续到莫霍面，电阻率高值区对应着在深反射地震剖面上无明显连续反射同相轴，尤其是在莫霍面之上呈现地震反射近似"空白区"。结合区域地热资料构建了研究区深部地热地质模型，对新区内深部地热机制进行了解释。该模型为"二元"生热模型，其热源包含两个部分，深部地幔热源和地壳放射性元素衰变生热。放射性元素衰变生热占地表热流的接近30%，而幔源热流在地表热流中的占比可达约70%。在牛驼镇下方，莫霍面以上，由于地幔热物质上涌造成下地壳上隆，幔源岩浆底侵作用于下地壳形成了局部热异常，该热异常具有低速高导的地球物理特征，认为是牛驼镇地热田和容城地热田的深部热源；以区域断裂为热通道，大地热流由深部向上传导、扩散到牛驼镇凸起和容城凸起顶部，对碳酸盐岩储水层进行加热，形成地热储层；上覆新近系沉积地层是良好的热盖层。     

A Block-Based Theoretical Model Suited to Gravitational Sliding

An original theoretical model has been devised to simulate mass flow over hill slopes due to gravitational sliding. The sliding mass is discretized into a sequence of contiguous blocks which are subjected to gravitational forces, to bottom friction and to surface resistance stresses that are generally negligible for subaerial flows, but are relevant for submarine slides. The blocks interact with each other while sliding down the hill flanks because of internal forces that dissipate mechanical energy and produce a momentum exchange between the individual blocks, yet conserving the total momentum of the mass. Internal forces are expressed in terms of interaction coefficients depending on the instantaneous distance between the block centers of mass, which is a measure of the deformation experienced by the blocks: the functional dependence includes three parameters, namely the interaction intensity ¯, the deformability parameter and the shape parameter , by means of which a wide range of interaction types can be fully accounted for. The time integration is performed numerically by solving the equations for the block velocities and positions at any time ti by means of the block accelerations at the previous time t_i-1, and by subsequently updating the block accelerations, which allows to proceed iteratively to the following times. The model has been tested against laboratory results available from literature and by means of several numerical experiments involving a simplified geometry both for the sliding body and the basal surface, with the purpose of clarifying the influence of the model parameters on the slide dynamics. The model improves the performance of the existing kinematic models for slides, moreover preserving an equivalent numerical simplicity. Future applications and possible improvements of this model are suggested.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

基于GIS城市洪水淹没模拟分析

城市洪水淹没范围动态模拟分析是城市防洪规划与防治决策的一个重要基础工作。本文主要研究洪水从翻堤口出发在地理空间蔓延、扩散、动态行进及确定淹没范围的数字模拟方法。为此,根据数学形态学及测地圆概念,研究设计了洪水扩散范围的“膨胀”模拟算法和淹没范围搜索算法,用于查询淹没通块中从翻堤处到任一点之间的淹没路径和ｔ时刻洪水扩散范围。     

哀牢山-红河剪切带的热年代学研究——多重扩散域模式的应用实例

洛弗拉等人(1989)把Dodson(1973)建立的单一扩散域模式扩展成矿物扩散域不连续分布的模式,即多重扩散域模式,与之对应具有某一封闭温度范围,因此利用单个K-长石矿物样品可以揭示更长时间段上的冷却过程。对云南哀牢山-红河剪切带的K-长石样品(FA-2-1)进行40Ar/39Ar同位素分析,利用多重扩散域模式,揭示在距今约20～19Ma问该带存在一快速冷却过程,这与西藏南部曲水地区快速抬升的时间是一致的     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

福建省第三产业发展的区域比较分析

运用可拓工程方法和主成分分析方法建立了福建省第三产业发展水平的综合评判物元模型。模型的输出结果表明：福建省各地区第三产业发展水平差距悬殊；综合关联度将第三产业发展水平划分成4个等级，单指标关联度则揭示了各地区第三产业发展的优、劣势。     

预测陕西关中地区需水量的改进GM(1,1)模型

需水量评价与预测是水资源规划与管理一项重要的研究内容。本文采用定额法计算需水量。通过对陕西关中地区需水量要素组成的分析,认为需水量要素属平稳时间序列,且一阶累加生成数据系列满足指数规律,符合灰色预测条件。为了提高GM(1,1)模型的精度,采用一种改变背景值的新方法,即:中心逼近方法。通过精度检验,中心逼近式GM(1,1)模型平均误差百分比、误差平方和分别比传统的GM(1,1)模型提高了0.147和1.579。用中心逼近式GM(1,1)模型预测需水量各要素近期和中长期值,经检验,其中非农业人口、农业人口、耕地面积、工业总产值预测模型满足一级模型,等级为良好。牲畜头数预测模型为不合格模型,分析其原因,主要是因为时间序列数据不完全符合指数规律。通过定额法计算,预计75%保证率时2005年和2010年关中地区工农业需水量分别占总需水量的92.21%和89.75%,生活需水量分别占7%和9.04%。