Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations     

M. Gerritsen  J. V. Lambers 《Computational Geosciences》2008,12(2):193-208

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.  相似文献   

基于贡献权重迭加法的滑坡风险区划   总被引：7，自引：1，他引：6  

乔建平  石莉莉王萌 《地质通报》2008,27(11):1787-1794

贡献权重迭加法是对滑坡本底因子和承灾体因子在滑坡发育中的贡献率进行统计后,通过贡献率均值化、规一化处理,利用权重转换模型计算出每一个因子内部的权重——自权重W和因子相互之间的权重——互权重W＇。将滑坡因子、承灾体因子和两者的自权重、互权重分别相乘迭加,得到滑坡危险度和易损度区划结果,再将2项结果相乘得出滑坡风险区划。  相似文献   

青藏高原西北缘祁漫塔格山中新世快速抬升 的磷灰石裂变径迹证据   总被引：1，自引：0，他引：1  

拜永山  任二峰  范桂兰  许长青 《地质通报》2008,27(7):1044-1048

对祁漫塔格山体不同海拔高度所取的9个磷灰石样品的裂变径迹分析结果表明,东昆仑西段中新世早中期为主要的隆升期且隆升速率较高,早期隆升速率为111m／Ma,晚期隆升速率为98m／Ma,总体隆升速率为100m／Ma。样品显示出磷灰石裂变径迹长度大致分2类,一类磷灰石裂变径迹长度为（12．21±10．20）-（13．75±0．30）μm,径迹长度分布图基本上为窄而对称的正态分布,反映具有快的剥露冷却速率,未受到后期热事件的干扰。另一类磷灰石裂变径迹长度为（11．88±0．33）～（13．32±0．27）μm,较前一类具有稍慢的剥露冷却速率,并且受到了后期热事件的干扰。  相似文献   

Actinic flux and photolysis frequency comparison computations using the model PHOTOGT     

Carsten Blindauer  Vladimir Rozanov  John Philip Burrows 《Journal of Atmospheric Chemistry》1996,24(1):1-21

The accurate radiative transfer model GOMETRAN, initially designed to yield radiances at TOA in the wavelength range 240–790 nm, has been extended to allow for the computation of actinic fluxes down to 175 nm and for the calculation of photolysis frequencies in the atmosphere. The capability of the extended model PHOTOGT (PHOTOGOMETRAN) is demonstrated in a number of successful comparison studies both with recent experiments (ground-based, balloonborne, airborne) and model calculations of radiances, actinic fluxes and photolysis frequencies in the stratosphere and troposphere. In an atmospheric case study, the impact of new quantum yield data for the O₃ » O₂+O(¹ d) photodissociation channel on the photolytic production of O(¹ d) atoms in the lower atmosphere has been quantified.  相似文献   

关于海平面上升及其环境效应   总被引：12，自引：0，他引：12  

陈梦熊 《地学前缘》1996,3(2):133-140

沿海地区地质环境复杂。由于国民经济的迅速发展,人类活动日益增多,对地质环境进一步造成严重的不利影响。特别是“温室效应”引发的海平面上升,使各种地质灾害更加激化。文章详细探讨了海平面上升的影响因素,理论海平面与相对海平面升降幅度的评估,以及海平面上升所造成的地质环境效应与相应的防治对策,特别强调控制沿海城市地面沉降,以减轻海平面相对上升造成的危害,具有特殊而重要的意义。  相似文献   

Rupture process of large earthquakes in the northern Mexico subduction zone     

Larry J. Ruff  Angus D. Miller 《Pure and Applied Geophysics》1994,142(1):101-172

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.  相似文献   

About the specification of erosion flux for soft stratified cohesive sediments     

Georges Chapalain  Y. Peter Sheng  André T. Temperville 《Mathematical Geology》1994,26(6):651-676

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.  相似文献   

Palaeolimnological studies of the eutrophication of volcanic Lake Albano (Central Italy)     

Andrea Lami  Frank Niessen  Piero Guilizzoni  Julieta Masaferro  Claudio A. Belis 《Journal of Paleolimnology》1994,10(3):181-197

We use palaeolimnological techniques to reconstruct the eutrophication history of a volcanic lake (Lake Albano, central Italy) over the past three centuries. The presence of annual varves down to the bottom of the core (c. 1700 A.D.) indicated the lack of bioturbation and likely long-term meromixis. Sedimentation rates were estimated by varve counts (calcite/diatom couplets), indicating a mean rate of 0.15 cm yr^–1. The reconstruction of eutrophication was traced using past populations of algal and photosynthetic bacteria (through their fossil pigment), and geochemistry, as well as fossil remains of chironomids. Phaeophorbidea and the red carotenoid astaxanthin were used to detect past zooplankton development.The first sign of trophic change related to human activities is datedc. 1870 A.D. From that period onward a sharp increase of authigenic CaCO₃, nitrogen, N:P ratio, and dinoxanthin, a characteristic carotenoid of Chrysophyceae and Dinophyceae, is observed.Chironomid analyses showed the near absence of a deep water fauna throughout the core length. The populations of chironomid larvae are restricted to oxygenated littoral zones. In fact, the few fossil remains found are primarily of littoral origin, representing shallow water midges that were transported to profundal waters. The reduction of total chironomid in the uppermost layers of the core is to be related to human land uses.  相似文献   

河北省唐山地区丰台—野鸡坨断裂第四纪活动性——来自¹⁴C和磁性地层年代学的证据          下载免费PDF全文

陈宏强  专少鹏  陈超  赵华平  杨瑞  何娇月  胥勤勉 《中国地质》2021,48(2):605-617

丰台—野鸡坨断裂为唐山地区主要断裂之一，西侧为鸦鸿桥凹陷，东侧为唐山凸起，断裂两侧第四系厚度之差巨大。本文依据该断层两侧钻孔对其第四纪以来活动性进行初步的探讨。通过对丰台—野鸡坨断裂上下两盘PZK14和PZK20孔磁性地层学研究，并结合钻孔岩石地层，及浅部光释光和¹⁴C测年结果，建立第四纪地层格架。结果表明：两孔底部“泥包砾”为新近纪沉积；PZK14孔下更新统底界埋深为387 m，中更新统底界埋深为114 m，上更新统底界埋深为71 m，全新统底界埋深为6 m；PZK20孔下更新统底界埋深为155 m，中更新统底界埋深为73 m，上更新统底界为36 m，无全新世地层。丰台—野鸡坨断裂活动在早更新世时表现为逐渐增强的特点，活动速率由早期的5.4 cm/ka增加到13.9 cm/ka。中更新世断裂活动基本处于停滞状态，活动速率为1.0 cm/ka。晚更新世以后，断裂重新活动，且更加剧烈，活动速率达到了54.5 cm/ka。  相似文献   

豫西后造山阶段存在变质核杂岩吗？   总被引：2，自引：1，他引：1  

曹高社  赵太平 《地质论评》1997,43(4):365-372

不少学者用变质核杂岩模式解释豫西造山阶段基底隆起和与之相伴的盆岭构造的构造格局。介一本区基底隆起的形态、机制、时限和地壳剥露层次均没有变质核杂岩特征与基底隆起相伴的断陷人舅地也不同于变质核杂岩模式中的半地堑。  相似文献