Yes, it is possible. However, you have to accept some errors caused by a simplification.
That means: You have to assume a simplified opamp open-loop gain expression A(s)=GBW/s=wt/s (wt=transit frequency).
This simplification would lead to a second-order lowpass response. But as you can see - your simulation reveals that the phase response exceeds -180 deg.
This is, of course, caused by the higher-order opamp model. However, it can be expected that the error will be relatively small.
The procedure is as follows:
* Calculate the second-order lowpass function of the circuit using A(s)=wt/s and identify the expression for the pole frequency wp (which contains wt=GBW)
* From the simulation (Bode diagram) find the pole frequency wp. Assuming a second-order response, this frequency wp can be found for phase=-90 deg.
* Compare both wp values (calculated vs. simulated) and solve for wt=GBW.
For determination of the transfer function I have used a symbolic simulation program.
After some manipulations I arrived at
H(s)=Vo/Vin=1/D(s) and D(s)=1+s*(K/wt)+[s*(K/wt)]^2 with K=1+R2/R1
Comparing the denominator D(s) with the classical second-order form D(s)=1+s/(Qp*wp)+[s/wp]^2
we immediately find: wp=wt/K or wt=GBW=wp*K.
From your simulation:
Pole frequency wp is approximately 400...420 kHz (find exact value at phase=-90 deg) .
Thus: wt=GBW=wp*K=wp*(1+R2/R1).
Hope it helps.