邹平王家庄铜矿床成矿地球化学及成因探讨

王家庄铜矿床的矿化脉石英中流体包裹体均一温度介于116 ~ 566 ℃之间,均值为 289 ℃;盐度介于7.2%～63.2% NaCleq,均值为21.1% NaCleq。流体的气相成分主要为H2O和CO2。在均一温度为240 ~ 440 ℃区间内,出现了富气相的两相水溶液包裹体、富液相的两相水溶液包裹体和含子晶的三相水溶液包裹体共存现象,以及加温后富气相包裹体均一到气相和同期富液相包裹体均一到液相的特征,表明成矿流体曾发生过沸腾作用;其中第一次发生于360 ~ 400 ℃,主要形成高温、高盐度含子晶的三相水溶液包裹体和高温、中盐度富液相的两相水溶液包裹体及高温、低盐度富气相的两相水溶液包裹体;第二次发生于240 ~ 320 ℃,主要形成高—中温、高盐度的含子晶的三相水溶液包裹体和高—中温、中盐度富液相的两相水溶液包裹体及高—中温、低盐度富气相的两相水溶液包裹体;之后主要形成富液相的两相水溶液包裹体,具有中低温和中盐度的特征。矿化脉石英中的δ18OH2O介于-1.14‰ ~ 1.79‰之间,均值为0.94‰;δD介于-63.70‰ ~ -56.50‰之间,均值为-59.8‰;说明王家庄铜矿床的成矿流体主要来源于岩浆,后期混入大气降水。矿石矿物的δ34S变化于-8.80‰ ~ -2.80‰之间,均值为-6.33‰。矿石矿物的n(206Pb)/n(204Pb)介于18.1684 ~ 18.3637之间,均值为18.2892;n(207Pb)/n(204Pb)介于155546 ~ 156342之间,均值为155777;n(208Pb)/n(204Pb)介于38.1286 ~ 38.4840之间,均值为38.2780。说明矿石具有贫重硫和富放射性成因铅的特征,硫、铅主要来源于深部,后期可能受到地壳物质或大气降水的混染。     

河南洛宁段河金矿流体包裹体研究和矿床成因

河南省洛宁县段河石英脉型金矿主要包括石寨沟和岭东两个矿区,分别由3～4条含金石英脉构成。矿化过程从早到晚包括石英-黄铁矿、石英-多金属硫化物和石英-碳酸盐等3个阶段．其中中阶段金矿化最强,次为早阶段。各阶段石英中流体包裹体以气液两相包裹体为主．次为纯液体包裹体。激光拉曼测试表明,气液两相包裹体的液相为H2O,气相主要为Ho和CO2混合、纯H2O,次为纯CO2;纯液体包裹体为纯H2O。石寨沟矿区包裹体均一温度从早到晚依次为240．9～315．9℃．188．7～304．5℃,137．3～259．3℃：流体盐度变化依次为（6．74～12．85）wt%NaCl．eq,（2．41～8．68）wt％NaCl．eq,（2．24-7．86）wt%NaCl．eq。岭东矿区均一温度从早到晚依次为303．7-343．1℃,251．8-325．4℃,305．7～355．0℃：流体盐度变化依次为（5．11～11．70）wt％NaCl．eq,（2．74-10．11）wt％NaCl．eq,（0．53-6．74）wt％NaCl．eq。两矿区主成矿期流体均为中温、低盐度,早阶段流体为改造热液和变质热液的混合体,含一定量CO2,且流体CO2含量和盐度从早到晚逐渐降低。石寨沟矿区包裹体均一温度逐渐降低,而岭东矿区包裹体均一温度先降后升,加之岭东矿区各阶段成矿温度均高于石寨沟矿区．表明成矿流体系统主要受岩浆热驱动,岭东矿区更靠近岩体,且在晚阶段又有脉动性的岩浆加热．段河金矿区南部存在隐伏岩体。     

无外加流体、350℃和差异应力条件下硫化物再活化实验研究

红透山块状硫化物矿石主要成分为黄铁矿、磁黄铁矿、黄铜矿、闪锌矿和石英、角闪石、黑云母等脉石矿物。将此矿石烘干后作为试料置于岩石三轴应力试验机,在13h内将轴压、围压和温度分步升至1276MPa、414MPa和350℃,然后在空气中自然冷却至室温。实验产物中黄铁矿、石英、长石等以脆性碎裂为主,而磁黄铁矿、黄铜矿、闪锌矿和云母等以塑性变形为主。再活化的黄铜矿、磁黄铁矿和少量闪锌矿呈脉和网脉穿插黄铁矿碎斑。实验结果表明,即使没有外来流体加入．构造动力作用所导致的流体包裹体破坏所释出的流体,就足以使成矿物质发生活化转移和再活化。     

Numerical Experiment on Two-Dimensional Line Thermal

The time evolution of a two-dimensional line thermal-a turbulent flow produced by an initial element with signifi-cant buoyancy released in a large water body, is numerically studied with the two-equation k - e model for turbulence closure. The numerical results show that the thermal is characterized by a vortex pair flow and a kidney shaped concentra-tion structure with double peak maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time around 10. There are two regions in the flow field of a line thermal: a mixing region where the concentration of tracer fluid is high and the flow is turbulent and rotational with a pair of vortex eyes, and an ambient region where the concentration is zero and the flow is potential and well-described by a model of doublet with strength very close to those given by early experimental and analytical studies. The added virtual mass coeffi-cient of the thermal motion is found to be approximat     

Numerical study of the motions and drift force of a floating OWC device

The motion and the drift force of a floating OWC (oscillating water column) wave energy device in regular waves are studied taking account of the oscillating surface-pressure due to the pressure drop across the duct of the air chamber. The potential problem inside the chamber is formulated by making use of the Green integral equation associated with the Rankine-type Green function while the outer problem with the Kelvin-type Green function. The added mass, wave damping and excitation coefficients as well as the motion and drift force of the OWC device are calculated for various values of parameter related to the pressure drop.     

Geophysical evidence of mud diapirism on the Mediterranean Ridge accretionary complex

Mud volcanoes, mud cones, and mud ridges have been identified on the inner portion of the crestal area, and possibly on the inner escarpment, of the Mediterranean Ridge accretionary complex. Four areas containing one or more mud diapirs have been investigated through bathymetric profiling, single channel seismic reflection profiling, heat flow measurements, and coring. A sequence of events is identified in the evolution of the mud diapirs: initially the expulsion on the seafloor of gasrich mud produces a seafloor depression outlined in the seismic record by downward dip of the host sediment reflectors towards the mud conduit; subsequent eruptions of fluid mud may create a flat topped mud volcano with step-like profile; finally, the intrusion of viscous mud produces a mud cone.The origin of the diapirs is deep within the Mediterranean Ridge. Although a minimum depth of about 400 m below the seafloor has been computed from the hydrostatic balance between the diapiric sediments and the host sediments, a maximum depth, suggested by geometric considerations, ranges between 5.3 and 7 km. The presence of thermogenic gas in the diapiric sediments suggests a better constrained origin depth of at least 2.2 km.The heat flow measured within the Olimpi mud diapir field and along a transect orthogonal to the diapiric field is low, ranging between 16 ± 5 and 41 ± 6 mW m^–2. Due to the presence of gas, the thermal conductivity of the diapiric sediments is lower than that of the host hemipelagic oozes (0.6–0.9 and 1.0–1.15 W m^–1 K^–1 respectively).We consider the distribution of mud diapirs to be controlled by the presence of tectonic features such as reverse faults or thrusts (inner escarpment) that develop where the thickness of the Late Miocene evaporites appears to be minimum. An upward migration through time of the position of the décollement within the stratigraphic column from the Upper Oligocene (diapiric sediments) to the Upper Miocene (present position) is identified.     

Validation of a finite difference method for the simulation of vortex-induced vibrations on a circular cylinder

The Karman Vortex Street generated by a circular cylinder is investigated by the numerical solution of the compressible Navier–Stokes equations in the incompressible Mach number range (Mach<0.3) using the Beam and Warming implicit scheme. The agreement with the fully incompressible projection method (Chorin, 1968) is fairly good while convergence time is very much better. The investigation suggests that the compressible Navier–Stokes equations may be used as an efficient alternative to study incompressible flows as well. Mach numbers just below 0.3 are enough to simulate incompressible flow behavior and at the same time do not cause numerical ill-conditioning in the solution. In addition, some relevant features of the vortices generated and carried by the wake of the cylinder could be fairly well captured.     

Numerical Solution of Lock-Release Gravity Current with Viscous Self-Similar Regime

Lock-release gravity currents with a viscous self-similar regime are simulated by use of the renormalization group(RNG) k - ε model for Reynolds-stress closure. Besides the turbulent regime with initially a slumping phase of a conslant current front speed and later an inviseid self-similar phase of front speed decreasing as t^-1/3(where t is the time measured from release), the viseous self-similar regime is satisfactorily reproduced with front speed decreasing as t^-4/5,consistent with well known experimental observations.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.