由等高线计算水域体积

介绍由等高线计算水域体积的过程,分析了水下地形的特点,对应注意的问题进行了讨论并给出了解决方案和算法。     

Comparison of theory and experiments for dispersion in homogeneous porous media

Modeling dispersion in homogeneous porous media with the convection–dispersion equation commonly requires computing effective transport coefficients. In this work, we investigate longitudinal and transverse dispersion coefficients arising from the method of volume averaging, for a variety of periodic, homogeneous porous media over a range of particle Péclet (Pe_p) numbers. Our objective is to validate the upscaled transverse dispersion coefficients and concentration profiles by comparison to experimental data reported in the literature, and to compare the upscaling approach to the more common approach of inverse modeling, which relies on fitting the dispersion coefficients to measured data. This work is unique in that the exact microscale geometry is available; thus, no simplifying assumptions regarding the geometry are required to predict the effective dispersion coefficients directly from theory. Transport of both an inert tracer and non-chemotactic bacteria is investigated for an experimental system that was designed to promote transverse dispersion. We highlight the occurrence of transverse dispersion coefficients that (1) depart from power-law behavior at relatively low Pe_p values and (2) are greater than their longitudinal counterparts for a specific range of Pe_p values. The upscaling theory provides values for the transverse dispersion coefficient that are within the 98% confidence interval of the values obtained from inverse modeling. The mean absolute error between experimental and upscaled concentration profiles was very similar to that between the experiments and inverse modeling. In all cases the mean absolute error did not exceed 12%. Overall, this work suggests that volume averaging can potentially be used as an alternative to inverse modeling for dispersion in homogeneous porous media.     

适于大规模数据的三维Kirchhoff 积分法体偏移实现方案

为适应实际生产中对大规模三维工区数据处理的效果及效率的要求,提出了按三维成像体输出成像结果的3D Kirchhoff积分法偏移实现方案.将地震数据按共偏移距道集形式排放,每个共偏移距数据的偏移类似于一个3D叠后Kirchhoff积分偏移,极大地降低了对计算机内存和局部盘及I/O通讯率的要求.每个地震道的成像(输出等时面)在由炮检点连线定义的旋转坐标系中进行,更好地考虑了偏移孔径计算及反假频处理.同时兼顾了超大规模地震数据PSTM成像处理中内存需求量、I/O通讯问题、并行处理方案及效率优化的细节问题.并行计算用偏移距号和每个共偏移距数据体中的线号作为一级和二级索引进行任务分解,更适应当前计算机集群中计算节点比较多的情况.最后考虑了在基本不影响效率的前提下的断点保护处理方案.理论及实际数据测试结果说明了该方案的可行性,与商业软件的对比验证了该方案的优越性.在此较完善的实现方案基础上,可以容易地把更优越的积分类偏移方法迅速推向实用化.     

Time series video analysis of bubble release processes at natural hydrocarbon seeps in the Northern Gulf of Mexico

This research quantifies the rate and volume of oil and gas released from two natural seep sites in the Gulf of Mexico: lease blocks GC600 (1200 m depth) and MC118 (850 m depth). Our objectives were to determine variability in release rates and bubble size at five individual vents and to investigate the effects of tidal fluctuations on bubble release. Observations with autonomous video cameras captured the formation of individual bubbles as they were released through partially exposed deposits of gas hydrate. Image processing techniques determined bubble type (oily, gaseous, and mixed: oily and gaseous), size distribution, release rate, and temporal variations (observation intervals ranged from 3 h to 26 d). A semi-automatic bubble counting algorithm was developed to analyze bubble count and release rates from video data. This method is suitable for discrete vents with small bubble streams commonly seen at seeps and is adaptable to multiple in situ set-ups. Two vents at GC600 (Birthday Candles 1 and Birthday Candles 2) were analyzed. They released oily bubbles with an average diameter of 5.0 mm at a rate of 4.7 bubbles s⁻¹, and 1.3 bubbles s⁻¹, respectively. Approximately 1 km away, within the GC600 seep site, two more vents (Mega Plume 1 and Mega Plume 2) were analyzed. These vents released a mixture of oily and gaseous bubbles with an average diameter of 3.9 mm at a rate of 49 bubbles s⁻¹, and 81 bubbles s⁻¹, respectively. The fifth vent at MC118 (Rudyville) released gaseous bubbles with an average diameter of 3.0 mm at a rate of 127 bubbles s⁻¹. Pressure records at Mega Plume and Rudyville showed a diurnal tidal cycle (24.5 h). Rudyville was the only vent that demonstrated any positive correlation (ρ = 0.60) to the 24.5 h diurnal tidal cycle. However, these observations were not conclusive regarding tidal effects on bubble release.     

Grain size,areal thickness distribution and controls on sedimentation of the 1991 Mount Pinatubo tephra layer in the South China Sea

The June 15, 1991 climactic eruption of Mt. Pinatubo produced an extensive, largely co-ignimbrite-derived airfall ash layer on Luzon Island and across the central South China Sea. The layer covers an area of ~4×10⁵ km² with a volume of 5.5 km³. Near the coast of Luzon, the deposit consists of two units: a normally graded basal ash bed, unimodal in grain size, and a finer-grained, internally structureless upper ash bed showing grain size bimodality. With increasing distance from the source, the coarse particle populations of the two units merge and migrate towards a near-constant fine population (~11 m); the distal region is covered by a fine-mode dominated, virtually ungraded single ash layer. The reversal of the winds from easterly directions at upper-tropospheric and stratospheric levels to westerly directions in the middle and lower troposphere indicates that both the coarse- and fine-mode components fell out from high-altitude eruption clouds. The high-velocity upper-level winds, however, would have transported fine-grained ash particles far beyond the South China Sea, which suggests that their settling was accelerated by aggregation. The boundary between the units thus marks a change from fallout of predominantly discrete pyroclasts to simultaneous fallout of aggregated fines and freely falling, coarse-grained particles. The particle populations composing the upper ash bed were almost completely removed from the proximal areas by the upper-level winds. At lower elevations, the counterclockwise circulation of a typhoon over the coastal area advected the ash south and eastward, producing a thickness maximum in the medial region (at about 160 km from source). The strong displacement of fines, possibly aided by wind turbulence, led to a break in bulk tephra thinning rates close to the coastline. In the distal region, outside the influence of the typhoon, southwest monsoonal winds caused a distinct lobe axis inflection and thickness asymmetry. Within this region, at about 420 km from source, fallout of particle aggregates created a second thickness maximum. Comparison of the field data with previous experimental observations and tephra flux records in the deep sea (Wiesner et al. 1995; Carey 1997; McCool 2002) implies that the transport of ash in the water column was largely determined by vertical density currents. Differences in the reaction of coarse and fine particles to turbulence in the descending plumes probably suppressed the segregation of fines but allowed the coarser pyroclasts to maintain their initial order of arrival at the sea surface. Considering typical fall rates of convective plumes, modifications of the initial fallout position of the particles by the South China Sea current system are on the order of only a few kilometers. The results suggest that convective sedimentation processes ensure the preservation of atmospheric particle transport directions, distances, and fallout modes in the deep sea.Editorial responsibility: R. CioniAn erratum to this article can be found at     

Molar volume and thermal expansion of glaucophane

The molar volume of glaucophane [Na₂Mg₃Al₂Si₈O₂₂(OH)₂] has been determined in this study by correcting synthetic glaucophane-rich amphiboles made in the system Na₂O–MgO–Al₂O₃–SiO₂–H₂O for very small deviations from ideal glaucophane composition using recent volume data on key amphibole components. The derived unit-cell volume for end-member glaucophane is 862.7±1.6 Å³, which gives a molar volume of 259.8±0.5 cm³/mol and a calculated density of 3.016±0.006 g/cm³. This value has been corroborated through an essentially independent method by correcting the volumes of natural sodic amphiboles reported in the literature for non-glaucophane components, particularly including calcium-rich components, to yield a value of 861.2±1.9 Å³. The unit-cell volume derived from the synthetic amphiboles, which is considered here to be more reliable, is somewhat smaller than that reported previously in the literature. A thermal expansion (α_V) at 298 K of 1.88±0.06×10^?5/K was derived from unit-cell volumes measured in the range of 25–500°C for a synthetic glaucophane sample, which is noticeably smaller than previously reported.     

北半球近地层典型区CO_2体积分数时空分布及成因

利用GEOS-Chem全球三维大气化学传输模式,分析了北半球近地层CO2体积分数的时空变化特征及其成因。2006—2010年的5 a的模拟结果表明:北半球中纬度近地层CO2体积分数存在着两个高值中心,即亚洲东部和北美东北部。在季节尺度上,亚洲东部CO2体积分数最大值出现在春季,而北美东北区域CO2体积分数最大值出现在冬季;而两个地区的CO2体积分数最低值都出现在夏季。在年际尺度上,两个区域CO2体积分数的年际变率增幅明显高于北半球其它区域,且CO2体积分数高值出现时间的年际差异较大。另外,模拟分析发现北半球森林、农田、草原典型区域,所对应的CO2体积分数具有不同的季节变化特点,它们的CO2季节内变幅依次减小。进一步分析发现3种不同典型区域的CO2体积分数与叶面积指数(LAI)季节变化,具有很好的负相关性。可见陆地生态系统作为碳汇,对近地层CO2体积分数的季节变化具有重要的作用。而温度和降水是影响LAI的最重要的两个气象因子,它们与CO2体积分数季节变化存在内在联系,模拟结果表明北半球大部分陆地近地层CO2体积分数与温度、降水呈现显著的负相关。     

气候变化对城市年径流总量控制率分区的影响

城市地区强降水发生频次和强度的增加容易诱发内涝现象,年径流总量控制率作为海绵城市的重要设计参数,更是直接受到降水变化的影响。以江苏省为例,利用全省70个国家级气象观测站1961—2019年最新的日降水量资料评估了气候变化对城市年径流总量控制率分区的影响。研究发现,有效降水的年代际变化十分明显,1991—2019年降水日数、降水量、降水强度均比其他时间段上更多、更强;太湖流域的设计雨量较小,连云港地区的设计雨量较大,南北差异随着控制率的提高而扩大,当控制率为85%时,全省设计雨量平均值为38.1 mm,最大值是最小值的1.7倍;气候变化对年径流总量控制率分区影响明显,江苏的苏南、江淮南部大部分地区的分区变大,导致全省IV区所占面积明显增加。不同等级降水的变化趋势是影响年径流控制率分区的关键因素,大雨以上的雨日、雨量在有效降水中占比增加,则分区变大。     

The Volume Calculation Method of Rock Collapses and Loess Landslides Triggered by the 1556 AD Huaxian M8 1/2 Earthquake in Shaanxi Province, China

Accurate volume calculation of each individual landslide triggered by strong historical earthquakes can help understand the characteristics of the typical earthquake-induced landslides, thus providing significant information for the modification of the focal parameters of historical earthquakes. In this study, we select one rock fall and three loess landslides triggered by the 1556 AD Huaxian M8 1/2 earthquake, compute their volumes using the low-altitude high-precision Unmanned Aerial Vehicle (UAV) photogrammetry and landslide profile restoration methods. The results show that:① the whole influencing area of the Huangjiagou Rock Fall is approximately 3.03×10⁵ m² and the area of the collapsed rock accumulated at the slope foot is 3.33×10⁴ m², accounting for approximately 10% of the entire influencing range. However, the estimated volume of the collapsed rock is only 0.699×10⁶ m³, indicating a rock fall with large influencing range but limited collapsed rock; ② the geological form of thethree loess landslides are preserved intactly, with volumes of 0.283×10⁸ m³, 0.074×10⁸ m³, and 0.377×10⁸ m³. These important geological hazard relics reflect the strong vibrations and severe casualties in the meizoseismal area; ③ loess landslides are the key reason of the serious death toll in the hilly-gully loess area. Our new method can be used to estimate the influencing area and the actual volume of each individual landslide, and rationally evaluate the role of earthquake landslides in the disaster. In addition, quantitative research on secondary disasters triggered by strong historical earthquakes is beneficial for understanding the surface process and focal parameters of the earthquakes.     

Redox and nonredox reactions of magnetite and hematite in rocks

Redox and nonredox reactions causing pseudomorphic replacement of hematite by magnetite and magnetite by hematite are compared.Pseudomorphic replacements resulting from redox reactions are known as martitization [replacement of magnetite by hematite due to oxidation; reaction (1)] and mushketovitization [replacement of hematite by magnetite due to reduction; reaction (2)]. These two replacements cause characteristic ore textures and volume changes (reaction (1): increase of 1.66%; reaction (2): decrease of 1.64%). These small volume changes are the reason that martitization and mushketovitization are widespread in many rocks under condition, however, that oxidizing or reducing fluids (solutions) are present.The same initial and end products may also be involved in nonredox reactions. Reaction (3) is the replacement of hematite by magnetite due to simple addition of Fe²⁺ atoms under basic conditions. This reaction causes an increase of the volume of 47.6%. Reaction (4), causing a volume decrease of 32.2%, is the replacement of magnetite by hematite due to leaching of Fe²⁺ atoms under acidic conditions. From these volume changes it is concluded that reaction (4) may occur in many rock types, whereas reaction (3) is restricted to unlithified sediments only. However, ore textures caused by nonredox reactions are unknown and therefore their occurrence in rocks is hypothetical.