Matlab求钻井作业开始之前所需准备的钻井液的总量

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