人流联系和经济联系视角下区域城市关联比较——基于手机信令数据和企业关联数据的研究

手机信令数据所代表的人流联系与国内外研究常用的企业关联数据在测度城市关联存在何种差异,尚缺乏明确的研究。本文以江西省北部地区为研究对象,分别以手机信令数据、企业关联数据测度城市关联,比较了两种城市关联与城市之间空间距离、人口规模的关系,比较了两种城市关联网络的层级、结构、腹地。主要研究发现：人流联系不仅集中在城市周边地区,而且集中于较高等级城市之间;经济联系则主要集中在较高等级城市之间。人流联系、经济联系受到城市人口规模影响均较小。人流联系显著地符合幂函数空间衰减规律。手机信令数据和企业关联数据是两种不同的城市关联数据,但可互为补充,较为全面地反映城市关联。     

Multiphase flow simulation through porous media with explicitly resolved fractures

Accurate simulation of multiphase flow in fractured porous media remains a challenge. An important problem is the representation of the discontinuous or near discontinuous behaviour of saturation in real geological formations. In the classical continuum approach, a refined mesh is required at the interface between fracture and porous media to capture the steep gradients in saturation and saturation‐dependent transport properties. This dramatically increases the computational load when large numbers of fractures are present in the numerical model. A discontinuous finite element method is reported here to model flow in fractured porous media. The governing multiphase porous media flow equations are solved in the adaptive mesh computational fluid dynamics code IC‐FERST on unstructured meshes. The method is based on a mixed control volume – discontinuous finite element formulation. This is combined with the P_N+1DG‐P_NDG element pair, which has discontinuous (order N+1) representation for velocity and discontinuous (order N) representation for pressure. A number of test cases are used to evaluate the method's ability to model fracture flow. The first is used to verify the performance of the element pair on structured and unstructured meshes of different resolution. Multiphase flow is then modelled in a range of idealised and simple fracture patterns. Solutions with sharp saturation fronts and computational economy in terms of mesh size are illustrated.     

Hybrid finite element–finite volume discretization of complex geologic structures and a new simulation workflow demonstrated on fractured rocks

The generation of computational meshes of complex geological objects is a challenge: their shape needs to be retained, resolution has to adapt to local detail, and variations of material properties in the objects have to be represented. Also mesh refinement and adaptation must be sufficient to resolve variations in the computed variable(s). Here, we present an unstructured hybrid finite element, node‐centred finite‐volume discretization suitable for solving fluid flow, reactive transport, and mechanical partial differential equations on a complex geometry with inhomogeneous material domains. We show that resulting meshes accurately capture free‐form material interfaces as defined by non‐uniform rational B‐spline curves and surfaces. The mesh discretization error is analysed for the elliptic pressure equation and an error metric is introduced to guide mesh refinement. Finite elements and finite volumes are represented in parametric space and integrations are conducted numerically. Subsequently, integral properties are mapped to physical space using Jacobian transformations. This method even retains its validity when the mesh is deformed. The resulting generic formulation is demonstrated for a transport calculation performed on a complex discrete fracture model.     

Disciplined Mobility and the Emotional Subject in Royal Dutch Lloyd's Early Twentieth Century Passenger Shipping Network

This paper examines the disciplined mobility and emotional geographies of “between‐deck” passengers in Royal Dutch Lloyd's early Twentieth Century passenger shipping network. Specifically, it is concerned with the ways in which the network was established and with the efforts made to maintain it. It is found that such a disciplinary network furthers the firm's goal of shipping healthy and productive bodies for corporate profits and that transhipment facility Lloyd Hotel in Amsterdam was integral to the performance and maintenance of such a transnational disciplinary network. The key consequence of such disciplined mobility was the creation of an emotional passenger‐migrant subject shaped in relation to the power of corporate, cultural and other authorities in maritime travel and migration. In identifying this historic network of disciplined mobility and its emotional subject, this paper seeks to reveal the emotional geographies relating to mobile subjectivities and the power relations associated with their historically significant travels.     

Fluid flow in shear zones: insights from the geometry and evolution of ore bodies at Renco gold mine,Zimbabwe

The geometry of mineral deposits can give insights into fluid flow in shear zones. Lode gold ore bodies at Renco Mine, in the Limpopo Belt, Zimbabwe, occur as siliceous breccias and mylonites within amphibolite facies shear zones that dip either gently or steeply. The two sets of ore bodies formed synchronously from hydrothermal fluids. The ore bodies are oblate, but have well‐defined long axes. Larger ore bodies are more oblate. High‐grade gold ore shoots have long axes that plunge down dip; this direction is perpendicular to the long axes of the low‐grade ore bodies. The centres of the high‐grade ore bodies align within the low‐grade ore bodies along strike in both gently and steeply dipping groups. The range of sizes and shapes of the ore bodies are interpreted as a growth sequence. Geometrical models are proposed for the gently and steeply dipping ore bodies, in which individual ore bodies grow with long axes plunging down dip, and merge to form larger, more oblate ore bodies. The models show that when three or more ore bodies coalesce, the long axis of the merged ore body is perpendicular to the component ore bodies, and that ore bodies in the deposit may have a range of shapes due to both growth of individual ore bodies, and their coalescence. The long axes of the high‐grade ore bodies are parallel to the shear directions of both the gently and steeply dipping dip slip shear zones, which were the directions of greatest permeability and fluid flow. The larger, lower grade bodies, which may have formed by coalescence, are elongate perpendicular to these directions.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

Heat transport at boiling,near‐critical conditions

Two‐phase flow and near‐critical phenomena are likely to enhance energy transport in high‐temperature hydrothermal systems. We present a series of two‐dimensional simulations of two‐phase flow of pure water at near‐critical conditions. The results show that at near‐critical conditions, two‐phase convection can be more efficient in transporting energy than single‐phase convection. The highest heat fluxes are attained when two‐phase heat‐pipes form near the bottom boundary, recharging the root of the upflow zone and thereby enabling the formation of broad upflow regions. When the system becomes more vapor‐dominated, it loses this ability, upflow zones become narrower and the energy efficiency drops to more moderate values.     

The significance of pockmarks to understanding fluid flow processes and geohazards

Underwater gas and liquid escape from the seafloor has long been treated as a mere curiosity. It was only after the advent of the side‐scan sonar and the subsequent discovery of pockmarks that the scale of fluid escape and the moon‐like terrain on parts of the ocean floor became generally known. Today, pockmarks ranging in size from the ‘unit pockmark’ (1–10 m wide, < 0.6 m deep) to the normal pockmark (10–700 m wide, up to 45 m deep) are known to occur in most seas, oceans, lakes and in many diverse geological settings. In addition to indicating areas of the seabed that are ‘hydraulically active’, pockmarks are known to occur on continental slopes with gas hydrates and in association with slides and slumps. However, possibly their potentially greatest significance is as an indicator of deep fluid pressure build‐up prior to earthquakes. Whereas only a few locations containing active (bubbling) pockmarks are known, those that become active a few days prior to major earthquakes may be important precursors that have been overlooked. Pockmark fields and individual pockmarks need to be instrumented with temperature and pressure sensors, and monitoring should continue over years. The scale of such research calls for a multinational project in several pockmark fields in various geological settings.     

Fluid flow and stability of the US continental slope offshore New Jersey from the Pleistocene to the present

We predict that portions of the New Jersey continental slope were unstable approximately 0.5 million years ago. This instability was caused by rapid sediment loading during a Pleistocene sea‐level lowstand and by flow focusing in underlying, permeable Miocene strata. The simulated instability is consistent with soft‐sediment deformation and small slumps in Pleistocene strata of the Hudson Apron. Stability of the New Jersey margin has increased since 0.3 Ma because sedimentation rate has decreased. Today, the modelled factor of safety (FS) for the upper slope is approximately 1.5 whereas in the lower slope it exceeds 3. We predict that sedimentation rate is a dominant factor on slope stability. When rapid and asymmetric loading of a highly permeable sedimentary layer occurs, the location of instability can shift seaward to regions where sedimentation rates are low. Stability calculations use pressures and effective stresses predicted by a coupled sedimentation‐fluid flow model. This hydrodynamic analysis demonstrates how the interplay of sedimentation and fluid migration affects the distribution, timing, and size of sedimentary failures.     

Pore pressure evolution and fluid flow during visco‐elastic single‐layer buckle folding

Pore pressure and fluid flow during the deformational history of geologic structures are directly influenced by tectonic deformation events. In this contribution, 2D plane strain finite element analysis is used to study the influence of different permeability distributions on the pore pressure field and associated flow regimes during the evolution of visco‐elastic single‐layer buckle folds. The buckling‐induced fluid flow regimes indicate that flow directions and, to a lesser degree, their magnitudes vary significantly throughout the deformation and as a function of the stratigraphic permeability distribution. The modelling results suggest that the volumetric strain and the permeability distribution significantly affect the resulting flow regime at different stages of fold development. For homogeneous permeability models (k > 10^?21 m²), low strain results in a mostly pervasive fluid flow regime and is in agreement with previous studies. For larger strain conditions, fluid focusing occurs in the buckling layer towards the top of the fold hinge. For low permeabilities (<10^?21 m²), local focused flow regimes inside the buckling layer emerge throughout the deformation history. For models featuring a low‐permeability layer embedded in a high‐permeability matrix or sandwiched between high‐permeability layers, focused flow regimes inside the folded layer result throughout the deformation history, but with significant differences in the flow vectors of the surrounding layers. Fluid flow vectors induced by the fold can result in different, even reversed, directions depending on the amount of strain. In summary, fluid flow regimes during single‐layer buckling can change from pervasive to focused and fluid flow vectors can be opposite at different strain levels, that is the flow vectors change significantly through time. Thus, a complete understanding of fluid flow regimes associated with single‐layer buckle folds requires consideration of the complete deformation history of the fold.