Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

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

河南省洛宁县段河石英脉型金矿主要包括石寨沟和岭东两个矿区,分别由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℃,然后在空气中自然冷却至室温。实验产物中黄铁矿、石英、长石等以脆性碎裂为主,而磁黄铁矿、黄铜矿、闪锌矿和云母等以塑性变形为主。再活化的黄铜矿、磁黄铁矿和少量闪锌矿呈脉和网脉穿插黄铁矿碎斑。实验结果表明,即使没有外来流体加入．构造动力作用所导致的流体包裹体破坏所释出的流体,就足以使成矿物质发生活化转移和再活化。     

预防油页岩地层钻进中压差卡钻的措施

分析了油页岩地层钻进中压差卡钻发生的原因,提出了预防压差卡钻的途径及处理方法,并通过工程实例来强调预防为主的重要性。     

气举反循环采矿方法单井试验及应用探讨

湖南道县一个矿区,铁锰矿体储量约2200万t。由于复杂的水文地质和工程地质条件,从20世纪70年代发现后,一直没有开采。笔者通过多年大直径工程井的实践经验,查阅国内外资料后,提出了钻井-气举反循环法试开采的方法,并在现场进行了单孔开采试验,获得了比较满意的效果。以该矿体Ⅱ号矿为例,介绍了钻井-气举反循环法采矿的适用条件、工艺过程以及试采情况,并就以后大面积开采提出了若干技术措施。     

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.