首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
A parametric two-phase, oil–water relative permeability/capillary pressure model for petroleum engineering and environmental applications is developed for porous media in which the smaller pores are strongly water-wet and the larger pores tend to be intermediate- or oil-wet. A saturation index, which can vary from 0 to 1, is used to distinguish those pores that are strongly water-wet from those that have intermediate- or oil-wet characteristics. The capillary pressure submodel is capable of describing main-drainage and hysteretic saturation-path saturations for positive and negative oil–water capillary pressures. At high oil–water capillary pressures, an asymptote is approached as the water saturation approaches the residual water saturation. At low oil–water capillary pressures (i.e. negative), another asymptote is approached as the oil saturation approaches the residual oil saturation. Hysteresis in capillary pressure relations, including water entrapment, is modeled. Relative permeabilities are predicted using parameters that describe main-drainage capillary pressure relations and accounting for how water and oil are distributed throughout the pore spaces of a porous medium with mixed wettability. The capillary pressure submodel is tested against published experimental data, and an example of how to use the relative permeability/capillary pressure model for a hypothetical saturation-path scenario involving several imbibition and drainage paths is given. Features of the model are also explained. Results suggest that the proposed model is capable of predicting relative permeability/capillary pressure characteristics of porous media mixed wettability.  相似文献   

2.
Three-phase flow is a key process occurring in subsurface reservoirs, for example, during $\text{ CO }_2$ sequestration and enhanced oil recovery techniques such as water alternating gas (WAG) injection. Predicting three-phase flow processes, for example, the increase in oil recovery during WAG, requires a sound understanding of the fundamental flow physics in water- to oil-wet rocks to derive physically robust flow functions, i.e. relative permeability and capillary pressure. In this study, we use pore-network modelling, a reliable and physically based simulation tool, to predict the flow functions. We have developed a new pore-scale network model for rocks with variable wettability, from water- to oil-wet. It comprises a constrained set of parameters that mimic the wetting state of a reservoir. Unlike other models, it combines three main features: (1) A novel thermodynamic criterion for formation and collapse of oil layers. The new model hence captures wetting film and layer flow of oil adequately, which affects the oil relative permeability at low oil saturation and leads to accurate prediction of residual oil. (2) Multiple displacement chains, where injection of one phase at the inlet triggers a chain of interface displacements throughout the network. This allows for the accurate modelling of the mobilisation of many disconnected phase clusters that arise, in particular, during higher order WAG floods. (3) The model takes realistic 3D pore-networks extracted from pore-space reconstruction methods and CT images as input, preserving both topology and pore shape of the sample. For water-wet systems, we have validated our model with available experimental data from core floods. For oil-wet systems, we validated our network model by comparing 2D network simulations with published data from WAG floods in oil-wet micromodels. This demonstrates the importance of film and layer flow for the continuity of the various phases during subsequent WAG cycles and for the residual oil saturations. A sensitivity analysis has been carried out with the full 3D model to predict three-phase relative permeabilities and residual oil saturations for WAG cycles under various wetting conditions with different flood end-points.  相似文献   

3.
We apply steady-state capillary-controlled upscaling in heterogeneous environments. A phase may fail to form a connected path across a given domain at capillary equilibrium. Moreover, even if a continuous saturation path exists, some regions of the domain may produce disconnected clusters that do not contribute to the overall connectivity of the system. In such cases, conventional upscaling processes might not be accurate since identification and removal of these isolated clusters are extremely important to the global connectivity of the system and the stability of the numerical solvers. In this study, we address the impact of percolation during capillary-controlled displacements in heterogeneous porous media and present a comprehensive investigation using random absolute permeability fields, for water-wet, oil-wet and mixed-wet systems, where J-function scaling is used to relate capillary pressure, porosity and absolute permeabilities in each grid cell. Important information is revealed about the average connectivity of the phases and trapping at the Darcy scale due to capillary forces. We show that in oil-wet and mixed-wet media, large-scale trapping of oil controlled by variations in local capillary pressure may be more significant than the local trapping, controlled by pore-scale displacement.  相似文献   

4.
The Effect of Wettability on Three-Phase Relative Permeability   总被引:3,自引:0,他引:3  
We study three-phase flow in water-wet, oil-wet, and fractionally-wet sandpacks. We use CT scanning to measure directly the oil and water relative permeabilites for three-phase gravity drainage. In an analogue experiment, we measure pressure gradients in the gas phase to determine the gas relative permeability. Thus we find all three relative permeabilities as a function of saturation. We find that the gas relative permeability is approximately half as much in a oil-wet medium than in an water-wet medium at the same gas saturation. The water relative permeability in the water-wet medium and the oil relative permeability in the oil-wet medium are similar. In the water-wet medium the oil relative permeability scales as k roS o 4 for S o>S or, where S or is the waterflood residual oil saturation. With octane as the oil phase, k roS o 2 for S o<S or, while with decane as the oil phase, k ro falls sharply for S o<S or. The water relative permeability in the oil-wet medium resembles the oil relative permeability in the water-wet medium for a non-spreading oil such as decane. These observations can be explained in terms of wetting, spreading, and the pore scale configurations of fluid.  相似文献   

5.

Three-phase flow in porous media is encountered in many applications including subsurface carbon dioxide storage, enhanced oil recovery, groundwater remediation and the design of microfluidic devices. However, the pore-scale physics that controls three-phase flow under capillary dominated conditions is still not fully understood. Recent advances in three-dimensional pore-scale imaging have provided new insights into three-phase flow. Based on these findings, this paper describes the key pore-scale processes that control flow and trapping in a three-phase system, namely wettability order, spreading and wetting layers, and double/multiple displacement events. We show that in a porous medium containing water, oil and gas, the behaviour is controlled by wettability, which can either be water-wet, weakly oil-wet or strongly oil-wet, and by gas–oil miscibility. We provide evidence that, for the same wettability state, the three-phase pore-scale events are different under near-miscible conditions—where the gas–oil interfacial tension is ≤?1 mN/m—compared to immiscible conditions. In a water-wet system, at immiscible conditions, water is the most-wetting phase residing in the corners of the pore space, gas is the most non-wetting phase occupying the centres, while oil is the intermediate-wet phase spreading in layers sandwiched between water and gas. This fluid configuration allows for double capillary trapping, which can result in more gas trapping than for two-phase flow. At near-miscible conditions, oil and gas appear to become neutrally wetting to each other, preventing oil from spreading in layers; instead, gas and oil compete to occupy the centre of the larger pores, while water remains connected in wetting layers in the corners. This allows for the rapid production of oil since it is no longer confined to movement in thin layers. In a weakly oil-wet system, at immiscible conditions, the wettability order is oil–water–gas, from most to least wetting, promoting capillary trapping of gas in the pore centres by oil and water during water-alternating-gas injection. This wettability order is altered under near-miscible conditions as gas becomes the intermediate-wet phase, spreading in layers between water in the centres and oil in the corners. This fluid configuration allows for a high oil recovery factor while restricting gas flow in the reservoir. Moreover, we show evidence of the predicted, but hitherto not reported, wettability order in strongly oil-wet systems at immiscible conditions, oil–gas–water, from most to least wetting. At these conditions, gas progresses through the pore space in disconnected clusters by double and multiple displacements; therefore, the injection of large amounts of water to disconnect the gas phase is unnecessary. We place the analysis in a practical context by discussing implications for carbon dioxide storage combined with enhanced oil recovery before suggesting topics for future work.

  相似文献   

6.
A look into the literature on the temperature dependency of oil and water relative permeabilities reveals contradictory reports. There are some publications reporting shifts in the water saturation range as well as variations in the relative permeability curves by temperature. On the other hand, some authors have blamed the experimental artifacts, viscous instabilities and fingering issues for these variations. We have performed core flooding experiments to further investigate this issue. Glass bead packs and sand packs were used as the porous media, and Athabasca bitumen with varying viscosities was displaced by hot water at differing temperatures. The unsteady-state method of relative permeability measurement was applied and the experimental data were history matched by a simulator that is tailor made to predict the relative permeabilities. The matches were obtained by varying the relative permeability correlation parameters. The results indicated that the initial water saturation has a direct relation with temperature, while residual oil saturation generally drops at higher temperatures. Although the water saturation range shifts, no direct and unique trend for either oil or water relative permeability is justified. The spread in relative permeabilities especially in the case of higher permeable cores suggests that viscous instabilities are present. As the same saturation shift happens by only changing the oil viscosity, the relative permeability variations with temperature can be attributed to oil to water viscosity ratio changes with temperature. Temperature dependency of relative permeabilities is more related to experimental artifacts, viscous fingering and viscosity changes than fundamental flow properties.  相似文献   

7.
Performance of a polymer flood process requires the knowledge of rheological behavior of the polymer solution and reservoir properties such as rock wettability. To provide a better understanding of effects of polymer chemistry and wettability on the performance of a polymer flood process, a comprehensive experimental study was conducted using a two-dimensional glass micromodel. A series of water and polymer flood processes were carried out at different polymer molecular weights, degrees of polymer hydrolysis, and polymer concentrations in both water-wet and oil-wet systems. Image processing technique was applied to analyze and compare microscopic and macroscopic displacement behaviors of polymer solution in each experiment. From micro-scale observations, the configuration of connate water film, polymer solution trapping, flow of continuous and discontinuous strings of polymer solution, piston-type displacement of oil, snap-off of polymer solution, distorted flow of polymer solution, emulsion formation, and microscopic pore-to-pore sweep of oil phase were observed and analyzed in the strongly oil-wet and water-wet media. Rheological experiments showed that a higher polymer molecular weight, degree of hydrolysis, and concentration result in a higher apparent viscosity for polymer solution and lower oil–polymer viscosity ratio. It is also shown that these parameters have different impacts on the oil recovery in different wettabilities. Moreover, a water-wet medium generally had higher recovery in contrast with an oil-wet medium. This experimental study illustrates the successful application of glass micromodel techniques for studying enhanced oil recovery (EOR) processes in five-spot pattern and provides a useful reference for understanding the displacement behaviors in a typical polymer flood process.  相似文献   

8.
The impact of fractional wettability on the production characteristics of a VAPEX process at the macroscale was investigated. Conventional VAPEX experiments were conducted in a 220 Darcy random packing of glass beads in a rectangular physical model and n-pentane was used to recover the Cold Lake bitumen from the oil-saturated model in the absence of connate water. The composition of oil-wet beads in the packed bed was altered from completely water-wet beads to completely oil-wet beads at different proportions of oil-wet beads mixed with water-wet beads. A substantial increase (about 40%) in the production rate of live oil was observed during the VAPEX process when the wettability of the porous packing was entirely oil-wet beads. A critical oil-wet fraction of 0.66 was found for the heterogeneous packing of water-wet and oil-wet beads of similar size distribution. Above this critical composition, the live oil production rate was not affected by further increase in the proportion of the oil-wet beads. It is believed that above this critical composition of the oil-wet beads, the crevice flow process is dominated by the continuity of higher conductivity live oil films between particles through the oil-wet regions. Below this critical composition, the live oil production rate increased linearly with the fraction of the oil-wet beads in the packing. The oil-wet regions favor the live oil drainage compared to that of the water-wet regions as they enhance the rate of imbibition of the live oil from the oil-filled pores to the vacated pores near the nominal VAPEX interface. These two factors enhance the live oil production rate during the VAPEX process. The solvent content of the live oil, the solvent-to-oil ratio (SOR), and the residual oil saturation did not correlate strongly with the proportion of the oil-wet beads in the packing. The average solvent content of the live oil and the residual oil saturation were measured to be 48% by weight and 7% by volume respectively.  相似文献   

9.
In this study, the main recovery mechanisms behind oil/water/gas interactions during the water-alternating-gas (WAG) injection process, in a network of matrix/fracture, were fundamentally investigated. A visual micromodel was utilized to provide insights into the potential applications of WAG process in fractured oil-wet media as well as the possibility of observing microscopic displacement behavior of fluids in the model. The model was made of an oil-wet facture/matrix network system, comprised of four matrix blocks surrounded with fractures. Different WAG injection scenarios, such as slug arrangements and the effects of fluid injection rates on oil recovery were studied. A new equation representing the capillary number, considering the fracture viscous force and matrix capillary force, was developed to make the experimental results more similar to a real field. In general, WAG tests performed in the fractured model showed a higher oil recovery factor compared with the results of gas and water injection tests at their optimum rates. The results showed that the presence of an oil film, in all cases, was the main reason for co-current drainage and double displacement of oil under applied driving forces. Furthermore, the formation of oil liquid bridges improved the recovery efficiency, which was greatly influenced by the size of fracture connecting the two matrix blocks; these connecting paths were more stable when there was initial water remaining in the media. Analyzing different recovery curves and microscopic view of the three phases in the transparent model showed that starting an injection mode with gas (followed by repeated small slugs of water and gas), could considerably improve oil recovery by pushing water into the matrix zone and increasing the total sweep efficiency.  相似文献   

10.
We present a mechanistic model of retrograde condensation processes in two- and three-dimensional capillary tube networks under gravitational forces. Condensate filling-emptying cycles in pore segments and gas connection–isolation cycles are included. With the pore-level distribution of gas and condensate in hand, we determine their corresponding relative permeabilities. Details of pore space and displacement are subsumed in pore conductances. Solving for the pressure field in each phase, we find a single effective conductance for each phase as a function of condensate saturation. Along with the effective conductance for the saturated network, the relative permeability for each phase is calculated. Our model porous media are two- and three-dimensional regular networks of pore segments with distributed size and square cross-section. With a Monte Carlo sampling we find the optimum network size to avoid size effects and then we investigate the effect of network dimensionality and pore size distribution on the relative permeabilities of gas and condensate.  相似文献   

11.
Displacement of a viscous fluid by a lower viscosity immiscible fluid (such as waterflood of a viscous oil) in a porous medium is unstable. The displacement front generates viscous fingers which lead to low oil recovery efficiency. These fingers are much smaller in width than typical reservoir simulation grid blocks, and capturing their effect in reservoir simulation is important. A dimensionless scaling group (viscous finger number) had been suggested in the past, which has a power-law relationship with the breakthrough recovery and cumulative recovery in unstable core floods. The relative permeability used in large grid block simulations had been modified to so-called pseudo-relative permeability on the basis of the dimensionless group, thus incorporating the effect of fingers in waterflood predictions. However, the previous proposed models were constructed from experiments in only water-wet rocks. This paper extends the recent viscous fingering models to oil-wet systems. Sandstone cores were treated to alter the wettability to oil-wet. Adverse viscosity water floods were performed in oil-wet cores. Viscosity ratio, velocity and diameter were varied. It is shown that the previously developed viscous finger number does not work for the oil-wet experiments. The correlating dimensionless number is modified for oil-wet systems; it is also different from the dimensionless group identified by Peters and Flock (Soc Petroleum Eng, 1981. doi: 10.2118/8371-PA) for oil-wet cores. A pseudo-relative permeability model has been developed for oil-wet cores. Corefloods have been matched by the new pseudo-relative permeability model to determine the model parameters. This pseudo-relative permeability model can be used in reservoir simulations of water and polymer floods in viscous oil-wet reservoirs.  相似文献   

12.
When regions of three-phase flow arise in an oil reservoir, each of the flow parameters, i.e. capillary pressures and relative permeabilities, are generally functions of two phase saturations and depend on the wettability state. The idea of this work is to generate consistent pore-scale based three-phase capillary pressures and relative permeabilities. These are then used as input to a 1-D continuum core- or reservoir-scale simulator. The pore-scale model comprises a bundle of cylindrical capillary tubes, which has a distribution of radii and a prescribed wettability state. Contrary to a full pore-network model, the bundle model allows us to obtain the flow functions for the saturations produced at the continuum-scale iteratively. Hence, the complex dependencies of relative permeability and capillary pressure on saturation are directly taken care of. Simulations of gas injection are performed for different initial water and oil saturations, with and without capillary pressures, to demonstrate how the wettability state, incorporated in the pore-scale based flow functions, affects the continuum-scale displacement patterns and saturation profiles. In general, wettability has a major impact on the displacements, even when capillary pressure is suppressed. Moreover, displacement paths produced at the pore-scale and at the continuum-scale models are similar, but they never completely coincide.  相似文献   

13.
Fractional wettability has been widely recognized in most of the oil reservoirs and it is a crucial factor that controls the fluid flow behavior in porous medium. The overall effect of the proportion of oil-wet grains on the fluid flow properties has been well discussed. However, recent studies found that the random distribution and coordination of oil-wet and water-wet grains could make multi-phase flow behaviors extremely complicated in such media. The multiphase flow mechanisms in fractional wettability media remains unclear. In this study, oil imbibition experiments were systematically conducted using glass cylinders packed with fractional-wet glass beads. To study the effect of fractional wettability on multiple-phase flow properties, samples with different oil-wet grain proportions were prepared, and fifteen repeated experiments were conducted for each oil-wet proportion. The experimental results showed that oil imbibition was largely dependent on but not strictly a function of the proportion of oil-wet grains in the medium. The imbibition behaviors of samples with the same fractional proportion could vary significantly, as some samples exhibited complete oil migration, while others did not. This probabilistic phenomenon is likely due to the random distribution of oil-wet and water-wet grains. A pore throat may behave as oil-wet or water-wet depending on the relative proportion of oil-wet grains the pore throat contains. When the grains that comprise the pore throat are dominated by oil-wet grains, the throat behaves as oil-wet, and vice versa. Only when these oil-wet pore throats are connected to form a complete oil-wet pathway throughout the medium can the oil continuously imbibe into the medium. Therefore, the extent of oil imbibition depends on the completeness of the oil-wet pathway, which is controlled by the proportion of oil-wet grains in the medium. The higher the proportion of oil-wet grains in the medium, the larger the number of oil-wet pore throats that can form; thus, the higher the possibility that those oil-wet pore throats can connect to form continuous oil-wet pathways.  相似文献   

14.

We predict waterflood displacement on a pore-by-pore basis using pore network modelling. The pore structure is captured by a high-resolution image. We then use an energy balance applied to images of the displacement to assign an average contact angle, and then modify the local pore-scale contact angles in the model about this mean to match the observed displacement sequence. Two waterflooding experiments on oil-wet rocks are analysed where the displacement sequence was imaged using time-resolved synchrotron imaging. In both cases the capillary pressure in the model matches the experimentally obtained values derived from the measured interfacial curvature. We then predict relative permeability for the full saturation range. Using the optimised contact angles distributed randomly in space has little effect on the predicted capillary pressures and relative permeabilities, indicating that spatial correlation in wettability is not significant in these oil-wet samples. The calibrated model can be used to predict properties outside the range of conditions considered in the experiment.

  相似文献   

15.
Functions of the relative phase permeabilities which take into account the tensor nature of the relations between the absolute and phase permeabilities are constructed on the basis of experimental data. Problems of the displacement of gas by water and of water by gas in an orthotropic pore space are studied. The relative phase permeabilities obtained are compared with functions that depend only on the local saturation. The effect of the phase compressibility and the anisotropy of the absolute permeability tensor is considered.  相似文献   

16.
A network model is established through the techniques of image reconstruction, a thinning algorithm, and pore–throat information extraction with the aid of an industrial microfocus CT scanning system. In order to characterize actual rock pore–throat structure, the established model is modified according to the matching of experimental factors such as porosity, permeability, and the relative permeability curve. On this basis, the impacts of wetting angle, pore radius, shape factor, pore–throat ratio, and coordination number as applied to microscopic remaining oil distribution after water flooding are discussed. For a partially wetting condition, the displacement result of a water-wet pore is somewhat better than that of an oil-wet pore as a whole, and the possibility of any remaining oil is relatively low. Taking the comprehensive effects of various factors into account, a prediction method of remaining oil distribution is presented through the use of fuzzy comprehensive evaluation. It is seen that this method can predict whether there is remaining oil or not in the pore space with satisfactory accuracy, which is above 75%. This method thus provides guidance for a better understanding of the microscopic causes of the remaining oil.  相似文献   

17.
We present results of high-pressure micromodel visualizations of pore-scale fluid distribution and displacement mechanisms during the recovery of residual oil by near-miscible hydrocarbon gas and SWAG (simultaneous water and gas) injection under conditions of very low gas–oil IFT (interfacial tension), negligible gravity forces and water-wet porous medium. We demonstrate that a significant amount of residual oil left behind after waterflooding can be recovered by both near-miscible gas and SWAG injection. In particular, we show that in both processes, the recovery of the contacted residual oil continues behind the main gas front and ultimately all of the oil that can be contacted by the gas will be recovered. This oil is recovered by a microscopic mechanism, which is strongly linked to the low IFT between the oil and gas and to the perfect spreading of the oil over water, both of which occur as the critical point of the gas–oil system is approached. Ultimate oil recovery by near-miscible SWAG injection was as high as near-miscible gas injection with SWAG injection using much less gas compared to gas injection. Comparison of the results of SWAG experiments with two different gas fractional flow values (SWAG ratio) of 0.5 and 0.2 shows that fractional flow of the near-miscible gas injected simultaneously with water is not a crucial factor for ultimate oil recovery. This makes SWAG injection an attractive IOR (improved oil recovery) process especially for reservoirs, where continuous and high-rate gas injection is not possible (e.g. due to supply constraint).  相似文献   

18.
In three-phase flow, the macroscopic constitutive relations of capillary pressure and relative permeability as functions of saturation depend in a complex manner on the underlying pore occupancies. These three-phase pore occupancies depend in turn on the interfacial tensions, the pore sizes and the degree of wettability of the pores, as characterised by the cosines of the oil–water contact angles. In this work, a quasi-probabilistic approach is developed to determine three-phase pore occupancies in media where the degree of wettability varies from pore to pore. Given a set of fluid and rock properties, a simple but novel graphical representation is given of the sizes and oil–water contact angles underlying three-phase occupancies for every allowed combination of capillary pressures. The actual phase occupancies are then computed using the contact angle probability density function. Since a completely accessible porous medium is studied, saturations, capillary pressures, and relative permeabilities are uniquely related to the pore occupancies. In empirical models of three-phase relative permeability it is of central importance whether a phase relative permeability depends only on its own saturation and how this relates to the corresponding two-phase relative permeability (if at all). The new graphical representation of pore sizes and wettabilities clearly distinguishes all three-phase pore occupancies with respect to these saturation-dependencies. Different types of saturation-dependencies may occur, which are shown to appear in ternary saturation diagrams of iso-relative permeability curves as well, thus guiding empirical approaches. However, for many saturation combinations three-phase and two-phase relative permeabilities can not be linked. In view of the latter, the present model has been used to demonstrate an approach for three-phase flow modelling on the basis of the underlying pore-scale processes, in which three-phase relative permeabilities are computed only along the actual flow paths. This process-based approach is used to predict an efficient strategy for oil recovery by simultaneous water-alternating-gas (SWAG) injection.  相似文献   

19.
A method has been developed for investigating the relative permeabilities of porous media for oil and for aqueous solutions of polymers; experimental equipment has been developed for determining the phase permeabilities by a stationary method. Investigations were made of the influence of polyacrylamide additives on the change in the relative permeabilities for the simultaneous flow of water and a nonpolar hydrocarbon liquid. It was established that addition of the polymer can lead to a simultaneous reduction in the relative permeability for the wetting liquid and an increase for the nonwetting liquid. The phase permeabilities were obtained for oil and water moving behind a fringe of polymer substance. It was established that the phase permeability for the water phase is a function of the saturation and the amount of sorbate. A cycle of experimental investigations was made into the influence of the rate of pumping and the concentration of the dissolved polymer on the change in the relative permeabilities.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 163–167, July–August, 1980.  相似文献   

20.
认识双重多孔介质中油水两相微观渗流机制是回答形成什么类型的裂隙网络可提高油藏采收率的关键. 微裂隙的分布可以提高多孔介质的绝对渗透率,但对于基质孔隙中的流体介质,微裂隙的存在会引起多孔介质中局部流体压力和流场的变化,导致局部流动以微裂隙流动为主,甚至出现窜流现象,降低驱油效率. 本文基于孔与裂隙双重网络模型,在网络进口设定两条平行等长且具有一定间隔的微裂隙,分析微裂隙的相对间隔(微裂隙之间距离/喉道长度)和微裂隙相对长度(微裂隙长度/喉道长度)对于微观渗流特征的影响. 结果表明:随微裂隙相对长度的增加,出现驱油效率逐渐降低,相对渗透率曲线中的油水共渗区水饱和度和等渗点增加,油水两相的共渗范围减小等现象;随着微裂隙之间相对间隔增大,周围越来越多的基质孔穴间的压力差减小,在毛管压力的限制下,驱替相绕过这些区域,而导致水窜现象.   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号