This paper presented empirical models of describing reaction rate vs. hydrochloric acid temperature and concentration by regressing experimental data. And this paper also introduced the dependent reaction heat model into the thermal non-equilibrium models and coupled with two-scale continuum model to obtain governing equations for describing wormholing under non-isothermal conditions. The governing equations were discretized by implicit difference method and solved by programing. The effects of temperature on wormholing have been investigated based on the simulation results on 2-D vision. A significant difference of the effluent temperature between the dependent reaction heat model in present work and the constant reaction heat model in available literatures was observed, especially in high injection rate and strong acid concentration. In addition, the tendencies of optimum injection capacity vs. injection temperature under isothermal and non-isothermal conditions were almost different. Finally, an optimum injection temperature was found by changing the injection temperature under non-isothermal conditions. 13Open Journal of Yangtze Gas and Oil the production of carbonate reservoirs were based on investigations of the initiation and propagation of wormholes or unstable dissolution fronts. Many mathematical models were developed to describe the wormhole patterns observed in experiment results [1] [2] [3] [4]: the capillary tube model [5] [6] [7]; the network model [8]; the averaged model [4] [9]; the fractal model [1] [2]. Pangaet al. (2005) coupled Darcy-scale dissolution model with the pore-scale model to summarize the governing equations which is called the two-scale continuum model [10]. In this model, the macroscopic properties such as mass transfer coefficients and dispersion properties were properly considered. Since the results calculated from the two-scale continuum model matched experiment results very well, most of the extension simulation studies were based on them. Kalia and Balakotaiah (2007) studied dissolution patterns in radial flow by extending the linear model to radial model [11]. Maheshwari et al. (2013) analyzed the 3-D dissolution patterns [12]. Smith et al. (2013) simulated core CO 2 flooding of low permeability carbonates coupled with X-ray computed microtomography (XCMT) [13]. The wormholing process is very complicated since there are many known factors influencing convection, dispersion and reaction properties, such as acid system, concentration of acid, porosity or permeability distribution, fracture or vug medium, injection conditions, acid flow types and temperature. In the physical heterogeneities, Kalia and Balakotaiah (2009) used a unique heterogeneity parameter to characterize porosity heterogeneity, and found that the dissolution depended on the orientation of fractures [14]; Liu et al. (2012) introduced a The designate parameters f γ , m Cl , B HCl , E and C HCl in Equation (8) are given by Alkattan et al. (1998) [28].
Pore Scale ModelThe dissolution process determined local qua...