3-D inversion of borehole-to-surface electrical data using a back-propagation neural network

The “fluid-flow tomography”, an advanced technique for geoelectrical survey based on the conventional mise-à-la-masse measurement, has been developed by Exploration Geophysics Laboratory at the Kyushu University. This technique is proposed to monitor fluid-flow behavior during water injection and production in a geothermal field. However data processing of this technique is very costly. In this light, this paper will discuss the solution to cost reduction by applying a neural network in the data processing. A case study in the Takigami geothermal field in Japan will be used to illustrate this. The achieved neural network in this case study is three-layered and feed-forward. The most successful learning algorithm in this network is the Resilient Propagation (RPROP). Consequently, the study advances the pragmatism of the “fluid-flow tomography” technique which can be widely used for geothermal fields. Accuracy of the solution is then verified by using root mean square (RMS) misfit error as an indicator.     

地磁转换函数方法之典型震例及数模对比

转换函数零线等值图在2003年新疆伽师-巴楚M_L7.0地震前密集交汇, 震后开始发散, 进一步证实了地磁转换函数方法对预测地震的有效性。 同时, 为适应全国地磁台站的记录数据由模拟向数字化的转变, 对模拟数据和数字化数据从时间序列和转换函数计算结果两方面进行了对比, 为下一步使用数字化数据计算转换函数进行地震预测研究打下了基础。     

非线性波浪波面追踪的一种新模式

基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。     

用多变量灰色预测模型模拟预测参考作物蒸散量的研究

当观测资料的数据量少而又存在多个相互影响或关联的变量时,常用的灰色预测模型GM（1,1）不能全面考虑多个变量。为此,采用自适应MGM（1,n）模型—多变量灰色预测模型,较好地解决了这一问题。针对一些地区气象数据较少甚至缺失的情况,以内蒙古正蓝旗的气象资料用Penman-Monteith计算的参考作物蒸散量（ET0）为研究对象,运用灰色系统理论建立MGM（1,3）模型,模拟预测参考作物蒸散量变化规律,并与GM（1,1）模型和BP神经网络模型比较,结果表明MGM（1,3）模型有较好的预测效果。     

二元样条最小二乘在地图投影数据加密中的应用

首先讨论了二元样条最小二乘拟合的原理和计算方法,然后给出了该方法在地图投影数据加密中的应用。结论表明,所述方法具有很好的精度和稳定性。     

新疆伽师地区地应力的GPS损伤模拟研究

采用几何损伤理论和方法推导了地壳存在不同方向和大小的断层时的有效应力的公式 ,讨论了利用几何损伤理论来分析和计算地壳中存在不连续断层的理论和方法 ,同时对中国新疆伽师地区的地应力变化用几何损伤理论和方法结合该地区的GPS资料进行了数值计算研究 ,通过比较计算结果与 1 997年 1月至 1 998年 1 2月发生在该地区的MS ≥ 5 .0地震活动可以看出 ,新疆南部的中强地震活动区主要在应力场变化激烈的地区 .     

保山井水位异常的数值模拟

利用1995年7月12日云南孟连Ms7.3地震和1996年2月3日云南丽江Ms7．0地震前保山井水位异常变化的原始资料,并根据介质体膨胀变化率和水位变化的耦合关系式,用有限差分法模拟了保山井水位的变化情况。通过拟合实际观测曲线,得出观测井附近含水层内平均应力变化率的演化情况。     

Impact of 1/8° to 1/64° resolution on Gulf Stream model–data comparisons in basin-scale subtropical Atlantic Ocean models

We investigate the impact of 1/8°, 1/16°, 1/32°, and 1/64° ocean model resolution on model–data comparisons for the Gulf Stream system mainly between the Florida Straits and the Grand Banks. This includes mean flow and variability, the Gulf Stream pathway, the associated nonlinear recirculation gyres, the large-scale C-shape of the subtropical gyre and the abyssal circulation. A nonlinear isopycnal, free surface model covering the Atlantic from 9°N to 47°N or 51°N, including the Caribbean and Gulf of Mexico, and a similar 1/16° global model are used. The models are forced by winds and by a global thermohaline component via ports in the model boundaries. When calculated using realistic wind forcing and Atlantic model boundaries, linear simulations with Munk western boundary layers and a Sverdrup interior show two unrealistic mean Gulf Stream pathways between Cape Hatteras and the Grand Banks, one proceeding due east from Cape Hatteras and a second one continuing northward along the western boundary until forced eastward by the regional northern boundary. The northern pathway is augmented when a linear version of the upper ocean global thermohaline contribution to the Gulf Stream is added as a Munk western boundary layer. A major change is required to obtain a realistic pathway in nonlinear models. Resolution of 1/8° is eddy-resolving but mainly gives a wiggly version of the linear model Gulf Stream pathway and weak abyssal flows except for the deep western boundary current (DWBC) forced by ports in the model boundaries. All of the higher resolution simulations show major improvement over the linear and 1/8° nonlinear simulations. Additional major improvement is seen with the increase from 1/16° to 1/32° resolution and modest improvement with a further increase to 1/64°. The improvements include (1) realistic separation of the Gulf Stream from the coast at Cape Hatteras and a realistic Gulf Stream pathway between Cape Hatteras and the Grand Banks based on comparisons with Gulf Stream pathways from satellite IR and from GEOSAT and TOPEX/Poseidon altimetry (but 1/32° resolution was required for robust results), (2) realistic eastern and western nonlinear recirculation gyres (which contribute to the large-scale C-shape of the subtropical gyre) based on comparisons with mean surface dynamic height from the generalized digital environmental model (GDEM) oceanic climatology and from the pattern and amplitude of sea surface height (SSH) variability surrounding the eastern gyre as seen in TOPEX/Poseidon altimetry, (3) realistic upper ocean and DWBC transports based on several types of measurements, (4) patterns and amplitude of SSH variability which are generally realistic compared to TOPEX/Poseidon altimetry, but which vary from simulation to simulation for specific features and which are most realistic overall in the 1/64° simulation, (5) a basin wide explosion in the number and strength of mesoscale eddies (with warm core rings (WCRs) north of the Gulf Stream, the regional eddy features best observed by satellite IR), (6) realistic statistics for WCRs north of the Gulf Stream based on comparison to IR analyses (low at 1/16° resolution and most realistic at 1/64° resolution for mean population and rings generated/year; realistic ring diameters at all resolutions), and (7) realistic patterns and amplitude of abyssal eddy kinetic energy (EKE) in comparison to historical measurements from current meters.     

初始应力各向异性土的弹塑性模型

剑桥模型沿球应力轴(p 轴)等向塑性体变硬化; 在日本广泛采用的关口 — 太田模型沿初始固结线( K0 线)不等向塑性体变硬化。三轴试验数据表明: 自 K0 状态向伸长方向剪切时, 前者方法计算的体积应变偏小, 而后者方法计算的体积应变偏大。 作者提出一种介于上述两者之间 、考虑初始应力各向异性(如 K0 固结)的不等向塑性体变硬化弹塑性模型。 为了使模型在三维应力下较好地反映土的强度和变形特性, 模型的剪切屈服准则使用 SMP 准则。模型的土性参数与剑桥模型一样, 其预测值与粘土实测值的比较表明, 提出的模型是简单合理的, 可望在实际工程计算中得到使用。     

Phytoplankton, nutrients and hydrography in the frontal zone between the Southwest Indian Subtropical gyre and the Southern Ocean

A survey was made of the Southwest Indian Ocean frontal region between 30 and 50°E containing the Agulhas Return, Subtropical and Subantarctic Fronts. From CTD, SeaSoar and extracted samples the distribution of nitrate, silicate and chlorophyll a is shown to be strongly linked to the front and water mass structure, varying zonally and meridionally. Surface chlorophyll a concentrations were low to the north and south leaving a band of elevated chlorophyll between the Subtropical and Subantarctic Fronts. The low concentration of chlorophyll a to the north, in Subtropical Water, was clearly due to nitrate limitation. Between the Subtropical and Subantarctic Fronts, where the chlorophyll a concentrations were highest, the surface layer showed silicate depletion limiting diatom growth. South of the Subantarctic Front there were deep extending, low concentrations of chlorophyll a, but despite plentiful supplies of macro-nutrients and a well-stratified surface layer, high concentrations of chlorophyll a were absent. Changes from west to east were associated with the meandering of the Southern Ocean Fronts, especially the Subtropical Front, and their strength and proximity to each other. Concentrations of chlorophyll a peaked where the Agulhas Return, Subtropical and Subantarctic Fronts were in close proximity. Combined frontal structures appear to have particularly pronounced vertical stability and are associated with enhanced upwelling of nutrients and leakage of nutrients across the front. Light levels are high within the shallow stable layer. Such conditions are clearly favourable for biological growth and support the development of larger-celled phytoplankton communities.