钻井优选参数的设计.doc

钻井作业（三）钻井优选参数的设计1. 变量说明c1牙齿磨损减慢系数D1钻头直径D2钻头直径Q1排量Q2排量c2牙齿磨损系数Cr钻机作业费Tr起下钻接单根时间Cb钻头成本te钻头与起下钻成本的折算时间W钻压Af地层研磨系数N转速b轴承工作系数CH水力净化系数Cp压差影响系数k比例系数M门限钻压j常量B1关于牙齿磨损量的一元二次方程的一次项系数a关于牙齿磨损量的一元二次方程的二次项系数B2代数值B3代数值c关于牙齿磨损量的一元二次方程的常数项系数X牙齿磨损量2. 编程Private Sub Command1_Click()c1 = 5D1 = 0.0146D2 = 6.44Q1 = 1.5Q2 = 0.00006525c2 = 3.68Cr = 350Tr = 6Cb = 850te = Cb / Cr + TrW = 304Af = 0.00123N = 100b = 7689784CH = 1Cp = 1k = 0.00755M = 9.8j = 0.68B1 = 1a = c1 / 2t = Text5.TextB2 = (D2 - D1 * W)B3 = Af * (Q1 * N + Q2 * (N 3) * tc = -(B3 / B2)If B1 2 - 4 * a * c = 0 ThenX = (-B1 + Sqr(B1 2 - 4 * a * c) / (2 * a)Bf = (c1 * X 2 / 2 + X) / (Af * (Q1 + Q2 * N 2) * b) * (D2 - D1 * W) * W 1.5)Text6.Text = BfText7.Text = XCost = (Cr * (te * Af * (Q1 * N + Q2 * N 3) / (D2 - D1 * W) + X + c1 * X 2) / 2) / (CH * Cp * k * (W - M) * N j * (c1 * X / c2 + (c2 - c1) * Log(1 + c2 * X) / (c2 2)Text8.Text = CostH = (Cb + Cr * (Tr + t) / CostText9.Text = HEnd IfEnd SubPrivate Sub Command2_Click()c2 = 3.68c1 = 5N = 100D1 = 0.0146D2 = 6.44Q1 = 1.5Q2 = 0.00006525c2 = 3.68Cr = 350Tr = 6Cb = 850te = Cb / Cr + TrW = 304Af = 0.00123N = 100b = 7689784CH = 1Cp = 1k = 0.00755M = 9.8j = 0.68D1 = 0.0146D2 = 6.44Q1 = 1.5Q2 = 0.00006525Cb = 850Cr = 350N = 100W = 304Tr = 6Af = 0.00123te = Cb / Cr + TrDoX1 = Text10.TextX2 = Text11.TextLoop While Sgn(c1 * X1 2 / 2 + (c1 / c2 - 1) * X1) - (c1 - c2) * (1 + c2 * X1) * Log(1 + c2 * X1) / (c2 2) - (Af * te * (Q1 * N + Q2 * N 3) / (D2 - D1 * W) = Sgn(c1 * X2 2 / 2 + (c1 / c2 - 1) * X2) - (c1 - c2) * (1 + c2 * X2) * Log(1 + c2 * X2) / (c2 2) - (Af * te * (Q1 * N + Q2 * N 3) / (D2 - D1 * W)i = 0Dof1 = c1 * X1 2 / 2 + (c1 / c2 - 1) * X1) - (c1 - c2) * (1 + c2 * X1) * Log(1 + c2 * X1) / (c2 2) - (Af * te * (Q1 * N + Q2 * N 3) / (D2 - D1 * W)f2 = c1 * X2 2 / 2 + (c1 / c2 - 1) * X2) - (c1 - c2) * (1 + c2 * X2) * Log(1 + c2 * X2) / (c2 2) - (Af * te * (Q1 * N + Q2 * N 3) / (D2 - D1 * W)X0 = (X1 + X2) / 2f0 = c1 * X0 2 / 2 + (c1 / c2 - 1) * X0) - (c1 - c2) * (1 + c2 * X0) * Log(1 + c2 * X0) / (c2 2) - (Af * te * (Q1 * N + Q2 * N 3) / (D2 - D1 * W)If Sgn(f0) Sgn(f1) ThenX2 = X0f2 = f0ElseX1 = X0f1 = f0End Ifi = i + 1Loop While Abs(f0) 0.00001 And i 100hf = X0Text1.Text = hft = (hf + (c1 * hf 2) / 2) * (D2 - D1 * W) / (Af * (Q1 * N + Q2 * N 3)Text2.Text = tCost = (Cr * (te * Af * (Q1 * N + Q2 * N 3) / (D2 - D1 * W) + hf + c1 * hf 2) / 2) / (CH * Cp * k * (W - M) * N j * (c1 * hf / c2 + (c2 - c1) * Log(1 + c2 * hf) / (c2 2)Text3.Text = CostH = (Cb + Cr * (Tr + t) / CostT