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...

730 lines
37 KiB
Raw Permalink Blame History 

  1. SUBROUTINE WellPressureDataTransfer

  2.     

  3.     !! This subroutine calculates pressure distribution in well

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

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

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

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

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

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

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

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

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

  13.     !! frictional pressure drop calculated by subroutine Sub_Press_Drop_Calculator

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

  15.     !! which causes acceleration or decceleration of flow

  16.     

  17.     USE FricPressDropVarsModule

  18.     use MudSystemModule

  19.     USE MudSystemVARIABLES

  20.     use SimulationVariables !@@@

  21.     USE GeoElements_FluidModule

  22.     USE Fluid_Flow_Startup_Vars

  23.     use KickVARIABLESModule

  24.     ! USE CMudPropertiesVariables

  25.     ! USE CBopStackVariables

  26.     ! USE sROP_Variables

  27.     ! USE CHOKEVARIABLES

  28.     ! use CChokeManifoldVariables

  29.     ! use SimulationVariables

  30.     ! USE CBitProblemsVariables

  31.     ! USE CChokeProblemsVariables

  32.     ! USE CStringConfigurationVariables

  33.     ! use CDataDisplayConsole

  34.     

  35.     IMPLICIT NONE

  36.     

  37.     INTEGER :: i , j

  38.     REAL :: InstantaneousTotalOpenChokeArea , VolumeDensityProduct

  39.     

  40.     ChokeIsClosing = .FALSE.

  41.     data%State%FricPressDrop%Kchoke = 0.0

  42.     data%State%FricPressDrop%KBit = 0.0

  43.     

  44.     

  45.     !!!!!!!!!!!!!!!! bit is present

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

  47.     BitCd = 0.98

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

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

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

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

  52.     IF (BitNozzleNum == 0) THEN

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

  54.     ELSE

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

  56.     END IF

  57.     

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

  59.     !FloatValveOpen = .TRUE.

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

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

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

  63.     !END IF

  64.     

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

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

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

  68.     

  69.     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  70.     

  71.     IF(ALLOCATED(GasPocketFlowEl))      DEALLOCATE(GasPocketFlowEl)

  72.     

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

  74.     ChokeKroneckerDelta = 0

  75. 

  76.     

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

  78.         

  79.         WellHeadOpen = .FALSE.

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

  81.     ELSE 

  82.         WellHeadOpen = .TRUE.

  83.     END IF

  84.     

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

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

  87.     

  88. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  89.     

  90.     

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

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

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

  94.     data%State%FricPressDrop%NoWellToChokeEl = 0

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

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

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

  98.     IF(ALLOCATED(FlowEl)) then

  99.         DEALLOCATE(FlowEl)

  100.     endif

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

  102.     

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

  104.     !!!!!!!!!!!!!!!!!!!!!!! Horizontal line

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

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

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

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

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

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

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

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

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

  114.         FlowEl(i)%Id = 0.0

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

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

  117.         

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

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

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

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

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

  123.             FlowEl(i)%MaterialType = 0

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

  125.             FlowEl(i)%MaterialType = 4

  126.         ELSE

  127.             FlowEl(i)%MaterialType = 1

  128.         END IF

  129.         

  130.         FlowEl(i)%dPdLgrav = 0.0

  131.         FlowEl(i)%StaticPressDiff = 0.0

  132. 

  133. 

  134.         

  135.     END DO

  136.     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  137.     

  138.     !!!!!!!!!!!!!!!!!!!!!!!!! String

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

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

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

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

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

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

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

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

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

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

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

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

  151.         FlowEl(i)%Id = 0.0

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

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

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

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

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

  157.         

  158.         

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

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

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

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

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

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

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

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

  167.             FlowEl(i)%MaterialType = 0

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

  169.             FlowEl(i)%MaterialType = 4

  170.         ELSE

  171.             FlowEl(i)%MaterialType = 1

  172.         END IF

  173.         

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

  175.             FlowEl(i)%dPdLgrav = 0.0

  176.         ELSE

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

  178.         END IF

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

  180.         

  181. 

  182. 

  183.         

  184.     END DO

  185.     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  186.     

  187.     !!!!!!!!!!!!!!!!!!!!!!!!!!!! Annulus

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

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

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

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

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

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

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

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

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

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

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

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

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

  201.         

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

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

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

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

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

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

  208.             FlowEl(i)%MaterialType = 0

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

  210.             FlowEl(i)%MaterialType = 4

  211.         ELSE

  212.             FlowEl(i)%MaterialType = 1

  213.         END IF

  214.         

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

  216.             FlowEl(i)%dPdLgrav = 0.0

  217.         ELSE

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

  219.         END IF

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

  221. 

  222. 

  223.         

  224.     END DO

  225.     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  226.     

  227.     !!!!!!!!!!!!!!! Well to choke manifold path

  228.     j = 1

  229.     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

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

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

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

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

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

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

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

  237.         FlowEl(i)%Id = 0.0

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

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

  240.             

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

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

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

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

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

  246.         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

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

  248.         ELSE

  249.             FlowEl(i)%MaterialType = 1

  250.         END IF

  251.         

  252.             

  253.         FlowEl(i)%dPdLgrav = 0.0

  254.         FlowEl(i)%StaticPressDiff = 0.0

  255.         

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

  257. 

  258.         j = j + 1

  259.             

  260.     END DO

  261.             

  262.     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  263.     

  264.     !!!!!!!!!!!!!!! Open Hole

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

  266.     j = 1

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

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

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

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

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

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

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

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

  275.         FlowEl(i)%Id = 0.0

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

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

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

  279.         

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

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

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

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

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

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

  286.             FlowEl(i)%MaterialType = 0

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

  288.             FlowEl(i)%MaterialType = 4

  289.         ELSE

  290.             FlowEl(i)%MaterialType = 1

  291.         END IF

  292.         

  293.         j = j + 1

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

  295.             FlowEl(i)%dPdLgrav = 0.0

  296.         ELSE

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

  298.         END IF

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

  300.     

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

  302. 

  303.     END DO

  304. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  305.     

  306. !!!!!!!!!!!!!!!!!!  Mud Element Viscosity Calculation

  307.     

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

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

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

  311.             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

  312.             

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

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

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

  316.             END DO

  317.         

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

  319.             CYCLE

  320.         END IF

  321.         

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

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

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

  325.     END DO

  326. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  327. 

  328.     

  329.     FlowEl(:)%FlowRate = 0.0

  330.     FlowEl(:)%StartPress = 0.0

  331.     FlowEl(:)%EndPress = 0.0

  332.     

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

  334.     

  335.     

  336.     

  337. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  338.     

  339. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  340.     

  341.     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

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

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

  344.     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

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

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

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

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

  349.     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

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

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

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

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

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

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

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

  357.         

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

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

  360.         data%State%FricPressDrop%TotalOpenChokeArea = InstantaneousTotalOpenChokeArea

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

  362.     ELSE

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

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

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

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

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

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

  369.         ELSE

  370.             data%State%FricPressDrop%TotalOpenChokeArea = InstantaneousTotalOpenChokeArea

  371.         END IF

  372.     END IF

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

  374.     

  375.         

  376.     

  377.     

  378.             

  379. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!            

  380.     

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

  382.     ChokeKroneckerDelta = 0

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

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

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

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

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

  388.         

  389.         WellHeadOpen = .FALSE.

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

  391.     ELSE

  392.         WellHeadOpen = .TRUE.

  393.     END IF

  394.     

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

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

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

  398. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  399. 

  400.     IF (NoGasPocket > 0) THEN

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

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

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

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

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

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

  407.                 CALL RemoveGasPocket(NoGasPocket)

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

  409.             END IF

  410.             

  411.         END IF

  412.         

  413.         

  414.         !NoGasPocket = SIZE(GasPocketWeight%Array)

  415.         IF (NoGasPocket > 0)         CALL GasPocketFlowElementTransformer

  416.     END IF

  417.     

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

  419.         OldChokeDensity = ChokeDensity

  420.         VolumeDensityProduct = 0.0

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

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

  423.         End Do

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

  425.                             

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

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

  428.     END IF

  429.     

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

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

  432.     END IF

  433.     

  434.     

  435.     

  436. END SUBROUTINE WellPressureDataTransfer

  437.     

  438.     

  439. SUBROUTINE GasPocketFlowElementTransformer

  440. 

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

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

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

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

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

  446. 

  447. 

  448.     USE FricPressDropVarsModule

  449.     USE MudSystemVARIABLES

  450.     use SimulationVariables

  451.     USE GeoElements_FluidModule

  452.     USE Fluid_Flow_Startup_Vars

  453.     use KickVARIABLESModule

  454.     USE CMudPropertiesVariables

  455.     USE CBopStackVariables

  456.     USE CDownHoleVariables

  457.     USE CError

  458.     

  459.     IMPLICIT NONE

  460.     

  461.     INTEGER :: i , j , k

  462.     REAL :: PressureCorrection

  463.     PressureCorrection = 1.0

  464.     

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

  466.     ALLOCATE(GasPocketFlowEl(NoGasPocket , 1))

  467.     GasPocketFlowEl(:,:) = 0

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

  469.     

  470.     i = 1

  471.     j = 1

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

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

  474.             IF (j > 1) THEN

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

  476.                     i = i + 1

  477.                     j = 1

  478.                 END IF

  479.             END IF

  480.             

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

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

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

  484.                 tempGasPocketFlowEl = GasPocketFlowEl

  485.                 DEALLOCATE(GasPocketFlowEl)

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

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

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

  489.                 DEALLOCATE(tempGasPocketFlowEl)

  490.             END IF

  491.             

  492.             GasPocketFlowEl(i , j) = k

  493.             j = j + 1

  494.         END IF

  495.     END DO

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

  497.     

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

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

  500.             IF (j > 1) THEN

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

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

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

  504.                     i = i + 1

  505.                     j = 1

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

  507.                 END IF

  508.             END IF

  509.             

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

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

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

  513.                 tempGasPocketFlowEl = GasPocketFlowEl

  514.                 DEALLOCATE(GasPocketFlowEl)

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

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

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

  518.                 DEALLOCATE(tempGasPocketFlowEl)

  519.             END IF

  520.             IF (i > NoGasPocket) THEN

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

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

  523.                 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

  524.             

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

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

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

  528.                 END DO

  529. 

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

  531.             END IF

  532. 

  533.             GasPocketFlowEl(i , j) = k

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

  535.             j = j + 1

  536.             

  537.         END IF

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

  539.     END DO

  540.     

  541. 

  542. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!    

  543.     

  544. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!    

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

  546.     

  547.     DO i = 1 , NoGasPocket

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

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

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

  551.         END DO

  552.     END DO

  553.     i = data%State%FricPressDrop%AnnulusLastEl

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

  555.     

  556.     

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

  558.     

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

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

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

  562.     

  563.     

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

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

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

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

  568.         END DO

  569.     END DO

  570.     

  571.     

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

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

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

  575.     

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

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

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

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

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

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

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

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

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

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

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

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

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

  589.             CALL RemoveGasPocket(NoGasPocket)

  590.             

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

  592.             

  593.         END IF

  594.         

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

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

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

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

  599.     

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

  601.         

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

  603.         

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

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

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

  607.         

  608.         

  609.         

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

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

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

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

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

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

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

  617.             END DO

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

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

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

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

  622.             j = NoGasPocket

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

  624.         END DO

  625.         

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

  627.         WRITE (*,*) 'After'

  628.         DO i = 1 , NoGasPocket

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

  630.         END DO

  631.         

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

  633.         

  634.     END IF

  635.     

  636.     

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

  638.         CALL RemoveGasPocket(1)

  639.         

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

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

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

  643.     END IF

  644.     

  645.     

  646.     

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

  648.     

  649.     

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

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

  652.         

  653.         

  654.     

  655.     END SUBROUTINE

  656.     

  657.     

  658. SUBROUTINE RemoveGasPocket(ilocal)

  659.     

  660.     use KickVARIABLESModule

  661.     USE MudSystemVARIABLES

  662. use SimulationVariables !@@@

  663.     

  664.     IMPLICIT NONE

  665.     

  666.     INTEGER :: ilocal

  667. 

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

  669.     

  670.         

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

  672.     call GasPocketOldPress%Remove(ilocal)

  673.     call GasPocketNewPress%Remove(ilocal)

  674.     call GasPocketOldTemp%Remove(ilocal)

  675.     call GasPocketNewTemp%Remove(ilocal)

  676.     call GasPocketOldVol%Remove(ilocal)

  677.     call GasPocketNewVol%Remove(ilocal)

  678.     call GasPocketDeltaVol%Remove(ilocal)

  679.     call GasPocketFlowInduced%Remove(ilocal)

  680.     call GasPocketModifiedVol%Remove(ilocal)

  681.     call GasPocketWeight%Remove(ilocal)

  682.     call GasPocketDensity%Remove(ilocal)

  683.     call GasPocketCompressibility%Remove(ilocal)

  684.     

  685.     NoGasPocket = NoGasPocket - 1

  686.      

  687.     IF (NoGasPocket > 0) THEN

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

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

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

  691.         

  692.         IF (ALLOCATED(GasPocketFlowEl)) THEN

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

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

  695.             DEALLOCATE(GasPocketFlowEl)

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

  697.             GasPocketFlowEl = tempGasPocketFlowEl

  698.             DEALLOCATE(tempGasPocketFlowEl)

  699.         END IF

  700.         

  701.     ELSE        ! NoGasPocket = 0

  702.         

  703.         

  704.         CALL GasPocketOldPress%Empty

  705.         CALL GasPocketNewPress%Empty

  706.         CALL GasPocketOldTemp%Empty

  707.         CALL GasPocketNewTemp%Empty

  708.         CALL GasPocketOldVol%Empty

  709.         CALL GasPocketNewVol%Empty

  710.         CALL GasPocketDeltaVol%Empty

  711.         CALL GasPocketFlowInduced%Empty

  712.         CALL GasPocketModifiedVol%Empty

  713.         CALL GasPocketWeight%Empty

  714.         CALL GasPocketDensity%Empty

  715.         

  716.         data%Equipments%DownHole%InfluxRate = 0.0

  717.         data%Equipments%DownHole%KickVolume = 0.0

  718.                 

  719.     END IF         

  720.     

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

  722.     

  723. END SUBROUTINE     

  724.                    

  725.                    

  726.                    

  727.                    

  728.                    

  729.                    

  730.