in the following description , numerous details are set forth to provide an understanding of the present invention . however , it will be understood by those skilled in the art that the present invention may be practiced without these details and that numerous variations or modifications from the described embodiments may be possible . in this respect , before explaining at least one embodiment of the invention in detail , it is to be understood that the invention is not limited in its application to the details of construction and to the arrangements of the components set forth in the following description or illustrated in the drawings . the invention is capable of other embodiments and of being practiced and carried out in various ways . also , it is to be understood that the phraseology and terminology employed herein are for the purpose of the description and should not be regarded as limiting . as used herein “ oilfield ” is a generic term including any hydrocarbon - bearing geologic formation , or formation thought to include hydrocarbons , including onshore and offshore . as used herein when discussing fluid flow , the terms “ divert ”, “ diverting ”, and “ diversion ” mean changing the direction , the location , the magnitude or all of these of all or a portion of a flowing fluid . a “ wellbore ” may be any type of well , including , but not limited to , a producing well , a non - producing well , an experimental well , and exploratory well , and the like . wellbores may be vertical , horizontal , some angle between vertical and horizontal , and combinations thereof , for example a vertical well with a non - vertical component . as mentioned previously , to increase the net permeability of a reservoir , it is common to perform a well stimulation treatment . a common stimulation technique consists of injecting an acid that reacts with and dissolves the formation damage or a portion of the formation thereby creating alternative flow paths for the hydrocarbons to migrate through the formation to the well . this technique known as acidizing ( or more generally as matrix stimulation ) may eventually be associated with fracturing if the injection rate and pressure is enough to induce the formation of a fracture in the reservoir . fluid placement is critical to the success of stimulation treatments . natural reservoirs are often heterogeneous ; the fluid will preferentially enter areas of higher permeability in lieu of entering areas where it is most needed . each additional volume of fluid follows the path of least resistance , and continues to invade in zones that have already been treated . therefore , it is difficult to place the treating fluids in severely damaged and lower permeability zones . in order to control placement of treating fluids , various techniques have been employed . mechanical techniques involve for instance the use of ball sealers and packers and of coiled tubing placement to specifically spot the fluid across the zone of interest . non - mechanical techniques typically make use of gelling agents as diverters for temporarily impairing the areas of higher permeability and increasing the proportion of the treating zone that goes into the areas of lower permeability . therefore , for evaluation and optimization of matrix treatments it is of interest to measure the placement of treating fluids . the present invention determines fluid placement in the reservoir by the measurement and interpretation of one or more of temperature , pressure , and flow rate of fluids injected into the wellbore and close to the fluid exit from an oilfield tubular , such as coiled tubing , using special diagnostic plots . methods in accordance with the invention may be used prior to , during and post treatment , and any combination thereof , including during all of these . using one or more methods within the invention prior to reservoir treatment will allow estimation of formation damage in each layer of the reservoir from measurements of injection of an inert fluid , such as brine , along some or all of the entire length of the wellbore . the bottomhole temperature data gathered during the injection test can be interpreted in real time by the method proposed and “ zones of interest ” can be identified . use of one or more methods within the present invention during the treatment will allow monitoring and optimization of the treatment in real time . the data may be transmitted to the surface ( such as , by a stream of optical signals ) and may be displayed on a computer screen , personal digital assistant , cellular phone , or other electronic device for real time interpretation . placement of fluids in the formation may be optimized in real time by the use of diversion agents such as foam , inflatable open hole packers , fibers , and the like , and combinations thereof , to divert the stimulation where desired to potential zones . for example , if one finds that a certain reservoir layer is not being treated the injection rate of the fluids or the diverter volume or type may be changed or adjusted to divert the treating fluids to that layer . post treatment use of one or more methods within the present invention will allow evaluation of the effectiveness of the treatment by monitoring the injection of an inert fluid ( such as brine used for post flush ) to evaluate the stimulation achieved in each zone . alternatively the entire data set may be recorded and analyzed post treatment ( such as when telemetry equipment is not available ). methods of the invention allow for monitoring fluid placement during matrix treatments by measuring temperature of the wellbore fluids at a fixed distance from the fluid injection point . the methods of the invention rely on gathering temperatures and / or pressures , and in certain methods using specialized diagnostic plots to estimate the placement of fluids . systems of the invention are exemplified in four embodiments illustrated in fig1 - 4 , wherein like numerals are employed to described like components unless otherwise noted . it should be pointed out that the system embodiments illustrated in fig1 - 4 are illustrative only , and not intended to be limiting in any way . fig1 illustrates embodiment 100 , including a tubular 2 inserted in a cased or uncased wellbore 3 in a formation 5 , tubular 2 comprising a section of tubing 4 having at least one fluid injection port 6 and at least one temperature sensor 8 attached at a known location on tubular section 4 . system 100 includes a pump 10 for injecting a fluid through tubular 2 , tubular section 4 , and the at least one fluid injection port 6 and into formation 5 . a unit 12 allows generating , in real time or at a later time , diagnostic plot curves of temperature derivative with respect to time and temperature derivative with respect to coiled tubing depth , both obtained at a known fixed distance from the fluid injection port . a communication link 7 connects temperature sensor 8 with unit 12 , and optionally other units not illustrated . communication link 7 may be fiber optic , hard wire , or wireless . a curve shape interpreting unit 14 allows for interpreting the curves generating by unit 12 to determine location of regions of a hydrocarbon - bearing reservoir exhibiting flow of the injected fluid , where the flow ranges from zero to a non - zero value . referring now to fig2 there is illustrated schematically another system embodiment 200 within the invention , comprising a tubular 2 inserted in a cased or uncased wellbore 3 in a formation 5 , tubular 2 comprising a section of tubing 4 having at least one fluid injection port 6 and at least one temperature sensor 8 placed at a known location on tubular section 4 . system 200 also includes a pump 10 for injecting a fluid through tubular 2 , tubular section 4 and the at least one fluid injection port 6 . system 200 includes a measuring unit 16 for measuring time of arrival of the injected fluid at temperature sensor 8 . a communication link 7 connects temperature sensor 8 with unit 16 , and optionally other units not illustrated . communication link 7 may be fiber optic , hard wire , or wireless . although communication link 7 is illustrated as traversing through tubular 2 and tubular section 4 , link 7 may traverse in the annulus between tubular 2 and wellbore or production casing 3 . fig3 illustrates schematically another system embodiment 300 within the invention , and includes a tubular 2 inserted in a cased or uncased wellbore 3 in a formation 5 , tubular 2 comprising a section of tubing 4 having at least one fluid injection port 6 and at least one sensor 8 placed at a known location on tubular section 4 . system 300 includes a first pump 10 a for injecting a first fluid through tubular 2 , tubular section 4 , and the at least one fluid injection port 6 , the first fluid having a first fluid property value , and a second pump 10 b for injecting a second fluid through an annulus between tubular 2 and the cased or uncased wellbore 3 , the second fluid having a second fluid property value that is different from the first fluid property value . system 300 includes a measuring unit 18 for measuring a differential between the first and second fluid property values . the first and second properties may be temperature , pressure , flow rate , conductance , or some other measurable parameter . a communication link 7 connects sensor 8 with unit 18 , and optionally other units not illustrated . communication link 7 may be fiber optic , hard wire , or wireless . although communication link 7 is illustrated as traversing through tubular 2 and tubular section 4 , link 7 may traverse in the annulus between tubular 2 and wellbore or production casing 3 . fig4 illustrates schematically a fourth system embodiment 400 within the invention , and comprises a prediction unit 20 for predicting temperature as a function of reservoir permeability distribution at one or more sensors placed at known locations on a tubular 2 injected into a cased or uncased wellbore 3 of a formation 3 . tubular 2 comprises a tubular section 4 having at least one fluid injection port 6 ; a pump 10 for injecting a fluid through tubular 2 , tubular section 4 , and the at least one fluid injection port 6 , and a measuring unit 22 for measuring actual temperatures at the one or more temperature sensors 8 attached to or integral with tubular section 4 . system 400 further includes a calculation unit 24 for calculating error between the predicted and the measured temperatures , and for minimizing the errors by iteratively adjusting the permeability distribution along the wellbore length . a communication link 7 connects sensor 8 with unit 18 , and optionally other units not illustrated . communication link 7 may be fiber optic , hard wire , or wireless . although communication link 7 is illustrated as traversing through tubular 2 and tubular section 4 , link 7 may traverse in the annulus between tubular 2 and wellbore or production casing 3 . systems of the invention include those wherein the temperature sensors may be selected from thermally active temperature sensors and thermally passive temperature sensors , and wherein the flow meters may be selected from flow meter spinners , electromagnetic flow meters , ph sensors , resistivity sensors , optical fluid sensors and radioactive and / or non - radioactive tracer sensors , such as dna or dye sensors . systems of the invention may include means for using this information in realtime to evaluate and change , if necessary , one or more parameters of the fluid diversion . means for using the information sensed may comprise command and control sub - systems located at the surface , at the tool , or both . systems of the invention may include downhole flow control devices and / or means for changing injection hydraulics in both the annulus and tubing injection ports at the surface . systems of the invention may comprise a plurality of sensors capable of detecting fluid flow out of the tubular , below the tubular and up the annulus between the tubular and the wellbore in realtime mode that may have programmable action both downhole and at the surface . this may be accomplished using one or more algorithms allow quick realtime interpretation of the downhole data , allowing changes to be made at surface or downhole for effective treatment . systems of the invention may comprise a controller for controlling fluid direction and / or shut off of flow from the surface . exemplary systems of the invention may include fluid handling sub - systems able to improve fluid diversion through command and control mechanisms . these sub - systems may allow controlled fluid mixing , or controlled changing of fluid properties . systems of the invention may comprise one or more downhole fluid flow control devices that may be employed to place a fluid in a prescribed location in the wellbore , change injection hydraulics in the annulus and / or tubular from the surface , and / or isolate a portion of the wellbore . the inventive systems may further include different combinations of sensors / measurements above and below , ( and may also be at ) a fluid injection port in the tubular to determine / verify diversion of the fluid . systems and methods of the invention may include surface / tool communication through one or more communication links , including but not limited to hard wire , optical fiber , radio , or microwave transmission . in exemplary embodiments , the sensor measurements , realtime data acquisition , interpretation software and command / control algorithms may be employed to ensure effective fluid diversion , for example , command and control may be performed via preprogrammed algorithms with just a signal sent to the surface that the command and control has taken place , the control performed via controlling placement of the injection fluid into the reservoir and wellbore . in other exemplary embodiments , the ability to make qualitative measurements that may be interpreted realtime during a pumping service on coiled tubing or jointed pipe is an advantage . systems and methods of the invention may include realtime indication of fluid movement ( diversion ) out the downhole end of the tubular , which may include down the completion , up the annulus , and in the reservoir . two or more flow meters , for example electromagnetic flow meters , or thermally active sensors that are spaced apart from the point of injection at the end of the tubular may be employed . other inventive methods and systems may comprise two identical diversion measurements spaced apart from each other and enough distance above the fluid injection port at the end or above the measurement devices , to measure the difference in the flow each sensor measures as compared to the known flow through the inside of the tubular ( as measured at the surface ). the inventive methods and systems may employ multiple sensors that are strategically positioned and take multiple measurements , and may be adapted for flow measurement in coiled tubing , drill pipe , or any other oilfield tubular . systems of the invention may be either moving or stationary while the operation is ongoing . treatment fluids , which may be liquid or gaseous , or combination thereof , and / or combinations of fluids and solids ( for example slurries ) may be used in stimulation methods , methods to provide conformance , methods to isolate a reservoir for enhanced production or isolation ( non - production ), or combination of these methods . data gathered may either be used in a “ program ” mode downhole ; alternatively , or in addition , surface data acquisition may be used to make real time “ action ” decisions for the operator to act on by means of surface and downhole parameter control . fiber optic telemetry may be used to relay information such as , but not limited to , pressure , temperature , casing collar location ( ccl ), and other information uphole . as described therein , due to the large id of a straddle tool , a measurement tool is placed inside the straddle tool housing . a hole is added to the bullnose , and a tube is run from below the lower seal to inside the measurement tool . the measurement tool may then measure treating pressure , bottomhole temperature , depth via casing collar location ( ccl ), or some other parameter , as well as pressure below the lower seal of the straddle , which may be measured in real - time . by measuring the pressure below the lower seal , the operator can tell if the lower seal is leaking , and also if there is cross - flow from one zone to another . this has the potential to change how the job is performed in real time and optimize the treatment . this data would be evaluated realtime to determine if another treatment of zone is necessary . the inventive methods and systems may be employed in any type of geologic formation , for example , but not limited to , reservoirs in carbonate and sandstone formations , and may be used to optimize the placement of treatment fluids , for example , to maximize wellbore coverage and diversion from high perm and water / gas zones , to maximize their injection rate ( such as to optimize damkohler numbers and fluid residence times in each layer ), and their compatibility ( such as ensuring correct sequence and optimal composition of fluids in each layer ). the interpretation method proposed in the invention is illustrated by the following examples . an acid stimulation treatment was performed in an openhole section of a horizontal well in a carbonate formation . the treatment objective was to remove drilling induced damage . by default , the injected treatment fluids take the path of least resistance and invade the regions that are more permeable than others . however , it was difficult to pin - point the regions where the fluids were being injected because the initial injectivity of the zones and how injectivity changes with time was not known . therefore , monitoring of fluid placement was performed for evaluation and optimization of the treatment . the plot in fig5 shows bottomhole temperature data obtained during the acid stimulation treatment . the bottom curve , shaped like an “ m ” depicts the coiled tubing depth , whereas the second curve shows the bottomhole temperature . a temperature sensor was located in the bottomhole assembly on the end of the ct . prior to start of the acid treatment , brine was pumped from the coiled tubing while running in the hole to the heel . during this phase of the treatment the well was open to the pit , where returns were monitored . during the main treatment the acid was continuously pumped with the ct moving up and down the lateral length at a rate of nearly 6 feet / min [ 1 . 83 m / min ] and the injection rate was constant at nearly 2 bbl / min [ 0 . 32 m 3 / min ]. at the start of the job ( left part of the plot ), one can see that as the pumping of acid began and the formation was exposed to the acid stimulation fluid , the bottomhole temperature started to decrease . however , the temperature increased as the ct traversed down the lateral section . this can mislead one into believing that the majority of the fluids invaded the heel of the lateral . therefore , though temperature is a useful measurement which may hold the key to solving the problem , its representation alone in graphical format was insufficient to draw any meaningful conclusions as to where the fluid invasion was actually taking place in the open - hole formation . the plot of fig6 represents the data of “ 1st acid ” treatment in context with entire job data shown in fig5 . a closer look at the data indicated that the rate of change of bottomhole temperature was not constant even though the ct was moving at a constant rate of 6 ft / min [ 1 . 83 m / min ] and the acid injection was taking place at a near constant rate of 2 bbl / min [ 0 . 32 m 3 / min ]. the bottomhole temperature sensor was placed a few feet before the distal tip of the ct , and thus a change in temperature ( increase or decrease ) was observed when the fluid came out of the ct , that had a different temperature than the surroundings . the injected acid fluid either passed over the sensor in a direction opposite to that of ct movement , or if the ct sensor entered a region which was invaded earlier , if the fluid flow was in the same direction as ct movement . the fact that the entire section of the well was completed as open hole meant that the fluid was free to take the path of least resistance ; in this case it seemed to be somewhat away from the heel towards the toe of the lateral . the initial rapid reduction in bottomhole temperature indicated that the bottomhole temperature sensor was moving into a “ cooler region ”; where most of the fluid had already invaded ahead of the sensor and had cooled the region down before the sensor reached that point . in short , this example showed that the initial fluid movement was mostly in the direction of ct movement . when the sensor reached the region marked “ i ” in fig6 , there was little change in the value of bottomhole temperature , which was indicated by “ flats ” in the bottomhole temperature profile . this arrested rate of change of bottomhole temperature indicated that the majority of the region marked under “ i ” was at an identical temperature ; the expanse of this region is easily seen from the difference in ct depth value from its curve . this lead to the first interpretation that sufficient fluid had penetrated this region to keep its temperature near constant over a length of nearly 75 ft [ 22 . 9 m ] from 6175 ft to 6250 ft [ 1882 m to 1905 m ]. in short , when looking at bottomhole temperature curves during acid stimulation treatment using point temperature measurement , one should try to identify “ flats ” or areas where bottomhole temperature shows little change with ct movement . in fig6 , a look at bottomhole temperature profile immediately after region i suggested that as the sensor moved away from the previously encountered “ colder ” region , it started experiencing slightly warmer temperatures ; the rate of change of bottomhole temperature had gained a positive value indicating a region where fluids may not have invaded . however , since the injection was continuously progressing , the fluid could go in the direction that offered lowest resistance . this may have been in the region that had been left “ behind ” the ct tip , the region “ ahead ” or both . for example , if there were no permeable zones after region i , the bottomhole temperature would have continued to increase , although now the direction of fluid flow would have been opposite to ct movement because there would have been fewer favorable zones ahead of the ct . in such cases , as the tip moved further away from the receptive zone that was left “ behind ”, higher annular friction may have been encountered for fluid which had to traverse the greater distance . this change in bottomhole pressure could have been detected by monitoring the bottomhole pressure curve , which may be plotted alongside . in this example though , the recurrence of bottomhole temperature “ flats ” ( rate of change of temperature close to zero ) indicated that there were other regions that had cooled down as a result of acid fluid invasion , and arrested the rate at which temperature was increasing before the sensor crossed those regions . fig6 shows an increase of nearly 2 ° f . [ 1 . 1 ° c .] to region ii and around 0 . 5 ° f . [ 0 . 28 ° c .] to the first part of region iii . note that the initial bottomhole temperature encountered in region iii was less than the preceding temperature , indicating a “ cool down ”. in this example , data for the acid job presented in example 1 is used to illustrate the use of temperature derivative plots for interpretation in accordance with a method of the invention . fig7 shows the temperature derivative curve ( lower curve ) which distinctly shows the regions where rate of change of bottomhole temperature was near zero . this provided a better indication of quantifying the extent of fluid taking regions , rather than getting an estimate from the bottomhole temperature curve alone . as is evident from a comparison of fig6 and 7 , the temperature derivative curve was able to “ split ” the larger region estimated between 6175 and 6250 ft [ 1882 m to 1905 m ] into several smaller regions . there were also a few other regions visible that were not clearly evident when using bottomhole temperature plot alone . therefore , the temperature derivative curve generated using t * dt / dt and d * dt / dd vs . time ( or any ( t + dt )/ dt where t = temperature , t = time , d = ct depth allowed much more accurate interpretation . smoothing of the curve , as is seen in plot of fig7 was performed by use of a standard , readily available algorithm . in this example data for the acid job presented in examples 1 and 2 was used to illustrate how fluid invasion can be quantified . the solid bars shown in fig7 represent the degree of fluid invasion across the various zones . based on the nature of the slope of the derivative in the identified “ zones ”, this method of the invention determined and assigned the degree of invasion of fluid and represented the same in graphic format ; fig7 shows them as “ bars ” of varying dimensions based on perceived effectiveness of stimulation . the method estimated the degree of invasion by taking into account the angle of separation from a base line of 0 degrees ; with the degree of invasion diminishing as the angle approaches 90 degrees . this example demonstrates a method to compute pre - job skin on - the - fly based on darcy &# 39 ; s equation . pre - stimulation treatment skin may be determined during the initial pass in this diagnostic method where an inert fluid is injected into the formation . some of the inputs required for the calculation , i . e ., pressure drop , rate of injection , height of pay ( or region of invasion ), volume factor , fluid viscosity , and the like are known values . unknowns are reservoir pressure and an estimated value of permeability , which may be obtained from the client . any change in the skin factor during the matrix acidizing treatment may then be computed with a better knowledge of fluid invasion profiles . in the acid job described in example 1 , locations in the reservoir section of the wellbore were visited multiple times ( fig5 ). this data may be used to create temperature history for various sections of the reservoir . the rate of change in temperature at any location may be correlated with fluid invasion in the zone . therefore , if the derivative and bottomhole temperature plots generated during various phases of the treatment are plotted together vs . depth along the wellbore , then the zones which show the most rapid change in temperature can be identified . although specific embodiments of the invention have been disclosed herein in some detail , this has been done solely for the purposes of describing various features and aspects of the invention , and is not intended to be limiting with respect to the scope of the invention . it is contemplated that various substitutions , alterations , and / or modifications , including but not limited to those implementation variations which may have been suggested herein , may be made to the disclosed embodiments without departing from the spirit and scope of the invention as defined by the appended claims which follow .