I use the below code for a kind of estimation. As you can see, it is so big and has long run time, I want to modify it to a short lines code and faster one. anyone can help me?
This code is written for estimating the static Elastic parameters of rock from dynamic one. there is many methods and formula for this. and i want to check most of them in one view. i think i can use a Data Frame which includes method name and formula but is i possible and fast?
methods=['Linear method of Eissa & Kazzi 1988 (E_s= 0.64E_D-0.32)',
'McCann & Entwise 1992 (E_s= 0.48E_d-3.26) it is for crystaline Rocks',
'McCann & Entwise 1992 (E_s= 0.69E_d-6.40) it is for all Rocks',
'the second linear formula of Eissa & Kazzi 1988 (E_s= 0.74E_d-0.82)',
'Christaras et al. 1994 (E_s= 1.05E_d-3.16)',
'Neville 1995 (E_s= 0.83E_d)',
'Brautigam et al. 1998 (E_s = 0.2807E_d)',
'Nur & Wang 1999 (E_s=1.153E_d-15.2)',
'Mokovciakova and Pandula 2003 (E_s = 0.7707E_d − 5854) ',
'Fahimifar and Soroush 2003 (Es = e(0.0477Ed))',
'Ameen et al. 2009 (Es = 0.541Ed + 12.852)',
'Brotons et al. 2014 (E_s= 0.86E_d-2.085) it is for crystaline Rocks',
'Brotons et al. 2016 (E_s= 0.932E_d-3.42)',
'Eissa and Kazzi 1988 nonlinear (log10(E_s)=0.77*log10(rho*E_d)+0.02)',
'Brotons et al 2016 nonlinear (log10(E_s)=0.976*log10(rho*E_d)-3.306) in this equation rho in kg/m3',
'Eissa & Kazzi 1988 in power law equation (E_s=(0.6019rho+0.42)*E_d^0.77)',
'Brotons et al. 2016 Power (E_s=(1.24rho+2.78)*E_d^0.96)',
'Brotons et al 2016 based on paper ((11.531*rho^-0.457)*E_d^1.251)',
'Oven 2003 (E_s=0.0158*E_d^2.74)',
'Horsrud 2001 (E_s=0.0428*E_d^2 + 0.2334*E_d)',
'Lacy 2007 for Sedimentary Rocks (E_s=0.018*E_d^2 + 0.422*E_d)',
'Davarpanah 2020 (E_s=0.89*E_d-5.26)',
'Davarpanah 2020 (E_s=0.164*E_d^1.373)',
'Davarpanah 2020 (log10(E_s)=1.37*log10(E_d*rho)-1.32)',
'NISOC Dynamic to Static Relationship (E_s=0.7E_d)',
'Fjear et al. 2008 for Sst, Shale and Lst',
'Afshari et al. 2010 (E_s= 0.4154*E_d - 1.0593)',
'Ezati PhD Thesis (E_s=0.0173E_d^2 + 0.9715E_d + 17.274)',
'Najibi et al. 2015 (Es = 0.014(Ed ^ 1.96))',
'Najibi et al. 2015 (Es = 0.169(Vp^3.324))',
'Shijie Shen 2024 (E_s = 1.0987E_d -13.519)',
'Lithology Base Najibi for Limestone, Vanheerden for Sst, Lacy for Shale and Davarpansh for else.']
Static_Elastic={}
for sheet_name in sheet_names:
dyn = dynamic_Elastic_Lithology[sheet_name].copy()
Well={}
for method in methods:
Well[method] = dyn.copy()
if method == 'Linear method of Eissa & Kazzi 1988 (E_s= 0.64E_D-0.32)':
Well[method]['St_Young_Modulus'] = ((0.64 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 0.32
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'McCann & Entwise 1992 (E_s= 0.48E_d-3.26) it is for crystaline Rocks':
Well[method]['St_Young_Modulus'] = ((0.48 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 3.26
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'McCann & Entwise 1992 (E_s= 0.69E_d-6.40) it is for all Rocks':
Well[method]['St_Young_Modulus'] = ((0.69 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 6.40
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'the second linear formula of Eissa & Kazzi 1988 (E_s= 0.74E_d-0.82)':
Well[method]['St_Young_Modulus'] = ((0.74 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 0.82
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Christaras et al. 1994 (E_s= 1.05E_d-3.16)':
Well[method]['St_Young_Modulus'] = ((1.05 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 3.16
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Neville 1995 (E_s= 0.83E_d)':
Well[method]['St_Young_Modulus'] = 0.83 * (dyn['Dy_Young_Modulus'] / (10 ** 9))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Brautigam et al. 1998 (E_s = 0.2807E_d)':
Well[method]['St_Young_Modulus'] = 0.2807 * (dyn['Dy_Young_Modulus'] / (10 ** 9))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Nur & Wang 1999 (E_s=1.153E_d-15.2)':
Well[method]['St_Young_Modulus'] = ((1.153 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 15.2
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Mokovciakova and Pandula 2003 (E_s = 0.7707E_d − 5854) ':
Well[method]['St_Young_Modulus'] = (0.7707 * (dyn['Dy_Young_Modulus']) / (10 ** 9)) - 5854
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Fahimifar and Soroush 2003 (Es = e(0.0477Ed))':
Well[method]['St_Young_Modulus'] = np.exp(0.0477 * (dyn['Dy_Young_Modulus']) / (10 ** 9))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Ameen et al. 2009 (Es = 0.541Ed + 12.852)':
Well[method]['St_Young_Modulus'] = (0.541 * (dyn['Dy_Young_Modulus']) / (10 ** 9)) + 12.852
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Brotons et al. 2014 (E_s= 0.86E_d-2.085) it is for crystaline Rocks':
Well[method]['St_Young_Modulus'] = ((0.86 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 2.085
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Brotons et al. 2016 (E_s= 0.932E_d-3.42)':
Well[method]['St_Young_Modulus'] = ((0.932 * dyn['Dy_Young_Modulus']) / (10 ** 9)) - 3.42
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Eissa and Kazzi 1988 nonlinear (log10(E_s)=0.77*log10(rho*E_d)+0.02)':
Well[method]['St_Young_Modulus'] = 10 ** ((np.log10((dyn['Dy_Young_Modulus'] * dyn['Density (Kg/m3)'])/(10 ** 12)) * (0.77)) + (0.02))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Brotons et al 2016 nonlinear (log10(E_s)=0.976*log10(rho*E_d)-3.306) in this equation rho in kg/m3':
Well[method]['St_Young_Modulus'] = 10 ** ((0.967 * (np.log10((dyn['Dy_Young_Modulus'] * dyn['Density (Kg/m3)'])/(10 ** 12)))) - 3.306)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Eissa & Kazzi 1988 in power law equation (E_s=(0.6019rho+0.42)*E_d^0.77)':
Well[method]['St_Young_Modulus'] = ((0.6019 * (dyn['Density (Kg/m3)']/1000)) + 0.42) * (dyn['Dy_Young_Modulus'] ** 0.77)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Brotons et al. 2016 Power (E_s=(1.24rho+2.78)*E_d^0.96)':
Well[method]['St_Young_Modulus'] = ((1.24 * (dyn['Density (Kg/m3)']/1000)) + 2.78) * (dyn['Dy_Young_Modulus'] ** 0.96)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Brotons et al 2016 based on paper ((11.531*rho^-0.457)*E_d^1.251)':
Well[method]['St_Young_Modulus'] = (11.531 * (dyn['Density (Kg/m3)'] ** (-0.457))) * (dyn['Dy_Young_Modulus'] ** 1.251)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Oven 2003 (E_s=0.0158*E_d^2.74)':
Well[method]['St_Young_Modulus'] = 0.0158 * (dyn['Dy_Young_Modulus'] ** 2.74)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Horsrud 2001 (E_s=0.0428*E_d^2 + 0.2334*E_d)':
Well[method]['St_Young_Modulus'] = ((0.0428 * dyn['Dy_Young_Modulus']) / (10 ** 9)) + (0.2334 * dyn['Dy_Young_Modulus'] / (10 ** 9))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Lacy 2007 for Sedimentary Rocks (E_s=0.018*E_d^2 + 0.422*E_d)':
Well[method]['St_Young_Modulus'] = ((0.018 * (((dyn['Dy_Young_Modulus'] / (10 ** 9)) * 0.145) ** 2)) + (0.422 * (dyn['Dy_Young_Modulus'] / (10 ** 9)) * 0.145)) * 6.8945
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Davarpanah 2020 (E_s=0.89*E_d-5.26)':
Well[method]['St_Young_Modulus'] = (0.89 * (dyn['Dy_Young_Modulus'] / (10 ** 9))) - 5.26
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Davarpanah 2020 (E_s=0.164*E_d^1.373)':
Well[method]['St_Young_Modulus'] = 0.164 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 1.373)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Davarpanah 2020 (log10(E_s)=1.37*log10(E_d*rho)-1.32)':
Well[method]['St_Young_Modulus'] = 10 ** ((1.37 * (np.log10((dyn['Dy_Young_Modulus'] * dyn['Density (Kg/m3)'])/(10 ** 12)))) - 1.32)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'NISOC Dynamic to Static Relationship (E_s=0.7E_d)':
Well[method]['St_Young_Modulus'] = 0.7 * dyn['Dy_Young_Modulus']
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Fjear et al. 2008 for Sst, Shale and Lst':
for i in range(len(dyn.index)):
if dyn.loc[i , 'Lithology'] == 'Limestone' or dyn.loc[i , 'Lithology'] == 'limestone' or dyn.loc[i , 'Lithology'] == 'Lst':
Well[method]['St_Young_Modulus'] = (0.18 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 2)) + (0.422 * (dyn['Dy_Young_Modulus'] / (10 ** 9)))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif dyn.loc[i , 'Lithology'] == 'Sandstone' or dyn.loc[i , 'Lithology'] == 'sandstone' or dyn.loc[i , 'Lithology'] == 'Sst':
Well[method]['St_Young_Modulus'] = (0.0293 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 2)) + (0.4533 * (dyn['Dy_Young_Modulus'] / (10 ** 9)))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif dyn.loc[i , 'Lithology'] == 'Shale' or dyn.loc[i , 'Lithology'] == 'shale':
Well[method]['St_Young_Modulus'] = (0.0428 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 2)) + (0.2334 * (dyn['Dy_Young_Modulus'] / (10 ** 9)))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Afshari et al. 2010 (E_s= 0.4154*E_d - 1.0593)':
Well[method]['St_Young_Modulus'] = (0.4154 * dyn['Dy_Young_Modulus']) - 1.0593
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Ezati PhD Thesis (E_s=0.0173E_d^2 + 0.9715E_d + 17.274)':
Well[method]['St_Young_Modulus'] = (0.0173 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 2)) - (0.9715 * (dyn['Dy_Young_Modulus'] / (10 ** 9))) + 17.274
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Najibi et al. 2015 (Es = 0.014(Ed ^ 1.96))':
Well[method]['St_Young_Modulus'] = 0.014 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 1.96)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Najibi et al. 2015 (Es = 0.169(Vp^3.324))':
Well[method]['St_Young_Modulus'] = 0.169 * ((dyn['Vp (m/s)'] / 1000) ** 3.324)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Shijie Shen 2024 (E_s = 1.0987E_d -13.519)':
Well[method]['St_Young_Modulus'] = (1.0987 * (dyn['Dy_Young_Modulus'] / (10 ** 9))) - 13.519
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif method == 'Lithology Base Najibi for Limestone, Vanheerden for Sst, Lacy for Shale and Davarpansh for else.':
for i in range(len(dyn.index)):
if dyn.loc[i , 'Lithology'] == 'Limestone' or dyn.loc[i , 'Lithology'] == 'limestone' or dyn.loc[i , 'Lithology'] == 'Lst':
Well[method]['St_Young_Modulus'] = 0.014 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 1.96)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif dyn.loc[i , 'Lithology'] == 'Sandstone' or dyn.loc[i , 'Lithology'] == 'sandstone' or dyn.loc[i , 'Lithology'] == 'Sst':
Well[method]['St_Young_Modulus'] = 0.97 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 1.38)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
elif dyn.loc[i , 'Lithology'] == 'Shale' or dyn.loc[i , 'Lithology'] == 'shale':
Well[method]['St_Young_Modulus'] = (0.0428 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 2)) + (0.2334 * (dyn['Dy_Young_Modulus'] / (10 ** 9)))
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
else:
Well[method]['St_Young_Modulus'] = 0.164 * ((dyn['Dy_Young_Modulus'] / (10 ** 9)) ** 1.373)
Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
#Well[method]['St_Shear_Modulus'] = Well[method]['St_Young_Modulus'] / (2 * (1 + Well[method]['Poission Ratio']))
#Well[method]['St_Bulk_Modulus'] = Well[method]['St_Young_Modulus'] / (3 * (1 - (2 * Well[method]['Poission Ratio'])))
Static_Elastic[sheet_name] = Well