function [Vz] = volume(H,dh,dn1,dn2,dn3,l1,l2,l3,Vc)
%钻井作业开始之前,所需准备的钻井液的总量,即Vz;
% 根据Vz=Vj+Vg+Vc+Vb求解,其中Vj为井筒钻井液量;
% Vg为地面高压循环管汇的钻井液量;Vc为循环罐内的钻井液量;
% Vb为备用钻井液量,它们的单位均为立方米;
% 求解: 输入需要的参数井深H,以及井眼直径dh,dh取米为单位输入;
% 输入各地面高压循环管汇的内径dn单位取毫米;确定各单位容积V(1,2,3),
% 输入需要的参数井深H,以及各管汇的长度l(1,2,3);
% 例:输入: H=3450;
% dh=0.2159;
% dn1=123;
% dn2=102;
% dn3=80;
% l1=40;
% l2=18;
% l3=14.5;
% Vc=160;
% 输出:Vz=539.56(立方米)。
%% 确定地面高压循环管汇的单位容积
A=xlsread('volume.xlsx'); %地面高压循环管汇的单位容积表
H=input('Enter the H=');
dh=input('Enter the dh=');
dn1=input('Enter the dn1=');
dn2=input('Enter the dn2=');
dn3=input('Enter the dn3=');
[m1,n1]=find(A==dn1);
V1=A(m1,n1+1);
[m2,n2]=find(A==dn2);
V2=A(m2,n2+2);
[m3,n3]=find(A==dn2);
V3=A(m3,n3+2);
%% 输入各参数并进行相关参数计算
l1=input('Enter the l1=');
l2=input('Enter the l2=');
l3=input('Enter the l3=');
Vc=input('Enter the Vc=');
Vj=1/4*pi*dh^2*H;
Vg=V1/1000*l1+V2/1000*l2+V3/1000*l3;
Vb=2*Vj;
Vz=Vj+Vg+Vc+Vb;
end
作者:TianSL