thanks Itag,
what i did, is
1-initialization match, by matching the RFT and water saturation for the first drilled well, to be sure that my contacts depths, PC curves , fluids densities and pressure at datum initially are correct.
2-i run the model under total fluid control mode, and i matched the pressure by some sensitivities in aquifer, transmissibilities and mainly in pore volume till i get acceptable RFT match for rest of wells.
3-i am working on Wcut match, previously i was changing the RP data manually, as i always has set of curves for each SATNUM region and i averag them to get one curve, i move this average up or down to match water production but from now i will do this by varying the corey exponent as vinomarky said it looks more scientific. and also changing the KRW endpoint is good idea also.

you mean grain size coarsening? coarsening upward or downward identified by many methods, simplest one by plotting gamma ray and its mirror image in opposite scale directions, if we have cone shape so there is coarsening towards the cone base, am i right?
dear Itag, i don't understand the part related to Pseudo functions and fractional flow match?
could you please give more explanations and upload some documentations related to these items?
thanks,