reaction due to acid rock reaction joule thomson effect additionally embodiments of methods according to the present disclosure presents at least one or more of the following features storing the computed fluid and temperature distribution in the rock from the previous time as initial condition for the current time step or storing the previous rate intake history or storing the evolution of rock properties such as permeability using an automated optimization algorithm for adjusting the flow rate distribution and or formation properties to match the predicted temperature profile with dts measurement and account for the overall temperature evolution during injection using supplemental shut in dts data to adjust flow rate distribution and match predicted temperature profile with dts measurement there is shown in fig3 a an example of the simulated temperature distribution in the wellbore and rock grids utilized for calculation purposes the lighter shaded areas indicate lower simulated temperatures while the darker shaded areas indicate higher simulated temperatures which represents an output 300 of the thermal model of the present disclosure there is shown in fig3 b a plurality of temperature profiles in the wellbore similar to the temperature profiles 202 204 206 and 208 of fig2 at spaced apart time intervals which is utilized to update the model as discussed hereinabove the overall workflow of the disclosed method is shown in fig4 wherein an embodiment of a method of the.