建立基于测井数据岩石力学解释的函数库

# -*- coding: utf-8 -*-
"""
Created on Thu Aug 29 05:15:01 2019
在大数据计算中，钻井测井数据一类必不可少的结构化数据。建立基于测井数据岩石力学解释的函数库，是一件很有意义的工作。这里只是做了一个开始，感兴趣的朋友，可以在此基础上进一步丰富和完善函数的数量就行了。
@author: nepu Liw
"""import math# =============================================================================
# 算法分析——钻井——测井解释 
# 建立基于测井数据的岩石力学解释的函数库
# =============================================================================
try: 
#================================================================# 从这里开始编写#================================================================# 第一个函数def hardness(AC=127):        '岩石硬度和声波时差---库车坳陷'result= 18260 * math.exp(-0.051 * AC)return round(result,8)
#================================================================def vpt2vst(Vpt=127):result =(11.44 * Vpt + 18.03) ** 0.5 - 5.686return round(result,6)
#================================================================def internal_friction_angle(C=10):#计算岩石的内摩擦角result = 36.545 - 0.4952 * Creturn round(result,6)
#================================================================        def Horizontal_geostress(V=0.2,Po=25,Pp=15):#假设两向水平应力相等，求水平应力result =V / (1 - V) * (Po - Pp)return round(result,6)
#================================================================
except:pass#================================================================#================================================================
# 函数测试 F5
#================================================================
print('Py=',hardness(100))
print('Vst=',vpt2vst())
print('Vst=',internal_friction_angle())
print('Horizontal_geostress=',Horizontal_geostress())计算结果：Py= 111.32659229
Vst= 32.666445
Vst= 31.593
Horizontal_geostress= 2.5

作者：LiW