长江口南支潮流的一个数值模拟

为了能够精确地拟合天然水域不规则的水陆边界,并能保持有限差分数值计算在矩形网格中进行,本文利用Thompson的数值网格生成技术,在长江口南支七丫口至横沙河段设计了一个椭圆型边界拟合坐标系,并在这个坐标系下建立了平面二维潮流数学模型。通过数值求解此数学模型,实现了计算区域内现状潮流的一个数值模拟,模拟结果表明,模式的设计是成功的。     

空间地理信息组织表达方法探讨

随着空间地理信息以及各种新型地图产品在各个领域的广泛应用,人们发现不同时期,应用需求对空间地理信息组织表达的要求是不同的,不同领域和部门对地理信息的组织表达方式也是不同的。但是,目前空间地理信息的共享和互操作、空间地理信息的多尺度分析与表达等,都要求空间地理信息的组织与表达要科学、合理和规范,否则难以满足测绘信息服务化及信息化测绘体系的建设。本文以交通类信息组织表达为例,详尽阐述当前空间地理信息组织表达的现状,分析其原因,并提出了解决这一问题的思路和方法,为今后信息化测绘体系建设、数据库数据组织和地理信息表达,提供理论和方法借鉴。     

Partial melting in the Inzie Head gneisses: the role of water and a petrogenetic grid in KFMASH applicable to anatectic pelitic migmatites

A sequence of prograde isograds is recognized within the Dalradian Inzie Head gneisses where pelitic compositions have undergone variable degrees of partial melting via incongruent melting reactions consuming biotite. Three leucosome types are identified. At the lowest grades, granitic leucosomes containing porphyroblasts of cordierite (CRD‐melt) are abundant. At intermediate grades, CRD‐melt mingles with garnetiferous leucosomes (GT‐melt). At the highest grades, CRD‐melt coexists with orthopyroxene‐bearing leucosomes (OPX‐melt), while garnet is conspicuously absent. The prograde metamorphic field gradient is constrained to pressures of 2–3 kbar below the CRD‐melt isograd, and no greater than 4.5 kbar at the highest grade around Inzie Head. A petrogenetic grid, calculated using thermocalc , is presented for the K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (KFMASH) system for the phases orthopyroxene, garnet, cordierite, biotite, sillimanite, H₂O and melt with quartz and K‐feldspar in excess. For the implied field gradient, the reaction sequence predicted by the grid is consistent with the successive prograde development of each leucosome type. Compatibility diagrams suggest that, as anatexis proceeded, bulk compositions may have been displaced towards higher MgO content by the removal of (relatively) ferroan granitic leucosome. An isobaric (P = 4 kbar) T–a_H2O diagram shows that premigmatization fluids must have been water‐rich (a_H2O > 0.85) and suggests that, following the formation of small volumes of CRD‐melt, the system became fluid‐absent and melting reactions buffered a_H2O to lower values as temperatures rose. GT‐ and OPX‐melt formed by fluid‐absent melting reactions, but a maximum of 7–11% CRD‐melt fraction can be generated under fluid‐absent conditions, much less than the large volumes observed in the field. There is strong evidence that the CRD‐melt leucosomes could not have been derived by buoyantly aided upwards migration from levels beneath the migmatites. Their formation therefore required a significant influx of H₂O‐rich fluid, but in a quantity insufficient to have exhausted the buffering capacity of the solid assemblage plus melt. Fluid : rock ratios cannot have exceeded 1 : 30. The fluid was channelled through a regionally extensive shear zone network following melt‐induced failure. Such an influx of fluid at such depths has obvious consequences for localized crustal magma production and possibly for cordierite‐bearing granitoids in general.     

自适应非线性有限元的多重网格法求解及应用

描述了非线性自适应有限元的多重网格法求解过程。 通过算例说明了方法的有效性以及自适应多重网格有限元用于开挖施工过程模拟分析的优越性。     

广州城市通达性综合评价及空间差异

广州市在2000 年提出“东进、西联、南拓、北优”的城市空间发展战略,2010 年初步形成“两心四城”城市空间格局,随后在发展战略中加入“中调”,以进一步优化城市发展布局,谋求在2020 年形成多中心、组团式、网络型的城市空间架构。城市空间发展战略调整有赖于城市交通基础设施建设的支撑,同时交通基础设施的完备程度与密度也是进行城市通达性差异研究的重要基础。近年来,广州的城市交通通达性发生了深刻变化,对城市空间结构布局产生了深远影响。在此背景下,借助地理信息系统,运用网格分区法对广州市交通基础设施进行区域划分和密度计算;进一步结合快速轨道交通（地铁）网络的时间通达性度量,对广州市的城区通达性进行综合评价,并分析其空间特征。研究结果表明,广州市城区的通达性呈现同心环状的空间分布格局,区域差异性较大,城市空间布局有待进一步完善和优化。在此基础上,基于地铁网络时间通达性和区域交通基础设施密度值进行区位等级划分,探讨区位因素与交通通达性的相互关系;并结合城市公共、商业服务设施和住宅区的空间特征差异,进一步甄别广州城市通达性空间格局特征和区域差异,从而为城市建设和空间发展提供科学参考。     

利用三维高斯射线束成像进行地震定位

常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低. 本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题. 通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取. 提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率. 利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明:基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题.     

与找矿信息量法结合的含矿网格单元蒙特卡洛矿产资源潜力评价方法的原理及应用

文章介绍了以含矿网格单元为统计单位、并与找矿信息量法结合的蒙特卡洛矿产资源潜力评价方法的原理和工作步骤,以及用该法进行丹池锡多金属成矿带锡矿矿产资源潜力评价工作案例,表明该法可以把评价工作区的矿产空间分布信息、成矿模式和找矿模型、成矿-找矿空间信息引入蒙特卡洛方法中,无需人类专家估计矿床个数,减少了工作步骤,提高了工作效率和自动化程度,减少了因人为主观性而造成的评价误差。     

栅格化属性精度损失的评估方法及其尺度效应分析——以四川省1:25 万土地覆被数据为例

选择在600 m~30 km 16 个尺度上,在ArcGIS 中利用常用的面积最大值法(Rule ofMaximum Area,RMA) 对2005 年四川省1:25 万土地覆被矢量数据进行栅格化,并采用两种属性精度损失评估方法：传统的常规分析方法和一种新的基于栅格单元分析方法,来对比分析在这两种评估方法下RMA栅格化的属性(这里是指面积) 精度损失随尺度的变化特征。结果表明：(1) 在同一尺度下采用基于栅格单元方法分析所得的研究区平均属性精度损失大于常规分析方法分析得到的平均属性精度损失,且二者之间的差异在1~10 km内很明显,当栅格单元大于10km时,两种方法得到的平均属性精度损失的差值稳定,且其随尺度的变化曲线趋于平行;(2) 基于栅格单元分析方法不仅能够准确地定量估计RMA栅格化的属性精度损失,而且能客观地反映属性精度损失的空间分布规律;(3) 对四川省1:25 万土地覆被数据进行面积最大值法(RMA)栅格化的适宜尺度域最好不要超过800 m,在该尺度域内数据工作量适宜,且RMA栅格化属性精度损失小于2.5%。     

Parabolic Wave Propagation Modelling in Orthogonal Coordinate Systems

This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild-slope equation in the orthogonal coordinate system. Two types of coordinate systems are studied: (a) a general form of orthogonal coordinate system and (b) the conformal system, a special form of orthogonal coordinate system. Two typical examples, namely, expanded breakwaters and a circular channel, are studied to validate the model. First, the examples are studied by use of the general orthogonal coordinates. Then the same examples are computed by use of the confonnal system. The computational results show that the confonnal coordinate system generally gives better predictions than the general orthogonal system. A numerical technique for generating the conformal grid is combined with the numerical model to improve the practicability of the model. The comparison between the result from the numerical grid system and that from the analytical grid system shows that reliable computational results can be obtained by use of the numerical confonnal grid system.     

网格嵌套技术对一次中尺度对流系统降水过程模拟的影响

利用非静力中尺度模式MM5对2002年7月22日12:00~23日12:00(世界时,下同)长江流域的一次梅雨锋暴雨过程进行数值模拟试验,主要讨论了网格嵌套技术对降水和中尺度对流系统的影响。结果表明:三重嵌套在D1,D2域选用积云参数化方案后,模拟的雨区收缩,虚假降水中心相对减少,降水强度及分布更接近观测值。在模式非线性动力、热力及湿物理过程共同驱动下,通过嵌套网格的双向相互作用,使可分辨云尺度的细网格域D3将其信息通过嵌套边界向次细网格域D2传递,然后再通过D2域边界向粗网格域D1域传递。同样,动力、热力反馈也会反向进行。结果将有助于改进各网格域的预报效果。但对D1网格域系统位置及其发展演变过程的影响相对小些;另外,通过双向多重嵌套,可提高模式预报区域的分辨率,特别是提高模式关键预报区域的分辨率,这也就有可能改进预报的中尺度物理场,使其能够较真实地描写大气实况。