珠江口盆地西江主洼油气差异分布机制

西江主洼是珠江口盆地一个低勘探程度洼陷，油气差异聚集特征明显，但其机理不清。本文在断陷盆地油气成藏理论指导下，利用研究区地质、地球物理和地球化学等资料，开展了烃源岩、断盖组合、储集体系和运聚模式等方面的研究。结果表明，裂陷期沉降、沉积中心有序迁移，造成主力烃源岩自东向西由文四段迁移至文三段、文一+二段，控制东、西部油气差异分布。恩平组区域性泥岩发育且晚期断裂缺乏，导致油气纵向上更易聚集于下构造层；区域泥岩减薄尖灭或晚期断裂切开盖层的区域，上构造层有一定油气分布。"源-汇"类型从宏观上控制储层优劣，影响下构造层油气富集程度；"仓储"运移是控制上构造层油气规模聚集的主要模式。下构造层古近系应作为重点部署方向，东部围绕文四段、西部围绕文三段、文一+二段烃源岩构成的含油气系统展开；上构造层勘探需关注珠海组，在隆起周边寻找具备"仓储"运移模式的有利区带。     

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

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用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.     

Research on temporal and spatial distribution, evolutionary character and mechanism of crustal deformation field before and after the Tangshan earthquake

Ｒｅｓｅａｒｃｈｏｎｔｅｍｐｏｒａｌａｎｄｓｐａｔｉａｌｄｉｓｔｒｉｂｕ┐ｔｉｏｎ,ｅｖｏｌｕｔｉｏｎａｒｙｃｈａｒａｃｔｅｒａｎｄｍｅｃｈａｎｉｓｍｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｂｅｆｏｒｅａｎｄａｆｔｅｒｔｈｅＴａｎｇｓｈａｎ...     

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.