多次透射公式的一种稳定实现措施

基于对双曲偏微分议程数值解稳定性的GKS（Gustafsson,Kreiss and Sundstrom）准则的物理解释，本文提出了在波动模拟中多次透射公式飘移失稳的一个简单措施，并用波源问题和散射问题的数值试验检验了建议的稳定措施的有效性。     

西藏盐湖卤水蒸发速率的实验与计算

本文对国内外水面蒸发速率的研究进行了综述,总结了西藏扎布耶盐湖Φ20cm蒸发皿淡水蒸发量与气温、降水、日照的相关关系,提出了改进的扩展彭曼公式法,用于较为准确地计算盐湖卤水蒸发速率,并以西藏扎布耶盐湖为例计算了盐湖卤水蒸发。该方法可以应用于盐湖湖面蒸发与水量均衡计算,也可以应用于盐湖开发中的盐田工艺设计计算与实际生产应用。     

Variability of tropical cyclone in high frequent occurrence regions over the western North Pacific

In this study, three high frequent occurrence regions of tropical cyclones(TCs), i.e., the northern South China Sea(the region S), the south Philippine Sea(the region P) and the region east of Taiwan Island(the region E), are defined with frequency of TC's occurrence at each grid for a 45-year period(1965–2009), where the frequency of occurrence(FO) of TCs is triple the mean value of the whole western North Pacific. Over the region S, there are decreasing trends in the FO of TCs, the number of TCs' tracks going though this region and the number of TCs' genesis in this region. Over the region P, the FO and tracks demonstrate decadal variation with periods of 10–12 year, while over the region E, a significant 4–5 years' oscillation appears in both FO and tracks. It is demonstrated that the differences of TCs' variation in these three different regions are mainly caused by the variation of the Western Pacific Subtropical High(WPSH) at different time scales. The westward shift of WPSH is responsible for the northwesterly anomaly over the region S which inhibits westward TC movement into the region S. On the decadal timescale, the WPSH stretches northwestward because of the anomalous anticyclone over the northwestern part of the region P, and steers more TCs reaching the region P in the greater FO years of the region P. The retreating of the WPSH on the interannual time scale is the main reason for the FO's oscillation over the region E.     

Massively parallel spatial point pattern analysis: Ripley’s K function accelerated using graphics processing units

This study presents a massively parallel spatial computing approach that uses general-purpose graphics processing units (GPUs) to accelerate Ripley’s K function for univariate spatial point pattern analysis. Ripley’s K function is a representative spatial point pattern analysis approach that allows for quantitatively evaluating the spatial dispersion characteristics of point patterns. However, considerable computation is often required when analyzing large spatial data using Ripley’s K function. In this study, we developed a massively parallel approach of Ripley’s K function for accelerating spatial point pattern analysis. GPUs serve as a massively parallel platform that is built on many-core architecture for speeding up Ripley’s K function. Variable-grained domain decomposition and thread-level synchronization based on shared memory are parallel strategies designed to exploit concurrency in the spatial algorithm of Ripley’s K function for efficient parallelization. Experimental results demonstrate that substantial acceleration is obtained for Ripley’s K function parallelized within GPU environments.     

INFLUENCE OF SERIES OF SQUARE GRIDS ON FRACTAL DIMENSIONS—A Case Study of Mountains of China＇s Mainland

MANDELBROT enunciated the uncertainty of the length of a coastline in his paper “How long is the coastline of Britain?” published in Science in 1967. The fractal concept was presented for the first time in that paper and has been applied to many fields ever since. Although fractal dimensions of lots of phenomena were calculated by the box-counting method, the quantitative influence of series of square grids on them is ignored. The issue is systematically discussed as a case study of the mountains of China‘s Mainland in this paper. And some significant conclusions are drawn as follows: 1) Although the fractal character objectively exists in the mountains of China‘s Mainland, and it does not vary with the changes of series of square grids, the fractal dimensions of the mountains of China‘s Mainland are different with these changes. 2) The fractal dimensions of the mountains of China‘s Mainlandvary with the average lengths of sides of series of square grids. The fractal dimension of the mountains of China‘s Mainland is the function of the average length of side of square grid. They conform to the formula D=f(r) (where D is the fractal dimension, and r is the average length of side of square grid). 3) Different dots of data collection can affect the fractal dimension of the mountains of China‘s Mainland. 4) The same range of length of side of square grid and dots of data collection can ensure the comparison of fractal dimensions of the mountains of China‘ s Mainland. The research is helpful to get the more understanding of fractal and fractal dimension, and ensure that the fractal studies would be scientific.     

有机酸对长石溶解度影响的热力学研究

为明确认识有机酸（草酸、醋酸）对长石溶解度及储层次生孔隙发育的影响，本文利用热力学方法，计算了不同温度、酸度下钾长石在含有醋酸根0.05mol/L、草酸0.0055mol/L的溶液中的溶解度及草酸和醋酸对溶液酸度的影响，计算的结果表明醒酸根与铝的络合作用只发生在偏酸性的溶液体系中，并且随温度升高要求酸度增加，实际的络合能力降低，醋酸根的络合作用对长石溶解度的增加作用极及有限，最高时不足一倍；但醋酸能有效地缓冲溶液酸度，从而提高长石溶解度。草酸根与铝的络合作用发生在弱碱性至酸性条件下，温度升高，络合能力增加，络合效应对长石溶解度，草酸根上铝的络合作用发生在弱碱性至酸性条件下，温度升高，络合能力增加，络合效应对长石溶出量的增加最多可达2个数量级以上，草酸的浓度太低无法通过缓冲溶液在较高的酸度下来促进长石的溶解，草酸和醒酸能否促进储层次生孔隙的发育取决于盆地流体的酸度及组成流体的缓冲体系。     

Bimodal magmatism of the Tucumã area,Carajás province: U-Pb geochronology,classification and processes

Geological mapping of the Tucumã area has enabled the identification of dike swarms intruded into an Archean basement. The disposition of these dikes is consistent with the well-defined NW-SE trending regional faults, and they can reach up to 20 km in length. They were divided into three main groups: (i) felsic dikes (70% of the dikes), composed exclusively of porphyritic rhyolite with euhedral phenocrysts of quartz and feldspars immersed in an aphyric felsite matrix; (ii) mafic dikes, with restricted occurrence, composed of basaltic andesite and subordinate basalt, with a mineralogical assembly consisting dominantly of plagioclase, clinopyroxene, orthopyroxene and olivine; and (iii) intermediate rocks, represented by andesite and dacite. Dacites are found in outcrops associated with felsic dikes, representing different degrees of hybridization or mixture of mafic and felsic magmas. This is evidenced by a large number of mafic enclaves in the felsic dikes and the frequent presence of embayment textures. SHRIMP U-Pb zircon dating of felsic dikes yielded an age of 1880.9 ± 3.3 Ma. The felsic dikes are peraluminous to slightly metaluminous and akin to A2, ferroan and reduced granites. The intermediate and mafic dikes are metaluminous and belong to the tholeiitic series. Geochemical modeling showed that mafic rocks evolved by pyroxene and plagioclase crystallization, while K-feldspar and biotite are the fractionate phases in felsic magma. A simple binary mixture model was used to determine the origin of intermediate rocks. It indicated that mixing 60% of rhyolite and 40% basaltic andesite melts could have generated the dacitic composition, while the andesite liquid could be produced by mixing of 60% and 40% basaltic andesite and rhyolite melts, respectively. The mixing of basaltic and andesitic magmas probably occurred during ascent and storage in the crust, where andesite dikes are likely produced by a more homogeneous mixture at high depths in the continental crust (mixing), while dacite dikes can be generated in the upper crust at a lower temperature, providing a less efficient mixing process (mingling). The affinities observed between the felsic to intermediate rocks of the Rio Maria and São Felix do Xingu areas and the bimodal magmatism of the Tucumã area reinforce the hypothesis that in the Paleoproterozoic the Carajás province was affected by processes involving thermal perturbations in the upper mantle, mafic underplating, and associated crustal extension or transtension. The 1.88 Ga fissure-controlled A-type magmatism of the Tucumã area was emplaced ∼1.0 to ∼0.65 Ga after stabilization of the Archean crust. Its origin is not related to subduction processes but to the disruption of the supercontinent at the end of the Paleoproterozoic.     

竖波钢板组合剪力墙震损修复研究

钢板组合剪力墙在强震作用下受损严重,震损修复问题较为突出。现设计一片以方钢管作为边缘约束构件的竖波钢板组合剪力墙,在最大弹塑性层间位移角时的震损情况进行分析:竖波钢板组合剪力墙边缘约束方钢管刚度较大,与墙体不匹配,最终发生方钢管从地梁拔动的脆性破坏模式。对竖波钢板组合剪力墙进行震损修复,得到墙趾可更换竖波钢板组合剪力墙,然后进行拟静力试验。试验结果表明:（1）在震损后的竖波钢板组合剪力墙的墙趾处安装阻尼器,其抗震性能基本得到恢复。（2）修复后的竖波钢板组合剪力墙承载力较修复前降低了23%,但延性和耗能能力显著提升。（3）改变了竖波钢板组合剪力墙的破坏模式,由脆性破坏转化为延性破坏。借助有限元软件详细讨论了内嵌竖向波形钢板厚度和波角、轴压比、阻尼器腹板数量对修复后竖波钢板组合剪力墙抗剪承载力的影响,结果表明:改变内嵌竖向波形钢板厚度对试件抗剪承载力影响较大,而改变内嵌竖向波形钢板波角、轴压比和阻尼器腹板数量对试件抗剪承载力影响较小。结合有限元算例,提出修复后竖波钢板组合剪力墙抗剪承载力的计算公式,可为工程实际提供参考。     

Yabuki&Matsu’ura方法在汶川M_W7.9地震反演中的应用

2008年5月12日,青藏高原东缘龙门山断裂带发生汶川MW7.9地震,该地震使得北川—映秀断裂、灌县—江油断裂发生了同震破裂。本文主要利用震后通过复测获得的GPS同震形变场,采用Yabuki&Matsu’ura反演计算方法和分段平面断层模型,反演了地震同震滑动分布。结果表明:映秀—北川主破裂带的断层错动,在映秀附近以逆冲滑动为主,而在北川以北,其走滑运动明显大于逆冲,这一结果与震后地质调查结果与通过地震波研究获得的断层破裂特征相一致;反演得到的最大滑动量达到9.3m和9.6m,分别对应于这次地震中地表破坏最为严重的北川和映秀地区;由所获得的滑动分布计算的地震矩为8.07×1020 N.m,对应的震级为MW7.9。研究结果初步显示,Yabuki&Matsu’ura反演方法可适用在内陆地震断层反演计算中。