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.     

Dynamical History of the Earth–Moon System

Differential equations describing the tidal evolution of the earth's rotation and of the lunar orbital motion are presented in a simple close form. The equations differ in form for orbits fixed to the terrestrial equator and for orbits with the nodes precessing along the ecliptic due to solar perturbations. Analytical considerations show that if the contemporary lunar orbit were equatorial the evolution would develop from an unstable geosynchronous orbit of the period about 4.42 h (in the past) to a stable geosynchronous orbit of the period about 44.8 days (in the future). It is also demonstrated that at the contemporary epoch the orbital plane of the fictitious equatorial moon would be unstable in the Liapunov's sense, being asymptotically stable at early stages of the evolution. Evolution of the currently near-ecliptical lunar orbit and of the terrestrial rotation is traced backward in time by numerical integration of the evolutional equations. It is confirmed that about 1.8 billion years ago a critical phase of the evolution took place when the equatorial inclination of the moon reached small values and the moon was in a near vicinity of the earth. Before the critical epoch t _cr two types of the evolution are possible, which at present cannot be unambiguously distinguished with the help of the purely dynamical considerations. In the scenario that seems to be the most realistic from the physical point of view, the evolution also has started from a geosynchronous equatorial lunar orbit of the period 4.19 h. At t < t _cr the lunar orbit has been fixed to the precessing terrestrial equator by strong perturbations from the earth's flattening and by tidal effects; at the critical epoch the solar perturbations begin to dominate and transfer the moon to its contemporary near-ecliptical orbit which evolves now to the stable geosynchronous state. Probably this scenario is in favour of the Darwin's hypothesis about originating the moon by its separation from the earth. Too much short time scale of the evolution in this model might be enlarged if the dissipative Q factor had somewhat larger values in the past than in the present epoch. Values of the length of day and the length of month, estimated from paleontological data, are confronted with the results of the developed model.     

天文(重力)水准的新手段及其计算方法

介绍当今确定天文大地垂偏差的新仪器－CCD自动天体测量仪和确定两点间高程异常差的新方法，该方法用单极坐标代替繁琐的双极坐标进行计算，讨论了天文重力水准的误差及垂线偏差非线性影响等问题，估计在不久的将来，用这一手段施测山区似大地水准面的精度可望达到厘米级。     

概论机场能见度的准确预报

通过对影响能见度的因素的分析，提出了一种能见度预报方法－综合分析法，极大地提高了能见度的预报准确率。     

High-pressure eclogites in northern Jiangsu – southern Shandong province, eastern China

Abstract Eclogites are distributed for more than 500 km along a major tectonic boundary between the Sino-Korean and Yangtze cratons in central and eastern China. These eclogites usually have high-P assemblages including omphacite + kyanite and/or coesite (or its pseudomorph), and form a high-P eclogite terrane. They occur as isolated lenses or blocks 10 cm to 300 m long in gneisses (Type I), serpentinized garnet peridotites (Type II) and marbles (Type III). Type I eclogites were formed by prograde metamorphism, and their primary metamorphic mineral assemblage consists mainly of garnet [pyrope (Prp) = 15–40 mol%], omphacite [jadeite (Jd) = 34–64 mol%], pargasitic amphibole, kyanite, phengitic muscovite, zoisite, an SiO₂ phase, apatite, rutile and zircon. Type II eclogites characteristically contain no SiO₂ phase, and are divided into prograde eclogites and mantle-derived eclogites. The prograde eclogites of Type II are petrographically similar to Type I eclogites. The mantle-derived eclogites have high MgO/(FeO + Fe₂O₃) and Cr₂O₃ compositions in bulk rock and minerals, and consist mainly of pyrope-rich garnet (Prp = 48–60 mol%), sodic augite (Jd = 10–27 mol%) and rutile. Type III eclogites have an unusual mineral assemblage of grossular-rich (Grs = 57 mol%) garnet + omphacite (Jd = 30–34 mol%) + pargasite + rutile. Pargasitic and taramitic amphiboles, calcic plagioclase (An₆₈), epidote, zoisite, K-feldspar and paragonite occur as inclusions in garnet and omphacite in the prograde eclogites. This suggests that the prograde eclogites were formed by recrystallization of epidote amphibolite and/or amphibolite facies rocks with near-isothermal compression reflecting crustal thickening during continent–continent collision of late Proterozoic age. Equilibrium conditions of the prograde eclogites range from P > 26 kbar and T= 500–750°C in the western part to P > 28 kbar and T= 810–880°C in the eastern part of the high-P eclogite terrane. The prograde eclogites in the eastern part are considered to have been derived from a deeper position than those in the western part. Subsequent reactions, manifested by (1) narrow rims of sodic plagioclase or paragonite on kyanite and (2) symplectites between omphacite and quartz are interpreted as an effect of near-isothermal decompression during the retrograde stage. The conditions at which symplectites re-equilibrated tend to increase from west (P < 10 kbar and T < 580°C) to east (P > 9 kbar and T > 680°C). Equilibrium temperatures of Type II mantle-derived eclogites and Type III eclogite are 730–750°C and 680°C, respectively.     

Late Quaternary landscapes in Central Australia: sedimentary history and palaeoecology of Puritjarra rock shelter

Puritjarra rock shelter provides a long record of late Quaternary vegetation in the Australian arid zone. Analysis of the sedimentary history of this rock shelter is combined with reanalysis of charcoal and phytolith records to provide a first‐order picture of changing landscapes in western Central Australia. These show a landscape responding to increasing aridity from 45 ka with deflation of clay‐rich red palaeosols (<45 ka) and sharp declines in grassland and other vegetation at 40–36 ka, and at the beginning of the Last Glacial Maximum (LGM) (24 ka). Vegetation in the catchment of the rock shelter recovered after 15 ka with expansion of both acacia woodland and spinifex grasslands, registering stronger summer rainfall in the interior of the continent. By 8.3 ka re‐vegetation of local palaeosols and dunes had choked off sediment supply to the rock shelter and the character of the sediments changed abruptly. Poaceae values peaked at 5.8 ka, suggesting the early–mid Holocene climatic optimum in Central Australia is bracketed between 8.3 and 5.8 ka. Local vegetation was disrupted in the late Holocene with a sharp decline in Poaceae at 3.8 ka, coinciding with an abrupt intensification of ENSO. Local grasslands recovered over the next two millennia and by 1.5 ka the modern vegetation appears to have become established. Copyright © 2009 John Wiley & Sons, Ltd.     

土地征用制度改革的探索

征地补偿安置是当前群众反映最强烈并迫切要求解决的问题。要解决该问题必须进行土地征用制度改革．由政府按照市场公平原则向农民集体征购土地，实行统一的土地征购制度。     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.