mahmood19227
/
SimulationCore2
1
0
Fork 0
Code Issues 0 Pull Requests 0 Releases 0 Wiki Activity
Simulation Core
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
69 Commits
2 Branches
134 MiB
 
 
 
 
 
 
Branch: master
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'master'
${ noResults }
SimulationCore2/FluidFlow/Well_Pressure_Data_Transfer...

732 lines
36 KiB
Raw Permalink Blame History 

  1. # 1 "/home/admin/SimulationCore2/FluidFlow/Well_Pressure_Data_Transfer.f90"

  2. SUBROUTINE WellPressureDataTransfer

  3.     

  4. !! This subroutine calculates pressure distribution in well

  5. !! We divide well hydraulic system to 4 zones. 1- pump to string zone; this zone is horizontal

  6. !! 2- string zone; this zone includes drill pipes and collars and is vertical or neraly vertical

  7. !! 3- annulus zone; this zone includes annulus space between string and openhole or casing

  8. !! 4- Openhole zone; this zone includes spaces under bit to bottom of the wellbore

  9. !! The last two zone includes annulus space and open hole, are influenced by kicks and have crucial role in

  10. !! computing pressure distribution because reference point (back pressure in well head open condition

  11. !! and bottomhole pressure in well head closed condition) are in one of them. because of this, pressure

  12. !! distribution calculates in two subroutines : PressureAnnAndOHDistribution and PressurePumptoBitDistribution

  13. !! pressure change through well consist of pressure change due to friction, change of elevation and acceleration

  14. !! frictional pressure drop calculated by subroutine Sub_Press_Drop_Calculator

  15. !! This subroutine includes hydrostatic pressure and later will include effect of corss-section area change

  16. !! which causes acceleration or decceleration of flow

  17.     

  18.     USE FricPressDropVarsModule

  19.     use MudSystemModule

  20.     USE MudSystemVARIABLES

  21.     use SimulationVariables !@@@

  22.     USE GeoElements_FluidModule

  23.     USE Fluid_Flow_Startup_Vars

  24.     use KickVARIABLESModule

  25. ! USE CMudPropertiesVariables

  26. ! USE CBopStackVariables

  27. ! USE sROP_Variables

  28. ! USE CHOKEVARIABLES

  29. ! use CChokeManifoldVariables

  30. ! use SimulationVariables

  31. ! USE CBitProblemsVariables

  32. ! USE CChokeProblemsVariables

  33. ! USE CStringConfigurationVariables

  34. ! use CDataDisplayConsole

  35.     

  36.     IMPLICIT NONE

  37.     

  38.     INTEGER :: i , j

  39.     REAL :: InstantaneousTotalOpenChokeArea , VolumeDensityProduct

  40.     

  41.     ChokeIsClosing = .FALSE.

  42.     data%State%FricPressDrop%Kchoke = 0.0

  43.     data%State%FricPressDrop%KBit = 0.0

  44.     

  45.     

  46. !!!!!!!!!!!!!!!! bit is present

  47.     BitTrue = data%Configuration%StringConfiguration%StringConfigurations(1)%ComponentType  ==  Bit_ComponentType           ! Bit_ComponentType = 0

  48.     BitCd = 0.98

  49.     BitNozzDia = 32.0 * data%Configuration%StringConfiguration%BitDefinition%BitNozzleSize       ! nozzle diameter in 1/32 in

  50.     BitNozzleArea = 7.6699E-4 * BitNozzDia**2

  51.     BitNozzleNum = data%Configuration%StringConfiguration%BitDefinition%BitNozzleNo - data%State%FricPressDrop%BitJetsPlugged * data%problems%BitProblems%PlugJetsCount

  52.     BitTotNozzArea = BitNozzleNum * BitNozzleArea + 0.5 * data%State%FricPressDrop%BitJetsWashedOut * data%problems%BitProblems%JetWashoutCount * BitNozzleArea

  53.     IF (BitNozzleNum == 0) THEN

  54.         data%State%FricPressDrop%BitTotallyPluged = .TRUE.

  55.     ELSE

  56.         data%State%FricPressDrop%BitTotallyPluged = .FALSE.

  57.     END IF

  58.     

  59.     data%State%FricPressDrop%FloatValveIn = data%Configuration%StringConfiguration%BitDefinition%FloatValve

  60. !FloatValveOpen = .TRUE.

  61. !IF (JetWashoutCount > 0 .OR. PlugJetsCount > 0) THEN

  62. !    WRITE (*,*) 'PlugJetsCount' , PlugJetsCount

  63. !    WRITE (*,*) 'JetWashoutCount', JetWashoutCount

  64. !END IF

  65.     

  66. !WRITE (*,*) 'BitNozzleNum', BitNozzleNum

  67. !WRITE (*,*) 'BitTotNozzArea', BitTotNozzArea

  68. !WRITE (*,*) 'Float Valve??', BitDefinition%FloatValve

  69.     

  70. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  71.     

  72.     IF(ALLOCATED(GasPocketFlowEl))      DEALLOCATE(GasPocketFlowEl)

  73.     

  74. !!!!!!!!!!!!!!!!!!!!  Well Head Condition ( Open or Closed )

  75.     ChokeKroneckerDelta = 0

  76. 

  77.     

  78.     IF (data%State%MudSystem%WellToPitsOpen == .FALSE. .AND. data%State%MudSystem%WellToChokeManifoldOpen == .FALSE.) THEN

  79.         

  80.         WellHeadOpen = .FALSE.

  81. !WRITE (*,*) ' Well Head is closed '

  82.     ELSE 

  83.         WellHeadOpen = .TRUE.

  84.     END IF

  85.     

  86.     IF (data%State%MudSystem%WellToPitsOpen == .FALSE. .AND. data%State%MudSystem%WellToChokeManifoldOpen)        ChokeKroneckerDelta = 1

  87. !WRITE (*,*) ' WelltoPitsOpen=' , WelltoPitsOpen, 'WellToChokeManifoldOpen= ' , WellToChokeManifoldOpen

  88.     

  89. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  90.     

  91.     

  92.     data%State%FricPressDrop%NoHorizontalEl = data%State%MudSystem%NoHorizontalMudElements

  93.     data%State%FricPressDrop%NoStringEl = data%State%MudSystem%NoStringMudElements

  94.     data%State%FricPressDrop%NoAnnulusEl = data%State%MudSystem%NoCasingMudElements

  95.     data%State%FricPressDrop%NoWellToChokeEl = 0

  96.     IF (ChokeKroneckerDelta == 1)       data%State%FricPressDrop%NoWellToChokeEl = data%State%MudSystem%ChokeLine_Density%Length()

  97.     data%State%FricPressDrop%NoOpenHoleEl = data%State%MudSystem%NoBottomHoleMudElements

  98.     data%State%FricPressDrop%NumbEl = data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + data%State%FricPressDrop%NoAnnulusEl + data%State%FricPressDrop%NoWellToChokeEl + data%State%FricPressDrop%NoOpenHoleEl

  99.     IF(ALLOCATED(FlowEl)) then

  100.         DEALLOCATE(FlowEl)

  101.     endif

  102.     ALLOCATE(FlowEl(data%State%FricPressDrop%NumbEl))

  103.     

  104. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  Transfering Data from module:MudSystemVARIABLES  to  module:FricPressDropVars to calculate Pressure distribution in flow path

  105. !!!!!!!!!!!!!!!!!!!!!!! Horizontal line

  106. !WRITE (*,*) 'H', NoHorizontalEl

  107.     DO i = 1 , data%State%FricPressDrop%NoHorizontalEl

  108. ! WRITE (*,*) 'H , i' , i

  109.         FlowEl(i)%StartX = data%State%MudSystem%Xstart_MudElement%Array(i) 

  110. !WRITE (*,*) 'StartX', FlowEl(i)%StartX

  111.         FlowEl(i)%EndX = data%State%MudSystem%Xend_MudElement%Array(i)

  112. !WRITE (*,*) 'Endx', FlowEl(i)%EndX

  113.         FlowEl(i)%StartTVD = 0.d0

  114.         FlowEl(i)%EndTVD = 0.d0

  115.         FlowEl(i)%Id = 0.0

  116.         FlowEl(i)%Od = data%State%MudSystem%PipeOD_MudElement%Array(i)

  117.         FlowEl(i)%density = data%State%MudSystem%Density_MudElement%Array(i)

  118.         

  119.         FlowEl(i)%Dhyd = FlowEl(i)%Od

  120.         FlowEl(i)%Length = ABS(FlowEl(i)%EndX - FlowEl(i)%StartX)

  121.         FlowEl(i)%DepthDiff = 0.d0

  122.         FlowEl(i)%Volume = PI / 4.0 * FlowEl(i)%Od**2 * REAL(FlowEl(i)%Length) / Convfttoinch**2

  123.         IF (data%State%MudSystem%MudType_MudElement%Array(i) == 0) THEN           ! = 0 for mud, = 1 for gas kick, = 4 for air

  124.             FlowEl(i)%MaterialType = 0

  125.         ELSE IF (data%State%MudSystem%MudType_MudElement%Array(i) == 4) THEN

  126.             FlowEl(i)%MaterialType = 4

  127.         ELSE

  128.             FlowEl(i)%MaterialType = 1

  129.         END IF

  130.         

  131.         FlowEl(i)%dPdLgrav = 0.0

  132.         FlowEl(i)%StaticPressDiff = 0.0

  133. 

  134. 

  135.         

  136.     END DO

  137. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  138.     

  139. !!!!!!!!!!!!!!!!!!!!!!!!! String

  140.     data%State%FricPressDrop%StringFirstEl = data%State%FricPressDrop%NoHorizontalEl + 1

  141.     data%State%FricPressDrop%StringLastEl = data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl

  142.     DO i = data%State%FricPressDrop%NoHorizontalEl + 1 , data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl

  143. ! WRITE (*,*) 'ST , i' , i

  144.         FlowEl(i)%StartX = data%State%MudSystem%Xstart_MudElement%Array(i)

  145. !WRITE (*,*) 'StartX', FlowEl(i)%StartX

  146.         FlowEl(i)%EndX = data%State%MudSystem%Xend_MudElement%Array(i)

  147. !WRITE (*,*) 'Endx', FlowEl(i)%EndX

  148.         FlowEl(i)%StartTVD = data%State%MudSystem%TVDstart_MudElement%Array(i)

  149. !WRITE (*,*) 'StartTVD', FlowEl(i)%StartTVD

  150.         FlowEl(i)%EndTVD = data%State%MudSystem%TVDend_MudElement%Array(i)

  151. !WRITE (*,*) 'EndTVD', FlowEl(i)%EndTVD

  152.         FlowEl(i)%Id = 0.0

  153. !WRITE (*,*) 'Id', FlowEl(i)%Id

  154.         FlowEl(i)%Od = data%State%MudSystem%PipeOD_MudElement%Array(i)

  155. !WRITE (*,*) 'Od', FlowEl(i)%Od

  156.         FlowEl(i)%Density = data%State%MudSystem%Density_MudElement%Array(i)

  157. !WRITE (*,*) 'Density', FlowEl(i)%Density

  158.         

  159.         

  160.         FlowEl(i)%Dhyd = FlowEl(i)%Od

  161.         FlowEl(i)%Length = ABS(FlowEl(i)%EndX - FlowEl(i)%StartX)

  162. !WRITE (*,*) 'Length', FlowEl(i)%Length

  163.         FlowEl(i)%DepthDiff = ABS(FlowEl(i)%StartTVD - FlowEl(i)%EndTVD)

  164. !WRITE (*,*) 'DepthDiff', FlowEl(i)%DepthDiff

  165.         FlowEl(i)%Area = PI / 4.0 * FlowEl(i)%Od**2 / Convfttoinch**2

  166.         FlowEl(i)%Volume = FlowEl(i)%Area * REAL(FlowEl(i)%Length)

  167.         IF (data%State%MudSystem%MudType_MudElement%Array(i) == 0) THEN           ! = 0 for mud, = 1 for gas kick, = 4 for air

  168.             FlowEl(i)%MaterialType = 0

  169.         ELSE IF (data%State%MudSystem%MudType_MudElement%Array(i) == 4) THEN

  170.             FlowEl(i)%MaterialType = 4

  171.         ELSE

  172.             FlowEl(i)%MaterialType = 1

  173.         END IF

  174.         

  175.         IF (FlowEl(i)%MaterialType == 1 .OR. FlowEl(i)%MaterialType == 4) THEN

  176.             FlowEl(i)%dPdLgrav = 0.0

  177.         ELSE

  178.             FlowEl(i)%dPdLgrav = 0.052 * FlowEl(i)%Density

  179.         END IF

  180.         FlowEl(i)%StaticPressDiff = FlowEl(i)%dPdLgrav * REAL(FlowEl(i)%DepthDiff)

  181.         

  182. 

  183. 

  184.         

  185.     END DO

  186. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  187.     

  188. !!!!!!!!!!!!!!!!!!!!!!!!!!!! Annulus

  189.     data%State%FricPressDrop%AnnulusFirstEl = data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + 1

  190.     data%State%FricPressDrop%AnnulusLastEl = data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + data%State%FricPressDrop%NoAnnulusEl

  191.     DO i = data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + 1 , data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + data%State%FricPressDrop%NoAnnulusEl

  192. ! WRITE (*,*) 'AN , i' , i

  193.         FlowEl(i)%StartX = data%State%MudSystem%Xstart_MudElement%Array(i)

  194. !WRITE (*,*) 'StartX', FlowEl(i)%StartX

  195.         FlowEl(i)%EndX = data%State%MudSystem%Xend_MudElement%Array(i)

  196. !WRITE (*,*) 'Endx', FlowEl(i)%EndX

  197.         FlowEl(i)%StartTVD = data%State%MudSystem%TVDstart_MudElement%Array(i)

  198.         FlowEl(i)%EndTVD = data%State%MudSystem%TVDend_MudElement%Array(i)

  199.         FlowEl(i)%Id = data%State%MudSystem%PipeID_MudElement%Array(i)

  200.         FlowEl(i)%Od = data%State%MudSystem%PipeOD_MudElement%Array(i)

  201.         FlowEl(i)%Density = data%State%MudSystem%Density_MudElement%Array(i)        

  202.         

  203.         FlowEl(i)%Dhyd = FlowEl(i)%Od - FlowEl(i)%Id

  204.         FlowEl(i)%Length = ABS(FlowEl(i)%EndX - FlowEl(i)%StartX)

  205.         FlowEl(i)%DepthDiff = ABS(FlowEl(i)%StartTVD - FlowEl(i)%EndTVD)

  206.         FlowEl(i)%Area = PI / 4.0 * (FlowEl(i)%Od**2 - FlowEl(i)%Id**2) / Convfttoinch**2

  207.         FlowEl(i)%Volume = FlowEl(i)%Area * REAL(FlowEl(i)%Length)

  208.         IF (data%State%MudSystem%MudType_MudElement%Array(i) == 0) THEN           ! = 0 for mud, = 1 for gas kick, = 4 for air

  209.             FlowEl(i)%MaterialType = 0

  210.         ELSE IF (data%State%MudSystem%MudType_MudElement%Array(i) == 4) THEN

  211.             FlowEl(i)%MaterialType = 4

  212.         ELSE

  213.             FlowEl(i)%MaterialType = 1

  214.         END IF

  215.         

  216.         IF (FlowEl(i)%MaterialType == 1 .OR. FlowEl(i)%MaterialType == 4) THEN

  217.             FlowEl(i)%dPdLgrav = 0.0

  218.         ELSE

  219.             FlowEl(i)%dPdLgrav = 0.052 * FlowEl(i)%density

  220.         END IF

  221.         FlowEl(i)%StaticPressDiff = FlowEl(i)%dPdLgrav * REAL(FlowEl(i)%DepthDiff)

  222. 

  223. 

  224.         

  225.     END DO

  226. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  227.     

  228. !!!!!!!!!!!!!!! Well to choke manifold path

  229.     j = 1

  230.     DO i = data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + data%State%FricPressDrop%NoAnnulusEl + 1 , data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + data%State%FricPressDrop%NoAnnulusEl + data%State%FricPressDrop%NoWelltoChokeEl

  231. ! WRITE (*,*) 'CH , i' , i

  232.         FlowEl(i)%StartX = data%State%MudSystem%ChokeLine_Mud_Backhead_X%Array(j)

  233. !WRITE (*,*) 'StartX', FlowEl(i)%StartX

  234.         FlowEl(i)%EndX = data%State%MudSystem%ChokeLine_Mud_Forehead_X%Array(j)

  235. !WRITE (*,*) 'Endx', FlowEl(i)%EndX

  236.         FlowEl(i)%StartTVD = 0.d0

  237.         FlowEl(i)%EndTVD = 0.d0

  238.         FlowEl(i)%Id = 0.0

  239.         FlowEl(i)%Od = data%Configuration%BopStack%ChokeLineId

  240.         FlowEl(i)%density = data%State%MudSystem%ChokeLine_Density%Array(j)

  241.             

  242.         FlowEl(i)%Dhyd = FlowEl(i)%Od

  243.         FlowEl(i)%Length = ABS(FlowEl(i)%EndX - FlowEl(i)%StartX)

  244.         FlowEl(i)%DepthDiff = 0.d0

  245.         FlowEl(i)%Area = PI / 4.0 * FlowEl(i)%Od**2 / Convfttoinch**2

  246.         FlowEl(i)%Volume = FlowEl(i)%Area * REAL(FlowEl(i)%Length)

  247.         IF (data%State%MudSystem%ChokeLine_MudOrKick%Array(j) == 0 .OR. data%State%MudSystem%ChokeLine_MudOrKick%Array(j) == 4) THEN      ! = 0 for mud, = 1 for gas kick, = 4 for air

  248.             FlowEl(i)%MaterialType = data%State%MudSystem%ChokeLine_MudOrKick%Array(j)

  249.         ELSE

  250.             FlowEl(i)%MaterialType = 1

  251.         END IF

  252.         

  253.             

  254.         FlowEl(i)%dPdLgrav = 0.0

  255.         FlowEl(i)%StaticPressDiff = 0.0

  256.         

  257. !WRITE (*,*) 'Ch density, length, Type' , i, FlowEl(i)%density , FlowEl(i)%length, FlowEl(i)%MaterialType

  258. 

  259.         j = j + 1

  260.             

  261.     END DO

  262.             

  263. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  264.     

  265. !!!!!!!!!!!!!!! Open Hole

  266.     data%State%FricPressDrop%OpenholeFirstEl = data%State%FricPressDrop%NoHorizontalEl + data%State%FricPressDrop%NoStringEl + data%State%FricPressDrop%NoAnnulusEl + data%State%FricPressDrop%NoWelltoChokeEl + 1

  267.     j = 1

  268.     DO i = data%State%FricPressDrop%OpenholeFirstEl , data%State%FricPressDrop%NumbEl ! = NoHorizontalEl + NoStringEl + NoCasingEl + NoOpenHoleEl

  269. ! WRITE (*,*) 'OP , i' , i

  270.         FlowEl(i)%StartX = data%State%MudSystem%Xstart_OpMudElement%Array(j)

  271. !WRITE (*,*) 'StartX', FlowEl(i)%StartX

  272.         FlowEl(i)%EndX = data%State%MudSystem%Xend_OpMudElement%Array(j)

  273. !WRITE (*,*) 'Endx', FlowEl(i)%EndX

  274.         FlowEl(i)%StartTVD = data%State%MudSystem%TVDstart_OpMudElement%Array(j)

  275.         FlowEl(i)%EndTVD = data%State%MudSystem%TVDend_OpMudElement%Array(j)

  276.         FlowEl(i)%Id = 0.0

  277.         FlowEl(i)%Od = data%State%MudSystem%PipeOD_OpMudElement%Array(j)

  278.         FlowEl(i)%density = data%State%MudSystem%Density_OpMudElement%Array(j)

  279.         FlowEl(i)%MaterialType = data%State%MudSystem%MudTypeOp_MudElement%Array(j)     ! = 0 for mud, = 1 for gas kick, = 4 for air

  280.         

  281.         FlowEl(i)%Dhyd = FlowEl(i)%Od

  282.         FlowEl(i)%Length = ABS(FlowEl(i)%EndX - FlowEl(i)%StartX)

  283.         FlowEl(i)%DepthDiff = ABS(FlowEl(i)%StartTVD - FlowEl(i)%EndTVD)

  284.         FlowEl(i)%Area = PI / 4.0 * FlowEl(i)%Od**2 / Convfttoinch**2        

  285.         FlowEl(i)%Volume = FlowEl(i)%Area * REAL(FlowEl(i)%Length)

  286.         IF (data%State%MudSystem%MudTypeOp_MudElement%Array(j) == 0) THEN           ! = 0 for mud, = 1 for gas kick, = 4 for air

  287.             FlowEl(i)%MaterialType = 0

  288.         ELSE IF (data%State%MudSystem%MudTypeOp_MudElement%Array(j) == 4) THEN

  289.             FlowEl(i)%MaterialType = 4

  290.         ELSE

  291.             FlowEl(i)%MaterialType = 1

  292.         END IF

  293.         

  294.         j = j + 1

  295.         IF (FlowEl(i)%MaterialType == 1) THEN

  296.             FlowEl(i)%dPdLgrav = 0.0

  297.         ELSE

  298.             FlowEl(i)%dPdLgrav = 0.052 * FlowEl(i)%density

  299.         END IF

  300.         FlowEl(i)%StaticPressDiff = FlowEl(i)%dPdLgrav * REAL(FlowEl(i)%DepthDiff)

  301.     

  302. !WRITE (*,*) 'O density, length, Type of ith element' , i, FlowEl(i)%density , FlowEl(i)%length, FlowEl(i)%MaterialType

  303. 

  304.     END DO

  305. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  306.     

  307. !!!!!!!!!!!!!!!!!!  Mud Element Viscosity Calculation

  308.     

  309.     DO i = 1 , data%State%FricPressDrop%NumbEl

  310.         IF (FlowEl(i)%MaterialType /= 0)    CYCLE

  311.         IF (FlowEl(i)%density < 8.33) THEN

  312.             WRITE (*,*) 'H, S, A, Ch, O', data%State%FricPressDrop%NoHorizontalEl , data%State%FricPressDrop%StringFirstEl , data%State%FricPressDrop%AnnulusFirstEl , data%State%FricPressDrop%NoWellToChokeEl , data%State%FricPressDrop%OpenHoleFirstEl

  313.             

  314.             WRITE (*,*) 'Element No, Volume(gal), Density(ppg), Type'

  315.             DO j = 1 , data%State%FricPressDrop%NumbEl

  316.                 WRITE (*,*) j, FlowEl(j)%Volume * ConvFt3toUSGal , FlowEl(j)%Density, FlowEl(j)%MaterialType

  317.             END DO

  318.         

  319.             CALL ErrorStop('Density below water density (8.33)', i)

  320.             CYCLE

  321.         END IF

  322.         

  323.         FlowEl(i)%Theta600 = Theta600Refrence + 11.0 * (MAX(FlowEl(i)%density , 8.33) - DensityRefrence)

  324.         FlowEl(i)%Theta300 = Theta300Refrence + 6.0 * (MAX(FlowEl(i)%density , 8.33) - DensityRefrence)

  325. !WRITE (*,*) 'i, Theta600, Theta600', i, FlowEl(i)%Theta600, FlowEl(i)%Theta300

  326.     END DO

  327. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  328. 

  329.     

  330.     FlowEl(:)%FlowRate = 0.0

  331.     FlowEl(:)%StartPress = 0.0

  332.     FlowEl(:)%EndPress = 0.0

  333.     

  334. !WRITE (*,*) 'NoGasPocket=' , SIZE(GasPocketWeight%Array)

  335.     

  336.     

  337.     

  338. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  339.     

  340. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  341.     

  342.     data%State%MudSystem%OpenArea32 = ((data%Equipments%ChokeManifold%LeftManualChoke * (1.0 - data%State%FricPressDrop%ManChoke1Plug * REAL(data%problems%ChokeProblems%ManualChoke1PluggedPercent)) / 100.0) + 0.5 * data%State%FricPressDrop%ManChoke1Washout) * data%State%Choke%ChokeAreaFullyOpen

  343.     data%State%MudSystem%OpenArea33 = data%State%CHOOKE(1)%AreaChokeFinal * Convfttoinch**2              !  (1.d0 - data%State%CHOOKE(1)%PercentClose) * ChokeAreaFullyOpen

  344.     data%State%MudSystem%OpenArea34 = data%State%CHOOKE(2)%AreaChokeFinal * Convfttoinch**2              !  (1.d0 - data%State%CHOOKE(2)%PercentClose) * ChokeAreaFullyOpen

  345.     data%State%MudSystem%OpenArea35 = ((data%Equipments%ChokeManifold%RightManualChoke * (1.0 - data%State%FricPressDrop%ManChoke2Plug * REAL(data%problems%ChokeProblems%ManualChoke2PluggedPercent)) / 100.0) + 0.5 * data%State%FricPressDrop%ManChoke2Washout) * data%State%Choke%ChokeAreaFullyOpen

  346. !write(*,*) 'OpenArea32=' , OpenArea32, active32, ManChoke1Plug, ManualChoke1PluggedPercent, ManChoke1Washout

  347. !write(*,*) 'OpenArea33=' , OpenArea33, active33

  348. !write(*,*) 'OpenArea34=' , OpenArea34, active34, !HydraulicChoke2PluggedPercent

  349. !write(*,*) 'OpenArea35=' , OpenArea35, active35, ManChoke2Plug, ManualChoke2PluggedPercent, ManChoke2Washout

  350.     InstantaneousTotalOpenChokeArea = data%State%MudSystem%OpenArea32 * active32 + data%State%MudSystem%OpenArea33 * active33 + data%State%MudSystem%OpenArea34 * active34 + data%State%MudSystem%OpenArea35 * active35 + data%State%FricPressDrop%ChokeBypassArea * active29

  351.     data%State%FricPressDrop%OldTotalOpenChokeArea = data%State%FricPressDrop%TotalOpenChokeArea

  352. !WRITE (*,*) 'Instantaneous / Old TotalOpenChokeArea', InstantaneousTotalOpenChokeArea, OldTotalOpenChokeArea

  353.     IF (data%State%FricPressDrop%OldTotalOpenChokeArea <= 0.01 * data%State%Choke%ChokeAreaFullyOpen)     data%State%FricPressDrop%OldTotalOpenChokeArea = 0.01 * data%State%Choke%ChokeAreaFullyOpen

  354.     IF (InstantaneousTotalOpenChokeArea <= 0.01 * data%State%Choke%ChokeAreaFullyOpen) THEN

  355.         data%State%MudSystem%WellToChokeManifoldOpen = .FALSE.

  356.         data%State%FricPressDrop%OldTotalOpenChokeArea = 0.01 * data%State%Choke%ChokeAreaFullyOpen

  357. !WRITE (*,*) ' Choke Controler Here 2'

  358.         

  359.     ELSE IF (InstantaneousTotalOpenChokeArea > 0.5 * data%State%Choke%ChokeAreaFullyOpen .OR. data%State%MudSystem%WellToPitsOpen) THEN

  360. ! mud flows through well to bell nipple, or choke is rather open

  361.         data%State%FricPressDrop%TotalOpenChokeArea = InstantaneousTotalOpenChokeArea

  362. !WRITE (*,*) ' Choke Controler Here 1'

  363.     ELSE

  364.         IF (InstantaneousTotalOpenChokeArea > data%State%FricPressDrop%OldTotalOpenChokeArea) THEN

  365.             data%State%FricPressDrop%TotalOpenChokeArea = MIN(1.1 * data%State%FricPressDrop%OldTotalOpenChokeArea , InstantaneousTotalOpenChokeArea)

  366.             WRITE (*,*) ' Choke is opening' !, TotalOpenChokeArea, OldTotalOpenChokeArea, InstantaneousTotalOpenChokeArea

  367.         ELSE IF (InstantaneousTotalOpenChokeArea < data%State%FricPressDrop%OldTotalOpenChokeArea) THEN

  368.             data%State%FricPressDrop%TotalOpenChokeArea = MAX(0.9 * data%State%FricPressDrop%OldTotalOpenChokeArea , InstantaneousTotalOpenChokeArea)

  369.             WRITE (*,*) ' Choke is closing' !, TotalOpenChokeArea, OldTotalOpenChokeArea, InstantaneousTotalOpenChokeArea

  370.         ELSE

  371.             data%State%FricPressDrop%TotalOpenChokeArea = InstantaneousTotalOpenChokeArea

  372.         END IF

  373.     END IF

  374. !WRITE (*,*) 'TotalOpenChokeArea=', TotalOpenChokeArea

  375.     

  376.         

  377.     

  378.     

  379.             

  380. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  381.     

  382. !!!!!!!!!!!!!!!!!!!!  Well Head Condition ( Open or Closed )

  383.     ChokeKroneckerDelta = 0

  384. !WRITE (*,*) ' WelltoPitsOpen = ', WelltoPitsOpen

  385. !WRITE (*,*) ' WellToChokeManifoldOpen = ', WellToChokeManifoldOpen

  386.     IF (data%State%MudSystem%WellToPitsOpen == .FALSE. .AND. data%State%MudSystem%WellToChokeManifoldOpen == .FALSE.) THEN

  387. !WRITE (*,*) ' WelltoPitsOpen = ', WelltoPitsOpen

  388. !WRITE (*,*) ' WellToChokeManifoldOpen = ', WellToChokeManifoldOpen

  389.         

  390.         WellHeadOpen = .FALSE.

  391. !WRITE (*,*) ' Well Head is closed '

  392.     ELSE

  393.         WellHeadOpen = .TRUE.

  394.     END IF

  395.     

  396.     IF (data%State%MudSystem%WellToPitsOpen == .FALSE. .AND. data%State%MudSystem%WellToChokeManifoldOpen)        ChokeKroneckerDelta = 1

  397. !WRITE (*,*) ' WelltoPitsOpen=' , WelltoPitsOpen, 'WellToChokeManifoldOpen= ' , WellToChokeManifoldOpen

  398. !WRITE (*,*) 'WellToChokeOpen , WellToChokeWasOpen', WellToChokeManifoldOpen, WellToChokeManifoldWasOpen

  399. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  400. 

  401.     IF (NoGasPocket > 0) THEN

  402.         IF (data%State%MudSystem%Choke_Kick_Saved_Volume_Final > 0.d0) THEN

  403. !WRITE (*,*) 'Choke_Kick_Saved_Volume_Final=', Choke_Kick_Saved_Volume_Final, (GasPocketModifiedVol%Array(NoGasPocket) * Convft3ToUSgal)

  404.             ExitMass = REAL(data%State%MudSystem%Choke_Kick_Saved_Volume_Final / (GasPocketModifiedVol%Array(NoGasPocket) * Convft3ToUSgal)) * GasPocketWeight%Array(NoGasPocket)

  405. !IF (NoGasPocket > 1)        ExitMass = ExitMass + (SUM(GasPocketDeltaVol%Array(1 : NoGasPocket - 1)) / GasPocketNewVol%Array(NoGasPocket)) * GasPocketWeight%Array(NoGasPocket)

  406. !WRITE (*,*) 'Total Exit mass=' , ExitMass !, DeltaVolumePipe, data%State%MudSystem%StringFlowRate * dt / ConvMinToSec

  407.             IF (FlowEl(data%State%FricPressDrop%OpenholeFirstEl - 1)%MaterialType /= 1) THEN

  408.                 CALL RemoveGasPocket(NoGasPocket)

  409.                 WRITE (*,*) 'Choke_Kick_Saved_Volume_Final=', data%State%MudSystem%Choke_Kick_Saved_Volume_Final

  410.             END IF

  411.             

  412.         END IF

  413.         

  414.         

  415. !NoGasPocket = SIZE(GasPocketWeight%Array)

  416.         IF (NoGasPocket > 0)         CALL GasPocketFlowElementTransformer

  417.     END IF

  418.     

  419.     IF (data%State%FricPressDrop%NoWelltoChokeEl > 0) THEN

  420.         OldChokeDensity = ChokeDensity

  421.         VolumeDensityProduct = 0.0

  422.         Do i = data%State%FricPressDrop%AnnulusLastEl + 1 , data%State%FricPressDrop%OpenholeFirstEl - 1

  423.             VolumeDensityProduct = VolumeDensityProduct + ((MAX(FlowEl(i)%Density , ChokeMinDensity)) * FlowEl(i)%Volume)

  424.         End Do

  425.         ChokeDensity = VolumeDensityProduct / SUM(FlowEl(data%State%FricPressDrop%AnnulusLastEl + 1 : data%State%FricPressDrop%OpenholeFirstEl - 1)%Volume)

  426.                             

  427. !WRITE (*,*) 'ChokeDensity=', ChokeDensity

  428.         data%State%FricPressDrop%Kchoke = (ChokeDensity / ((2.0 * 89158.0) * (0.26 * 0.61 * data%State%FricPressDrop%TotalOpenChokeArea)**2)) * 4.d0    ! *4.d0: seyyed gofte

  429.     END IF

  430.     

  431.     IF (data%State%FricPressDrop%BitTotallyPluged == .FALSE.) THEN

  432.         data%State%FricPressDrop%KBit = FlowEl(data%State%FricPressDrop%StringLastEl)%density / 12042.0 / BitCd**2 / BitTotNozzArea**2

  433.     END IF

  434.     

  435.     

  436.     

  437. END SUBROUTINE WellPressureDataTransfer

  438.     

  439.     

  440. SUBROUTINE GasPocketFlowElementTransformer

  441. 

  442. !! This Subroutine makes relationship between fluid flow elements and gas pocket elements

  443. !! The main variable of this subroutine is GasPocketFlowEl

  444. !! For each gas pocket, GasPocketFlowEl has a row, which columns are numbers of flow elements which are contains that gas pocket

  445. !! For example, if we have two gas pockets and they are in (7,8) and (12,3,4) flow elements respectively,

  446. !! GasPocketFlowEl(1,:) = [12 , 3 , 4] and GasPocketFlowEl(2,:) = [7 , 8 , 0]

  447. 

  448. 

  449.     USE FricPressDropVarsModule

  450.     USE MudSystemVARIABLES

  451.     use SimulationVariables

  452.     USE GeoElements_FluidModule

  453.     USE Fluid_Flow_Startup_Vars

  454.     use KickVARIABLESModule

  455.     USE CMudPropertiesVariables

  456.     USE CBopStackVariables

  457.     USE CDownHoleVariables

  458.     USE CError

  459.     

  460.     IMPLICIT NONE

  461.     

  462.     INTEGER :: i , j , k

  463.     REAL :: PressureCorrection

  464.     PressureCorrection = 1.0

  465.     

  466.     if (print_log) print *, "NoGasPocket = ",NoGasPocket

  467.     ALLOCATE(GasPocketFlowEl(NoGasPocket , 1))

  468.     GasPocketFlowEl(:,:) = 0

  469. !WRITE (*,*) 'NoGasPocket=' , NoGasPocket

  470.     

  471.     i = 1

  472.     j = 1

  473.     DO k = data%State%FricPressDrop%OpenholeFirstEl , data%State%FricPressDrop%NumbEl

  474.         IF (FlowEl(k)%MaterialType == 1) THEN

  475.             IF (j > 1) THEN

  476.                 IF (GasPocketFlowEl(i , j - 1) /= k - 1) THEN

  477.                     i = i + 1

  478.                     j = 1

  479.                 END IF

  480.             END IF

  481.             

  482.             IF (SIZE(GasPocketFlowEl , dim = 2) < j) THEN

  483.                 if (print_log) print *, "GasPocketFlowEl = ",GasPocketFlowEl

  484.                 ALLOCATE(tempGasPocketFlowEl(SIZE(GasPocketFlowEl , dim = 1) , SIZE(GasPocketFlowEl , dim = 2)))

  485.                 tempGasPocketFlowEl = GasPocketFlowEl

  486.                 DEALLOCATE(GasPocketFlowEl)

  487.                 ALLOCATE(GasPocketFlowEl(SIZE(tempGasPocketFlowEl , dim = 1) , SIZE(tempGasPocketFlowEl , dim = 2) + 1))

  488.                 GasPocketFlowEl(: , 1 : SIZE(tempGasPocketFlowEl , dim = 2)) = tempGasPocketFlowEl

  489.                 GasPocketFlowEl(: , SIZE(GasPocketFlowEl , dim = 2)) = 0

  490.                 DEALLOCATE(tempGasPocketFlowEl)

  491.             END IF

  492.             

  493.             GasPocketFlowEl(i , j) = k

  494.             j = j + 1

  495.         END IF

  496.     END DO

  497. !WRITE (*,*) 'GasPocketFlowEl=' , GasPocketFlowEl

  498.     

  499.     DO k = data%State%FricPressDrop%AnnulusFirstEl , data%State%FricPressDrop%OpenholeFirstEl - 1

  500.         IF (FlowEl(k)%MaterialType == 1) THEN

  501.             IF (j > 1) THEN

  502.                 IF (k == data%State%FricPressDrop%AnnulusFirstEl  .AND.  GasPocketFlowEl(i , j - 1) == data%State%FricPressDrop%NumbEl) THEN

  503. !WRITE (*,*) 'Kick is Around Bit'

  504.                 ELSE IF (GasPocketFlowEl(i , j - 1) /= k - 1) THEN

  505.                     i = i + 1

  506.                     j = 1

  507. !WRITE (*,*) 'i, j, k', i, j, k, FlowEl(k)%MaterialType

  508.                 END IF

  509.             END IF

  510.             

  511.             IF (j > SIZE(GasPocketFlowEl , dim = 2)) THEN

  512.                 if (print_log) print *, "GasPocketFlowEl = ",GasPocketFlowEl

  513.                 ALLOCATE(tempGasPocketFlowEl(SIZE(GasPocketFlowEl , dim = 1) , SIZE(GasPocketFlowEl , dim = 2)))

  514.                 tempGasPocketFlowEl = GasPocketFlowEl

  515.                 DEALLOCATE(GasPocketFlowEl)

  516.                 ALLOCATE(GasPocketFlowEl(SIZE(tempGasPocketFlowEl , dim = 1) , SIZE(tempGasPocketFlowEl , dim = 2) + 1))

  517.                 GasPocketFlowEl(: , 1 : SIZE(tempGasPocketFlowEl , dim = 2)) = tempGasPocketFlowEl

  518.                 GasPocketFlowEl(: , SIZE(GasPocketFlowEl , dim = 2)) = 0

  519.                 DEALLOCATE(tempGasPocketFlowEl)

  520.             END IF

  521.             IF (i > NoGasPocket) THEN

  522.                 WRITE (*,*) 'GasPocketFlowEl', GasPocketFlowEl

  523.                 WRITE (*,*) 'i , j , k', i, j, k

  524.                 WRITE (*,*) 'H, S, A, Ch, O', data%State%FricPressDrop%NoHorizontalEl , data%State%FricPressDrop%StringFirstEl , data%State%FricPressDrop%AnnulusFirstEl , data%State%FricPressDrop%NoWellToChokeEl , data%State%FricPressDrop%OpenHoleFirstEl

  525.             

  526.                 WRITE (*,*) 'Element No, Volume(gal), Density(ppg),         Type'

  527.                 DO j = 1 , data%State%FricPressDrop%NumbEl

  528.                     WRITE (*,*) j, FlowEl(j)%Volume * ConvFt3toUSGal , FlowEl(j)%Density, FlowEl(j)%MaterialType

  529.                 END DO

  530. 

  531.                 CALL ErrorStop('Error in calculating GasPocketFlowEl')

  532.             END IF

  533. 

  534.             GasPocketFlowEl(i , j) = k

  535. !WRITE (*,*) 'GasPocketFlowEl=' , GasPocketFlowEl

  536.             j = j + 1

  537.             

  538.         END IF

  539. !IF (i > NoGasPocket)        CALL ErrorStop('Error in GasPocketFlowEl' , i)

  540.     END DO

  541.     

  542. 

  543. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  544.     

  545. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  546.     GasPocketModifiedVol%Array(:) = 0.d0

  547.     

  548.     DO i = 1 , NoGasPocket

  549. !WRITE(*,*) 'GasPocketFlowEl Data' , i , GasPocketFlowEl(: , i)

  550.         DO j = 1 , SIZE(GasPocketFlowEl , dim = 2)

  551.             IF (GasPocketFlowEl(i , j) > 0)     GasPocketModifiedVol%Array(i) = GasPocketModifiedVol%Array(i) + FlowEl(GasPocketFlowEl(i , j))%Volume

  552.         END DO

  553.     END DO

  554.     i = data%State%FricPressDrop%AnnulusLastEl

  555.     j = data%State%FricPressDrop%OpenholeFirstEl - 1

  556.     

  557.     

  558.     GasPocketOldVol%Array(:) = GasPocketModifiedVol%Array(:)

  559.     

  560.     GasPocketOldPress%Array(:) = GasPocketNewPress%Array(:)

  561. !WRITE (*,*) ' GasPocketweight%Array(1)', GasPocketweight%Array(1)

  562.     IF (GasPocketWeight%Array(1) < 0.0)     CALL ErrorStop(' Error in Calculating Mass ' , KickmdotBCoef)

  563.     

  564.     

  565.     GasPocketDensity%Array(:) = (GasPocketweight%Array(:) / GasPocketModifiedVol%Array(:)) / convft3toUSgal     ! [lbm/ft^3 to ppg]

  566.     DO i = 1 , SIZE(GasPocketFlowEl , dim = 1)

  567.         DO j = 1 , SIZE(GasPocketFlowEl , dim = 2)

  568.             IF (GasPocketFlowEl(i , j) > 0)     FlowEl(GasPocketFlowEl(i , j))%Density = GasPocketDensity%Array(i)

  569.         END DO

  570.     END DO

  571.     

  572.     

  573. !WRITE (*,*) 'Kick density (ppg)=' , GasPocketDensity(1)

  574.     data%Equipments%DownHole%InfluxRate = MAX(((KickmdotACoef * (KickmdotBCoef - GasPocketNewPress%Array(1))) / GasPocketDensity%Array(1) * ConvMinToSec) , 0.0)

  575. !WRITE (*,*) ' InfluxRate (gpm) =', InfluxRate

  576.     

  577.     i = data%State%FricPressDrop%OpenholeFirstEl - 1

  578.     IF (data%State%MudSystem%Choke_Kick_Saved_Volume_Final > 0.d0 .AND. WellHeadOpen) THEN

  579.         WRITE (*,*) 'Kick is exiting through the choke'

  580. !ExitMass = GasPocketDensity%Array(NoGasPocket) * (SUM(GasPocketFlowInduced%Array(:)) + (DeltaVolumePipe * ConvMinToSec / dt) + data%State%MudSystem%StringFlowRate) / ConvMinToSec * dt

  581.         GasPocketWeight%Array(NoGasPocket) = GasPocketweight%Array(NoGasPocket) - ExitMass

  582.         IF (GasPocketWeight%Array(NoGasPocket) > 0.0) THEN

  583.             GasPocketOldPress%Array(NoGasPocket) = GasPocketWeight%Array(NoGasPocket) * data%State%GasType(KickGasType)%GasConstant * &

  584.                                                     GasPocketCompressibility%Array(NoGasPocket) * GasPocketNewTemp%Array(NoGasPocket) / GasPocketModifiedVol%Array(NoGasPocket)

  585.             GasPocketDensity%Array(NoGasPocket) = (GasPocketweight%Array(NoGasPocket) / GasPocketModifiedVol%Array(NoGasPocket)) / convft3toUSgal

  586.             GasPocketNewVol%Array(:) = GasPocketOldVol%Array(:)

  587.         ELSE    ! gas pocket is escaped from the well completely

  588.             WRITE (*,*) ' Last Pocket Removed'

  589.             WRITE (*,*) ' GasPocketFlowEl', GasPocketFlowEl

  590.             CALL RemoveGasPocket(NoGasPocket)

  591.             

  592. !WRITE (*,*) ' GasPocketFlowEl', GasPocketFlowEl

  593.             

  594.         END IF

  595.         

  596.     ELSE IF (NoGasPocket == 1 .OR. WellHeadOpen) THEN        ! kick is not last element of choke line or well head is closed

  597.             GasPocketNewPress%Array(:) = GasPocketNewPress%Array(:) * REAL(GasPocketNewVol%Array(:) / GasPocketOldVol%Array(:))

  598.             GasPocketOldPress%Array(:) = GasPocketNewPress%Array(:)

  599.             GasPocketNewVol%Array(:) = GasPocketOldVol%Array(:)

  600.     

  601.     ELSE    ! Many gas pocket is in the well and well head is closed

  602.         

  603. 103     FORMAT (2I, 4X, (F9.3), 5X, (F9.3))

  604.         

  605.         GasPocketNewPress%Array(NoGasPocket) = GasPocketOldPress%Array(NoGasPocket)

  606.         GasPocketNewVol%Array(NoGasPocket) = GasPocketCompressibility%Array(NoGasPocket) * data%State%GasType(KickGasType)%GasConstant * &

  607.                             GasPocketNewTemp%Array(NoGasPocket) * GasPocketWeight%Array(NoGasPocket) / GasPocketNewPress%Array(NoGasPocket)

  608.         

  609.         

  610.         

  611. !WRITE (*,*) 'Mid'

  612.         DO WHILE (ABS(PressureCorrection * 10.0) > KickConvergenceTolerance)

  613.             DO j = NoGasPocket - 1 , 1 , -1

  614.                 CALL KickFunctionsCalculator(GasPocketNewPress%Array(j) , j , 5)

  615.                 GasPocketNewVol%Array(j) = GasPocketCompressibility%Array(j) * data%State%GasType(KickGasType)%GasConstant * &

  616.                             GasPocketNewTemp%Array(j) * GasPocketWeight%Array(j) / GasPocketNewPress%Array(j)

  617. !WRITE (*,*) j, REAL(GasPocketNewPress%Array(j)), REAL(GasPocketNewVol%Array(j) * Convft3toUSgal)

  618.             END DO

  619.             PressureCorrection = (SUM(GasPocketNewVol%Array(:)) - SUM(GasPocketOldVol%Array(:))) / SUM(GasPocketNewVol%Array(:) / GasPocketNewPress%Array(:))

  620.             GasPocketNewPress%Array(NoGasPocket) = GasPocketNewPress%Array(NoGasPocket) + PressureCorrection

  621.             GasPocketNewVol%Array(NoGasPocket) = GasPocketCompressibility%Array(NoGasPocket) * data%State%GasType(KickGasType)%GasConstant * &

  622.                             GasPocketNewTemp%Array(NoGasPocket) * GasPocketWeight%Array(NoGasPocket) / GasPocketNewPress%Array(NoGasPocket)

  623.             j = NoGasPocket

  624. !WRITE (*,*) j, REAL(GasPocketNewPress%Array(j)), REAL(GasPocketNewVol%Array(j) * Convft3toUSgal)

  625.         END DO

  626.         

  627. !WRITE (*,*) 'Pressure Correction = ', PressureCorrection

  628.         WRITE (*,*) 'After'

  629.         DO i = 1 , NoGasPocket

  630.             WRITE (*,*) i, REAL(GasPocketNewPress%Array(i)), REAL((GasPocketNewVol%Array(i) * Convft3toUSgal))

  631.         END DO

  632.         

  633.         GasPocketOldPress%Array(:) = GasPocketNewPress%Array(:)

  634.         

  635.     END IF

  636.     

  637.     

  638.     IF (GasPocketFlowEl(NoGasPocket , 1) == 0) THEN

  639.         CALL RemoveGasPocket(1)

  640.         

  641. !WRITE (*,*) 'GasPocketFlowEl', GasPocketFlowEl

  642.         WRITE (*,*) 'First Pocket Removed'

  643. !WRITE (*,*) 'GasPocketFlowEl', GasPocketFlowEl

  644.     END IF

  645.     

  646.     

  647.     

  648.     data%Equipments%DownHole%KickVolume = INT(SUM(GasPocketOldVol%Array(:)) * convft3toUSgal / 42. * 10.0) / 10.0

  649.     

  650.     

  651. !WRITE (*,*) ' Gas Kick Pressure (psi) = ' , GasPocketOldPress(1) , INT((GasPocketNewVol(1) / GasPocketOldVol(1)) * 1000.d0) / 1000.d0

  652. !WRITE (*,*) ' b) GasPocketNewVol (ft^3) = ' , GasPocketOldVol(1)

  653.         

  654.         

  655.     

  656.     END SUBROUTINE

  657.     

  658.     

  659. SUBROUTINE RemoveGasPocket(ilocal)

  660.     

  661.     use KickVARIABLESModule

  662.     USE MudSystemVARIABLES

  663. use SimulationVariables !@@@

  664.     

  665.     IMPLICIT NONE

  666.     

  667.     INTEGER :: ilocal

  668. 

  669.     DEALLOCATE(KickJacobian , OldKickJacobian , KickVandPFunction , KickUnknownVector , KickCorrectionVector)

  670.     

  671.         

  672.     WRITE (*,*) 'Gas Pocket Will be Removed', ilocal

  673.     call GasPocketOldPress%Remove(ilocal)

  674.     call GasPocketNewPress%Remove(ilocal)

  675.     call GasPocketOldTemp%Remove(ilocal)

  676.     call GasPocketNewTemp%Remove(ilocal)

  677.     call GasPocketOldVol%Remove(ilocal)

  678.     call GasPocketNewVol%Remove(ilocal)

  679.     call GasPocketDeltaVol%Remove(ilocal)

  680.     call GasPocketFlowInduced%Remove(ilocal)

  681.     call GasPocketModifiedVol%Remove(ilocal)

  682.     call GasPocketWeight%Remove(ilocal)

  683.     call GasPocketDensity%Remove(ilocal)

  684.     call GasPocketCompressibility%Remove(ilocal)

  685.     

  686.     NoGasPocket = NoGasPocket - 1

  687.      

  688.     IF (NoGasPocket > 0) THEN

  689.         if (print_log) print *, "2 * NoGasPocket , 2 * NoGasPocket = ",2 * NoGasPocket , 2 * NoGasPocket

  690.         ALLOCATE(KickJacobian(2 * NoGasPocket , 2 * NoGasPocket) , OldKickJacobian(2 * NoGasPocket , 2 * NoGasPocket) , KickVandPFunction(2 * NoGasPocket))

  691.         ALLOCATE(KickUnknownVector(2 * NoGasPocket) , KickCorrectionVector(2 * NoGasPocket))

  692.         

  693.         IF (ALLOCATED(GasPocketFlowEl)) THEN

  694.             ALLOCATE(tempGasPocketFlowEl(NoGasPocket , SIZE(GasPocketFlowEl , dim = 2)))

  695.             tempGasPocketFlowEl = GasPocketFlowEl(1 : NoGasPocket , :)

  696.             DEALLOCATE(GasPocketFlowEl)

  697.             ALLOCATE(GasPocketFlowEl(SIZE(tempGasPocketFlowEl , dim = 1) , SIZE(tempGasPocketFlowEl , dim = 2)))

  698.             GasPocketFlowEl = tempGasPocketFlowEl

  699.             DEALLOCATE(tempGasPocketFlowEl)

  700.         END IF

  701.         

  702.     ELSE        ! NoGasPocket = 0

  703.         

  704.         

  705.         CALL GasPocketOldPress%Empty

  706.         CALL GasPocketNewPress%Empty

  707.         CALL GasPocketOldTemp%Empty

  708.         CALL GasPocketNewTemp%Empty

  709.         CALL GasPocketOldVol%Empty

  710.         CALL GasPocketNewVol%Empty

  711.         CALL GasPocketDeltaVol%Empty

  712.         CALL GasPocketFlowInduced%Empty

  713.         CALL GasPocketModifiedVol%Empty

  714.         CALL GasPocketWeight%Empty

  715.         CALL GasPocketDensity%Empty

  716.         

  717.         data%Equipments%DownHole%InfluxRate = 0.0

  718.         data%Equipments%DownHole%KickVolume = 0.0

  719.                 

  720.     END IF         

  721.     

  722. !WRITE (*,*) 'Gas Pocket Removed', NoGasPocket + 1, ALLOCATED(GasPocketWeight%Array)

  723.     

  724. END SUBROUTINE     

  725.                    

  726.                    

  727.                    

  728.                    

  729.                    

  730.                    

  731.