A parallelizable method for two‐phase flows in naturally‐fractured reservoirs     

Jim Douglas Jr.  Felipe Pereira  Li‐Ming Yeh 《Computational Geosciences》1997,1(3-4):333-368

A parallelizable, semi‐implicit numerical method is proposed for the study of naturally‐fractured reservoir systems. It has proved to be computationally efficient in producing accurate numerical solutions for the dual‐porosity model for immiscible, two‐phase flow in such reservoirs. The method combines hybridized mixed finite elements, a new version of the modified method of characteristics, a sophisticated operator‐splitting procedure for separating the pressure calculation in the fractures from that of the saturation, another operator splitting to handle the interaction of the matrix blocks and the fractures, and domain decomposition iterative procedures for both the pressure and the saturation. It permits moderately long time steps for the pressure and the saturation in the fractures and matrix blocks by using short, inexpensive microsteps to treat the transport portion of the saturation equation in the fractures. This paper is devoted to the formulation of the method and a discussion of numerical results for five‐spot and vertical cross‐section examples.  相似文献   

北京/上海-香港客源潜力研究^*   总被引：2，自引：3，他引：2  

金凤君  叶嘉安 《地理研究》1997,16(4):7-14

文章对内地主要城市与香港间的客源潜力进行了研究,包括居民的赴港愿望、交通方式选择,以及上述决策与收入、职业、费用、旅行时间、安全舒适性等因子的关系。分析得出:内地居民的赴港愿望和交通方式选择具有明显特征,以上海、北京为源地或中转地的居民,作为未来访港的客源,具有较大潜力,并在交通方式选择上有较强倾向性,偏好直达列车、直达飞机等便捷、快速的交通方式,且被访者的职业、收入水平等对其旅行方式选择有重要影响;对北京、上海近期与远期赴港客源潜力作了测算。  相似文献   

Effect of Near-fault Terrain upon Dislocation Modeling     

B.-S. Huang  Y. T. Yeh 《Pure and Applied Geophysics》1997,150(1):1-18

—Coseismic surface deformation provides important information needed to determine source rupture geometry and slip distribution as well as to estimate seismic moment. In this study, numerical experiments were designed to analyze and classify how free-surface topography affects surface deformation. The investigation was performed by 3-D finite element modeling. Results of this study show that crustal deformation induced by near-fault terrain is significant and can be measured with present geodetic survey techniques. The characteristics of the terrain effects show that a hill structure produces more crustal deformation than a half-space model, and that the crustal deformation of a basin structure is less than that of the half-space model. The topographic correction is in the order of five percent of the fault dislocation. On the basis of the relationship between fault offset and earthquake magnitude, it is suggested that the terrain effects on the coseismic crustal deformation of shallow earthquakes with a magnitude greater than 5.6 should be considered as one of the major errors in coseismic deformation modeling which ignored the surface topography on the order of 300 meters.  相似文献   

Linear spatial interpolation: Analysis with an application to San Joaquin Valley     

H. A. Loaiciga  R. H. Shumway  W. W. -G. Yeh 《Stochastic Environmental Research and Risk Assessment (SERRA)》1988,2(2):113-136

The properties of linear spatial interpolators of single realizations and trend components of regionalized variables are examined in this work. In the case of the single realization estimator explicit and exact expressions for the weighting vector and the variances of estimator and estimation error were obtained from a closed-form expression for the inverse of the Lagrangian matrix. The properties of the trend estimator followed directly from the Gauss-Markoff theorem. It was shown that the single realization estimator can be decomposed into two mutually orthogonal random functions of the data, one of which is the trend estimator. The implementation of liear spatial estimation was illustrated with three different methods, i.e., full information maximum likelihood (FIML), restricted maximum likelihood (RML), and Rao's minimum norm invariant quadratic unbiased estimation (MINQUE) for the single realization case and via generalized least squares (GLS) for the trend. The case study involved large correlation length-scale in the covariance of specific yield producing a nested covariance structure that was nearly positive semidefinite. The sensitivity of model parameters, i.e., drift and variance components (local and structured) to the correlation length-scale, choice of covariance model (i.e., exponential and spherical), and estimation method was examined. the same type of sensitivity analysis was conducted for the spatial interpolators. It is interesting that for this case study, characterized by a large correlation length-scale of about 50 mi (80 km), both parameter estimates and linear spatial interpolators were rather insensitive to the choice of covariance model and estimation method within the range of credible values obtained for the correlation length-scale, i.e., 40–60 mi (64–96 km), with alternative estimates falling within ±5% of each other.  相似文献   

Numerical Identification of Parameters in Leaky Aquifers   总被引：2，自引：0，他引：2  

Hund-Der Yeh  Hung-Yuang Han 《Ground water》1989,27(5):655-663

  相似文献   

Attenuation relations of Arias intensity based on the Chi-Chi Taiwan earthquake data   总被引：1，自引：0，他引：1  

H. Hwang  C. K. Lin  Y. T. Yeh  S. N. Cheng  K. C. Chen 《Soil Dynamics and Earthquake Engineering》2004,24(7):509-517

This paper presents a derivation of attenuation relations of Arias intensity for various site conditions based on the strong-motion data recorded in the 1999 Chi-Chi Taiwan earthquake. The data are from the mainshock and three large aftershocks with stations in the footwall area and in the area away from the fault. At each station, Arias intensity is computed from two horizontal acceleration time histories. The Arias intensity data are separated into four groups according to site classes assigned to recording stations. For each site class, the attenuation relation of Arias intensity is derived using a two-step regression analysis method. The attenuation relations established in this study can be used to estimate Arias intensity from a rupture of a thrust fault for sites in the footwall area or in the area away from the fault.  相似文献   

长时间序列Radarsat图像的神经网络模拟及土地覆盖变化的快速检测     

钱峻屏  黎夏  叶嘉安  艾彬  刘凯  陈晓越 《遥感学报》2007,11(6):931-940

土地覆盖的短期时空变化模式研究,对土地覆盖的快速、动态监测具有重要意义,也是遥感研究的新热点。本文利用2000—2001年的时间序列Radarsat图像,采用功率谱分析方法,对土地覆盖的短期时—空变化的周期特征进行了分析,由此建立了基于时间序列影像分析的神经网络预测模型,从植被主要生长季节的时间序列雷达卫星影像获取训练样本,对研究区域的典型土地覆盖的短期动态变化过程进行了学习。学习后的模型能够利用多个时间序列的Radarsat影像对下一时刻的影像进行模拟,并进一步检测变化。在模拟结果基础上,定义相对变化距离函数和检测门限,对模拟影像及实际影像中的变化区域进行了检测。检测精度范围在66.67%(农村居民点)—91.67%(水体)之间,平均检测精度为81.66%。由于时间序列信号的引入,神经网络模型能够较好地获取土地覆盖的短期动态变化信息。  相似文献   

Experimental verification of a wireless sensing and control system for structural control using MR dampers     

Chin‐Hsiung Loh  Jerome P. Lynch  Kung‐Chun Lu  Yang Wang  Chia‐Ming Chang  Pei‐Yang Lin  Ting‐Hei Yeh 《地震工程与结构动力学》2007,36(10):1303-1328

The performance aspects of a wireless ‘active’ sensor, including the reliability of the wireless communication channel for real‐time data delivery and its application to feedback structural control, are explored in this study. First, the control of magnetorheological (MR) dampers using wireless sensors is examined. Second, the application of the MR‐damper to actively control a half‐scale three‐storey steel building excited at its base by shaking table is studied using a wireless control system assembled from wireless active sensors. With an MR damper installed on each floor (three dampers total), structural responses during seismic excitation are measured by the system's wireless active sensors and wirelessly communicated to each other; upon receipt of response data, the wireless sensor interfaced to each MR damper calculates a desired control action using an LQG controller implemented in the wireless sensor's computational core. In this system, the wireless active sensor is responsible for the reception of response data, determination of optimal control forces, and the issuing of command signals to the MR damper. Various control solutions are formulated in this study and embedded in the wireless control system including centralized and decentralized control algorithms. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Source of low frequency modulation of ENSO amplitude in a CGCM     

Byung-Kwon Moon  Sang-Wook Yeh  Boris Dewitte  Jong-Ghap Jhun  In-Sik Kang 《Climate Dynamics》2007,29(1):101-111

We study the relationship between changes in equatorial stratification and low frequency El Niño/Southern Oscillation (ENSO) amplitude modulation in a coupled general circulation model (CGCM) that uses an anomaly coupling strategy to prevent climate drifts in the mean state. The stratification is intensified at upper levels in the western and central equatorial Pacific during periods of high ENSO amplitude. Furthermore, changes in equatorial stratification are connected with subsurface temperature anomalies originating from the central south tropical Pacific. The correlation analysis of ocean temperature anomalies against an index for the ENSO modulation supports the hypothesis of the existence of an oceanic “tunnel” that connects the south tropical Pacific to the equatorial wave guide. Further analysis of the wind stress projection coefficient onto the oceanic baroclinic modes suggests that the low frequency modulation of ENSO amplitude is associated with a significant contribution of higher-order modes in the western and central equatorial Pacific. In the light of these results, we suggest that, in the CGCM, change in the baroclinic mode energy distribution associated with low frequency ENSO amplitude modulation have its source in the central south tropical Pacific.  相似文献   

Interaction Between A Fast Plasmoid and the Solar Wind     

Wang  S.  Wang  X.Y.  Wu  C.S.  Li  Y.  Chao  J.K.  Yeh  T. 《Solar physics》2001,202(2):385-393

In this note a kinetic interaction process between a fast plasmoid ejected by the Sun, which represents another form of CME, and the background solar wind in the corona is discussed. We consider a system which consists of the plasmoid ions moving faster than the solar wind. We are interested in the time evolution of the ion distribution functions due to wave–particle interactions. Simulation results show that both perpendicular and parallel temperatures of the solar wind ions increase when the relative velocity between the plasmoid and the solar wind is sufficiently greater than the Alfvén velocity of the plasmoid ions. We suggest that this process is significant for the heating and acceleration of the solar wind in the low-heliographic latitude regions near the Sun.  相似文献