CN105510966B - A kind of pore media SVEL computational methods and device containing non-newtonian fluid - Google Patents

A kind of pore media SVEL computational methods and device containing non-newtonian fluid Download PDF

Info

Publication number
CN105510966B
CN105510966B CN201511019018.7A CN201511019018A CN105510966B CN 105510966 B CN105510966 B CN 105510966B CN 201511019018 A CN201511019018 A CN 201511019018A CN 105510966 B CN105510966 B CN 105510966B
Authority
CN
China
Prior art keywords
mrow
fluid
mfrac
newtonian
pore media
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201511019018.7A
Other languages
Chinese (zh)
Other versions
CN105510966A (en
Inventor
曹宏
孙卫涛
卢明辉
杨志芳
晏信飞
李晓明
未晛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201511019018.7A priority Critical patent/CN105510966B/en
Publication of CN105510966A publication Critical patent/CN105510966A/en
Application granted granted Critical
Publication of CN105510966B publication Critical patent/CN105510966B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The embodiment of the present application, which discloses a kind of the pore media SVEL computational methods containing non-newtonian fluid and device, methods described, to be included：Establish pore media microchannel model, according to the mean flow rate of shear stress and fluid between the dissipation energy inside flow model export non-newtonian fluid, Fractional Derivative Maxwell fluid solid, according to the shear stress and the mean flow rate, stream solid phase is calculated to kinergety dissipation factor；Based on the solid relative motion Dissipated energy of the stream, the pore media wave equation for including non-newtonian flow bulk effect is established；Plane wave analysis is carried out to the pore media wave equation of the foundation, obtains including the dispersion relation of non-newtonian flow bulk effect, calculates velocity of longitudinal wave frequency dispersion and decay.The pore media SVEL computational methods and device containing non-newtonian fluid that the embodiment of the present application provides, can relatively accurately predict the frequency dispersion of SVEL and decay in unconventionaloil pool reservoir.

Description

A kind of pore media SVEL computational methods and device containing non-newtonian fluid
Technical field
The application is related to earthquake rock physicses technical field, more particularly to a kind of pore media sound wave containing non-newtonian fluid Speed calculation method and device.
Background technology
Fine and close oily geologic characteristics are complicated, and hydrocarbon reservoir rock includes the heterogeneities such as mineral grain and solid organic matter Phase, anisotropism is strong, and capillary pressure is big.Connected system, pore structure are shouted in major developmental micro/nano level hole in this rocks Complicated, poor connectivity, pore diameter range is low (being generally less than 10%) in nanometer, micron dimension, and porosity, and permeability is generally less than 0.1mD.These factors cause the complex fluid permeability behavior in reservoir pore space often not follow the constitutive relation of Newtonian fluid, Brought challenges to conventional meanses such as seismic explorations.
Fine and close oil oil is light, based on light oil, good fluidity.In conventional gas and oil exploration research work in general window gap Crude oil be Newtonian fluid, however, the characteristic of fluid changes in the pipeline of micro-nano pore size.Due to the flowing space Limitation, the fluid physical mechanics behavior in compact reservoir pore media can be more much more complex than perfect fluid, Newtonian fluid etc., table Reveal the characteristic of non-newtonian fluid.For example, the viscosity of fluid is the amount of a change in micro-nano pipeline, and conventional Newtonian liquid Viscosity is constant in body constitutive equation.When caliber reduces, arranging close to the regular molecular of solid occurs in fluid in micro-nano hole, It is the complex state material between ideal solid and perfect fluid to show now fluid, and its viscosity has changed, therefore Conventional Newtonian fluid has not been suitable for compact reservoir situation, more appropriate using non newtonian constitutive model.
Meet linear relationship between the shear stress and shear strain rate of Newtonian fluid, linear coefficient is fluid viscosity system Number.In non-newtonian fluid, linear relationship is not present between shear stress and shear strain rate, therefore can not also define a perseverance Fixed constant viscosity.According to the variation relation between shear stress and shear strain rate, non-newtonian fluid can be divided into Polytype, wherein mainly including：STF (dilatant fluid), fluid viscosity coefficient with the increase of shearing rate and Increase；Shear-thinning fluid (pseudoplastic fluid), fluid viscosity coefficient reduce with the increase of shearing rate；Bingham's plasticity stream Body, fluid need to reach certain yield stress before flowing, after fluid starts flowing, fluid viscosity coefficient and shearing Rate meets linear relationship.
Propagation of the sound wave in complicated pore media is used to distribution and the property of Underground fluid.When saturation in hole During fluid, the speed and amplitude that can cause sound wave change, and produce frequency dispersion and attenuation effect.Biot is given full of Newtonian liquid The quantitative analysis of the pore media Elastic Wave field of body, and predict the presence of Slow P-wave.One kind etc. is given in Biot models The rock matrix of effect saturated with fluid velocity of wave Forecasting Methodology under global mobility status, this is the idealized model of actual conditions, and The influence that the micro-nano hole spatial liquid flow of compact reservoir is brought is not considered.
Tsiklauri and Beresnev introduces Maxwell's viscoelastic fluid in Biot classical theories, provides rotation Ripple (rotational wave) and dilatational wave (dilatational wave) are in the pore media full of Maxwell fluid The velocity dispersion of propagation and decay.The stress-strain relation of perfect elastic body meets Hooke's law, i.e., The stress of preferable viscoid and the relation of shearing rate meet Newton's law, i.e.,Fine and close oily reservoir fluid seepage flow Passage it is narrow (micro-nano magnitude), fluid has shown the state of the transition between fluid-solid under micro-nano-scale, Therefore in addition to thinking pore-fluid and can bear shearing, it is no longer Newtonian fluid to be additionally considered that fluid, and between perfect elastic body and reason Think the state between viscoid.But due to the complexity of oil gas molecular structure in fine and close pore structure itself, it is caused to be difficult to Described with the classical non-Newtonian model of textbook coideal, i.e. constitutive relationshipMiddle strain rate is closed In the derivative exponent number α of time may not be a positive integer.
Inventor has found that at least there are the following problems in the prior art：Using viscous in Biot pore media wave equations Property Newtonian fluid model, this hypothesis can not meet the application of non-homogeneous complicated pore media.In addition, introduce Maxwell fluid Biot theory hypothesis holes in be full of a kind of Utopian viscoelastic fluid, and the constitutive relation of real fluid be it is unknown, This idealization is unsatisfactory for it is assumed that can not the accurate complete wave field communication process described in compact reservoir rock.
The content of the invention
The purpose of the embodiment of the present application be to provide a kind of pore media SVEL computational methods containing non-newtonian fluid and Device, relatively accurately to predict the frequency dispersion of SVEL and decay in unconventionaloil pool reservoir.
In order to solve the above technical problems, the embodiment of the present application provides a kind of pore media SVEL containing non-newtonian fluid What computational methods and device were realized in：
A kind of pore media SVEL computational methods containing non-newtonian fluid, including：
Pore media microchannel model is established, the Dissipated energy inside non-newtonian fluid is exported according to the flow model Amount, the mean flow rate of shear stress and fluid between Fractional Derivative Maxwell fluid-solid, according to the shear stress and The mean flow rate, stream solid phase is calculated to kinergety dissipation factor；
Based on energy dissipation caused by the solid relative motion of stream, the pore media fluctuation side for including non-newtonian flow bulk effect is established Journey；
Plane wave analysis is carried out to the pore media wave equation of the foundation, obtains including the frequency of non-newtonian flow bulk effect The relation of dissipating, determines velocity of longitudinal wave frequency dispersion and decay.
It is described according to the shear stress and the mean flow rate in preferred scheme, calculate the energy dissipation coefficient tool Body includes：
According to fluid continuity equation, momentum conservation equation and fractional order Maxwell's non-newtonian fluid constitutive relation, obtain Space velocity distribution and wall shear stress to pipeline flowing；
It is distributed according to the space velocity and calculates non-newtonian fluid pipeline flowing mean flow rate；
The ratio relation of mean flow rate and wall shear stress is flowed by non-newtonian fluid pipeline again, obtains energy dissipation Coefficient.
In preferred scheme, energy dissipation caused by the solid relative motion of stream is calculated by following formula：
In formula, D represents energy dissipation caused by the solid relative motion of stream；B represents stream solid phase to kinergety dissipation factor； F (κ) is high frequency correction factor；U, U represent the displacement of solid and fluid respectively；The displacement of solid and fluid is represented respectively To the speed of the first derivative of time, i.e. solid and fluid.
In preferred scheme,
The fluid continuity equation is：
The momentum conservation equation is：
The constitutive relation of the fractional order Maxwell non-newtonian fluid is：
Wherein, vfFluid velocity is represented, p represents pressure, and τ represents fluid pipe walls shear stress, and η represents fluid viscosity, μ tables Show modulus of shearing；λ=η/μ, represent the fluid Relaxation Characteristics time；ρfThe density of pore-fluid is represented, α represents stress to the time The exponent number of Fractional Derivative, β represent exponent number of the strain to time score order derivative；
Fractional DerivativeIt is defined as：
Wherein, m=1,2,3 ..., and 0≤m-1≤α≤m.
In preferred scheme, the VELOCITY DISTRIBUTION that pipeline flowing is calculated, it is calculated according to following formula：
Wherein,
V represents the speed of pipeline flowing；P represents pressure, and η represents fluid viscosity；ρfFluid density is represented, r represents cylindricality Radius in coordinate；U represents the displacement of solid；λ=η/μ, represent the fluid Relaxation Characteristics time；ω represents wave field frequency, and α is represented Stress is to the exponent number of the Fractional Derivative of time, and β expressions strain is to the exponent number of time score order derivative, μ expression hole non newtonians The modulus of shearing of fluid.In preferred scheme, the calculating mean flow rate obtains according to following formula：
In formula,Represent desired mean flow rate；A is runner radius.
In preferred scheme, the fluid pipe walls shear stress, it is calculated using following formula：
Wherein, τ represents fluid pipe walls shear stress, and η represents fluid viscosity, and V represents the speed of pipeline flowing, and r represents post Radius in shape coordinate；ν=η/ρf, ρfRepresent fluid density；λ=η/μ, represent fluid The Relaxation Characteristics time；A is runner radius；ω represents wave field frequency, and α represents stress to the exponent number of the Fractional Derivative of time, β tables Show exponent number of the strain to time score order derivative, μ represents the modulus of shearing of hole non-newtonian fluid.
In preferred scheme, the fluid dissipation energy isWherein, τ represents shear stress, and γ represents that shearing should Become,Represent first derivatives of the shear strain γ to the time；Represent porosity.
It is described that plane wave analysis is carried out to pore media wave equation in preferred scheme, obtain imitating comprising non-newtonian fluid The dispersion relation answered, including：The wave equation is transformed into frequency wavenumber domain, obtains the dispersion equation of compressional wave.
A kind of pore media SVEL computing device containing non-newtonian fluid, including：Solid phase is flowed to dissipate to kinergety Coefficients calculation block, pore media wave equation establish module and plane wave analysis module；Wherein,
The stream solid phase is to kinergety dissipation factor computing module, for establishing pore media microchannel model, According between the dissipation energy inside flow model export non-newtonian fluid, Fractional Derivative Maxwell fluid-solid The mean flow rate of shear stress and fluid, according to the shear stress and the mean flow rate, stream solid phase is calculated to kinergety Dissipation factor；
The pore media wave equation establishes module, for based on energy dissipation caused by the solid relative motion of stream, establishing Pore media wave equation comprising non-newtonian flow bulk effect；
The plane wave analysis module, for carrying out plane wave analysis to the pore media wave equation of the foundation, obtain To the dispersion relation for including non-newtonian flow bulk effect, velocity of longitudinal wave frequency dispersion and decay are calculated.
The technical scheme provided from above the embodiment of the present application, the embodiment of the present application provide containing non-newtonian fluid Pore media SVEL computational methods and device, the pore media wave equation model for including non-newtonian fluid parameter is established, The influence of pore character size, non-newtonian fluid constitutive relation feature to SVEL can be reflected.Simultaneously, it is contemplated that non newtonian Fractional Derivative relation in fluid constitutive relation, by the change to fractional order, it is specific that model is no longer limited to certain Non-Newtonian models are idealized, this method is can apply to the SVEL comprising different type complex fluid effect and predicts Situation.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this Some embodiments described in application, for those of ordinary skill in the art, do not paying the premise of creative labor Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is the flow chart of the pore media SVEL computational methods one embodiment of the application containing non-newtonian fluid；
Fig. 2 is the schematic diagram of the embodiment of the present application mesopore medium microchannel model；
Fig. 3 is the fast p wave interval velocity dispersion curve that non-newtonian fluid is considered in the embodiment of the present application；
Fig. 4 is the fast P attenuation curves that non-newtonian fluid is considered in the embodiment of the present application；
Fig. 5 is the module map of the pore media SVEL computing device one embodiment of the application containing non-newtonian fluid.
Embodiment
The embodiment of the present application provides a kind of pore media SVEL computational methods and device containing non-newtonian fluid.
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area The every other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection Scope.
Fig. 1 is the flow chart of the pore media SVEL computational methods one embodiment of the application containing non-newtonian fluid. As shown in figure 1, methods described can include：
S101：Pore media microchannel model is established, is exported according to the flow model inside non-newtonian fluid The mean flow rate of shear stress and fluid between dissipation energy, Fractional Derivative Maxwell fluid-solid, according to the shearing Stress and the mean flow rate, stream solid phase is calculated to kinergety dissipation factor.
Fig. 2 is the schematic diagram of the embodiment of the present application mesopore medium microchannel model.As shown in Figure 2, it can be assumed that Flowing is non-newtonian fluid in the microchannel, and the non-newtonian fluid is in the microchannel internal flow, and stream solid phase is to fortune Dynamic to produce energy dissipation, energy dissipation caused by the solid relative motion of stream can be calculated by following formula (1)：
In formula (1), D represents energy dissipation caused by the solid relative motion of stream；B represents stream solid phase to kinergety dissipative system Number；F (κ) is high frequency correction factor；U, U represent the displacement of solid and fluid respectively, and τ represents shear stress, and γ represents that shearing should Become,Represent porosity.The speed of the first derivative, i.e. solid and fluid of the displacement versus time of solid and fluid is represented respectively Degree.
It is described according to the shear stress and the mean flow rate, calculating the energy dissipation coefficient can specifically include： According to fluid continuity equation, momentum conservation equation and fractional order Maxwell's non-newtonian fluid constitutive relation, pipeline stream is obtained Dynamic space velocity distribution and wall shear stress；It is distributed according to the space velocity and calculates the flowing of non-newtonian fluid pipeline averagely Flow velocity；The ratio relation of mean flow rate and wall shear stress is flowed by non-newtonian fluid pipeline again, obtains energy dissipation system Number.
Specifically, calculating process can be as follows：
The constitutive relation of the fluid continuity equation, momentum conservation equation and fractional order Maxwell's non-newtonian fluid can So that such as following formula (2) (3) (4) are shown respectively：
Wherein, vfFluid velocity is represented, p represents pressure, and τ represents fluid pipe walls shear stress, and η represents fluid viscosity, μ tables Show modulus of shearing；λ=η/μ, represent the fluid Relaxation Characteristics time.ρfThe density of fluid is represented, α represents fraction of the stress to the time The exponent number of order derivative, β represent exponent number of the strain to time score order derivative；
Fractional DerivativeIt is defined as
Wherein, m=1,2,3 ..., and 0≤m-1≤α≤m.
It can be obtained from formula (3) and (4)：
Fluid is V with respect to the speed of solid,
U and u represents the displacement of fluid and solid along runner axis direction, v respectivelyfWithFluid and solid are represented respectively Along the speed of runner axis direction；
So, formula (6) can be written as following form：
Wherein,ν=η/ρf,
If the physical quantity in formula meets the form e of fluctuation-iωt, then, it can be obtained under cylindrical coordinates：
Wherein,ω represents wave field frequency, and α represents exponent number of the stress to the Fractional Derivative of time, β represents exponent number of the strain to time score order derivative, and μ represents modulus of shearing.
Above-mentioned formula (8) is Bessel equation, and its solution is：
Wherein,
According to fluid on tube wall without slip boundary condition, it may be determined that the constant C in formula (9), therefore formula (9) Can further it obtain：
After runner VELOCITY DISTRIBUTION is calculated using formula (10), the mean flow rate on cross section can be calculated：
In formula (11),A is runner radius, represents pore character size.
For Maxwell's non-newtonian fluid, fluid pipe walls shear stress can be calculated from current meter using following formula：
So, the ratio of shear stress and mean flow rate can be：
IntroduceAbove formula (13) can be simplified to obtain：
So,
In formula (15),It can be expressed asAndThat is characteristic time of viscous effect and non- The ratio of Newtonian fluid Relaxation Characteristics time, the parameter can determine that fluid is to be in non newtonian dissipation section or elastic region Between.
Work as DeWhen larger, pore-fluid is closer to viscous dissipation fluid (Newtonian fluid) feature；Work as DeWhen smaller, hole stream Body is more obvious the characteristics of similar solid elastic closer to viscoelastic fluid (non-newtonian fluid) feature.
Here can be according to pore character size a, fluid viscosity η and density pf, the basic parameter such as wave field frequency calculates Go out F (κ), then provide the dissipation factor of non-newtonian fluid, for subsequently establishing wave equation and predicting SVEL.
S102：Based on energy dissipation caused by the solid relative motion of stream, the pore media for including non-newtonian flow bulk effect is established Wave equation.
Specifically according to pore media kinetic energy, potential energy, fluid dissipation energy and the solid relative motion dissipation energy of stream, bag is established The pore media wave equation of the bulk effect containing non-newtonian flow
Flowing of the fluid inside pipeline causes two parts energy loss：As caused by flowing solid relative motion dissipation energy and Fluid dissipation energy caused by internal fluid shear action.
Wherein, dissipation energy caused by internal fluid shear action isWherein, τ represents shear stress, and γ is represented Shear strain,First derivatives of the shear strain γ to the time is represented,Represent porosity.
Therefore, calculating the formula (1) of the dissipation of non-newtonian fluid pore media can be expressed as follows：
Wherein b=b0F (κ) Gu be stream-relative motion dissipation factor, b0It is Newtonian fluid energy dissipation coefficient, F (κ) is high Frequency correction factor, u, U are the displacement of solid and fluid respectively, and τ is shear stress, and γ is shear strain, and φ is porosity.
Based on Hamiton's principle, the single hole medium wave equation for meeting Darcy's law for exporting non-newtonian fluid is as follows：
Wherein,α order derivatives, 0 ＜ α ＜ 1 are asked time t in expression.And in formula (17),
In above-mentioned formula (18), quality coefficient ρ11,ρ12,ρ22It is identical with Biot theories.
S103：Plane wave analysis is carried out to the pore media wave equation of the foundation, obtains imitating comprising non-newtonian fluid The dispersion relation answered, calculate velocity of longitudinal wave frequency dispersion and decay.
It is described that plane wave analysis is carried out to the pore media wave equation, obtain including the frequency dispersion of non-newtonian flow bulk effect Relation, including：The wave equation is transformed into frequency wavenumber domain, obtains the dispersion equation of compressional wave.Asked for by above-mentioned Complex wave number, the required predictor formula for being used to describe decay and the p-and s-wave velocity to dissipate and inverse quality factor can be obtained.
Consider isotropism pore media, elastic wave is propagated consistent along all directions.In the case of plane longitudinal wave propagation, ripple Dynamic non trivial solution is set to following form：
In formula (19), k is wave number, and w is circular frequency, and remaining each displacement component is zero.Wave equation is substituted into, i.e., In formula (17), following form can be obtained：
Introduce following mark：
A=c11c22-c12c21, B=c12d21+c21d12-c11d22-c22d11, C=d11d22-d12d21 (21)
Set up because formula (20) is permanent, the determinant of coefficient is necessary for zero, therefore has：
According to being defined as below, you can calculating speed frequency dispersion and attenuation curve.
Formula (23) represents velocity dispersion formula；In formula (24), Q represents quality factor, and 1/Q represents inverse quality factor； In formula (23) and (24), Im and Re represent imaginary part and real part respectively.
The velocity dispersion formula and attenuation curve can be used for description decay and the p-and s-wave velocity to dissipate.
Pore media SVEL computational methods of the application containing non-newtonian fluid are introduced with reference to practical application.
According to the geologic feature of fine and close oily reservoir, it have chosen one group of petrophysical parameter and calculated, the Biot with classics Theory analyzes.Calculating parameter is as shown in table 1.
Fractional order Maxwell fluid model is employed in calculating.Fig. 3 is that non-newtonian flow is considered in the embodiment of the present application The fast p wave interval velocity dispersion curve of body.In Fig. 3, abscissa represents frequency, and ordinate represents speed.From the figure 3, it may be seen that work as fractional order When derivative takes 1, fluid is classical Maxwell's non-Newtonian models, and p wave interval velocity, which has, significantly to move to right.This shows for Biot For model, the presence of non-newtonian fluid makes velocity dispersion curve be moved to high frequency section.0.8 and 0.6 are taken in Fractional Derivative α When, dispersion curve is further to high-frequency mobile, but mobile range is little.Fig. 4 is that non newtonian is considered in the embodiment of the present application The fast P attenuation curves of fluid.In Fig. 4, abscissa represents frequency, and ordinate represents inverse quality factor.From fig. 4, it can be seen that adopt During with classical Maxwell's non-Newtonian models, the attenuation peak of p wave interval velocity is substantially to high-frequency mobile, when the rank of Fractional Derivative When number takes 0.8 and 0.6, attenuation peak decreases, and moves right simultaneously, but amplitude all very littles.It is worth noting that, when consideration During fluid nonNewtonian percolation, to high-frequency mobile, this oozes Biot model velocities dispersion curve with increase fluid viscosity and reduction fluid Saturating rate effect is identical, and the mobility for being equal to fluid in pore media is deteriorated.Meanwhile with the reduction of Fractional Derivative exponent number, Also it is equal to fluid mobility variation, causes frequency dispersion attenuation curve to high-frequency mobile.This shows the non-newtonian flow that the application proposes Body wave equation can reflect pore-fluid feature, directly be built between wavefield velocity and the constitutive relation parameter of non-newtonian fluid Contact has been stood, the characteristic parameter of pore-fluid can be obtained by measuring the change of velocity of wave, is had a clear superiority compared with other method.
The single pore media non-Newtonian models petrophysical parameter of table 1
In the pore media SVEL computational methods containing non-newtonian fluid that above-described embodiment provides, foundation includes non-ox The pore media wave equation model of fluid parameter, can reflect pore character size, non-newtonian fluid constitutive relation feature Influence to SVEL.The method that the application provides considers the Fractional Derivative relation in non-newtonian fluid constitutive relation, By the change to fractional order, model is no longer limited to certain specific idealization non-Newtonian models, makes this method can With applied to the SVEL prediction case for including different type complex fluid effect.
Pore media SVEL computing device of the application containing non-newtonian fluid is described below.Fig. 5 is the application containing non- The module map of pore media SVEL computing device one embodiment of Newtonian fluid.As shown in figure 4, described device can wrap Include：Stream solid phase establishes module 502 and plane wavelength-division to kinergety dissipation factor computing module 501, pore media wave equation Analyse module 503.Wherein,
The stream solid phase is to kinergety dissipation factor computing module 501, for establishing pore media microchannel mould Type, according to the flow model export non-newtonian fluid inside dissipation energy, Fractional Derivative Maxwell fluid-solid it Between the mean flow rate of shear stress and fluid, according to the shear stress and the mean flow rate, calculate the solid relative motion energy of stream Measure dissipation factor.
The pore media wave equation establishes module 502, for based on energy dissipation caused by the solid relative motion of stream, building The vertical pore media wave equation for including non-newtonian flow bulk effect.
The plane wave analysis module 503, for carrying out plane wave analysis to the pore media wave equation of the foundation, Obtain including the dispersion relation of non-newtonian flow bulk effect, calculate velocity of longitudinal wave frequency dispersion and decay.
The pore media SVEL computing device and the present processes containing non-newtonian fluid that above-described embodiment provides Embodiment is corresponding, it is possible to achieve the application embodiment of the method and the technique effect for reaching the application embodiment of the method.
In the 1990s, the improvement for a technology can clearly distinguish be on hardware improvement (for example, Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So And as the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit. Designer nearly all obtains corresponding hardware circuit by the way that improved method flow is programmed into hardware circuit.Cause This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, PLD (Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate Array, FPGA)) it is exactly such a integrated circuit, its logic function is determined by user to device programming.By designer Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, without asking chip maker to design and make Special IC chip 2.Moreover, nowadays, substitution manually makes IC chip, and this programming is also used instead mostly " logic compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development Seemingly, and the source code before compiling also handy specific programming language is write, this is referred to as hardware description language (Hardware Description Language, HDL), and HDL is also not only a kind of, but have many kinds, such as ABEL (Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL (Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language) etc., VHDL (Very-High-Speed are most generally used at present Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art It will be apparent to the skilled artisan that only need method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages In, it is possible to it is readily available the hardware circuit for realizing the logical method flow.
Controller can be implemented in any suitable manner, for example, controller can take such as microprocessor or processing Device and storage can by the computer of the computer readable program code (such as software or firmware) of (micro-) computing device Read medium, gate, switch, application specific integrated circuit (Application Specific Integrated Circuit, ASIC), the form of programmable logic controller (PLC) and embedded microcontroller, the example of controller include but is not limited to following microcontroller Device：ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320, are deposited Memory controller is also implemented as a part for the control logic of memory.
It is also known in the art that in addition to realizing controller in a manner of pure computer readable program code, it is complete Entirely can by by method and step carry out programming in logic come controller with gate, switch, application specific integrated circuit, may be programmed The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind Hardware component, and it is used to realize that the device of various functions can also to be considered as the structure in hardware component to what is included in it.Or Even, it not only can be able to will be the software module of implementation method for realizing that the device of various functions is considered as but also can be Hardware Subdivision Structure in part.
System, device, module or the unit that above-described embodiment illustrates, it can specifically be realized by computer chip or entity, Or realized by the product with certain function.
For convenience of description, it is divided into various units during description apparatus above with function to describe respectively.Certainly, this is being implemented The function of each unit can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application On the part that is contributed in other words to prior art can be embodied in the form of software product, in a typical configuration In, computing device includes one or more processors (CPU), input/output interface, network interface and internal memory.The computer is soft Part product can include some instructions make it that a computer equipment (can be personal computer, server, or network Equipment etc.) perform method described in some parts of each embodiment of the application or embodiment.The computer software product can To be stored in internal memory, internal memory may include the volatile memory in computer-readable medium, random access memory (RAM) and/or the form such as Nonvolatile memory, such as read-only storage (ROM) or flash memory (flash RAM).Internal memory is computer The example of computer-readable recording medium.Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by Any method or technique come realize information store.Information can be computer-readable instruction, data structure, the module of program or its His data.The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM (SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only storage (ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc are read-only Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic cassette tape, tape magnetic rigid disk storage or Other magnetic storage apparatus or any other non-transmission medium, the information that can be accessed by a computing device available for storage.According to Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation Data-signal and carrier wave.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment Divide mutually referring to what each embodiment stressed is the difference with other embodiment.It is real especially for system For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as：Personal computer, clothes Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer including any of the above system or equipment DCE etc..
The application can be described in the general context of computer executable instructions, such as program Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by Task is performed and connected remote processing devices by communication network.In a distributed computing environment, program module can be with In the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's Spirit.

Claims (10)

  1. A kind of 1. pore media SVEL computational methods containing non-newtonian fluid, it is characterised in that including：
    Establish pore media microchannel model, according to the flow model export non-newtonian fluid inside dissipation energy, The mean flow rate of shear stress and fluid between Fractional Derivative Maxwell fluid-solid, according to the shear stress and institute Mean flow rate is stated, calculates stream solid phase to kinergety dissipation factor；
    Based on energy dissipation caused by the solid relative motion of stream, the pore media wave equation for including non-newtonian flow bulk effect is established；
    Plane wave analysis is carried out to the pore media wave equation of the foundation, the frequency dispersion comprising non-newtonian flow bulk effect is obtained and closes System, determines velocity of longitudinal wave frequency dispersion and decay.
  2. A kind of 2. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 1, it is characterised in that It is described according to the shear stress and the mean flow rate, calculate the energy dissipation coefficient and specifically include：
    According to fluid continuity equation, momentum conservation equation and fractional order Maxwell's non-newtonian fluid constitutive relation, pipe is obtained The space velocity distribution of road flowing and wall shear stress；
    It is distributed according to the space velocity and calculates non-newtonian fluid pipeline flowing mean flow rate；
    The ratio relation of mean flow rate and wall shear stress is flowed by non-newtonian fluid pipeline again, obtains energy dissipation system Number.
  3. A kind of 3. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 2, it is characterised in that Energy dissipation caused by the solid relative motion of stream is calculated by following formula：
    <mrow> <mi>D</mi> <mo>=</mo> <mfrac> <mrow> <mi>b</mi> <mi>F</mi> <mrow> <mo>(</mo> <mi>&amp;kappa;</mi> <mo>)</mo> </mrow> </mrow> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <mover> <mi>u</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>-</mo> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mover> <mi>u</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>-</mo> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>)</mo> </mrow> </mrow>
    In formula, D represents energy dissipation caused by the solid relative motion of stream；B represents stream solid phase to kinergety dissipation factor；F(κ) It is high frequency correction factor；U, U represent the displacement of solid and fluid respectively；Respectively during the displacement pair of expression solid and fluid Between first derivative, i.e. solid and fluid speed.
  4. A kind of 4. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 2, it is characterised in that
    The fluid continuity equation is：▽·vf=0；
    The momentum conservation equation is：
    The constitutive relation of the fractional order Maxwell non-newtonian fluid is：
    Wherein, vfFluid velocity is represented, p represents pressure, and τ represents fluid pipe walls shear stress, and η represents fluid viscosity, and μ represents to cut Shear modulu；λ=η/μ, represent the fluid Relaxation Characteristics time；ρfThe density of fluid is represented, α represents that stress is led to the fractional order of time Several exponent numbers, β represent exponent number of the strain to time score order derivative；
    Fractional DerivativeIt is defined as：
    <mrow> <msubsup> <mi>D</mi> <mi>t</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msup> <mi>d</mi> <mi>&amp;alpha;</mi> </msup> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <mi>dt</mi> <mi>&amp;alpha;</mi> </msup> </mrow> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>t</mi> </msubsup> <mfrac> <mrow> <msup> <mi>f</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> </mrow> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mrow> <mi>&amp;alpha;</mi> <mo>+</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mfrac> <mi>d</mi> <mi>&amp;tau;</mi> </mrow>
    Wherein, m=1,2,3 ..., and 0≤m-1≤α≤m.
  5. A kind of 5. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 2, it is characterised in that The VELOCITY DISTRIBUTION that pipeline flowing is calculated, is calculated according to following formula：
    <mrow> <mi>V</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mi>X</mi> <mrow> <mi>i</mi> <mi>&amp;omega;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <msub> <mi>J</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mi>r</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>J</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mi>a</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mi>A</mi> <mrow> <msup> <mi>k</mi> <mn>2</mn> </msup> <mi>v</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <msub> <mi>J</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mi>r</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>J</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mi>a</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>)</mo> </mrow> <mi>X</mi> </mrow>
    Wherein,
    V represents the speed of pipeline flowing；P represents pressure, and η represents fluid viscosity；ρfFluid density is represented, r is represented in cylindricality coordinate Radius；U represents the displacement of solid；λ=η/μ, represent the fluid Relaxation Characteristics time；ω represents wave field frequency, and α represents stress pair The exponent number of the Fractional Derivative of time, β represent exponent number of the strain to time score order derivative, and μ represents modulus of shearing, and a represents stream Road radius.
  6. A kind of 6. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 5, it is characterised in that The calculating mean flow rate obtains according to following formula：
    <mrow> <mover> <mi>V</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mi>&amp;pi;a</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>a</mi> </msubsup> <mn>2</mn> <mi>r</mi> <mi>V</mi> <mi>d</mi> <mi>r</mi> <mo>=</mo> <mi>K</mi> <mrow> <mo>(</mo> <mi>&amp;omega;</mi> <mo>)</mo> </mrow> <mi>X</mi> </mrow>
    In formula,Represent desired mean flow rate；A is runner radius.
  7. A kind of 7. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 2, it is characterised in that The fluid pipe walls shear stress, is calculated using following formula：
    <mrow> <mi>T</mi> <mo>=</mo> <mfrac> <mi>&amp;eta;</mi> <mi>A</mi> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>V</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>r</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>r</mi> <mo>=</mo> <mi>a</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mi>&amp;eta;</mi> <mrow> <mi>k</mi> <mi>v</mi> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msub> <mi>J</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mi>a</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>J</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mi>a</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mi>X</mi> </mrow>
    Wherein, τ represents fluid pipe walls shear stress, and η represents fluid viscosity, and V represents the speed of pipeline flowing, and r represents cylindricality coordinate In radius；ν=η/ρf, ρfRepresent fluid density； λ=η/μ, represent the fluid Relaxation Characteristics time；A is runner radius；ω represents wave field frequency, and α represents fraction of the stress to the time The exponent number of order derivative, β represent exponent number of the strain to time score order derivative, and μ represents modulus of shearing, and u represents the displacement of solid；p Represent pressure.
  8. A kind of 8. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 1, it is characterised in that The fluid dissipation energy isWherein, τ represents shear stress, and γ represents shear strain,Represent γ pairs of shear strain The first derivative of time；Represent porosity.
  9. A kind of 9. pore media SVEL computational methods containing non-newtonian fluid as claimed in claim 1, it is characterised in that It is described that plane wave analysis is carried out to pore media wave equation, obtain including the dispersion relation of non-newtonian flow bulk effect, including：Will The wave equation is transformed into frequency wavenumber domain, obtains the dispersion equation of compressional wave.
  10. A kind of 10. pore media SVEL computing device containing non-newtonian fluid, it is characterised in that including：Solid phase is flowed to fortune Energy dissipation factor computing module, pore media wave equation establish module and plane wave analysis module；Wherein,
    The stream solid phase is to kinergety dissipation factor computing module, for establishing pore media microchannel model, according to Sheared between dissipation energy, Fractional Derivative Maxwell fluid-solid inside the flow model export non-newtonian fluid The mean flow rate of stress and fluid, according to the shear stress and the mean flow rate, calculate stream solid phase and kinergety is dissipated Coefficient；
    The pore media wave equation establishes module, for being included based on energy dissipation, foundation caused by the solid relative motion of stream The pore media wave equation of non-newtonian flow bulk effect；
    The plane wave analysis module, for carrying out plane wave analysis to the pore media wave equation of the foundation, wrapped The dispersion relation of the bulk effect containing non-newtonian flow, calculate velocity of longitudinal wave frequency dispersion and decay.
CN201511019018.7A 2015-12-30 2015-12-30 A kind of pore media SVEL computational methods and device containing non-newtonian fluid Active CN105510966B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511019018.7A CN105510966B (en) 2015-12-30 2015-12-30 A kind of pore media SVEL computational methods and device containing non-newtonian fluid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511019018.7A CN105510966B (en) 2015-12-30 2015-12-30 A kind of pore media SVEL computational methods and device containing non-newtonian fluid

Publications (2)

Publication Number Publication Date
CN105510966A CN105510966A (en) 2016-04-20
CN105510966B true CN105510966B (en) 2018-02-02

Family

ID=55719074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511019018.7A Active CN105510966B (en) 2015-12-30 2015-12-30 A kind of pore media SVEL computational methods and device containing non-newtonian fluid

Country Status (1)

Country Link
CN (1) CN105510966B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107704423B (en) * 2017-08-21 2020-09-29 重庆工商大学 User-defined function-based dual-field coupling analysis method
CN110703322B (en) * 2019-10-10 2020-11-13 清华大学 Wave propagation processing method, device and equipment
CN110879249B (en) * 2019-11-29 2022-06-03 南昌航空大学 System device for measuring limit shear modulus of biological Newtonian fluid
CN111025396B (en) * 2020-01-06 2021-11-05 中国石油化工股份有限公司 Oil reservoir physical property parameter earthquake prediction method based on artificial intelligence algorithm
CN115270651B (en) * 2022-06-20 2024-03-15 北京科技大学 Monocular video-oriented non-Newtonian fluid simulation reconstruction method
CN115774286B (en) * 2022-12-06 2023-09-05 上海电子信息职业技术学院 Longitudinal wave speed prediction method for pore medium containing viscoelastic fluid and application of longitudinal wave speed prediction method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412323A (en) * 2013-07-05 2013-11-27 中国石油天然气股份有限公司 Longitudinal wave speed prediction method of rock rich in solid organic matter
CN103984027A (en) * 2014-03-28 2014-08-13 清华大学 Rock longitudinal wave speed prediction method based on ellipsoid double porosity model

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2510685B (en) * 2012-12-12 2015-11-18 Schlumberger Holdings Measurement of fluid properties

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412323A (en) * 2013-07-05 2013-11-27 中国石油天然气股份有限公司 Longitudinal wave speed prediction method of rock rich in solid organic matter
CN103984027A (en) * 2014-03-28 2014-08-13 清华大学 Rock longitudinal wave speed prediction method based on ellipsoid double porosity model

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A time-domain numerical method for Biot-JKD poroelastic waves in 2D heterogeneous media;Emilie Blanc et al.;《The Journal of the Acoustical Society of America》;20131231;第1-17页 *
UNSTEADY FLOW OF VISCOELASTIC FLUID DUE TO IMPULSIVE MOTION OF PLATE;AMIR KHAN et al.;《ASIAN JOURNAL OF MATHEMATICS AND APPLICATIONS》;20141231;第1-8页 *
非饱和岩石中的纵波频散与衰减：双重孔隙介质波传播方程;巴晶等;《地球物理学报》;20120131;第55卷(第1期);第219-229页 *

Also Published As

Publication number Publication date
CN105510966A (en) 2016-04-20

Similar Documents

Publication Publication Date Title
CN105510966B (en) A kind of pore media SVEL computational methods and device containing non-newtonian fluid
Mangeney‐Castelnau et al. Numerical modeling of avalanches based on Saint Venant equations using a kinetic scheme
Geiger et al. Black-oil simulations for three-component, three-phase flow in fractured porous media
Cueto-Felgueroso et al. A phase-field model of two-phase Hele-Shaw flow
Iassonov et al. A model for enhanced fluid percolation in porous media by application of low‐frequency elastic waves
Castro‐Orgaz et al. Nonhydrostatic granular flow over 3‐D terrain: New Boussinesq‐type gravity waves?
Malek-Mohammadi et al. New methodology for laboratory generation of solitary waves
Moreno-Casas et al. Computation of the Basset force: recent advances and environmental flow applications
CN107525553B (en) A kind of method and device of determining heterogeneous fluid composition flow rate
Holmedal et al. Wave-induced steady streaming, mass transport and net sediment transport in rough turbulent ocean bottom boundary layers
CN110703322B (en) Wave propagation processing method, device and equipment
Dorrell et al. The structure of the deposit produced by sedimentation of polydisperse suspensions
Helland et al. Footprints of spontaneous fluid redistribution on capillary pressure in porous rock
Deng et al. Dynamics and dislodgment from pore constrictions of a trapped nonwetting droplet stimulated by seismic waves
Nienhuis et al. Wavelength selection and symmetry breaking in orbital wave ripples
Zhao et al. Spontaneous imbibition in a square tube with corner films: Theoretical model and numerical simulation
Minatti et al. SPH numerical approach in modelling 2D muddy debris flow
Jing et al. Runout scaling and deposit morphology of rapid mudflows
Vahab et al. An X-FEM formulation for the optimized graded proppant injection into hydro-fractures within saturated porous media
Afzal et al. Three‐dimensional streaming in the seabed boundary layer beneath propagating waves with an angle of attack on the current
CN107367763A (en) Determine the method and device of rock stratum Biot coefficients
Kelly et al. The structure and entrainment characteristics of partially confined gravity currents
Liang et al. Study of the dilatancy/contraction mechanism of landslide fluidization behavior using an initially saturated granular column collapse simulation
Rozhko et al. Contact line friction and surface tension effects on seismic attenuation and effective bulk moduli in rock with a partially saturated crack
Maleki et al. Tracking fluid interfaces in primary cementing of surface casing

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant