Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

高原大地形及海洋热力强迫对东亚及北半球环流和气候变化影响的观测研究

本文总结了“青藏高原大地形及西太平洋暖池势力强迫对东亚及全球气候变化的影响”专题五年来的主要研究工作。其中包括青世故高原东部大气热源的时间演变特征,夏半年高原热源异常对我国降水和北半球环流的影响;西沙海温变化特征及其与我南方降水的关系,北太平洋海温主因子特征及其与华南前汛期降水的变化。     

板块绝对运动及地幔热对流

本文以板块绝对运动AM1-2模型为边界条件探讨了不同的瑞利数下地幔热对流模型.结果表明,瑞利数小于10000（529.41）时,地幔对流呈现以板块驱动图式,运动的极型场和环型场由板块运动激发,两种场占有差不多相同的功率.当瑞利数增加到接近或略超过最低临界值时（约1.5倍）,对流呈现出复杂状态:1.板块运动速率小于下伏地幔对流速率;2.区域性的双层对流环出现;3.对流谱成分发生变化;4.环型场仅在地幔很浅的区域中起作用,而在地幔深部对流图式中影响很小.     

Subsurface thermal and hydrothermal characterization based on geothermal resources map, drill hole thermal gradients, Curie point depths and hypocenter distribution — Examples of Tohoku district, Japan

Subsurface thermal structure in Tohoku district are characterized by existing data such as geothermal resources maps, drill hole thermal gradients, Curie point depths and hypocenters distribution maps. The collected data are registered in a database system, then, compared in plan view, cross-section and bird's-eye pictures. The comparison indicates that subsurface temperatures extrapolated from drill hole thermal gradients are generally concordant to the Curie point depth, assumed to be 650 °C. Tohoku district is generally divided into 5 type areas; fore arc lowland, fore arc mountain country, Quaternary volcanic terrain, back arc lowland and back arc mountain country. The surface thermal manifestations in Quaternary volcanic terrain are mainly controlled by the magma chambers as heat sources, while, surface thermal features such as hot springs in non-volcanic areas are controlled by degrees of heat flows, and hydrothermal flows in permeable Cenozoic formations and along permeable fault zones.     

Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

爆破夯实技术在水下地基处理中的应用

以某杂货码头和集装箱重力式码头水下地基处理工程为例，介绍了其地基处理的爆夯试验过程，得出了一次性爆破夯实厚层块石抛石基床的爆夯参数，并对其施工过程进行了监测。总结了类似工程施工的经验与教训。     

我国入境旅游和经济增长关系分析

文章分析了入境旅游发展的区域差异，并用Granger因果检验来考察入境旅游和各区域经济增长之间是否存在相应的因果关系，Granger因果检验的结果表明在东部地区，入境旅游和经济增长存在单向因果关系。东部地区的入境旅游对区域的经济增长产生显著影响。但是全国和中西部地区的入境旅游与经济增长之间不存在显著的因果关系。类似地，东部地区入境旅游增长是其第三产业增长的Granger原因。但是全国和中西部地区的入境旅游和相应区域第三产业的增长并不存在显著的因果关系。在区域层面上。只有少数区域支持入境旅游是区域GDP增长的Granger原因，它们是东部地区的北京、广东、上海、天津、福建、江苏、浙江、海南和广西；中部地区的黑龙江和湖南以及西部地区的云南和陕西。在这些省区市。入境旅游对区域经济增长产生显著的影响。但总体而言，入境旅游在我国还没有足够的发展。     

近几十年我国极端气温变化特征分区方法探讨

采用聚类统计检验分析和旋转主分量分析相结合确定中心站的方法,利用我国多年极端气温资料,对我国最高和最低气温年际变化型态进行区划。结果表明,这两种方法结合可以互相补充,使分区结果更具客观性。中国极端高温和极端低温年际变化分别可划为12和11个不同类型的区域,分别计算了各区域第一主成分的方差贡献率以及各区域之间的两两相关系数,检验证明分区是合理的。     

高新技术产业集聚的城市化响应与反馈机制研究

高新技术产业集聚与城市化存在互动机制。一方面,高新技术产业结构效应影响了城市化动力机制,集聚效应则促进了城市外部经济的发展,使城市化表现出不同的模式特征和地域特征;另一方面,城市化促进了高新技术发展所需要素资源的聚集与流动,为高新技术产业集聚提供人力资本保障,还为高新技术企业提供创新的源动力。文章以苏州高新区的高新技术产业发展和城市化互动关系佐证了这一理论,并认为建立二者之间切实可行的调控模式是当务之急。     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.