{"id":"1bd5f875-3191-4b7e-bb3c-c31fce57942c","arxiv_id":"1908.10495","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A fractional-derivative Maxwellian fluid is inserted into Biot's poroelastic wave equations, producing P- and S-wave dispersion and attenuation curves with viscoelastic resonance peaks.","lead":"A model that extends Biot's theory of waves in fluid-filled rocks to a stretchy, time-memory (fractional Maxwell) saturating fluid predicts velocity oscillations and attenuation peaks that ordinary viscous fluids do not show. The work matters for seismic interpretation in heavy-oil and polymer-flooded reservoirs, where non-Newtonian pore fluids are common.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pore radius for sandstone wave predictions is unstated: with a=1 mm vs. a≈0.94 μm from κ0, the Deborah number changes by ~10^6, so the predicted resonance peaks may not apply to the stated medium.","rationale":"The reader's verdict is CONDITIONAL, and I agree that the paper should not be accepted without revision. My concern is narrower and sharper than the reader's weakest assumption: it is not primarily that real pore networks deviate from a single capillary, but that even within the single-capillary model the paper never fixes the capillary radius for the sandstone wave calculations. The two plausible radii (1 mm illustration vs. 0.94 μm from Table 1's permeability) imply Deborah numbers differing by six orders and dimensionless frequencies differing by three orders at the same physical frequency, so the predicted attenuation peaks may not survive a recomputation with the permeability-consistent radius. This is a checkable, internal consistency issue rather than a matter of theoretical taste. The model deserves credit for recovering Biot's Newtonian limit and for reproducing the peak locations in the 25 mm tube experiment; those results support the tube-flow part of the paper. However, the porous-medium wave predictions, which are the main contribution, cannot be evaluated without the missing radius specification and a recalculation. Therefore, the verdict should remain CONDITIONAL, with the additional requirement that the authors disclose the radius and Deborah number used and recompute the sandstone curves with the permeability-consistent value.","tokens_in":15647,"tokens_out":16101,"duration_ms":164311,"concrete_test":"Ask the authors to report the exact channel radius a and Deborah number used in Eqs. (13)-(16) for Figures 9-11, then recompute the P- and S-wave attenuation curves with a obtained from κ0=1.1e-13 m² via κ0=a²/8 (a≈0.94 μm, De≈1.2e11). If the multiple attenuation peaks in Figs. 9b-11b persist at the stated frequencies, the central claim survives; if the peaks shift, broaden, or vanish, the resonance-peak claim is an artifact of the illustrative 1 mm radius and the verdict should move to REJECT.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central numerical claim—multiple attenuation peaks and the dissipative-to-elastic transition in a sandstone saturated by CPyCl/NaSal (Figs. 9-11)—depends on the dimensionless frequency χ=a√(ω/ν) and the Deborah number De=νλ/a² that enter F(χ) in Eq. (12). Section IV.A fixes a=1.0 mm for the illustrative F(χ) curves in Fig. 2, giving De≈1.1e5 for CPyCl/NaSal (η=60 Pa·s, ρ=1050 kg/m³, λ=1.9 s). For the French Vosgian sandstone, Table 1 gives κ0=1.1e-13 m², so the Hagen-Poiseuille relation κ0=a²/8 used in Section III yields a≈0.94 μm and De≈1.2e11—six orders of magnitude larger. The paper never states which radius is used when F(χ) is inserted into the wave equations (Eqs. 13-16) for Figs. 9-11. At a fixed frequency, χ scales as a: with a=0.94 μm, χ≈0.31 at 1 GHz, whereas with a=1 mm, χ≈332 at 1 GHz. Thus the two choices probe entirely different branches of F(χ): the illustrative radius places the plotted band in the resonant/elastic regime, while the permeability-consistent radius keeps χ<1 up to about 10 GHz. If the authors used a=1 mm, the implied static permeability is 1.25e-7 m², contradicting Table 1 by six orders. The experimental tube in Fig. 5 uses a=25 mm (De≈1.7e2), yet a third regime. The wave-propagation predictions are therefore unreproducible as written, and the claimed multiplicity of attenuation peaks may be an artifact of an unstated, permeability-inconsistent radius.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a viscous dissipation function F(χ) for a fractional-derivative Maxwellian fluid flowing in a circular tube, incorporates it into Biot's poroelastic wave equations, and obtains analytical expressions for P- and S-wave dispersion and attenuation. It shows that the Newtonian limit recovers Biot's classical results, and it applies the model to sandstone saturated with brine and with the viscoelastic solution CPyCl/NaSal. The predicted tube-center fluid velocities are compared with laboratory data from Castrejon-Pita et al. (2003), and the paper argues that viscoelastic fluid behavior produces multiple attenuation peaks and a dissipative-to-elastic transition that are absent for Newtonian saturants.","tokens_in":16032,"tokens_out":6980,"duration_ms":76556,"significance":"If the results hold, this is a useful extension of Biot theory to non-Newtonian saturants: the derivation of F(χ) from the fractional Maxwell constitutive equation is self-contained, and the check that the Newtonian limit reproduces Biot's F(χ) is a genuine strength. The paper also offers a concrete mechanism for velocity oscillations in viscoelastic-fluid-saturated rocks. However, the sandstone wave predictions depend on an unspecified and potentially inconsistent pore radius, and the experimental agreement in Figure 5 is obtained only after tuning the fractional order α to 1.03; these issues currently limit the strength of the central claims.","major_comments":[{"comment":"The pore radius used in the sandstone wave predictions in Figures 9-11 is never stated. Section IV.A fixes a=1.0 mm to illustrate F(χ), while Section IV.C uses French Vosgian sandstone with κ=1.1e-13 m² and, in Figure 3, invokes the Hagen-Poiseuille relation κ0=a²/8, which implies a≈0.94 μm. The two radii differ by three orders of magnitude, changing De=νλ/a² from about 1.1e5 to about 1.2e11 and changing χ=a√(ω/ν) by a factor of 10^3 at fixed frequency. With the permeability-consistent radius, χ remains below unity up to about 10 GHz, so the resonant/elastic regime displayed in Figures 2 and 9-11 may not be probed at all. The multiple attenuation peaks in Figures 9-11 may therefore be artifacts of an unstated, permeability-inconsistent radius. Please state the radius used and recompute the sandstone results using a derived from κ0=a²/8, or justify an alternative pore-geometry relation.","section":"Section IV.A, IV.C, Eq. (12), Table 1"},{"comment":"The claimed consistency with laboratory observations is weakened because the fractional derivative order α is adjusted from 1.0 to 1.03 to make the predicted peaks lower and closer to the experimental data. This is a calibration step rather than an independent prediction, so the abstract's statement that 'the predicted fluid velocities are consistent with the laboratory observations' overstates the evidence. Please report the fitting procedure, the misfit before and after adjustment, and a sensitivity analysis over α and β, and reframe the comparison as a calibration plus qualitative validation rather than a direct prediction.","section":"Section IV.B, Figure 5"},{"comment":"The wave equations insert the friction term ηϕ²F(χ)/κ using the static permeability κ from Table 1, while F(χ) is evaluated with a pore radius a that is not specified for Figures 9-11. If the radius inside F(χ) is not the one implied by κ through κ0=a²/8, the friction term mixes incompatible length scales. This is not just a presentation issue: because De and χ change by orders of magnitude with a, the qualitative behavior of the predicted dispersion and attenuation curves depends directly on this choice. The manuscript must make the value of a explicit for every figure and ensure it is consistent with the stated permeability.","section":"Section III, Appendix, Eqs. (13)-(16)"}],"minor_comments":[{"comment":"Typo: 'pock parameters' should be 'rock parameters'.","section":"Section IV.C"},{"comment":"The fractional derivative in Eq. (4) is defined as a Riemann-Liouville derivative; please state explicitly whether the constitutive equation uses this definition or the Caputo definition, since initial-condition treatment differs.","section":"Section II, Eq. (4)"},{"comment":"The text mentions 'α=1 and β=1.5' while the caption says fractional order pairs (1.0,1.0) and (1.0,1.5); please make the notation for the order pair consistent throughout.","section":"Section IV.A, Figure 2 caption"},{"comment":"Figure 5 would be much more informative with error bars or individual experimental points; the current comparison is only qualitative, and the claim 'peak locations are consistent' is not quantified.","section":"Section IV.B"},{"comment":"The statement that the fdMaxwell model 'would be valid even at extremely high frequency' is stronger than the evidence: the model still ignores scattering and other high-frequency effects that limit Biot theory, and the large characteristic frequency alone does not establish validity.","section":"Section III"}],"recommendation":"major_revision","confidential_remarks":"The central derivation appears sound and the Newtonian-limit check is a good sign, but the sandstone numerical claims need to be recomputed with a stated, permeability-consistent pore radius. If the multiple attenuation peaks disappear at the correct radius, the paper's main novelty is reduced to a theoretical extension whose experimental support is limited to a single tube-flow comparison with tuned parameters."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely new: it generalizes Biot's viscous dissipation function to a fractional-derivative Maxwell fluid, deriving F(χ) from the constitutive equation with Bessel-function analysis, and it recovers the Newtonian limit. That part is solid.\n\nThe reason I'd be cautious: the numerical section that claims multiple attenuation peaks in sandstone saturated by CPyCl/NaSal (Figs. 9-11) never states which pore radius was used. In Sec. IV.A, a=1.0 mm is used to illustrate F(χ), giving De≈1.1e5. But the sandstone in Table 1 has κ0=1.1e-13 m², and the paper itself gives the Hagen-Poiseuille relation κ0=a²/8, so a≈0.94 μm and De≈1.2e11. Those two choices put you on completely different branches of F(χ): with the small radius, χ stays below 1 up to ~10 GHz, while with a=1 mm, χ is hundreds at 1 GHz. The paper doesn't say which was used. If it was 1 mm, the implied permeability contradicts Table 1 by six orders; if it was the small radius, the resonance peaks plotted may not appear in that frequency band. Either way, Figs. 9-11 are not reproducible as written. This is the load-bearing soft spot.\n\nSecond soft spot: the abstract says the predicted fluid velocities are 'consistent with' the lab data, but the actual match required tuning α from 1.0 to 1.03, and the text admits the peak amplitudes are not perfect. That's a fit, not an independent confirmation. Should be disclosed as such.\n\nThird: the claim that the fdMaxwell model is valid at 'extremely high frequency' because the characteristic frequency grows exponentially is overreach. The scattering limit of Biot's theory is a physical grain-size effect; redefining the characteristic frequency doesn't remove it.\n\nThe derivation itself is the contribution, and the experimental tube-velocity peak locations are a point in its favor. The paper deserves a serious referee, but as written it needs major revision: state the radius used, rerun the sandstone predictions with a permeability-consistent radius and report sensitivity, and temper the experimental and high-frequency claims. No code or data is provided, which makes the numerical curves even harder to check.","headline":"The fractional-Maxwell extension of Biot's dissipation function is a real derivation, but the sandstone wave predictions are unreproducible as written because the pore radius is never reconciled with Table 1's permeability.","tokens_in":16581,"tokens_out":3012,"would_cite":false,"duration_ms":30106,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By inserting a fractional-derivative Maxwell constitutive law into Biot's poroelastic equations, this paper claims that viscoelastic pore fluids create multiple attenuation peaks and velocity oscillations that a Newtonian fluid cannot…","keywords":["wave propagation","poroelastic media","fractional derivative Maxwell model","viscoelastic fluid","seismic attenuation","dynamic permeability","Deborah number","Biot theory"],"falsifier":"Measure P-wave phase velocity and attenuation in a rock of known porosity and permeability saturated with a viscoelastic fluid whose relaxation time and viscosity are measured independently, sweeping frequency from tens of hertz to tens of megahertz. The central claim fails if the attenuation curve shows one smooth Biot-like peak rather than multiple peaks at lower frequencies, or if the measured centerline velocity peaks in an oscillating tube occur at frequencies that do not match the model's predicted locations.","tokens_in":15409,"feed_emoji":"🌊","tokens_out":9918,"duration_ms":86992,"temperature":0.7,"pith_summary":"The paper attempts to show that the standard assumption of a Newtonian pore fluid is too restrictive for rocks saturated with viscoelastic fluids, and that replacing it with a fractional derivative Maxwell model changes the predicted wave behavior in a way experiments support. The authors derive a frequency-dependent frictional dissipation function $F(\\chi)$ from the fractional Maxwell constitutive relation and insert it into Biot's poroelastic wave equations, yielding analytical dispersion and attenuation expressions for P-, slow P-, and S-waves. The central result is that a Maxwellian saturant produces pulsatile velocity enhancement and multiple attenuation peaks, while a Newtonian saturant produces a single smooth peak; the fractional derivative orders control the transition from dissipative to elastic fluid-solid coupling. In the Newtonian limit the equations recover Biot's theory, and for a Maxwellian fluid the predicted velocity peaks in an oscillating tube align with laboratory observations. If right, the work matters because seismic and ultrasonic measurements in reservoirs containing heavy oil or polymer solutions would carry a viscoelastic fingerprint, and the frequency location of attenuation peaks would report on fluid rheology rather than just viscosity.","feed_headline":"Maxwellian fluid leaves multiple peaks in rock wave attenuation","feed_subtitle":"Heavy-oil pore fluids should ring at resonance frequencies where Newtonian theory predicts one smooth peak.","key_machinery":"The load-bearing object is the viscous dissipation function $F(\\chi)$ of Eq. (12), a parameter-free correction to the friction term in Biot's equations. It comes from solving oscillating flow of the fractional Maxwell fluid in a circular tube with no-slip walls, and it depends on the dimensionless frequency $\\chi=a\\sqrt{\\omega/\\nu}$, the Deborah number $De=\\lambda/\\lambda_\\nu$, and the fractional derivative orders $\\alpha$ and $\\beta$. Replacing the Newtonian viscous term by $\\eta\\phi^2 F(\\chi)/\\kappa$ in the Biot equations makes the rheology of the saturating fluid control the wave field: $F(\\chi)$ determines the dispersion and attenuation of the P-, slow P-, and S-waves, sets the frequency-dependent characteristic frequency, and encodes the transition from a dissipative to an elastic coupling regime. It is the single point where the constitutive model enters the wave equations.","core_discovery":"The paper's central claim is that a poroelastic medium saturated by a fractional derivative Maxwellian fluid supports elastic waves whose dispersion and attenuation are qualitatively different from Biot's Newtonian prediction, and that this difference is controlled by a frequency-dependent frictional dissipation function derived from the fractional Maxwell constitutive law. In the Newtonian limit the equations reduce to Biot's theory; for a viscoelastic fluid such as CPyCl/NaSal they produce pulsatile velocity enhancement and multiple attenuation peaks at resonance frequencies, with the fractional derivative orders determining whether the fluid-solid coupling is in a dissipative or an elastic regime. The same formulation predicts centerline fluid velocities in an oscillating tube whose peak locations match laboratory measurements for a Maxwellian fluid, which the paper offers as evidence that viscoelastic fluid effects can account for velocity oscillations observed in experiments. If correct, the result means that seismic and ultrasonic wave trains in rocks saturated by heavy oil or polymer solutions carry a measurable viscoelastic fingerprint that a Newtonian model cannot reproduce.","pith_inferences":["The same frequency-dependent friction function could be used to invert field seismic or sonic-log data for a reservoir fluid's relaxation time and fractional order, since the predicted peak locations are sensitive to those parameters; the paper does not build that inversion.","The single-tube no-slip geometry is likely to understate how much real pore networks and surface roughness smear the resonance peaks, so the sharpest test of the mechanism would compare model predictions with waves measured in rocks of independently characterized pore geometry.","A practical extension, not pursued in the paper, is to test the low-attenuation slow P-wave prediction directly in a laboratory sandstone saturated with a wormlike micelle solution.","If the fractional derivative orders can be tied to the fractal geometry of the pore network, wave dispersion could become a probe of pore structure rather than only fluid rheology; that link is suggested by work the paper cites but not established here."],"forward_implications":["When the pore fluid is viscoelastic, attenuation curves for P- and S-waves develop multiple peaks and the main transition shifts to lower frequencies than Biot's theory predicts for the same viscosity.","Increasing the fractional derivative order $\\beta$ smooths the resonance peaks and moves the attenuation peak to intermediate frequencies, so the model spans the range from Newtonian to ideal Maxwell behavior.","For Newtonian pore fluids the new equations reduce to Biot's results, so existing interpretations remain valid where the fluid is truly Newtonian.","The slow P-wave in a Maxwellian-fluid-saturated rock has much lower attenuation at low frequencies, so it should be observable where it is usually lost in Newtonian saturation.","The characteristic frequency below which the poroelastic model is valid becomes frequency dependent and grows exponentially with frequency for the Maxwellian fluid, extending the model's range beyond Biot's high-frequency limit."],"supporting_citations":[{"why":"Supplies the poroelastic wave equations and the high-frequency friction correction that this paper extends to a fractional Maxwell fluid.","marker":"Biot 1956a, b"},{"why":"Introduces the fractional Maxwell constitutive relaxation and retardation functions used as the fluid stress-strain relation.","marker":"Friedrich 1991"},{"why":"Provides the generalized fractional viscoelastic model formulation from which the constitutive equation is taken.","marker":"Schiessel et al. 1995"},{"why":"Establishes the resonant enhancement of viscoelastic tube flow and the Deborah-number scaling that underpins the predicted peaks.","marker":"Rio, Haro, and Whitaker 1998"},{"why":"Predicts resonant peaks in the dynamic response of a viscoelastic fluid in a vibrating tube, the behavior the wave model carries into porous media.","marker":"Tsiklauri and Beresnev 2001"},{"why":"Supplies the experimental velocity data for glycerol and CPyCl/NaSal used to validate the predicted peak locations.","marker":"Castrejon-Pita et al. 2003"},{"why":"Provides the French Vosgian sandstone parameters used in the numerical dispersion and attenuation calculations.","marker":"Bacri and Salin 1986"},{"why":"Treats elastic waves in a Maxwell-fluid-saturated porous medium, the baseline that this fractional model generalizes.","marker":"Tsiklauri and Beresnev 2003"}],"fun_headline_variants":["Maxwell fluid triggers extra wave attenuation peaks in rocks","Viscoelastic pore fluids add resonance peaks to wave attenuation","Maxwell model explains velocity oscillations in porous media","Fractional Maxwell fluid shifts seismic wave attenuation to multiple peaks","Poroelastic rock with Maxwell fluid shows extra attenuation peaks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the friction felt by waves in a real rock is adequately represented by the solution for a single straight cylindrical pore with no-slip walls, under a long-wave approximation, once that friction is inserted into the Biot equations.","fun_headline_variants_meta":{"raw":{"variants":["Maxwell fluid triggers extra wave attenuation peaks in rocks","Viscoelastic pore fluids add resonance peaks to wave attenuation","Maxwell model explains velocity oscillations in porous media","Fractional Maxwell fluid shifts seismic wave attenuation to multiple peaks","Poroelastic rock with Maxwell fluid shows extra attenuation peaks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000734,"raw_usage":{"total_tokens":3268,"prompt_tokens":920,"completion_tokens":2348,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":536,"completion_tokens_details":{"reasoning_tokens":2284}},"tokens_in":536,"tokens_out":2348,"duration_ms":18008,"temperature":1.0,"reasoning_tokens":2284,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:41:54.218345+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure P-wave phase velocity and attenuation in a rock of known porosity and permeability saturated with a viscoelastic fluid whose relaxation time and viscosity are measured independently, sweeping frequency from tens of hertz to tens of megahertz. The central claim fails if the attenuation curve shows one smooth Biot-like peak rather than multiple peaks at lower frequencies, or if the measured centerline velocity peaks in an oscillating tube occur at frequencies that do not match the model's predicted locations.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the French Vosgian sandstone parameters used in the numerical dispersion and attenuation calculations."}],"review_version":1}